FLED: a full-length eccDNA detector for long-reads sequencing data.

Reconstructing the full-length sequence of extrachromosomal circular DNA (eccDNA) from short sequencing reads has proved challenging given the similarity of eccDNAs and their corresponding linear DNAs. Previous sequencing methods were unable to achieve high-throughput detection of full-length eccDNAs. Herein, a novel algorithm was developed, called Full-Length eccDNA Detection (FLED), to reconstruct the sequence of eccDNAs based on the strategy that combined rolling circle amplification and nanopore long-reads sequencing technology. Seven human epithelial and cancer cell line samples were analyzed by FLED and over 5000 full-length eccDNAs were identified per sample. The structures of identified eccDNAs were validated by both Polymerase Chain Reaction (PCR) and Sanger sequencing. Compared to other published nanopore-based eccDNA detectors, FLED exhibited higher sensitivity. In cancer cell lines, the genes overlapped with eccDNA regions were enriched in cancer-related pathways and cis-regulatory elements can be predicted in the upstream or downstream of intact genes on eccDNA molecules, and the expressions of these cancer-related genes were dysregulated in tumor cell lines, indicating the regulatory potency of eccDNAs in biological processes. The proposed method takes advantage of nanopore long reads and enables unbiased reconstruction of full-length eccDNA sequences. FLED is implemented using Python3 which is freely available on GitHub (https://github.com/FuyuLi/FLED).

Circular RNAs in the human brain are tailored to neuron identity and neuropsychiatric disease.

Little is known about circular RNAs (circRNAs) in specific brain cells and human neuropsychiatric disease. Here, we systematically identify over 11,039 circRNAs expressed in vulnerable dopamine and pyramidal neurons laser-captured from 190 human brains and non-neuronal cells using ultra-deep, total RNA sequencing. 1526 and 3308 circRNAs are custom-tailored to the cell identity of dopamine and pyramidal neurons and enriched in synapse pathways. 29% of Parkinson's and 12% of Alzheimer's disease-associated genes produced validated circRNAs. circDNAJC6, which is transcribed from a juvenile-onset Parkinson's gene, is already dysregulated during prodromal, onset stages of common Parkinson's disease neuropathology. Globally, addiction-associated genes preferentially produce circRNAs in dopamine neurons, autism-associated genes in pyramidal neurons, and cancers in non-neuronal cells. This study shows that circular RNAs in the human brain are tailored to neuron identity and implicate circRNA-regulated synaptic specialization in neuropsychiatric diseases.

Circular RNAs in the human brain are tailored to neuron identity and neuropsychiatric disease.

Little is known about circular RNAs (circRNAs) in specific brain cells and human neuropsychiatric disease. Here, we systematically identified over 11,039 circRNAs expressed in vulnerable dopamine and pyramidal neurons laser-captured from 190 human brains and non-neuronal cells using ultra-deep, total RNA sequencing. 1,526 and 3,308 circRNAs were custom-tailored to the cell identity of dopamine and pyramidal neurons and enriched in synapse pathways. 88% of Parkinson's and 80% of Alzheimer's disease-associated genes produced circRNAs. circDNAJC6, produced from a juvenile-onset Parkinson's gene, was already dysregulated during prodromal, onset stages of common Parkinson's disease neuropathology. Globally, addiction-associated genes preferentially produced circRNAs in dopamine neurons, autism-associated genes in pyramidal neurons, and cancers in non-neuronal cells. This study shows that circular RNAs in the human brain are tailored to neuron identity and implicate circRNA- regulated synaptic specialization in neuropsychiatric diseases.

Genome-wide association meta-analysis identifies 29 new acne susceptibility loci.

Acne vulgaris is a highly heritable skin disorder that primarily impacts facial skin. Severely inflamed lesions may leave permanent scars that have been associated with long-term psychosocial consequences. Here, we perform a GWAS meta-analysis comprising 20,165 individuals with acne from nine independent European ancestry cohorts. We identify 29 novel genome-wide significant loci and replicate 14 of the 17 previously identified risk loci, bringing the total number of reported acne risk loci to 46. Using fine-mapping and eQTL colocalisation approaches, we identify putative causal genes at several acne susceptibility loci that have previously been implicated in Mendelian hair and skin disorders, including pustular psoriasis. We identify shared genetic aetiology between acne, hormone levels, hormone-sensitive cancers and psychiatric traits. Finally, we show that a polygenic risk score calculated from our results explains up to 5.6% of the variance in acne liability in an independent cohort.

The genetic source tracking of human urinary exosomes.

The genetic origins of nanoscale extracellular vesicles in our body fluids remains unclear. Here, we perform a tracking analysis of urinary exosomes via RNA sequencing, revealing that urine exosomes mostly express tissue-specific genes for the bladder and have close cell-genetic relationships to the endothelial cell, basal cell, monocyte, and dendritic cell. Tracking the differentially expressed genes of cancers and corresponding enrichment analysis show urine exosomes are intensively involved in immune activities, indicating that they may be harnessed as reliable biomarkers of noninvasive liquid biopsy in cancer genomic diagnostics and precision medicine.

Understanding the effect of smoking and drinking behavior on Parkinson's disease risk: a Mendelian randomization study.

Previous observational studies have identified correlations between Parkinson's disease (PD) risk and lifestyle factors. However, whether or not those associations are causal remains unclear. To infer causality between PD risk and smoking or alcohol intake, we conducted a two-sample Mendelian randomization study using genome-wide association study summary statistics from the GWAS & Sequencing Consortium of Alcohol and Nicotine use study (1.2 million participants) and the latest meta-analysis from the International Parkinson's Disease Genomics Consortium (37,688 PD cases and 18,618 proxy-cases). We performed sensitivity analyses, including testing for pleiotropy with MR-Egger and MR-PRESSO, and multivariable MR modeling to account for the genetic effects of competing substance use traits on PD risk. Our results revealed causal associations of alcohol intake (OR 0.79; 95% CI 0.65-0.96; p = 0.021) and smoking continuation (which compares current vs. former smokers) (OR 0.64; 95% CI 0.46-0.89; p = 0.008) with lower PD risk. Multivariable MR analyses showed that the causal association between drinks per week and PD is unlikely due to confounding by smoking behavior. Finally, frailty analyses suggested that the causal effects of both alcohol intake and smoking continuation on PD risk estimated from MR analysis are not explained by the presence of survival bias alone. Our findings support the role of smoking as a protective factor against PD, but only when comparing current vs. former smokers. Similarly, increased alcohol intake had a protective effect over PD risk, with the alcohol dehydrogenase 1B (ADH1B) locus as a potential candidate for further investigation of the mechanisms underlying this association.

High-resolution analysis of differential gene expression during skeletal muscle atrophy and programmed cell death.

Skeletal muscles can undergo atrophy and/or programmed cell death (PCD) during development or in response to a wide range of insults, including immobility, cachexia, and spinal cord injury. However, the protracted nature of atrophy and the presence of multiple cell types within the tissue complicate molecular analyses. One model that does not suffer from these limitations is the intersegmental muscle (ISM) of the tobacco hawkmoth Manduca sexta. Three days before the adult eclosion (emergence) at the end of metamorphosis, the ISMs initiate a nonpathological program of atrophy that results in a 40% loss of mass. The ISMs then generate the eclosion behavior and initiate a nonapoptotic PCD during the next 30 h. We have performed a comprehensive transcriptomics analysis of all mRNAs and microRNAs throughout ISM development to better understand the molecular mechanisms that mediate atrophy and death. Atrophy involves enhanced protein catabolism and reduced expression of the genes involved in respiration, adhesion, and the contractile apparatus. In contrast, PCD involves the induction of numerous proteases, DNA methylases, membrane transporters, ribosomes, and anaerobic metabolism. These changes in gene expression are largely repressed when insects are injected with the insect steroid hormone 20-hydroxyecdysone, which delays death. The expression of the death-associated proteins may be greatly enhanced by reductions in specific microRNAs that function to repress translation. This study not only provides fundamental new insights into basic developmental processes, it may also represent a powerful resource for identifying potential diagnostic markers and molecular targets for therapeutic intervention.

ß2-Adrenoreceptor is a regulator of the α-synuclein gene driving risk of Parkinson's disease.

Copy number mutations implicate excess production of α-synuclein as a possibly causative factor in Parkinson's disease (PD). Using an unbiased screen targeting endogenous gene expression, we discovered that the ß2-adrenoreceptor (ß2AR) is a regulator of the α-synuclein gene (SNCA). ß2AR ligands modulate SNCA transcription through histone 3 lysine 27 acetylation of its promoter and enhancers. Over 11 years of follow-up in 4 million Norwegians, the ß2AR agonist salbutamol, a brain-penetrant asthma medication, was associated with reduced risk of developing PD (rate ratio, 0.66; 95% confidence interval, 0.58 to 0.76). Conversely, a ß2AR antagonist correlated with increased risk. ß2AR activation protected model mice and patient-derived cells. Thus, ß2AR is linked to transcription of α-synuclein and risk of PD in a ligand-specific fashion and constitutes a potential target for therapies.

The Role of H3K4me3 in Transcriptional Regulation Is Altered in Huntington's Disease.

Huntington's disease (HD) is an autosomal-dominant neurodegenerative disorder resulting from expansion of CAG repeats in the Huntingtin (HTT) gene. Previous studies have shown mutant HTT can alter expression of genes associated with dysregulated epigenetic modifications. One of the most widely studied chromatin modifications is trimethylated lysine 4 of histone 3 (H3K4me3). Here, we conducted the first comprehensive study of H3K4me3 ChIP-sequencing in neuronal chromatin from the prefrontal cortex of six HD cases and six non-neurologic controls, and its association with gene expression measured by RNA-sequencing. We detected 2,830 differentially enriched H3K4me3 peaks between HD and controls, with 55% of them down-regulated in HD. Although H3K4me3 signals are expected to be associated with mRNA levels, we found an unexpected discordance between altered H3K4me3 peaks and mRNA levels. Gene ontology (GO) term enrichment analysis of the genes with differential H3K4me3 peaks, revealed statistically significantly enriched GO terms only in the genes with down-regulated signals in HD. The most frequently implicated biological process terms are organ morphogenesis and positive regulation of gene expression. More than 9,000 H3K4me3 peaks were located not near any recognized transcription start sites and approximately 36% of these "distal" peaks co-localized to known enhancer sites. Six transcription factors and chromatin remodelers are differentially enriched in HD H3K4me3 distal peaks, including EZH2 and SUZ12, two core subunits of the polycomb repressive complex 2 (PRC2). Moreover, PRC2 repressive state was significantly depleted in HD-enriched peaks, suggesting the epigenetic role of PRC2 inhibition associated with up-regulated H3K4me3 in Huntington's disease. In summary, our study provides new insights into transcriptional dysregulation of Huntington's disease by analyzing the differentiation of H3K4me3 enrichment.

Epigenetic dysregulation of hairy and enhancer of split 4 (HES4) is associated with striatal degeneration in postmortem Huntington brains.

To investigate epigenetic contributions to Huntington's disease (HD) pathogenesis, we carried out genome-wide mapping of the transcriptional mark, trimethyl-histone H3-lysine 4 (H3K4me3) in neuronal nuclei extracted from prefrontal cortex of HD cases and controls using chromatin immunoprecipitation followed by deep-sequencing. Neuron-specific mapping of the genome-wide distribution of H3K4me3 revealed 136 differentially enriched loci associated with genes implicated in neuronal development and neurodegeneration, including GPR3, TMEM106B, PDIA6 and the Notch signaling genes hairy and enhancer of split 4 (HES4) and JAGGED2, supporting the view that the neuronal epigenome is affected in HD. Importantly, loss of H3K4me3 at CpG-rich sequences on the HES4 promoter was associated with excessive DNA methylation, reduced binding of nuclear proteins to the methylated region and altered expression of HES4 and HES4 targeted genes MASH1 and P21 involved in striatal development. Moreover, hypermethylation of HES4 promoter sequences was strikingly correlated with measures of striatal degeneration and age-of-onset in a cohort of 25 HD brains (r = 0.56, P = 0.006). Lastly, shRNA knockdown of HES4 in human neuroblastoma cells altered MASH1 and P21 mRNA expression and markedly increased mutated HTT-induced aggregates and cell death. These findings, taken together, suggest that epigenetic dysregulation of HES4 could play a critical role in modifying HD disease pathogenesis and severity.

MicroRNAs located in the Hox gene clusters are implicated in huntington's disease pathogenesis.

Transcriptional dysregulation has long been recognized as central to the pathogenesis of Huntington's disease (HD). MicroRNAs (miRNAs) represent a major system of post-transcriptional regulation, by either preventing translational initiation or by targeting transcripts for storage or for degradation. Using next-generation miRNA sequencing in prefrontal cortex (Brodmann Area 9) of twelve HD and nine controls, we identified five miRNAs (miR-10b-5p, miR-196a-5p, miR-196b-5p, miR-615-3p and miR-1247-5p) up-regulated in HD at genome-wide significance (FDR q-value<0.05). Three of these, miR-196a-5p, miR-196b-5p and miR-615-3p, were expressed at near zero levels in control brains. Expression was verified for all five miRNAs using reverse transcription quantitative PCR and all but miR-1247-5p were replicated in an independent sample (8HD/8C). Ectopic miR-10b-5p expression in PC12 HTT-Q73 cells increased survival by MTT assay and cell viability staining suggesting increased expression may be a protective response. All of the miRNAs but miR-1247-5p are located in intergenic regions of Hox clusters. Total mRNA sequencing in the same samples identified fifteen of 55 genes within the Hox cluster gene regions as differentially expressed in HD, and the Hox genes immediately adjacent to the four Hox cluster miRNAs as up-regulated. Pathway analysis of mRNA targets of these miRNAs implicated functions for neuronal differentiation, neurite outgrowth, cell death and survival. In regression models among the HD brains, huntingtin CAG repeat size, onset age and age at death were independently found to be inversely related to miR-10b-5p levels. CAG repeat size and onset age were independently inversely related to miR-196a-5p, onset age was inversely related to miR-196b-5p and age at death was inversely related to miR-615-3p expression. These results suggest these Hox-related miRNAs may be involved in neuroprotective response in HD. Recently, miRNAs have shown promise as biomarkers for human diseases and given their relationship to disease expression, these miRNAs are biomarker candidates in HD.

An ancient transcription factor initiates the burst of piRNA production during early meiosis in mouse testes.

Animal germ cells produce PIWI-interacting RNAs (piRNAs), small silencing RNAs that suppress transposons and enable gamete maturation. Mammalian transposon-silencing piRNAs accumulate early in spermatogenesis, whereas pachytene piRNAs are produced later during postnatal spermatogenesis and account for >95% of all piRNAs in the adult mouse testis. Mutants defective for pachytene piRNA pathway proteins fail to produce mature sperm, but neither the piRNA precursor transcripts nor the trigger for pachytene piRNA production is known. Here, we show that the transcription factor A-MYB initiates pachytene piRNA production. A-MYB drives transcription of both pachytene piRNA precursor RNAs and the mRNAs for core piRNA biogenesis factors including MIWI, the protein through which pachytene piRNAs function. A-MYB regulation of piRNA pathway proteins and piRNA genes creates a coherent feedforward loop that ensures the robust accumulation of pachytene piRNAs. This regulatory circuit, which can be detected in rooster testes, likely predates the divergence of birds and mammals.

Mbd3/NURD complex regulates expression of 5-hydroxymethylcytosine marked genes in embryonic stem cells.

Numerous chromatin regulators are required for embryonic stem (ES) cell self-renewal and pluripotency, but few have been studied in detail. Here, we examine the roles of several chromatin regulators whose loss affects the pluripotent state of ES cells. We find that Mbd3 and Brg1 antagonistically regulate a common set of genes by regulating promoter nucleosome occupancy. Furthermore, both Mbd3 and Brg1 play key roles in the biology of 5-hydroxymethylcytosine (5hmC): Mbd3 colocalizes with Tet1 and 5hmC in vivo, Mbd3 knockdown preferentially affects expression of 5hmC-marked genes, Mbd3 localization is Tet1-dependent, and Mbd3 preferentially binds to 5hmC relative to 5-methylcytosine in vitro. Finally, both Mbd3 and Brg1 are themselves required for normal levels of 5hmC in vivo. Together, our results identify an effector for 5hmC, and reveal that control of gene expression by antagonistic chromatin regulators is a surprisingly common regulatory strategy in ES cells.

Transcriptional features of genomic regulatory blocks.

BACKGROUND: Genomic regulatory blocks (GRBs) are chromosomal regions spanned by highly conserved non-coding elements (HCNEs), most of which serve as regulatory inputs of one target gene in the region. The target genes are most often transcription factors involved in embryonic development and differentiation. GRBs often contain extensive gene deserts, as well as additional 'bystander' genes intertwined with HCNEs but whose expression and function are unrelated to those of the target gene. The tight regulation of target genes, complex arrangement of regulatory inputs, and the differential responsiveness of genes in the region call for the examination of fundamental rules governing transcriptional activity in GRBs. Here we use extensive CAGE tag mapping of transcription start sites across different human tissues and differentiation stages combined with expression data and a number of sequence and epigenetic features to discover these rules and patterns. RESULTS: We show evidence that GRB target genes have properties that set them apart from their bystanders as well as other genes in the genome: longer CpG islands, a higher number and wider spacing of alternative transcription start sites, and a distinct composition of transcription factor binding sites in their core/proximal promoters. Target gene expression correlates with the acetylation state of HCNEs in the region. Additionally, target gene promoters have a distinct combination of activating and repressing histone modifications in mouse embryonic stem cell lines. CONCLUSIONS: GRB targets are genes with a number of unique features that are the likely cause of their ability to respond to regulatory inputs from very long distances.