Sequencing and de novo assembly of 150 genomes from Denmark as a population reference.

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 traits. Genetic variation is identified mainly by mapping short reads to the reference genome or by performing local assembly. 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 technology. 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.

circZNF827 nucleates a transcription inhibitory complex to balance neuronal differentiation.

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.

Dynamic Copy Number Evolution of X- and Y-Linked Ampliconic Genes in Human Populations.

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.

Transcriptome analysis of the adult human Klinefelter testis and cellularity-matched controls reveals disturbed differentiation of Sertoli- and Leydig cells.

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.

Distinct differences in global gene expression profiles in non-implanted blastocysts and blastocysts resulting in live birth.

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.

Striated muscle fiber size, composition, and capillary density in diabetes in relation to neuropathy and muscle strength.

OBJECTIVE: Diabetic polyneuropathy (DPN) leads to progressive loss of muscle strength in the lower extremities due to muscular atrophy. Changes in vascularization occur in diabetic striated muscle; however, the relationship between these changes and DPN is as yet unexplored. The aim of the present study was to evaluate histologic properties and capillarization of diabetic skeletal muscle in relation to DPN and muscle strength. METHODS: Twenty type 1 and 20 type 2 diabetic (T1D and T2D, respectively) patients underwent biopsy of the gastrocnemic muscle, isokinetic dynamometry at the ankle, electrophysiological studies, clinical examination, and quantitative sensory examinations. Muscle biopsies were stained immunohistochemically and muscle fiber diameter, fiber type distribution, and capillary density determined. Twenty control subjects were also included in the study. RESULTS: No relationship was found between muscle fiber diameter, muscle fiber type distribution, or capillary density and degree of neuropathy or muscle strength for either patient group. Muscle fiber diameter and the proportion of Type II fibers were greater for T1D patients than both T2D patients and controls. The T2D patients had fewer capillaries per muscle fiber than T1D patients and controls. CONCLUSIONS: Striated muscle fiber size, muscle fiber distribution, and vascularization are unrelated to DPN and muscle strength in diabetes.