1.全基因组分析确定68个与眼压相关的新位点，提高原发性开角型青光眼的风险预测

Genome-wide analyses identify 68 new loci associated with intraocular pressure and improve risk prediction for primary open-angle glaucoma (Nature Genetics)

Abstract

Glaucoma is the leading cause of irreversible blindness globally. Despite its gravity, the disease is frequently undiagnosed in the community. Raised intraocular pressure (IOP) is the most important risk factor for primary open-angle glaucoma (POAG). Here we present a meta-analysis of 139,555 European participants, which identified 112 genomic loci associated with IOP, 68 of which are novel. These loci suggest a strong role for angiopoietin-receptor tyrosine kinase signaling, lipid metabolism, mitochondrial function and developmental processes underlying risk for elevated IOP. In addition, 48 of these loci were nominally associated with glaucoma in an independent cohort, 14 of which were significant at a Bonferroni-corrected threshold. Regression-based glaucoma-prediction models had an area under the receiver operating characteristic curve (AUROC) of 0.76 in US NEIGHBORHOOD study participants and 0.74 in independent glaucoma cases from the UK Biobank. Genetic-prediction models for POAG offer an opportunity to target screening and timely therapy to individuals most at risk.

2. 全基因组meta关联分析强调光诱导信号作为屈光不正的驱动因素

Genome-wide association meta-analysis highlights light-induced signaling as a driver for refractive error(Nature Genetics)

Abstract

Refractive errors, including myopia, are the most frequent eye disorders worldwide and an increasingly common cause of blindness. This genome-wide association meta-analysis in 160,420 participants and replication in 95,505 participants increased the number of established independent signals from 37 to 161 and showed high genetic correlation between Europeans and Asians (>0.78). Expression experiments and comprehensive in silico analyses identified retinal cell physiology and light processing as prominent mechanisms, and also identified functional contributions to refractive-error development in all cell types of the neurosensory retina, retinal pigment epithelium, vascular endothelium and extracellular matrix. Newly identified genes implicate novel mechanisms such as rod-and-cone bipolar synaptic neurotransmission, anterior-segment morphology and angiogenesis. Thirty-one loci resided in or near regions transcribing small RNAs, thus suggesting a role for post-transcriptional regulation. Our results support the notion that refractive errors are caused by a light-dependent retina-to-sclera signaling cascade and delineate potential pathobiological molecular drivers.

3.癌症基因组中亚克隆群体的拷贝数分析和推断工具Sclust

Copy-number analysis and inference of subclonal populations in cancer genomes using Sclust(Nature Protocols)

Abstract

The genomes of cancer cells constantly change during pathogenesis. This evolutionary process can lead to the emergence of drug-resistant mutations in subclonal populations, which can hinder therapeutic intervention in patients. Data derived from massively parallel sequencing can be used to infer these subclonal populations using tumor-specific point mutations. The accurate determination of copy-number changes and tumor impurity is necessary to reliably infer subclonal populations by mutational clustering. This protocol describes how to use Sclust, a copy-number analysis method with a recently developed mutational clustering approach. In a series of simulations and comparisons with alternative methods, we have previously shown that Sclust accurately determines copy-number states and subclonal populations. Performance tests show that the method is computationally efficient, with copy-number analysis and mutational clustering taking <10 min. Sclust is designed such that even non-experts in computational biology or bioinformatics with basic knowledge of the Linux/Unix command-line syntax should be able to carry out analyses of subclonal populations.

4.重建恐龙祖先基因组实现鸟类和非鸟类恐龙的染色体进化追踪

Reconstruction of the diapsid ancestral genome permits chromosome evolution tracing in avian and non-avian dinosaurs(Nature Communications)

Abstract

Genomic organisation of extinct lineages can be inferred from extant chromosome-level genome assemblies. Here, we apply bioinformatic and molecular cytogenetic approaches to determine the genomic structure of the diapsid common ancestor. We then infer the events that likely occurred along this lineage from theropod dinosaurs through to modern birds. Our results suggest that most elements of a typical ‘avian-like’ karyotype (40 chromosome pairs, including 30 microchromosomes) were in place before the divergence of turtles from birds ~255 mya. This genome organisation therefore predates the emergence of early dinosaurs and pterosaurs and the evolution of flight. Remaining largely unchanged interchromosomally through the dinosaur–theropod route that led to modern birds, intrachromosomal changes nonetheless reveal evolutionary breakpoint regions enriched for genes with ontology terms related to chromatin organisation and transcription. This genomic structure therefore appears highly stable yet contributes to a large degree of phenotypic diversity, as well as underpinning adaptive responses to major environmental disruptions via intrachromosomal repatterning.

5.全基因组测序揭示与NSCLC患者炎性微环境相关的基因组特征

Whole-genome sequencing reveals genomic signatures associated with the inflammatory microenvironments in Chinese NSCLC patients(Nature Communications)

Abstract

Chinese lung cancer patients have distinct epidemiologic and genomic features, highlighting the presence of specific etiologic mechanisms other than smoking. Here, we present a comprehensive genomic landscape of 149 non-small cell lung cancer (NSCLC) cases and identify 15 potential driver genes. We reveal that Chinese patients are specially characterized by not only highly clustered EGFR mutations but a mutational signature (MS3, 33.7%), that is associated with inflammatory tumor-infiltrating B lymphocytes (P = 0.001). The EGFR mutation rate is significantly increased with the proportion of the MS3 signature (P = 9.37 × 10−5). TCGA data confirm that the infiltrating B lymphocyte abundance is significantly higher in the EGFR-mutated patients (P = 0.007). Additionally, MS3-high patients carry a higher contribution of distant chromosomal rearrangements >1 Mb (P = 1.35 × 10−7), some of which result in fusions involving genes with important functions (i.e.,ALK and RET). Thus, inflammatory infiltration may contribute to the accumulation of EGFR mutations, especially in never-smokers.

6.3800年生鼠疫菌的分析提示黑死病在青铜时代的起源

Analysis of 3800-year-old Yersinia pestisgenomes suggests Bronze Age origin for bubonic plague(Nature Communications)

Abstract

The origin of Yersinia pestis and the early stages of its evolution are fundamental subjects of investigation given its high virulence and mortality that resulted from past pandemics. Although the earliest evidence of Y. pestis infections in humans has been identified in Late Neolithic/Bronze Age Eurasia (LNBA 5000–3500y BP), these strains lack key genetic components required for flea adaptation, thus making their mode of transmission and disease presentation in humans unclear. Here, we reconstruct ancient Y. pestis genomes from individuals associated with the Late Bronze Age period (~3800 BP) in the Samara region of modern-day Russia. We show clear distinctions between our new strains and the LNBA lineage, and suggest that the full ability for flea-mediated transmission causing bubonic plague evolved more than 1000 years earlier than previously suggested. Finally, we propose that several Y. pestis lineages were established during the Bronze Age, some of which persist to the present day.

7.RNA M6A甲基化参与小鼠小脑后生发育的调控

RNA m6A methylation participates in regulation of postnatal development of the mouse cerebellum(Genome Biology)

Abstract

Background

N6-methyladenosine (m6A) is an important epitranscriptomic mark with high abundance in the brain. Recently, it has been found to be involved in the regulation of memory formation and mammalian cortical neurogenesis. However, while it is now established that m6A methylation occurs in a spatially restricted manner, its functions in specific brain regions still await elucidation.

Results

We identify widespread and dynamic RNA m6A methylation in the developing mouse cerebellum and further uncover distinct features of continuous and temporal-specific m6A methylation across the four postnatal developmental processes. Temporal-specific m6A peaks from P7 to P60 exhibit remarkable changes in their distribution patterns along the mRNA transcripts. We also show spatiotemporal-specific expression of m6A writers METTL3, METTL14, and WTAP and erasers ALKBH5 and FTO in the mouse cerebellum. Ectopic expression of METTL3 mediated by lentivirus infection leads to disorganized structure of both Purkinje and glial cells. In addition, under hypobaric hypoxia exposure, Alkbh5-deletion causes abnormal cell proliferation and differentiation in the cerebellum through disturbing the balance of RNA m6A methylation in different cell fate determination genes. Notably, nuclear export of the hypermethylated RNAs is enhanced in the cerebellum of Alkbh5-deficient mice exposed to hypobaric hypoxia.

Conclusions

Together, our findings provide strong evidence that RNA m6A methylation is controlled in a precise spatiotemporal manner and participates in the regulation of postnatal development of the mouse cerebellum.

8.单细胞RNA测序的UMI计数模型及差异表达分析

UMI-count modeling and differential expression analysis for single-cell RNA sequencing(Genome Biology)

Abstract

Read counting and unique molecular identifier (UMI) counting are the principal gene expression quantification schemes used in single-cell RNA-sequencing (scRNA-seq) analysis. By using multiple scRNA-seq datasets, we reveal distinct distribution differences between these schemes and conclude that the negative binomial model is a good approximation for UMI counts, even in heterogeneous populations. We further propose a novel differential expression analysis algorithm based on a negative binomial model with independent dispersions in each group (NBID). Our results show that this properly controls the FDR and achieves better power for UMI counts when compared to other recently developed packages for scRNA-seq analysis.

9.BGP:用分支高斯过程识别单细胞数据的基因分支动态

BGP: identifying gene-specific branching dynamics from single-cell data with a branching Gaussian process(Genome Biology)

Abstract

High-throughput single-cell gene expression experiments can be used to uncover branching dynamics in cell populations undergoing differentiation through pseudotime methods. We develop the branching Gaussian process (BGP), a non-parametric model that is able to identify branching dynamics for individual genes and provide an estimate of branching times for each gene with an associated credible region. We demonstrate the effectiveness of our method on simulated data, a single-cell RNA-seq haematopoiesis study and mouse embryonic stem cells generated using droplet barcoding. The method is robust to high levels of technical variation and dropout, which are common in single-cell data.

10.单分子实时测序技术在真核基因组中定位和鉴定N6-甲基腺嘌呤

Mapping and characterizing N6-methyladenine in eukaryotic genomes using single-molecule real-time sequencing(Genome Research)

Abstract

N6-Methyladenine (m6dA) has been discovered as a novel form of DNA methylation prevalent in eukaryotes; however, methods for high-resolution mapping of m6dA events are still lacking. Single-molecule real-time (SMRT) sequencing has enabled the detection of m6dA events at single-nucleotide resolution in prokaryotic genomes, but its application to detecting m6dA in eukaryotic genomes has not been rigorously examined. Herein, we identified unique characteristics of eukaryotic m6dA methylomes that fundamentally differ from those of prokaryotes. Based on these differences, we describe the first approach for mapping m6dA events using SMRT sequencing specifically designed for the study of eukaryotic genomes and provide appropriate strategies for designing experiments and carrying out sequencing in future studies. We apply the novel approach to study two eukaryotic genomes. For green algae, we construct the first complete genome-wide map of m6dA at single-nucleotide and single-molecule resolution. For human lymphoblastoid cells (hLCLs), it was necessary to integrate SMRT sequencing data with independent sequencing data. The joint analyses suggest putative m6dA events are enriched in the promoters of young full-length LINE-1 elements (L1s), but call for validation by additional methods. These analyses demonstrate a general method for rigorous mapping and characterization of m6dA events in eukaryotic genomes

11.多倍体植株的基因保留、分化和亚基因组差异

Gene retention, fractionation and subgenome differences in polyploid plants(Nature Plants)

Abstract

All natural plant species are evolved from ancient polyploids. Polyloidization plays an important role in plant genome evolution, species divergence and crop domestication. We review how the pattern of polyploidy within the plant phylogenetic tree has engendered hypotheses involving mass extinctions, lag-times following polyploidy, and epochs of asexuality. Polyploidization has happened repeatedly in plant evolution and, we conclude, is important for crop domestication. Once duplicated, the effect of purifying selection on any one duplicated gene is relaxed, permitting duplicate gene and regulatory element loss (fractionation). We review the general topic of fractionation, and how some gene categories are retained more than others. Several explanations, including neofunctionalization, subfunctionalization and gene product dosage balance, have been shown to influence gene content over time. For allopolyploids, genetic differences between parental lines immediately manifest as subgenome dominance in the wide-hybrid, and persist and propagate for tens of millions of years. While epigenetic modifications are certainly involved in genome dominance, it has been difficult to determine which came first, the chromatin marks being measured or gene expression. Data support the conclusion that genome dominance and heterosis are antagonistic and mechanically entangled; both happen immediately in the synthetic wide-cross hybrid. Also operating in this hybrid are mechanisms of ‘paralogue interference’. We present a foundation model to explain gene expression and vigour in a wide hybrid/new allotetraploid. This Review concludes that some mechanisms operate immediately at the wide-hybrid, and other mechanisms begin their operations later. Direct interaction of new paralogous genes, as measured using high-resolution chromatin conformation capture, should inform future research and single cell transcriptome sequencing should help achieve specificity while studying gene sub- and neo-functionalization.

12.WEGO 2：一个用于分析和绘制GO注释的Web工具，2018更新

WEGO 2.0: a web tool for analyzing and plotting GO annotations, 2018 update(Nucleic Acids Research)

Abstract

WEGO (Web Gene Ontology Annotation Plot), created in 2006, is a simple but useful tool for visualizing, comparing and plotting GO (Gene Ontology) annotation results. Owing largely to the rapid development of high-throughput sequencing and the increasing acceptance of GO, WEGO has benefitted from outstanding performance regarding the number of users and citations in recent years, which motivated us to update to version 2.0. WEGO uses the GO annotation results as input. Based on GO’s standardized DAG (Directed Acyclic Graph) structured vocabulary system, the number of genes corresponding to each GO ID is calculated and shown in a graphical format. WEGO 2.0 updates have targeted four aspects, aiming to provide a more efficient and up-to-date approach for comparative genomic analyses. First, the number of input files, previously limited to three, is now unlimited, allowing WEGO to analyze multiple datasets. Also added in this version are the reference datasets of nine model species that can be adopted as baselines in genomic comparative analyses. Furthermore, in the analyzing processes each Chi-square test is carried out for multiple datasets instead of every two samples. At last, WEGO 2.0 provides an additional output graph along with the traditional WEGO histogram, displaying the sorted P-values of GO terms and indicating their significant differences. At the same time, WEGO 2.0 features an entirely new user interface. WEGO is available for free at http://wego.genomics.org.cn.

生信人交流群WeGAP－欢迎加入

欢迎关注生信人

一个物种一个家

1对1 VIP培训服务

TCGA | 小工具 | 数据库 |组装| 注释 |   基因家族  |  Pvalue

基因预测  |bestorf |  sci | NAR | 在线工具 | 生存分析 | 热图

 生信不死 | 初学者 | circRNA | 一箭画心| 十二生肖 | circos

 舞台|基因组 | 黄金测序 | 套路 | 杂谈组装 |  进化 | 测序简史