A major part of the human Y chromosome consists of palindromes with multiple copies of genes primarily expressed in testis, many of which have been claimed to affect male fertility. Here we examine copy number variation in these palindromes based on whole genome sequence data from 11,527 Icelandic men. Using a subset of 7947 men grouped into 1449 patrilineal genealogies, we infer 57 large scale de novo copy number mutations affecting palindrome 1. This corresponds to a mutation rate of 2.34 × 10 −3 mutations per meiosis, which is 4.1 times larger than our phylogenetic estimate of the mutation rate (5.72 × 10 −4 ), suggesting that de novo mutations on the Y are lost faster than expected under neutral evolution. Although simulations indicate a selection coefficient of 1.8% against non-reference copy number carriers, we do not observe differences in fertility among sequenced men associated with their copy number genotype, but we lack statistical power to detect differences resulting from weak negative selection. We also perform association testing of a diverse set of 341 traits to palindromic copy number without any significant associations. We conclude that large-scale palindrome copy number variation on the Y chromosome has little impact on human phenotype diversity.

Circular RNAs are important for many cellular processes but their mechanisms of action remain poorly understood. Here, we map circRNA inventories of mouse embryonic stem cells, neuronal progenitor cells and differentiated neurons and identify hundreds of highly expressed circRNAs. By screening several candidate circRNAs for a potential function in neuronal differentiation, we find that circZNF827 represses expression of key neuronal markers, suggesting that this molecule negatively regulates neuronal differentiation. Among 760 tested genes linked to known neuronal pathways, knockdown of circZNF827 deregulates expression of numerous genes including nerve growth factor receptor (NGFR), which becomes transcriptionally upregulated to enhance NGF signaling. We identify a circZNF827-nucleated transcription-repressive complex containing hnRNP-K/L proteins and show that knockdown of these factors strongly augments NGFR regulation. Finally, we show that the ZNF827 protein is part of the mRNP complex, suggesting a functional co-evolution of a circRNA and the protein encoded by its linear pre-mRNA host.

Incoming viruses challenge the cell with diverse foreign molecules, which need to be sensed quickly to initiate immune responses and to remove the viral components. In this study, we investigate the cellular requirements for sensing and degradation of incoming viral DNA and capsids during herpes simplex virus type 1 (HSV-1) infections. Using click chemistry labeling of the viral genome, we found that HSV-1 DNA was released from a subset of capsids into the cytosol early in infection. By next-generation sequencing of cyclic GMP-AMP (cGAMP) synthase (cGAS)-bound DNA from HSV-1-infected cells, we show that HSV-1 DNA was bound by the cytosolic DNA sensor cGAS. Activation of cGAS enzymatic activity by viral DNA did not require proteasomal activity, indicating that viral DNA release into the cytosol is not proteasome-dependent. However, induction of interferon (IFN)-β expression was blocked by inhibition of the proteasome, suggesting a contribution of the proteasome to IFN-β induction through the cGAS-stimulator of interferon genes pathway. Viral DNA was cleared from the cytosol within few hours, in a manner dependent on TREX1 and a cGAS-dependent process. Capsid material in the cytoplasm was also degraded rapidly. This was partially blocked by treatment with a proteasome inhibitor. This treatment led to accumulation of DNA-containing viral capsids near the nucleus and reduced nuclear entry of viral DNA. Thus, cells infected with HSV-1 use a panel of mechanisms to eliminate viral DNA and capsids. This represents a barrier for establishment of infection and potentially enables the host to gear the IFN-β response to a level required for antiviral defense without causing immunopathology.

The most common human sex chromosomal disorder is Klinefelter syndrome (KS; 47,XXY). Adult patients with KS display a diverse phenotype but are nearly always infertile, due to testicular degeneration at puberty. To identify mechanisms causing the selective destruction of the seminiferous epithelium, we performed RNA-sequencing of 24 fixed paraffin-embedded testicular tissue samples. Analysis of informative transcriptomes revealed 235 differentially expressed transcripts (DETs) in the adult KS testis showing enrichment of long non-coding RNAs, but surprisingly not of X-chromosomal transcripts. Comparison to 46,XY samples with complete spermatogenesis and Sertoli cell-only-syndrome allowed prediction of the cellular origin of 71 of the DETs. DACH2 and FAM9A were validated by immunohistochemistry and found to mark apparently undifferentiated somatic cell populations in the KS testes. Moreover, transcriptomes from fetal, pre-pubertal, and adult KS testes showed a limited overlap, indicating that different mechanisms are likely to operate at each developmental stage. Based on our data, we propose that testicular degeneration in men with KS is a consequence of germ cells loss initiated during early development in combination with disturbed maturation of Sertoli- and Leydig cells.

Ampliconic genes are multicopy genes often located on sex chromosomes and enriched for testis-expressed genes. Here, Lucotte et al. developed new bioinformatic approaches to investigate the ampliconic gene copy number and their coding... Ampliconic genes are multicopy, with the majority found on sex chromosomes and enriched for testis-expressed genes. While ampliconic genes have been associated with the emergence of hybrid incompatibilities, we know little about their copy number distribution and their turnover in human populations. Here, we explore the evolution of human X- and Y-linked ampliconic genes by investigating copy number variation (CNV) and coding variation between populations using the Simons Genome Diversity Project. We develop a method to assess CNVs using the read depth on modified X and Y chromosome targets containing only one repetition of each ampliconic gene. Our results reveal extensive standing variation in copy number both within and between human populations for several ampliconic genes. For the Y chromosome, we can infer multiple independent amplifications and losses of these gene copies even within closely related Y haplogroups, that diversified < 50,000 years ago. Moreover, X- and Y-linked ampliconic genes seem to have a faster amplification dynamic than autosomal multicopy genes. Looking at expression data from another study, we also find that X- and Y-linked ampliconic genes with extensive CNV are significantly more expressed than genes with no CNV during meiotic sex chromosome inactivation (for both X and Y) and postmeiotic sex chromosome repression (for the Y chromosome only). While we cannot rule out that the XY-linked ampliconic genes are evolving neutrally, this study gives insights into the distribution of copy number within human populations and demonstrates an extremely fast turnover in copy number of these regions.

In humans, the most common sex chromosomal disorder is Klinefelter syndrome (KS), caused by the presence of one or more extra X-chromosomes. KS patients display a varying adult phenotype but usually present with azoospermia due to testicular degeneration, which accelerates at puberty. The timing of the germ cell loss and whether it is caused by dysgenetic fetal development of the testes is not known. We investigated eight fetal KS testes and found a marked reduction in MAGE-A4-positive pre-spermatogonia compared with testes from 15 age-matched controls, indicating a failure of the gonocytes to differentiate into pre-spermatogonia. Transcriptome analysis by RNA-sequencing of formalin-fixed, paraffin-embedded testes originating from four fetal KS and five age-matched controls revealed 211 differentially expressed transcripts in the fetal KS testis. We found a significant enrichment of upregulated X-chromosomal transcripts and validated the expression of the pseudoautosomal region 1 (PAR1) gene, AKAP17A. Moreover, we found enrichment of long non-coding RNAs in the KS testes (e.g. LINC01569 and RP11-485F13.1). In conclusion, our data indicate that the testicular phenotype observed among adult men with KS is initiated already in fetal life by failure of the gonocyte differentiation into pre-spermatogonia, which could be due to aberrant expression of long non-coding RNAs.

A report of high-depth, short-read sequencing and de novo assemblies for 150 individuals from 50 parent–offspring trios as part of establishing a population reference genome for the GenomeDenmark project. Genome sequencing of individuals across a population is an important component of precision medicine initiatives, and can be used to characterize genetic variation and in association mapping for diseases and complex traits. Mikkel Schierup and colleagues report efforts to establish a population reference genome for the Danish population as part of the GenomeDenmark project. The authors report high-depth short-read sequencing and de novo assemblies for 150 individuals from 50 parent–offspring trios. They demonstrate that this approach provides similar quality metrics to long-read approaches and helps to resolve structural variation and complex genomic regions. This provides a cost-effective way to establish a population reference genome that will be useful for association mapping and precision medicine initiatives. Hundreds of thousands of human genomes are now being sequenced to characterize genetic variation and use this information to augment association mapping studies of complex disorders and other phenotypic traits1,2,3,4. Genetic variation is identified mainly by mapping short reads to the reference genome or by performing local assembly2,5,6,7. However, these approaches are biased against discovery of structural variants and variation in the more complex parts of the genome. Hence, large-scale de novo assembly is needed. Here we show that it is possible to construct excellent de novo assemblies from high-coverage sequencing with mate-pair libraries extending up to 20 kilobases. We report de novo assemblies of 150 individuals (50 trios) from the GenomeDenmark project. The quality of these assemblies is similar to those obtained using the more expensive long-read technology4,8,9,10,11,12,13. We use the assemblies to identify a rich set of structural variants including many novel insertions and demonstrate how this variant catalogue enables further deciphering of known association mapping signals. We leverage the assemblies to provide 100 completely resolved major histocompatibility complex haplotypes and to resolve major parts of the Y chromosome. Our study provides a regional reference genome that we expect will improve the power of future association mapping studies and hence pave the way for precision medicine initiatives, which now are being launched in many countries including Denmark.

Genes in the major histocompatibility complex (MHC, also known as HLA) play a critical role in the immune response and variation within the extended 4-Mb region shows association with major risks of many diseases. Yet, deciphering the underlying causes of these associations is difficult because the MHC is the most polymorphic region of the genome with a complex linkage disequilibrium structure. Here, we reconstruct full MHC haplotypes from de novo assembled trios without relying on a reference genome and perform evolutionary analyses. We report 100 full MHC haplotypes and call a large set of structural variants in the regions for future use in imputation with GWAS data. We also present the first complete analysis of the recombination landscape in the entire region and show how balancing selection at classical genes have linked effects on the frequency of variants throughout the region.

Results from animal models points towards the existence of a gene expression profile that is distinguishably different in viable embryos compared with non-viable embryos. Knowledge of human embryo transcripts is however limited, in particular with regard to how gene expression is related to clinical outcome. The purpose of the present study was therefore to determine the global gene expression profiles of human blastocysts. Next Generation Sequencing was used to identify genes that were differentially expressed in non-implanted embryos and embryos resulting in live birth. Three trophectoderm biopsies were obtained from morphologically high quality blastocysts resulting in live birth and three biopsies were obtained from non-implanting blastocysts of a comparable morphology. Total RNA was extracted from all samples followed by complete transcriptome sequencing. Using a set of filtering criteria, we obtained a list of 181 genes that were differentially expressed between trophectoderm biopsies from embryos resulting in either live birth or no implantation (negative hCG), respectively. We found that 37 of the 181 genes displayed significantly differential expression (p<0.05), e.g. EFNB1, CYTL1 and TEX26 and TESK1, MSL1 and EVI5 in trophectoderm biopsies associated with live birth and non-implanting, respectively. Out of the 181 genes, almost 80% (145 genes) were up-regulated in biopsies from un-implanted embryos, whereas only 20% (36 genes) showed an up-regulation in the samples from embryos resulting in live birth. Our findings suggest the presence of molecular differences visually undetectable between implanted and non-implanted embryos, and represent a proof of principle study.

The mammalian target of rapamycin complex 1 (mTORC1) is a critical regulator of protein synthesis. The best studied targets of mTORC1 in translation are the eukaryotic initiation factor-binding protein 1 (4E-BP1) and ribosomal protein S6 kinase 1 (S6K1). In this study, we identify the La-related protein 1 (LARP1) as a key novel target of mTORC1 with a fundamental role in terminal oligopyrimidine (TOP) mRNA translation. Recent genome-wide studies indicate that TOP and TOP-like mRNAs compose a large portion of the mTORC1 translatome, but the mechanism by which mTORC1 controls TOP mRNA translation is incompletely understood. Here, we report that LARP1 functions as a key repressor of TOP mRNA translation downstream of mTORC1. Our data show the following: (i) LARP1 associates with mTORC1 via RAPTOR; (ii) LARP1 interacts with TOP mRNAs in an mTORC1-dependent manner; (iii) LARP1 binds the 5'TOP motif to repress TOP mRNA translation; and (iv) LARP1 competes with the eukaryotic initiation factor (eIF) 4G for TOP mRNA binding. Importantly, from a drug resistance standpoint, our data also show that reducing LARP1 protein levels by RNA interference attenuates the inhibitory effect of rapamycin, Torin1, and amino acid deprivation on TOP mRNA translation. Collectively, our findings demonstrate that LARP1 functions as an important repressor of TOP mRNA translation downstream of mTORC1.