The UCSC Genome Browser database: 2017 update

From Gene to Whole Cell: Modeling, Visualization, and Analysis

Proteogenomics combines proteomic and genetic data to gain new insights in molecular mechanisms. Here, we extend this approach toward structural biology from a tool perspective. The chapter starts with tools which can be used to explore genetic information and then enrich those with proteomic data. Based on the corresponding identifiers, three-dimensional structures of proteins are identified and used to embed them in their molecular environment, here the surrounding membrane. This membrane is then mapped onto the surface of an interpretative three-dimensional cell model. Then, the embedded protein and the cell environment are associated with a metabolic pathway, again based on the identifiers provided by biomedical databases. Accompanying the different chapters, related work is discussed which can alternatively be used. Finally, an outlook toward immersive analytics is given.

Chromosome-level genome assembly of the morabine grasshopper Vandiemenella viatica19

Morabine grasshoppers in the Vandiemenella viatica species group, which show karyotype diversity, have been studied for their ecological distribution and speciation in relation to their genetic and chromosomal diversity. They are good models for studying sex chromosome evolution as "old" and newly emerged sex chromosomes co-exist within the group. Here we present a reference genome for the viatica19 chromosomal race, that possesses the ancestral karyotype within the group. Using PacBio HiFi and Hi-C sequencing, we generated a chromosome-level assembly of 4.09 Gb in span, scaffold N50 of 429 Mb, and complete BUSCO score of 98.1%, containing 10 pseudo-chromosomes. We provide Illumina datasets of males and females, used to identify the X chromosome. The assembly contains 19,034 predicted protein-coding genes, and a total of 75.21% of repetitive DNA sequences. By leveraging HiFi reads, we mapped the genome-wide distribution of methylated bases (5mC and 6 mA). This comprehensive assembly offers a robust reference for morabine grasshoppers and supports further research into speciation and sex chromosome diversification within the group and its related species.

Single-cell analysis of chromatin accessibility in the adult mouse brain

Recent advances in single-cell technologies have led to the discovery of thousands of brain cell types; however, our understanding of the gene regulatory programs in these cell types is far from complete1-4. Here we report a comprehensive atlas of candidate cis-regulatory DNA elements (cCREs) in the adult mouse brain, generated by analysing chromatin accessibility in 2.3 million individual brain cells from 117 anatomical dissections. The atlas includes approximately 1 million cCREs and their chromatin accessibility across 1,482 distinct brain cell populations, adding over 446,000 cCREs to the most recent such annotation in the mouse genome. The mouse brain cCREs are moderately conserved in the human brain. The mouse-specific cCREs-specifically, those identified from a subset of cortical excitatory neurons-are strongly enriched for transposable elements, suggesting a potential role for transposable elements in the emergence of new regulatory programs and neuronal diversity. Finally, we infer the gene regulatory networks in over 260 subclasses of mouse brain cells and develop deep-learning models to predict the activities of gene regulatory elements in different brain cell types from the DNA sequence alone. Our results provide a resource for the analysis of cell-type-specific gene regulation programs in both mouse and human brains.

Evolution of neuronal cell classes and types in the vertebrate retina

The basic plan of the retina is conserved across vertebrates, yet species differ profoundly in their visual needs1. Retinal cell types may have evolved to accommodate these varied needs, but this has not been systematically studied. Here we generated and integrated single-cell transcriptomic atlases of the retina from 17 species: humans, two non-human primates, four rodents, three ungulates, opossum, ferret, tree shrew, a bird, a reptile, a teleost fish and a lamprey. We found high molecular conservation of the six retinal cell classes (photoreceptors, horizontal cells, bipolar cells, amacrine cells, retinal ganglion cells (RGCs) and Müller glia), with transcriptomic variation across species related to evolutionary distance. Major subclasses were also conserved, whereas variation among cell types within classes or subclasses was more pronounced. However, an integrative analysis revealed that numerous cell types are shared across species, based on conserved gene expression programmes that are likely to trace back to an early ancestral vertebrate. The degree of variation among cell types increased from the outer retina (photoreceptors) to the inner retina (RGCs), suggesting that evolution acts preferentially to shape the retinal output. Finally, we identified rodent orthologues of midget RGCs, which comprise more than 80% of RGCs in the human retina, subserve high-acuity vision, and were previously believed to be restricted to primates2. By contrast, the mouse orthologues have large receptive fields and comprise around 2% of mouse RGCs. Projections of both primate and mouse orthologous types are overrepresented in the thalamus, which supplies the primary visual cortex. We suggest that midget RGCs are not primate innovations, but are descendants of evolutionarily ancient types that decreased in size and increased in number as primates evolved, thereby facilitating high visual acuity and increased cortical processing of visual information.

GenArk: towards a million UCSC genome browsers

Interactive graphical genome browsers are essential tools in genomics, but they do not contain all the recent genome assemblies. We create Genome Archive (GenArk) collection of UCSC Genome Browsers from NCBI assemblies. Built on our established track hub system, this enables fast visualization of annotations. Assemblies come with gene models, repeat masks, BLAT, and in silico PCR. Users can add annotations via track hubs and custom tracks. We can bulk-import third-party resources, demonstrated with TOGA and Ensembl gene models for hundreds of assemblies.Three thousand two hundred sixty-nine GenArk assemblies are listed at https://hgdownload.soe.ucsc.edu/hubs/ and can be searched for on the Genome Browser gateway page.

ABE8e Corrects Pax6-Aniridic Variant in Humanized Mouse ESCs and via LNPs in Ex Vivo Cortical Neurons

INTRODUCTION

Aniridia is a rare congenital vision-loss disease caused by heterozygous variants in the PAX6 gene. There is no vision-saving therapy, but one exciting approach is to use CRISPR/Cas9 to permanently correct the causal genomic variants. Preclinical studies to develop such a therapy in animal models face the challenge of showing efficacy when binding human DNA. Thus, we hypothesized that a CRISPR gene therapy can be developed and optimized in humanized mouse embryonic stem cells (ESCs) that will be able to distinguish between an aniridia patient variant and nonvariant chromosome and lay the foundation for human therapy.

METHODS

To answer the challenge of binding human DNA, we proposed the "CRISPR Humanized Minimally Mouse Models" (CHuMMMs) strategy. Thus, we minimally humanized Pax6 exon 9, the location of the most common aniridia variant c.718C > T. We generated and characterized a nonvariant CHuMMMs mouse, and a CHuMMMs cell-based disease model, in which we tested five CRISPR enzymes for therapeutic efficacy. We then delivered the therapy via lipid nanoparticles (LNPs) to alter a second variant in ex vivo cortical primary neurons.

RESULTS

We successfully established a nonvariant CHuMMMs mouse and three novel CHuMMMs aniridia cell lines. We showed that humanization did not disrupt Pax6 function in vivo, as the mouse showed no ocular phenotype. We developed and optimized a CRISPR therapeutic strategy for aniridia in the in vitro system, and found that the base editor, ABE8e, had the highest correction of the patient variant at 76.8%. In the ex vivo system, the LNP-encapsulated ABE8e ribonucleoprotein (RNP) complex altered the second patient variant and rescued 24.8% Pax6 protein expression.

CONCLUSION

We demonstrated the usefulness of the CHuMMMs approach, and showed the first genomic editing by ABE8e encapsulated as an LNP-RNP. Furthermore, we laid the foundation for translation of the proposed CRISPR therapy to preclinical mouse studies and eventually patients with aniridia.

ncOrtho: efficient and reliable identification of miRNA orthologs

MicroRNAs (miRNAs) are post-transcriptional regulators that finetune gene expression via translational repression or degradation of their target mRNAs. Despite their functional relevance, frameworks for the scalable and accurate detection of miRNA orthologs are missing. Consequently, there is still no comprehensive picture of how miRNAs and their associated regulatory networks have evolved. Here we present ncOrtho, a synteny informed pipeline for the targeted search of miRNA orthologs in unannotated genome sequences. ncOrtho matches miRNA annotations from multi-tissue transcriptomes in precision, while scaling to the analysis of hundreds of custom-selected species. The presence-absence pattern of orthologs to 266 human miRNA families across 402 vertebrate species reveals four bursts of miRNA acquisition, of which the most recent event occurred in the last common ancestor of higher primates. miRNA families are rarely modified or lost, but notable exceptions for both events exist. miRNA co-ortholog numbers faithfully indicate lineage-specific whole genome duplications, and miRNAs are powerful markers for phylogenomic analyses. Their exceptionally low genetic diversity makes them suitable to resolve clades where the phylogenetic signal is blurred by incomplete lineage sorting of ancestral alleles. In summary, ncOrtho allows to routinely consider miRNAs in evolutionary analyses that were thus far reserved to protein-coding genes.

Dysregulation of a lncRNA within the TNFRSF10A locus activates cell death pathways

TNFRSF10A (tumor necrosis factor receptor superfamily member 10A) encodes a cell surface receptor protein involved in apoptotic, necroptotic, and inflammatory pathways. Dysregulation of TNFRSF10A has been implicated in sensitization to apoptosis and to the development of multiple diseases, yet little is known of the AC100861.1 long noncoding RNA (lncRNA) that lies head-to-head with TNFRSF10A. Given its genomic positioning, we sought to investigate the function of AC100861.1, focusing on its potential relationship with TNFRSF10A and the role it may play in death receptor signaling. Using knockdown and overexpression strategies, we probed cell viability and examined transcript and protein-level changes in key genes involved in apoptosis, necroptosis, and inflammation. Decreased cell viability was observed upon TNFRSF10A overexpression, regardless of whether the cells were subjected to the chemical stressor tunicamycin. Similarly, overexpression of AC100861.1 led to increased cell death, with a further increase observed under conditions of cellular stress. Knockdown of TNFRSF10A increased cell death only when the cells were stressed, and AC100861.1 knockdown exhibited no effect on cell death. Neither knockdown nor overexpression of either of these genes greatly affected the expression of the other. Manipulating AC100861.1, however, led to marked changes in the expression of genes involved in necroptosis and inflammatory cell-signaling pathways. Additionally, RNA fluorescence in situ hybridization (RNA-FISH) revealed that the AC100861.1 transcript is localized primarily to the cytoplasm. Together, these data suggest that AC100861.1 may have a role in regulating necroptotic and inflammatory signaling pathways and that this function is separate from changes in TNFRSF10A expression. Given the importance of this genomic locus for cell survival, these data provide insight into the function of a poorly understood lncRNA with potential implications regarding disease pathology and treatment.

Transposons contribute to the acquisition of cell type-specific cis-elements in the brain

Mammalian brains have evolved in stages over a long history to acquire higher functions. Recently, several transposable element (TE) families have been shown to evolve into cis-regulatory elements of brain-specific genes. However, it is not fully understood how TEs are important for gene regulatory networks. Here, we performed a single-cell level analysis using public data of scATAC-seq to discover TE-derived cis-elements that are important for specific cell types. Our results suggest that DNA elements derived from TEs, MER130 and MamRep434, can function as transcription factor-binding sites based on their internal motifs for Neurod2 and Lhx2, respectively, especially in glutamatergic neuronal progenitors. Furthermore, MER130- and MamRep434-derived cis-elements were amplified in the ancestors of Amniota and Eutheria, respectively. These results suggest that the acquisition of cis-elements with TEs occurred in different stages during evolution and may contribute to the acquisition of different functions or morphologies in the brain.

Integrated analysis of single-cell chromatin state and transcriptome identified common vulnerability despite glioblastoma heterogeneity

In 2021, the World Health Organization reclassified glioblastoma, the most common form of adult brain cancer, into isocitrate dehydrogenase (IDH)-wild-type glioblastomas and grade IV IDH mutant (G4 IDHm) astrocytomas. For both tumor types, intratumoral heterogeneity is a key contributor to therapeutic failure. To better define this heterogeneity, genome-wide chromatin accessibility and transcription profiles of clinical samples of glioblastomas and G4 IDHm astrocytomas were analyzed at single-cell resolution. These profiles afforded resolution of intratumoral genetic heterogeneity, including delineation of cell-to-cell variations in distinct cell states, focal gene amplifications, as well as extrachromosomal circular DNAs. Despite differences in IDH mutation status and significant intratumoral heterogeneity, the profiled tumor cells shared a common chromatin structure defined by open regions enriched for nuclear factor 1 transcription factors (NFIA and NFIB). Silencing of NFIA or NFIB suppressed in vitro and in vivo growths of patient-derived glioblastomas and G4 IDHm astrocytoma models. These findings suggest that despite distinct genotypes and cell states, glioblastoma/G4 astrocytoma cells share dependency on core transcriptional programs, yielding an attractive platform for addressing therapeutic challenges associated with intratumoral heterogeneity.

The landscape of tolerated genetic variation in humans and primates

Personalized genome sequencing has revealed millions of genetic differences between individuals, but our understanding of their clinical relevance remains largely incomplete. To systematically decipher the effects of human genetic variants, we obtained whole genome sequencing data for 809 individuals from 233 primate species, and identified 4.3 million common protein-altering variants with orthologs in human. We show that these variants can be inferred to have non-deleterious effects in human based on their presence at high allele frequencies in other primate populations. We use this resource to classify 6% of all possible human protein-altering variants as likely benign and impute the pathogenicity of the remaining 94% of variants with deep learning, achieving state-of-the-art accuracy for diagnosing pathogenic variants in patients with genetic diseases.

One Sentence Summary

Deep learning classifier trained on 4.3 million common primate missense variants predicts variant pathogenicity in humans.

GenArk: Towards a million UCSC Genome Browsers

Interactive graphical genome browsers are essential tools for biologists working with DNA sequences. Although tens of thousands of new genome assemblies have become available over the last decade, accessibility is limited by the work involved in manually creating browsers and curating annotations. The results can push the limits of data storage infrastructure. To facilitate managing this increasing number of genome assemblies, we created the Genome Archive (GenArk) collection of UCSC Genome Browsers from assemblies hosted at NCBI(1). Built on our established assembly hub system, this collection enables fast, on-demand visualization of chromosome regions without requiring a database server. Available annotations include gene models, some mapped through whole-genome alignments, repeat masks, GC content, and others. We also modified our popular BLAT(2) aligner and in-silico PCR to support a large number of genomes using limited RAM. Users can upload additional annotations themselves via track hubs(3) and custom tracks. We can import more annotations in bulk from third-party resources, demonstrated here with TOGA(4) gene models. 2,430 GenArk assemblies are listed at https://hgdownload.soe.ucsc.edu/hubs/ and can be found by searching on the main UCSC gateway page. We will continue to add human high-quality assemblies and for other organisms, we are looking forward to receiving requests from the research community for ever more browsers and whole-genome alignments via http://genome.ucsc.edu/assemblyRequest.html.

Comparison of the loci associated with HbA1c and blood glucose levels identified by a genome-wide association study in the Japanese population

Aims

Hemoglobin A1c (HbA1c) levels are widely employed to diagnose diabetes. However, estimates of the heritability of HbA1c and glucose levels are different. Therefore, we explored HbA1c- and blood glucose-associated loci in a non-diabetic Japanese population.

Methods

We conducted a two-stage genome-wide association study (GWAS) on variants associated with HbA1c and blood glucose levels in a Japanese population. In the initial stage, data of 4911 participants of the Japan Multi-Institutional Collaborative Cohort (J-MICC) were subjected to discovery analysis. In the second stage, two datasets from the Tohoku Medical Megabank project, with 8175 and 40,519 participants, were used for the replication study. Association of the imputed variants with HbA1c and blood glucose levels was determined via linear regression analyses adjusted for age, sex, body mass index (BMI), smoking, and genetic principal components (PC1-PC10). Moreover, we performed a BMI-stratified GWAS on HbA1c levels in the J-MICC. The discovery analysis and BMI-stratified GWAS results were validated with re-analyses of normalized HbA1c levels adjusted for site in addition to the above, and blood glucose adjusted for fasting time as an additional covariate.

Results

Genetic variants associated with HbA1c levels were identified in KCNQ1 and TMC6. None of the genetic variants associated with blood glucose levels in the discovery analysis were replicated. Association of rs2299620 in KCNQ1 with HbA1c levels showed heterogeneity between individuals with BMI ≥ 25 kg/m2 and BMI < 25 kg/m2.

Conclusions

The variant rs2299620 in KCNQ1 might affect HbA1c levels differentially based on BMI grouping in the Japanese population.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13340-023-00618-0.

Exploring microbial functional biodiversity at the protein family level-From metagenomic sequence reads to annotated protein clusters

Metagenomics has enabled accessing the genetic repertoire of natural microbial communities. Metagenome shotgun sequencing has become the method of choice for studying and classifying microorganisms from various environments. To this end, several methods have been developed to process and analyze the sequence data from raw reads to end-products such as predicted protein sequences or families. In this article, we provide a thorough review to simplify such processes and discuss the alternative methodologies that can be followed in order to explore biodiversity at the protein family level. We provide details for analysis tools and we comment on their scalability as well as their advantages and disadvantages. Finally, we report the available data repositories and recommend various approaches for protein family annotation related to phylogenetic distribution, structure prediction and metadata enrichment.

Gene expression variability across cells and species shapes the relationship between renal resident macrophages and infiltrated macrophages

BACKGROUND

Two main subclasses of macrophages are found in almost all solid tissues: embryo-derived resident tissue macrophages and bone marrow-derived infiltrated macrophages. These macrophage subtypes show transcriptional and functional divergence, and the programs that have shaped the evolution of renal macrophages and related signaling pathways remain poorly understood. To clarify these processes, we performed data analysis based on single-cell transcriptional profiling of renal tissue-resident and infiltrated macrophages in human, mouse and rat.

RESULTS

In this study, we (i) characterized the transcriptional divergence among species and (ii) illustrated variability in expression among cells of each subtype and (iii) compared the gene regulation network and (iv) ligand-receptor pairs in human and mouse. Using single-cell transcriptomics, we mapped the promoter architecture during homeostasis.

CONCLUSIONS

Transcriptionally divergent genes, such as the differentially TF-encoding genes expressed in resident and infiltrated macrophages across the three species, vary among cells and include distinct promoter structures. The gene regulatory network in infiltrated macrophages shows comparatively better species-wide consistency than resident macrophages. The conserved transcriptional gene regulatory network in infiltrated macrophages among species is uniquely enriched in pathways related to kinases, and TFs associated with largely conserved regulons among species are uniquely enriched in kinase-related pathways.

WINNER: A network biology tool for biomolecular characterization and prioritization

BACKGROUND AND CONTRIBUTION

In network biology, molecular functions can be characterized by network-based inference, or "guilt-by-associations." PageRank-like tools have been applied in the study of biomolecular interaction networks to obtain further the relative significance of all molecules in the network. However, there is a great deal of inherent noise in widely accessible data sets for gene-to-gene associations or protein-protein interactions. How to develop robust tests to expand, filter, and rank molecular entities in disease-specific networks remains an ad hoc data analysis process.

RESULTS

We describe a new biomolecular characterization and prioritization tool called Weighted In-Network Node Expansion and Ranking (WINNER). It takes the input of any molecular interaction network data and generates an optionally expanded network with all the nodes ranked according to their relevance to one another in the network. To help users assess the robustness of results, WINNER provides two different types of statistics. The first type is a node-expansion p-value, which helps evaluate the statistical significance of adding "non-seed" molecules to the original biomolecular interaction network consisting of "seed" molecules and molecular interactions. The second type is a node-ranking p-value, which helps evaluate the relative statistical significance of the contribution of each node to the overall network architecture. We validated the robustness of WINNER in ranking top molecules by spiking noises in several network permutation experiments. We have found that node degree-preservation randomization of the gene network produced normally distributed ranking scores, which outperform those made with other gene network randomization techniques. Furthermore, we validated that a more significant proportion of the WINNER-ranked genes was associated with disease biology than existing methods such as PageRank. We demonstrated the performance of WINNER with a few case studies, including Alzheimer's disease, breast cancer, myocardial infarctions, and Triple negative breast cancer (TNBC). In all these case studies, the expanded and top-ranked genes identified by WINNER reveal disease biology more significantly than those identified by other gene prioritizing software tools, including Ingenuity Pathway Analysis (IPA) and DiAMOND.

CONCLUSION

WINNER ranking strongly correlates to other ranking methods when the network covers sufficient node and edge information, indicating a high network quality. WINNER users can use this new tool to robustly evaluate a list of candidate genes, proteins, or metabolites produced from high-throughput biology experiments, as long as there is available gene/protein/metabolic network information.

Systemic lupus erythematosus variants modulate the function of an enhancer upstream of TNFAIP3

TNFAIP3/A20 is a prominent autoimmune disease risk locus that is correlated with hypomorphic TNFAIP3 expression and exhibits complex chromatin architecture with over 30 predicted enhancers. This study aimed to functionally characterize an enhancer ∼55 kb upstream of the TNFAIP3 promoter marked by the systemic lupus erythematosus (SLE) risk haplotype index SNP, rs10499197. Allele effects of rs10499197, rs58905141, and rs9494868 were tested by EMSA and/or luciferase reporter assays in immune cell types. Co-immunoprecipitation, ChIP-qPCR, and 3C-qPCR were performed on patient-derived EBV B cells homozygous for the non-risk or SLE risk TNFAIP3 haplotype to assess haplotype-specific effects on transcription factor binding and chromatin regulation at the TNFAIP3 locus. This study found that the TNFAIP3 locus has a complex chromatin regulatory network that spans ∼1M bp from the promoter region of IL20RA to the 3' untranslated region of TNFAIP3. Functional dissection of the enhancer demonstrated co-dependency of the RelA/p65 and CEBPB binding motifs that, together, increase IL20RA and IFNGR1 expression and decreased TNFAIP3 expression in the context of the TNFAIP3 SLE risk haplotype through dynamic long-range interactions up- and downstream. Examination of SNPs in linkage disequilibrium (D' = 1.0) with rs10499197 identified rs9494868 as a functional SNP with risk allele-specific increase in nuclear factor binding and enhancer activation in vitro. In summary, this study demonstrates that SNPs carried on the ∼109 kb SLE risk haplotype facilitate hypermorphic IL20RA and IFNGR1 expression, while suppressing TNFAIP3 expression, adding to the mechanistic potency of this critically important locus in autoimmune disease pathology.

HPV genotyping by L1 amplicon sequencing of archived invasive cervical cancer samples: a pilot study

BACKGROUND

Human papillomavirus (HPV) is the primary cause of invasive cervical cancer (ICC). The prevalence of various HPV genotypes, ranging from oncogenically low- to high-risk, may be influenced by geographic and demographic factors, which could have critical implications for the screening and prevention of HPV infection and ICC incidence. However, many technical factors may influence the identification of high-risk genotypes associated with ICC in different populations.

METHODS

We used high-throughput sequencing of a single amplicon within the HPV L1 gene to assess the influence of patient age, race/ethnicity, histological subtype, sample type, collection date, experimental factors, and computational parameters on the prevalence of HPV genotypes detected in archived DNA (n = 34), frozen tissue (n = 44), and formalin-fixed paraffin-embedded (FFPE) tissue (n = 57) samples collected in the Los Angeles metropolitan area.

RESULTS

We found that the percentage of off-target human reads and the concentration of DNA amplified from each sample varied by HPV genotype and by archive type. After accounting for the percentage of human reads and excluding samples with especially low levels of amplified DNA, the HPV prevalence was 95% across all ICC samples: HPV16 was the most common genotype (in 56% of all ICC samples), followed by HPV18 (in 21%). Depending upon the genotyping parameters, the prevalence of HPV58 varied up to twofold in our cohort. In archived DNA and frozen tissue samples, we detected previously established differences in HPV16 and HPV18 frequencies based on histological subtype, but we could not reproduce those findings using our FFPE samples.

CONCLUSIONS

In this pilot study, we demonstrate that sample collection, preparation, and analysis methods can influence the detection of certain HPV genotypes and must be carefully considered when drawing any biological conclusions based on HPV genotyping data from ICC samples.

A genome-wide association study on adherence to low-carbohydrate diets in Japanese

BACKGROUND/OBJECTIVES

Low-carbohydrate diets (LCD) are useful for weight reduction, and 50-55% carbohydrate consumption is associated with minimal risk. Genetic differences were related to nutritional consumption, food preferences, and dietary patterns, but whether particular genetic differences in individuals influence LCD adherence is unknown.

SUBJECTS/METHODS

We conducted a GWAS on adherence to LCD utilizing 14,076 participants from the Japan Multi-Institutional Collaborative Cohort study. We used a previously validated semiquantitative food frequency questionnaire to estimate food consumption. Association of the imputed variants with the LCD score by Halton et al. we used linear regression analysis adjusting for sex, age, total dietary energy consumption, and components 1 to 10 by principal component analysis. We repeated the analysis with adjustment for alcohol consumption (g/day) in addition to the above-described variables.

RESULTS

Men and women combined analysis without adjustment for alcohol consumption; we found 395 variants on chromosome 12 associated with the LCD score having P values <5 × 10^-8. A conditional analysis with the addition of the dosage data of rs671 on chromosome 12 as a covariate, P values for all 395 SNPs on chromosome 12 turned out to be insignificant. In the analysis with additional adjustment for alcohol consumption, we did not identify any SNPs associated with the LCD score.

CONCLUSION

We found rs671 was inversely associated with adherence to LCD, but that was strongly confounded by alcohol consumption.

Downregulation of CRTC1 Is Involved in CUMS-Induced Depression-Like Behavior in the Hippocampus and Its RNA Sequencing Analysis

Chronic stress is an important risk factor for mood disorders including depression. The decreased level of CREB (cAMP-responsive element binding)-regulated transcription coactivator 1 (CRTC1) expression in hippocampus may be involved in depression-like behavior in some stress-induced depression models. But the mechanism of CRTC1 in mediating depression-like behavior remains unknown. In this study, chronic unpredictable mild stress (CUMS)-treated mice showed depression-like behavior accompanied by the downregulation of CRTC1 in the hippocampus. Adeno-associated virus (AAV)-CRTC1-mediated overexpression of CRTC1 in the hippocampus by stereotactic brain injection could significantly prevent depression-like behavior in CUMS-treated mice. The above data reveal that the downregulation of hippocampal CRTC1 expression participates in CUMS-induced depression-like behavior. In order to explore the key targets regulated by CRTC1, AAV-mediated CRTC1 short hairpin (shRNA) was constructed to achieve knockdown of CRTC1 in the hippocampus, and then the hippocampi were collected for RNA-sequencing (RNA-seq). The RNA-seq data show that upregulated genes were enriched in stress and immune system-associated GO terms and pathways such as response to stress and external stimulus and regulation of immune response and that downregulated genes were enriched in neural activity such as synaptic transmission and cognitive behavior. We further provided RT-qPCR data that the inflammation-related factors including Gpr84, Tlr2, Lyz2, and Icam1 were significantly upregulated in the hippocampus of both CUMS- and CRTC1 shRNA-induced models, some of them were also validated in protein levels by Western blotting. We propose a hypothesis that CUMS induces downregulation of CRTC1, which might lead to depression-like behavior via neuroinflammation pathway. This study provides new explanation for the inflammatory hypothesis of depression and some clues for exploring the molecular mechanism of CRTC1 regulation.

Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin

BACKGROUND

Evolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation of many enhancers, distal regulatory elements that control spatial and temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high.

RESULTS

We present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that can predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Next, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss of regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity have a stronger tendency than the general population to have predicted lineage-specific open chromatin.

CONCLUSION

The framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements.

Whole-genome sequencing reveals that variants in the Interleukin 18 Receptor Accessory Protein 3'UTR protect against ALS

The noncoding genome is substantially larger than the protein-coding genome but has been largely unexplored by genetic association studies. Here, we performed region-based rare variant association analysis of >25,000 variants in untranslated regions of 6,139 amyotrophic lateral sclerosis (ALS) whole genomes and the whole genomes of 70,403 non-ALS controls. We identified interleukin-18 receptor accessory protein (IL18RAP) 3' untranslated region (3'UTR) variants as significantly enriched in non-ALS genomes and associated with a fivefold reduced risk of developing ALS, and this was replicated in an independent cohort. These variants in the IL18RAP 3'UTR reduce mRNA stability and the binding of double-stranded RNA (dsRNA)-binding proteins. Finally, the variants of the IL18RAP 3'UTR confer a survival advantage for motor neurons because they dampen neurotoxicity of human induced pluripotent stem cell (iPSC)-derived microglia bearing an ALS-associated expansion in C9orf72, and this depends on NF-κB signaling. This study reveals genetic variants that protect against ALS by reducing neuroinflammation and emphasizes the importance of noncoding genetic association studies.

In vivo temporal resolution of acute promyelocytic leukemia progression reveals a role of Klf4 in suppressing early leukemic transformation

In this study, Mas et al. used primary hematopoietic stem and progenitor cells (HSPCs) and leukemic blasts that express the fusion protein PML-RARα as a paradigm to temporally dissect the dynamic changes in the epigenome, transcriptome, and genome architecture induced during oncogenic transformation. Their multiomics-integrated analysis identified Klf4 as an early down-regulated gene in PML-RARα-driven leukemogenesis, and they characterized the dynamic alterations in the Klf4 cis-regulatory network during APL progression and demonstrated that ectopic Klf4 overexpression can suppress self-renewal and reverse the differentiation block induced by PML-RARα.

Genome organization plays a pivotal role in transcription, but how transcription factors (TFs) rewire the structure of the genome to initiate and maintain the programs that lead to oncogenic transformation remains poorly understood. Acute promyelocytic leukemia (APL) is a fatal subtype of leukemia driven by a chromosomal translocation between the promyelocytic leukemia (PML) and retinoic acid receptor α (RARα) genes. We used primary hematopoietic stem and progenitor cells (HSPCs) and leukemic blasts that express the fusion protein PML-RARα as a paradigm to temporally dissect the dynamic changes in the epigenome, transcriptome, and genome architecture induced during oncogenic transformation. We found that PML-RARα initiates a continuum of topologic alterations, including switches from A to B compartments, transcriptional repression, loss of active histone marks, and gain of repressive histone marks. Our multiomics-integrated analysis identifies Klf4 as an early down-regulated gene in PML-RARα-driven leukemogenesis. Furthermore, we characterized the dynamic alterations in the Klf4 cis-regulatory network during APL progression and demonstrated that ectopic Klf4 overexpression can suppress self-renewal and reverse the differentiation block induced by PML-RARα. Our study provides a comprehensive in vivo temporal dissection of the epigenomic and topological reprogramming induced by an oncogenic TF and illustrates how topological architecture can be used to identify new drivers of malignant transformation.

Naïve-like pluripotency to pave the way for saving the northern white rhinoceros from extinction

The northern white rhinoceros (NWR) is probably the earth's most endangered mammal. To rescue the functionally extinct species, we aim to employ induced pluripotent stem cells (iPSCs) to generate gametes and subsequently embryos in vitro. To elucidate the regulation of pluripotency and differentiation of NWR PSCs, we generated iPSCs from a deceased NWR female using episomal reprogramming, and observed surprising similarities to human PSCs. NWR iPSCs exhibit a broad differentiation potency into the three germ layers and trophoblast, and acquire a naïve-like state of pluripotency, which is pivotal to differentiate PSCs into primordial germ cells (PGCs). Naïve culturing conditions induced a similar expression profile of pluripotency related genes in NWR iPSCs and human ESCs. Furthermore, naïve-like NWR iPSCs displayed increased expression of naïve and PGC marker genes, and a higher integration propensity into developing mouse embryos. As the conversion process was aided by ectopic BCL2 expression, and we observed integration of reprogramming factors, the NWR iPSCs presented here are unsuitable for gamete production. However, the gained insights into the developmental potential of both primed and naïve-like NWR iPSCs are fundamental for in future PGC-specification in order to rescue the species from extinction using cryopreserved somatic cells.

A comprehensive bioinformatics analysis to identify a candidate prognostic biomarker for ovarian cancer

BACKGROUND

This study aimed to investigate prognostic genes in ovarian cancer (OC) and to explore their potential underlying biological mechanisms through a comprehensive bioinformatics analysis.

METHODS

Common differentially expressed genes (DEGs) in 3 OC datasets from the Gene Expression Omnibus (GEO) (GSE26712, GSE18520, and GSE14407) were screened out. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed by Metascape. The protein-protein interaction (PPI) network of the DEGs was constructed using the STRING database. The prognostic value of DEGs were determined using the Kaplan-Meier plotter. The ONCOMINE and Human Protein Atlas databases were used to verify the expression levels of prognostic genes in OC. Genomic analysis of prognostic genes were also investigated by cBio Cancer Genomics Portal (cBioPortal) database, UCSC Xena browser and UALCAN. Gene set enrichment analysis (GSEA) was used to predict the possible pathways and biological processes of the prognostic genes.

RESULTS

Integration of the 3 datasets have found 879 common DEGs. A high expression of structural maintenance of chromosomes protein 4 (SMC4) was revealed in the Kaplan-Meier plotter analysis to be meaningful for the prognosis of OC and was verified at both the mRNA and protein levels. The results from cBioPortal showed that SMC4 alterations accounted for 7 to 18% of genetic alterations in OC, and the majority alterations were copy number amplifications. Finally, the GSEA results showed that samples with SMC4 overexpression were mainly enriched in the cell cycle, spliceosome, ubiquitin mediated proteolysis, and adherens junctions.

CONCLUSIONS

High SMC4 expression is linked with a poor prognosis in patients with OC and might serve as a prognostic biomarker for the disease.

Ferroptosis Patterns Correlate with Immune Microenvironment Characterization in Gastric Cancer

Objective

We aimed to build a ferroptosis-based classifier to characterize the molecular features of gastric cancers (GC) and investigate the relationship between different ferroptosis patterns and GC tumor microenvironment (TME).

Methods

Based on the genomic and clinical information from TCGA portal and GEO database, non-negative matrix factorization (NMF) was used to identify ferroptosis subtypes in GC patients. In order to estimate the ferroptosis levels, we established ferroptosis subtype score (FSS) to quantify ferroptosis patterns and ferroptosis potential index (FPI) by principal component analysis (PCA). The correlations of different ferroptosis patterns with TME cell-infiltrating characteristics (including immune cell infiltration, immune checkpoints expression levels, tumor mutational burden (TMB) and immunotherapy response) were systematically analyzed.

Results

Two ferroptosis subtypes, C1 (with lower FSS) and C2 (with higher FSS), were determined. C2 displayed a significantly lower FPI than C1. Besides, C2 was associated with diffuse subtype while C1 with intestinal subtype. As for TME characteristics, C2 was in accordance with the immune-excluded phenotype as it showed more active immune and stromal activities but lower TMB, less probability of immunotherapy response and poorer prognosis. C1 was linked to immune-inflamed phenotype as it had lower stromal activities but increased neoantigen load, enhanced response to immunotherapy and relatively better prognosis.

Conclusion

The systematic assessment of ferroptosis patterns and ferroptosis levels presented in our study implied that ferroptosis serves as an important factor in the formation of TME, which may expand the understanding of TME and provide a novel perspective for the development of targeted immunotherapeutic strategies for GC patients.

An atlas of gene regulatory elements in adult mouse cerebrum

The mammalian cerebrum performs high-level sensory perception, motor control and cognitive functions through highly specialized cortical and subcortical structures1. Recent surveys of mouse and human brains with single-cell transcriptomics2-6 and high-throughput imaging technologies7,8 have uncovered hundreds of neural cell types distributed in different brain regions, but the transcriptional regulatory programs that are responsible for the unique identity and function of each cell type remain unknown. Here we probe the accessible chromatin in more than 800,000 individual nuclei from 45 regions that span the adult mouse isocortex, olfactory bulb, hippocampus and cerebral nuclei, and use the resulting data to map the state of 491,818 candidate cis-regulatory DNA elements in 160 distinct cell types. We find high specificity of spatial distribution for not only excitatory neurons, but also most classes of inhibitory neurons and a subset of glial cell types. We characterize the gene regulatory sequences associated with the regional specificity within these cell types. We further link a considerable fraction of the cis-regulatory elements to putative target genes expressed in diverse cerebral cell types and predict transcriptional regulators that are involved in a broad spectrum of molecular and cellular pathways in different neuronal and glial cell populations. Our results provide a foundation for comprehensive analysis of gene regulatory programs of the mammalian brain and assist in the interpretation of noncoding risk variants associated with various neurological diseases and traits in humans.

Paths and timings of the peopling of Polynesia inferred from genomic networks

Polynesia was settled in a series of extraordinary voyages across an ocean spanning one third of the Earth1, but the sequences of islands settled remain unknown and their timings disputed. Currently, several centuries separate the dates suggested by different archaeological surveys2-4. Here, using genome-wide data from merely 430 modern individuals from 21 key Pacific island populations and novel ancestry-specific computational analyses, we unravel the detailed genetic history of this vast, dispersed island network. Our reconstruction of the branching Polynesian migration sequence reveals a serial founder expansion, characterized by directional loss of variants, that originated in Samoa and spread first through the Cook Islands (Rarotonga), then to the Society (Tōtaiete mā) Islands (11th century), the western Austral (Tuha'a Pae) Islands and Tuāmotu Archipelago (12th century), and finally to the widely separated, but genetically connected, megalithic statue-building cultures of the Marquesas (Te Henua 'Enana) Islands in the north, Raivavae in the south, and Easter Island (Rapa Nui), the easternmost of the Polynesian islands, settled in approximately AD 1200 via Mangareva.

LncRNA KCNQ1OT1 activated by c-Myc promotes cell proliferation via interacting with FUS to stabilize MAP3K1 in acute promyelocytic leukemia

Uncontrolled proliferation is the hallmark of cancer cells. Previous studies mainly focused on the role of protein-coding genes in cancer cell proliferation. Emerging evidence showed that long non-coding RNAs (lncRNAs) also play critical roles in cancer cell proliferation and growth. LncRNA KCNQ1OT1 is found to contribute to carcinogenesis, but its role in acute promyelocytic leukemia (APL) is unclear. In this study, by analyzing data from Gene Expression Omnibus, The Cancer Genome Atlas database and our clinical samples, we found that KCNQ1OT1 was selectively highly expressed in APL. Functional assays demonstrated that knockdown of KCNQ1OT1 reduced APL cell proliferation and increased apoptosis. Further evidence showed that KCNQ1OT1 was mainly located in the cytoplasm of APL patient-derived NB4 cells and APL patient bone marrow samples. Mechanistically, KCNQ1OT1 bound to RNA binding protein FUS, and silencing either KCNQ1OT1 or FUS reduced the expression level and stability of MAP3K1 mRNA. Whereas KCNQ1OT1 and FUS did not affect each other. Importantly, knockdown of MAP3K1 impaired APL cell proliferation. Finally, c-Myc transactivated KCNQ1OT1 in APL cells through binding to its promoter while knockdown of c-Myc decreased KCNQ1OT1 expression. Our results not only revealed that c-Myc transactivated KCNQ1OT1 and upregulated KCNQ1OT1 promoted APL cell proliferation, but also demonstrated that KCNQ1OT1 bound to FUS to synergistically stabilize MAP3K1 mRNA, thus facilitating APL cell proliferation. This study established a previously unidentified role of KCNQ1OT1 in the development of APL, and KCNQ1OT1 may serve as a potential therapeutic target for APL.

PacBio Iso-Seq Improves the Rainbow Trout Genome Annotation and Identifies Alternative Splicing Associated With Economically Important Phenotypes

Rainbow trout is an important model organism that has received concerted international efforts to study the transcriptome. For this purpose, short-read sequencing has been primarily used over the past decade. However, these sequences are too short of resolving the transcriptome complexity. This study reported a first full-length transcriptome assembly of the rainbow trout using single-molecule long-read isoform sequencing (Iso-Seq). Extensive computational approaches were used to refine and validate the reconstructed transcriptome. The study identified 10,640 high-confidence transcripts not previously annotated, in addition to 1,479 isoforms not mapped to the current Swanson reference genome. Most of the identified lncRNAs were non-coding variants of coding transcripts. The majority of genes had multiple transcript isoforms (average ∼3 isoforms/locus). Intron retention (IR) and exon skipping (ES) accounted for 56% of alternative splicing (AS) events. Iso-Seq improved the reference genome annotation, which allowed identification of characteristic AS associated with fish growth, muscle accretion, disease resistance, stress response, and fish migration. For instance, an ES in GVIN1 gene existed in fish susceptible to bacterial cold-water disease (BCWD). Besides, under five stress conditions, there was a commonly regulated exon in prolyl 4-hydroxylase subunit alpha-2 (P4HA2) gene. The reconstructed gene models and their posttranscriptional processing in rainbow trout provide invaluable resources that could be further used for future genetics and genomics studies. Additionally, the study identified characteristic transcription events associated with economically important phenotypes, which could be applied in selective breeding.

Bioinformatics integrated analysis to investigate candidate biomarkers and associated metabolites in osteosarcoma

BACKGROUND

This study hoped to explore the potential biomarkers and associated metabolites during osteosarcoma (OS) progression based on bioinformatics integrated analysis.

METHODS

Gene expression profiles of GSE28424, including 19 human OS cell lines (OS group) and 4 human normal long bone tissue samples (control group), were downloaded. The differentially expressed genes (DEGs) in OS vs. control were investigated. The enrichment investigation was performed based on DEGs, followed by protein-protein interaction network analysis. Then, the feature genes associated with OS were explored, followed by survival analysis to reveal prognostic genes. The qRT-PCR assay was performed to test the expression of these genes. Finally, the OS-associated metabolites and disease-metabolic network were further investigated.

RESULTS

Totally, 357 DEGs were revealed between the OS vs. control groups. These DEGs, such as CXCL12, were mainly involved in functions like leukocyte migration. Then, totally, 38 feature genes were explored, of which 8 genes showed significant associations with the survival of patients. High expression of CXCL12, CEBPA, SPARCL1, CAT, TUBA1A, and ALDH1A1 was associated with longer survival time, while high expression of CFLAR and STC2 was associated with poor survival. Finally, a disease-metabolic network was constructed with 25 nodes including two disease-associated metabolites cyclophosphamide and bisphenol A (BPA). BPA showed interactions with multiple prognosis-related genes, such as CXCL12 and STC2.

CONCLUSION

We identified 8 prognosis-related genes in OS. CXCL12 might participate in OS progression via leukocyte migration function. BPA might be an important metabolite interacting with multiple prognosis-related genes.

Prime editing in mice reveals the essentiality of a single base in driving tissue-specific gene expression

BACKGROUND

Most single nucleotide variants (SNVs) occur in noncoding sequence where millions of transcription factor binding sites (TFBS) reside. Here, a comparative analysis of CRISPR-mediated homology-directed repair (HDR) versus the recently reported prime editing 2 (PE2) system was carried out in mice over a TFBS called a CArG box in the Tspan2 promoter.

RESULTS

Quantitative RT-PCR showed loss of Tspan2 mRNA in aorta and bladder, but not heart or brain, of mice homozygous for an HDR-mediated three base pair substitution in the Tspan2 CArG box. Using the same protospacer, mice homozygous for a PE2-mediated single-base substitution in the Tspan2 CArG box displayed similar cell-specific loss of Tspan2 mRNA; expression of an overlapping long noncoding RNA was also nearly abolished in aorta and bladder. Immuno-RNA fluorescence in situ hybridization validated loss of Tspan2 in vascular smooth muscle cells of HDR and PE2 CArG box mutant mice. Targeted sequencing demonstrated variable frequencies of on-target editing in all PE2 and HDR founders. However, whereas no on-target indels were detected in any of the PE2 founders, all HDR founders showed varying levels of on-target indels. Off-target analysis by targeted sequencing revealed mutations in many HDR founders, but none in PE2 founders.

CONCLUSIONS

PE2 directs high-fidelity editing of a single base in a TFBS leading to cell-specific loss in expression of an mRNA/long noncoding RNA gene pair. The PE2 platform expands the genome editing toolbox for modeling and correcting relevant noncoding SNVs in the mouse.

COMSUC: A web server for the identification of consensus molecular subtypes of cancer based on multiple methods and multi-omics data

Extensive amounts of multi-omics data and multiple cancer subtyping methods have been developed rapidly, and generate discrepant clustering results, which poses challenges for cancer molecular subtype research. Thus, the development of methods for the identification of cancer consensus molecular subtypes is essential. The lack of intuitive and easy-to-use analytical tools has posed a barrier. Here, we report on the development of the COnsensus Molecular SUbtype of Cancer (COMSUC) web server. With COMSUC, users can explore consensus molecular subtypes of more than 30 cancers based on eight clustering methods, five types of omics data from public reference datasets or users' private data, and three consensus clustering methods. The web server provides interactive and modifiable visualization, and publishable output of analysis results. Researchers can also exchange consensus subtype results with collaborators via project IDs. COMSUC is now publicly and freely available with no login requirement at http://comsuc.bioinforai.tech/ (IP address: http://59.110.25.27/). For a video summary of this web server, see S1 Video and S1 File.

Conserved regulatory logic at accessible and inaccessible chromatin during the acute inflammatory response in mammals

The regulatory elements controlling gene expression during acute inflammation are not fully elucidated. Here we report the identification of a set of NF-κB-bound elements and common chromatin landscapes underlying the acute inflammatory response across cell-types and mammalian species. Using primary vascular endothelial cells (human/mouse/bovine) treated with the pro-inflammatory cytokine, Tumor Necrosis Factor-α, we identify extensive (~30%) conserved orthologous binding of NF-κB to accessible, as well as nucleosome-occluded chromatin. Regions with the highest NF-κB occupancy pre-stimulation show dramatic increases in NF-κB binding and chromatin accessibility post-stimulation. These 'pre-bound' regions are typically conserved (~56%), contain multiple NF-κB motifs, are utilized by diverse cell types, and overlap rare non-coding mutations and common genetic variation associated with both inflammatory and cardiovascular phenotypes. Genetic ablation of conserved, 'pre-bound' NF-κB regions within the super-enhancer associated with the chemokine-encoding CCL2 gene and elsewhere supports the functional relevance of these elements.

Multiple microarray analyses identify key genes associated with the development of Non-Small Cell Lung Cancer from Chronic Obstructive Pulmonary Disease

Introduction: Chronic obstructive pulmonary disease (COPD) is an independent risk factor of non-small cell lung cancer (NSCLC). This study aimed to analyze the key genes and potential molecular mechanisms that are involved in the development from COPD to NSCLC.

Methods: Expression profiles of COPD and NSCLC in GSE106899, GSE12472, and GSE12428 were downloaded from the Gene Expression Omnibus (GEO) database, followed by identification of the differentially expressed genes (DEGs) between COPD and NSCLC. Based on the identified DEGs, functional pathway enrichment and lung carcinogenesis-related networks analyses were performed and further visualized with Cytoscape software. Then, principal component analysis (PCA), cluster analysis, and support vector machines (SVM) verified the ability of the top modular genes to distinguish COPD from NSCLC. Additionally, the corrections between these key genes and clinical staging of NSCLC were studied using the UALCAN and HPA websites. Finally, a prognostic risk model was constructed based on multivariate Cox regression analysis. Kaplan-Meier survival curves of the top modular genes on the training and verification sets were generated.

Results: A total of 2350, 1914, and 1850 DEGs were obtained from GSE106899, GSE12472, and GSE12428 datasets, respectively. Following analysis of protein-protein interaction networks, the identified modular gene signatures containing H2AFX, MCM2, MCM3, MCM7, POLD1, and RPA1 were identified as markers for discrimination between COPD and NSCLC. The modular gene signatures were mainly enriched in the processes of DNA replication, cell cycle, mismatch repair, and others. Besides, the expression levels of these genes were significantly higher in NSCLC than in COPD, which was further verified by the immunohistochemistry. In addition, the high expression levels of H2AFX, MCM2, MCM7, and POLD1 correlate with poor prognosis of lung adenocarcinoma (LUAD). The Cox regression prognostic risk model showed the similar results and the predictive ability of this model is independent of other clinical variables.

Conclusions: This study revealed several key modules that closely relate to NSCLC with underlying disease COPD, which provide a deeper understanding of the potential mechanisms underlying the malignant development from COPD to NSCLC. This study provides valuable prognostic factors in high-risk lung cancer patients with COPD.

A genome-wide association study in Japanese identified one variant associated with a preference for a Japanese dietary pattern

BACKGROUND/OBJECTIVES

Individual eating habits may be influenced by genetic factors, in addition to environmental factors. Previous studies suggested that adherence to Japanese food patterns was associated with a decreased risk of all-cause and cardiovascular disease mortality. We conducted a genome-wide association study (GWAS) in a Japanese population to find genetic variations that affect adherence to a Japanese food pattern.

SUBJECTS/METHODS

We analyzed GWAS data using 14,079 participants from the Japan Multi-Institutional Collaborative Cohort study. We made a Japanese food score based on six food groups. Association of the imputed variants with the Japanese food score was performed by linear regression analysis with adjustments for age, sex, total energy intake, alcohol intake (g/day), and principal components 1-10 omitting variants in the major histocompatibility region.

RESULTS

We found one SNP in the 14q11.2 locus that was significantly associated with the Japanese food score with P values <5 × 10^-8. Functional annotation revealed that the expression levels of two genes (BCL2L2, SLC22A17) were significantly inversely associated with this SNP. These genes are known to be related to olfaction and obesity.

CONCLUSION

We found a new SNP that was associated with the Japanese food score in a Japanese population. This SNP is inversely associated with genes link to olfaction and obesity.

Prevalent intron retention fine-tunes gene expression and contributes to cellular senescence

Intron retention (IR) is the least well‐understood alternative splicing type in animals, and its prevalence and function in physiological and pathological processes have long been underestimated. Cellular senescence contributes to individual aging and age‐related diseases and can also serve as an important cancer prevention mechanism. Dynamic IR events have been observed in senescence models and aged tissues; however, whether and how IR impacts senescence remain unclear. Through analyzing polyA⁺ RNA‐seq data from human replicative senescence models, we found IR was prevalent and dynamically regulated during senescence and IR changes negatively correlated with expression alteration of corresponding genes. We discovered that knocking down (KD) splicing factor U2AF1, which showed higher binding density to retained introns and decreased expression during senescence, led to senescence‐associated phenotypes and global IR changes. Intriguingly, U2AF1‐KD‐induced IR changes also negatively correlated with gene expression. Furthermore, we demonstrated that U2AF1‐mediated IR of specific gene (CPNE1 as an example) contributed to cellular senescence. Decreased expression of U2AF1, higher IR of CPNE1, and reduced expression of CPNE1 were also discovered in dermal fibroblasts with age. We discovered prevalent IR could fine‐tune gene expression and contribute to senescence‐associated phenotypes, largely extending the biological significance of IR.

Molecular consequences of fetal alcohol exposure on amniotic exosomal miRNAs with functional implications for stem cell potency and differentiation

Alcohol (ethanol, EtOH) consumption during pregnancy can result in fetal alcohol spectrum disorders (FASDs), which are characterized by prenatal and postnatal growth restriction and craniofacial dysmorphology. Recently, cell-derived extracellular vesicles, including exosomes and microvesicles containing several species of RNAs (exRNAs), have emerged as a mechanism of cell-to-cell communication. However, EtOH's effects on the biogenesis and function of non-coding exRNAs during fetal development have not been explored. Therefore, we studied the effects of maternal EtOH exposure on the composition of exosomal RNAs in the amniotic fluid (AF) using rat fetal alcohol exposure (FAE) model. Through RNA-Seq analysis we identified and verified AF exosomal miRNAs with differential expression levels specifically associated with maternal EtOH exposure. Uptake of purified FAE AF exosomes by rBMSCs resulted in significant alteration of molecular markers associated with osteogenic differentiation of rBMSCs. We also determined putative functional roles for AF exosomal miRNAs (miR-199a-3p, miR-214-3p and let-7g) that are dysregulated by FAE in osteogenic differentiation of rBMSCs. Our results demonstrate that FAE alters AF exosomal miRNAs and that exosomal transfer of dysregulated miRNAs has significant molecular effects on stem cell regulation and differentiation. Our results further suggest the usefulness of assessing molecular alterations in AF exRNAs to study the mechanisms of FAE teratogenesis that should be further investigated by using an in vivo model.

Single-cell multiomic profiling of human lungs reveals cell-type-specific and age-dynamic control of SARS-CoV2 host genes

Respiratory failure associated with COVID-19 has placed focus on the lungs. Here, we present single-nucleus accessible chromatin profiles of 90,980 nuclei and matched single-nucleus transcriptomes of 46,500 nuclei in non-diseased lungs from donors of ~30 weeks gestation,~3 years and ~30 years. We mapped candidate cis-regulatory elements (cCREs) and linked them to putative target genes. We identified distal cCREs with age-increased activity linked to SARS-CoV-2 host entry gene TMPRSS2 in alveolar type 2 cells, which had immune regulatory signatures and harbored variants associated with respiratory traits. At the 3p21.31 COVID-19 risk locus, a candidate variant overlapped a distal cCRE linked to SLC6A20, a gene expressed in alveolar cells and with known functional association with the SARS-CoV-2 receptor ACE2. Our findings provide insight into regulatory logic underlying genes implicated in COVID-19 in individual lung cell types across age. More broadly, these datasets will facilitate interpretation of risk loci for lung diseases.

PINCER: improved CRISPR/Cas9 screening by efficient cleavage at conserved residues

CRISPR/Cas9 functional genomic screens have emerged as essential tools in drug target discovery. However, the sensitivity of available genome-wide CRISPR libraries is impaired by guides which inefficiently abrogate gene function. While Cas9 cleavage efficiency optimization and essential domain targeting have been developed as independent guide design rationales, no library has yet combined these into a single cohesive strategy to knock out gene function. Here, in a massive reanalysis of CRISPR tiling data using the most comprehensive feature database assembled, we determine which features of guides and their targets best predict activity and how to best combine them into a single guide design algorithm. We present the ProteIN ConsERvation (PINCER) genome-wide CRISPR library, which for the first time combines enzymatic efficiency optimization with conserved length protein region targeting, and also incorporates domains, coding sequence position, U6 termination (TTT), restriction sites, polymorphisms and specificity. Finally, we demonstrate superior performance of the PINCER library compared to alternative genome-wide CRISPR libraries in head-to-head validation. PINCER is available for individual gene knockout and genome-wide screening for both the human and mouse genomes.

Molecular characteristics associated with ferroptosis in hepatocellular carcinoma progression

The aim of this study was to investigate the genes associated with ferroptosis and the progression of hepatocellular carcinoma (HCC). The RNA sequencing data of erastin-induced ferroptosis in HCC cells were downloaded from the Sequence Read Archive database with accession number SRP119173. The microarray dataset GSE89377 of HCC progression was downloaded from the Gene Expression Omnibus database. The ferroptosis-related genes were screened by differential analysis and HCC progression-related genes were screened by cluster analysis using Mfuzz. Then, the genes associated with ferroptosis and HCC progression were screened by Venn analysis, followed by functional enrichment, protein-protein interaction (PPI) analysis, and transcription factor (TF) prediction. Finally, survival analysis was performed using data from the Cancer Genome Atlas database. A total of 33 upregulated and 52 downregulated genes associated with HCC progression and ferroptosis were obtained, and these genes were significantly involved in the negative regulation of ERK1 and ERK2 cascades; the NAD biosynthetic process; alanine, aspartate, and glutamate metabolism; and other pathways. The PPI network contained 52 genes and 78 interactions, of which, cell division cycle 20 (CDC20) and heat shock protein family B (small) member 1 (HSPB1) were hub genes found in higher degrees. Among the 85 genes associated with HCC progression and ferroptosis, two TFs (activating TF 3 (ATF3) and HLF) were predicted, with HSPB1 targeted by ATF3. In addition, 26 genes that were found to be significantly correlated with the overall survival of HCC patients were screened, including CDC20 and thyroid hormone receptor interactor 13. Several genes associated with HCC progression and ferroptosis were screened based on a comprehensive bioinformatics analysis. These genes played roles in HCC progression and ferroptosis via the negative regulation of the ERK1 and ERK2 cascades; the NAD biosynthetic process; and alanine, aspartate, and glutamate metabolism. ATF3 and HSPB1 played important roles in HCC progression and ferroptosis, with HSPB1 possibly regulated by ATF3.

Macrophage development and activation involve coordinated intron retention in key inflammatory regulators

Monocytes and macrophages are essential components of the innate immune system. Herein, we report that intron retention (IR) plays an important role in the development and function of these cells. Using Illumina mRNA sequencing, Nanopore direct cDNA sequencing and proteomics analysis, we identify IR events that affect the expression of key genes/proteins involved in macrophage development and function. We demonstrate that decreased IR in nuclear-detained mRNA is coupled with increased expression of genes encoding regulators of macrophage transcription, phagocytosis and inflammatory signalling, including ID2, IRF7, ENG and LAT. We further show that this dynamic IR program persists during the polarisation of resting macrophages into activated macrophages. In the presence of proinflammatory stimuli, intron-retaining CXCL2 and NFKBIZ transcripts are rapidly spliced, enabling timely expression of these key inflammatory regulators by macrophages. Our study provides novel insights into the molecular factors controlling vital regulators of the innate immune response.

Hyperconserved Elements in Human 5'UTRs Shape Essential Post-transcriptional Regulatory Networks

Post-transcriptional regulation (PTR) of gene expression is a powerful determinant of cellular phenotypes. The 5′ and 3′ untranslated regions of the mRNA (UTRs) mediate this role through sequence and secondary structure elements bound by RNA-binding proteins (RBPs) and non-coding RNAs. While functional regions in the 3′UTRs have been extensively studied, the 5′UTRs are still relatively uncharacterized. To fill this gap, we used a computational approach exploiting phylogenetic conservation to identify hyperconserved elements in human 5′UTRs (5′HCEs). Our assumption was that 5′HCEs would represent evolutionarily stable and hence important PTR sites. We identified over 5000 5′HCEs occurring in 10% of human protein-coding genes. These sequence elements are rather short and mostly found in narrowly-spaced clusters. 5′HCEs-containing genes are enriched in essential cellular functions and include 20% of all homeotic genes. Homeotic genes are essential transcriptional regulators, driving body plan and neuromuscular development. However, the role of PTR in their expression is mostly unknown. By integrating computational and experimental approaches we identified RBMX as the initiator RBP of a post-transcriptional cascade regulating many homeotic genes. This work thus establishes 5′HCEs as mediators of essential post-transcriptional regulatory networks.

Discovery and Functional Prediction of Long Non-Coding RNAs Common to Ischemic Stroke and Myocardial Infarction

Objective

Ischemic stroke and myocardial infarction are 2 of the leading causes of mortality. Both conditions are caused by arterial occlusion, resulting in ischemic necrosis of the cells in the cortex and heart. Long non-coding RNAs (lncRNAs) are a group of non-coding RNAs longer than 200 nucleotides without protein-coding potential. Thousands of lncRNAs have been identified but their involvement in ischemic stroke and myocardial infarction has not been studied extensively. Therefore, this study aimed to identify the role of lncRNAs, particularly those that are commonly altered in these two ischemic injuries.

Methods

We combined diverse RNA sequencing data obtained from public databases and performed extensive bioinformatics analyses to determine reliable lncRNAs commonly identified from these datasets. Using sequence analysis, we also detected the lncRNAs that may act as microRNA (miRNA) regulators.

Results

We found several altered lncRNAs that were common in ischemic stroke and myocardial infarction models. Some of these lncRNAs, including zinc finger NFX1-type containing 1 antisense RNA 1 and small nucleolar RNA host gene 1, were previously reported to be involved in the pathogenesis of each of these models. Interestingly, several lncRNAs had binding sites for miRNAs that were previously reported to be involved in the hypoxic response, suggesting the possible role of these lncRNAs as regulators in ischemic responses.

Conclusion

The lncRNAs identified in this study will be useful in determining the regulatory networks in ischemic stroke and myocardial infarction and in identifying potential specific markers for each of these ischemic diseases.

Impact of diet and genes on murine autoimmune pancreatitis

The impact of environmental factors, such as diet, and the genetic basis of autoimmune pancreatitis (AIP) are largely unknown. Here, we used an experimental murine AIP model to identify the contribution of diet to AIP development, as well as to fine-map AIP-associated genes in outbred mice prone to develop the disease. For this purpose, we fed mice of an autoimmune-prone intercross line (AIL) three different diets (control, calorie-reduced and western diet) for 6 months, at which point the mice were genotyped and phenotyped for AIP. Overall, 269 out of 734 mice (36.6%) developed AIP with signs of parenchymal destruction, equally affecting mice of both sexes. AIP prevalence and severity were reduced by approximately 50% in mice held under caloric restriction compared to those fed control or western diet. We identified a quantitative trait locus (QTL) on chromosome 4 to be associated with AIP, which is located within a previously reported QTL. This association does not change when considering diet or sex as an additional variable for the mapping. Using whole-genome sequences of the AIL founder strains, we resolved this QTL to a single candidate gene, namely Map3k7. Expression of Map3k7 was largely restricted to islet cells as well as lymphocytes found in the exocrine pancreas of mice with AIP. Our studies suggest a major impact of diet on AIP. Furthermore, we identify Map3k7 as a novel susceptibility gene for experimental AIP. Both findings warrant clinical translation.

Insights into variation in meiosis from 31,228 human sperm genomes

Meiosis, although essential for reproduction, is also variable and error-prone: rates of chromosome crossover vary among gametes, between the sexes, and among humans of the same sex, and chromosome missegregation leads to abnormal chromosome numbers (aneuploidy)1-8. To study diverse meiotic outcomes and how they covary across chromosomes, gametes and humans, we developed Sperm-seq, a way of simultaneously analysing the genomes of thousands of individual sperm. Here we analyse the genomes of 31,228 human gametes from 20 sperm donors, identifying 813,122 crossovers and 787 aneuploid chromosomes. Sperm donors had aneuploidy rates ranging from 0.01 to 0.05 aneuploidies per gamete; crossovers partially protected chromosomes from nondisjunction at the meiosis I cell division. Some chromosomes and donors underwent more-frequent nondisjunction during meiosis I, and others showed more meiosis II segregation failures. Sperm genomes also manifested many genomic anomalies that could not be explained by simple nondisjunction. Diverse recombination phenotypes-from crossover rates to crossover location and separation, a measure of crossover interference-covaried strongly across individuals and cells. Our results can be incorporated with earlier observations into a unified model in which a core mechanism, the variable physical compaction of meiotic chromosomes, generates interindividual and cell-to-cell variation in diverse meiotic phenotypes.

Native American gene flow into Polynesia predating Easter Island settlement

The possibility of voyaging contact between prehistoric Polynesian and Native American populations has long intrigued researchers. Proponents have pointed to the existence of New World crops, such as the sweet potato and bottle gourd, in the Polynesian archaeological record, but nowhere else outside the pre-Columbian Americas1-6, while critics have argued that these botanical dispersals need not have been human mediated7. The Norwegian explorer Thor Heyerdahl controversially suggested that prehistoric South American populations had an important role in the settlement of east Polynesia and particularly of Easter Island (Rapa Nui)2. Several limited molecular genetic studies have reached opposing conclusions, and the possibility continues to be as hotly contested today as it was when first suggested8-12. Here we analyse genome-wide variation in individuals from islands across Polynesia for signs of Native American admixture, analysing 807 individuals from 17 island populations and 15 Pacific coast Native American groups. We find conclusive evidence for prehistoric contact of Polynesian individuals with Native American individuals (around AD 1200) contemporaneous with the settlement of remote Oceania13-15. Our analyses suggest strongly that a single contact event occurred in eastern Polynesia, before the settlement of Rapa Nui, between Polynesian individuals and a Native American group most closely related to the indigenous inhabitants of present-day Colombia.

Delving into Eukaryotic Origins of Replication Using DNA Structural Features

DNA replication in eukaryotes is an intricate process, which is precisely synchronized by a set of regulatory proteins, and the replication fork emanates from discrete sites on chromatin called origins of replication (Oris). These spots are considered as the gateway to chromosomal replication and are stereotyped by sequence motifs. The cognate sequences are noticeable in a small group of entire origin regions or totally absent across different metazoans. Alternatively, the use of DNA secondary structural features can provide additional information compared to the primary sequence. In this article, we report the trends in DNA sequence-based structural properties of origin sequences in nine eukaryotic systems representing different families of life. Biologically relevant DNA secondary structural properties, namely, stability, propeller twist, flexibility, and minor groove shape were studied in the sequences flanking replication start sites. Results indicate that Oris in yeasts show lower stability, more rigidity, and narrow minor groove preferences compared to genomic sequences surrounding them. Yeast Oris also show preference for A-tracts and the promoter element TATA box in the vicinity of replication start sites. On the contrary, Drosophila melanogaster, humans, and Arabidopsis thaliana do not have such features in their Oris, and instead, they show high preponderance of G-rich sequence motifs such as putative G-quadruplexes or i-motifs and CpG islands. Our extensive study applies the DNA structural feature computation to delve into origins of replication across organisms ranging from yeasts to mammals and including a plant. Insights from this study would be significant in understanding origin architecture and help in designing new algorithms for predicting DNA trans-acting factor recognition events.

Proteasome subunit PSMC3 variants cause neurosensory syndrome combining deafness and cataract due to proteotoxic stress

The ubiquitin–proteasome system degrades ubiquitin‐modified proteins to maintain protein homeostasis and to control signalling. Whole‐genome sequencing of patients with severe deafness and early‐onset cataracts as part of a neurological, sensorial and cutaneous novel syndrome identified a unique deep intronic homozygous variant in the PSMC3 gene, encoding the proteasome ATPase subunit Rpt5, which lead to the transcription of a cryptic exon. The proteasome content and activity in patient's fibroblasts was however unaffected. Nevertheless, patient's cells exhibited impaired protein homeostasis characterized by accumulation of ubiquitinated proteins suggesting severe proteotoxic stress. Indeed, the TCF11/Nrf1 transcriptional pathway allowing proteasome recovery after proteasome inhibition is permanently activated in the patient's fibroblasts. Upon chemical proteasome inhibition, this pathway was however impaired in patient's cells, which were unable to compensate for proteotoxic stress although a higher proteasome content and activity. Zebrafish modelling for knockout in PSMC3 remarkably reproduced the human phenotype with inner ear development anomalies as well as cataracts, suggesting that Rpt5 plays a major role in inner ear, lens and central nervous system development.