The genetic background to bipolar disorder (BPD) has been attributed to different genetic and genomic risk factors. In the present study we hypothesized that inherited copy number variations (CNVs) contribute to susceptibility of BPD. We screened 637 BP-pedigrees from the NIMH Genetic Initiative and gave priority to 46 pedigrees. In this subsample we performed parametric and non-parametric genome-wide linkage analyses using ~21,000 SNP-markers. We developed an algorithm to test for linkage restricted to regions with CNVs that are shared within and across families. For the combined CNV and linkage analysis, one region on 19q13 survived correction for multiple comparisons and replicates a previous BPD risk locus. The shared CNV map to the pregnancy-specific glycoprotein (PSG) gene, a gene-family not previously implicated in BPD etiology. Two SNPs in the shared CNV are likely transcription factor binding sites and are linked to expression of an F-box binding gene, a key regulator of neuronal pathways suggested to be involved in BPD etiology. Our CNV-weighted linkage approach identifies a risk locus for BPD on 19q13 and forms a useful tool to future studies to unravel part of the genetic vulnerability to BPD.

Tumor-associated macrophages (TAMs) are attractive targets for immunotherapy. Recently, studies in animal models showed that treatment with an anti-TAM antibody directed against the scavenger receptor MARCO resulted in suppression of tumor growth and metastatic dissemination. Here we investigated the expression of MARCO in relation to other macrophage markers and immune pathways in a non-small cell lung cancer (NSCLC) cohort (n = 352). MARCO, CD68, CD163, MSR1 and programmed death ligand-1 (PD-L1) were analyzed by immunohistochemistry and immunofluorescence, and associations to other immune cells and regulatory pathways were studied in a subset of cases (n = 199) with available RNA-seq data. We observed a large variation in macrophage density between cases and a strong correlation between CD68 and CD163, suggesting that the majority of TAMs present in NSCLC exhibit a protumor phenotype. Correlation to clinical data only showed a weak trend toward worse survival for patients with high macrophage infiltration. Interestingly, MARCO was expressed on a distinct subpopulation of TAMs, which tended to aggregate in close proximity to tumor cell nests. On the transcriptomic level, we found a positive association between MARCO gene expression and general immune response pathways including strong links to immunosuppressive TAMs, T-cell infiltration and immune checkpoint molecules. Indeed, a higher macrophage infiltration was seen in tumors expressing PD-L1, and macrophages residing within tumor cell nests co-expressed MARCO and PD-L1. Thus, MARCO is a potential new immune target for anti-TAM treatment in a subset of NSCLC patients, possibly in combination with available immune checkpoint inhibitors.

The etiology of hereditary ataxia syndromes is heterogeneous, and the mechanisms underlying these disorders are often unknown. Here, we utilized exome sequencing in two siblings with progressive ataxia and muscular weakness and identified a novel homozygous splice mutation (c.3020-1G > A) in neurofascin ( NFASC). In RNA extracted from fibroblasts, we showed that the mutation resulted in inframe skipping of exon 26, with a deprived expression of the full-length transcript that corresponds to NFASC isoform NF186. To further investigate the disease mechanisms, we reprogrammed fibroblasts from one affected sibling to induced pluripotent stem cells, directed them to neuroepithelial stem cells and finally differentiated to neurons. In early neurogenesis, differentiating cells with selective depletion of the NF186 isoform showed significantly reduced neurite outgrowth as well as fewer emerging neurites. Furthermore, whole-cell patch-clamp recordings of patient-derived neuronal cells revealed a lower threshold for openings, indicating altered Na+ channel kinetics, suggesting a lower threshold for openings as compared to neuronal cells without the NFASC mutation. Taken together, our results suggest that loss of the full-length NFASC isoform NF186 causes perturbed neurogenesis and impaired neuronal biophysical properties resulting in a novel early-onset autosomal recessive ataxia syndrome.

The consequences of the Neolithic transition in Europe--one of the most important cultural changes in human prehistory--is a subject of great interest. However, its effect on prehistoric and modern-day people in Iberia, the westernmost frontier of the European continent, remains unresolved. We present, to our knowledge, the first genome-wide sequence data from eight human remains, dated to between 5,500 and 3,500 years before present, excavated in the El Portalón cave at Sierra de Atapuerca, Spain. We show that these individuals emerged from the same ancestral gene pool as early farmers in other parts of Europe, suggesting that migration was the dominant mode of transferring farming practices throughout western Eurasia. In contrast to central and northern early European farmers, the Chalcolithic El Portalón individuals additionally mixed with local southwestern hunter-gatherers. The proportion of hunter-gatherer-related admixture into early farmers also increased over the course of two millennia. The Chalcolithic El Portalón individuals showed greatest genetic affinity to modern-day Basques, who have long been considered linguistic and genetic isolates linked to the Mesolithic whereas all other European early farmers show greater genetic similarity to modern-day Sardinians. These genetic links suggest that Basques and their language may be linked with the spread of agriculture during the Neolithic. Furthermore, all modern-day Iberian groups except the Basques display distinct admixture with Caucasus/Central Asian and North African groups, possibly related to historical migration events. The El Portalón genomes uncover important pieces of the demographic history of Iberia and Europe and reveal how prehistoric groups relate to modern-day people.

Distinct types of dorsal root ganglion sensory neurons may have unique contributions to chronic pain. Identification of primate sensory neuron types is critical for understanding the cellular origin and heritability of chronic pain. However, molecular insights into the primate sensory neurons are missing. Here we classify non-human primate dorsal root ganglion sensory neurons based on their transcriptome and map human pain heritability to neuronal types. First, we identified cell correlates between two major datasets for mouse sensory neuron types. Machine learning exposes an overall cross-species conservation of somatosensory neurons between primate and mouse, although with differences at individual gene level, highlighting the importance of primate data for clinical translation. We map genomic loci associated with chronic pain in human onto primate sensory neuron types to identify the cellular origin of chronic pain. Genome-wide associations for chronic pain converge on two different neuronal types distributed between pain disorders that display different genetic susceptibilities, suggesting both unique and shared mechanisms between different pain conditions.

Cancer is driven by various events leading to cell differentiation and disease progression. Molecular tools are powerful approaches for describing how and why these events occur. With the growing field of next-generation DNA sequencing, there is an increasing need for high-quality nucleic acids derived from human cells and tissues-a prerequisite for successful cell profiling. Although advances in RNA preservation have been made, some of the largest biobanks still do not employ RNA blood preservation as standard because of limitations in low blood-input volume and RNA stability over the whole gene body. Therefore, we have developed a robust protocol for blood preservation and long-term storage while maintaining RNA integrity. Furthermore, we explored the possibility of using the protocol for preserving rare cell samples, such as circulating tumor cells. The results of our study confirmed that gene expression was not impacted by the preservation procedure (r(2) > 0.88) or by long-term storage (r(2) = 0.95), with RNA integrity number values averaging over 8. Similarly, cell surface antigens were still available for antibody selection (r(2) = 0.95). Lastly, data mining for fusion events showed that it was possible to detect rare tumor cells among a background of other cells present in blood irrespective of fixation. Thus, the developed protocol would be suitable for rare blood cell preservation followed by RNA sequencing analysis.

The Atlantic herring is a model species for exploring the genetic basis for ecological adaptation, due to its huge population size and extremely low genetic differentiation at selectively neutral loci. However, such studies have so far been hampered because of a highly fragmented genome assembly. Here, we deliver a chromosome-level genome assembly based on a hybrid approach combining a de novo Pacific Biosciences (PacBio) assembly with Hi-C-supported scaffolding. The assembly comprises 26 autosomes with sizes ranging from 12.4 to 33.1 Mb and a total size, in chromosomes, of 726 Mb, which has been corroborated by a high-resolution linkage map. A comparison between the herring genome assembly with other high-quality assemblies from bony fishes revealed few inter-chromosomal but frequent intra-chromosomal rearrangements. The improved assembly facilitates analysis of previously intractable large-scale structural variation, allowing, for example, the detection of a 7.8-Mb inversion on Chromosome 12 underlying ecological adaptation. This supergene shows strong genetic differentiation between populations. The chromosome-based assembly also markedly improves the interpretation of previously detected signals of selection, allowing us to reveal hundreds of independent loci associated with ecological adaptation.