Ensembl variation resources

Full methylation of H3K27 by PRC2 is dispensable for initial embryoid body formation but required to maintain differentiated cell identity

Polycomb repressive complex 2 (PRC2) catalyzes methylation of histone H3 on lysine 27 and is required for normal development of complex eukaryotes. The nature of that requirement is not clear. H3K27me3 is associated with repressed genes, but the modification is not sufficient to induce repression and, in some instances, is not required. We blocked full methylation of H3K27 with both a small molecule inhibitor, GSK343, and by introducing a point mutation into EZH2, the catalytic subunit of PRC2, in the mouse CJ7 cell line. Cells with substantively decreased H3K27 methylation differentiate into embryoid bodies, which contrasts with EZH2 null cells. PRC2 targets had varied requirements for H3K27me3, with a subset that maintained normal levels of repression in the absence of methylation. The primary cellular phenotype of blocked H3K27 methylation was an inability of altered cells to maintain a differentiated state when challenged. This phenotype was determined by H3K27 methylation in embryonic stem cells through the first 4 days of differentiation. Full H3K27 methylation therefore was not necessary for formation of differentiated cell states during embryoid body formation but was required to maintain a stable differentiated state.

Structure analysis of deleterious nsSNPs in human PALB2 protein for functional inference

Partner and Localizer of BRCA2 or PALB2 is a typical tumor suppressor protein, that responds to DNA double stranded breaks through homologous recombination repair. Heterozygous mutations in PALB2 are known to contribute to the susceptibility of breast and ovarian cancer. However, there is no comprehensive study characterizing the structural and functional impacts of SNPs located in the PALB2 gene. Therefore, it is of interest to document a comprehensive analysis of coding and non-coding SNPs located at the PALB2 loci using in silico tools. The data for 1455 non-synonymous SNPs (nsSNPs) located in the PALB2 loci were retrieved from the dbSNP database. Comprehensive characterization of the SNPs using a combination of in silico tools such as SIFT, PROVEAN, PolyPhen, PANTHER, PhD-SNP, Pmut, MutPred 2.0 and SNAP-2, identified 28 functionally important SNPs. Among these, 16 nsSNPs were further selected for structural analysis using conservation profile and protein stability. The most deleterious nsSNPs were documented within the WD40 domain of PALB2. A general outline of the structural consequences of each variant was developed using the HOPE project data. These 16 mutant structures were further modelled using SWISS Model and three most damaging mutant models (rs78179744, rs180177123 and rs45525135) were identified. The non-coding SNPs in the 3' UTR region of the PALB2 gene were analyzed for altered miRNA target sites. The comprehensive characterization of the coding and non-coding SNPs in the PALB2 locus has provided a list of damaging SNPs with potential disease association. Further validation through genetic association study will reveal their clinical significance.

Genome-wide RNAi screen for regulators of UPRmt in Caenorhabditis elegans mutants with defects in mitochondrial fusion

Mitochondrial dynamics plays an important role in mitochondrial quality control and the adaptation of metabolic activity in response to environmental changes. The disruption of mitochondrial dynamics has detrimental consequences for mitochondrial and cellular homeostasis and leads to the activation of the mitochondrial unfolded protein response (UPRmt), a quality control mechanism that adjusts cellular metabolism and restores homeostasis. To identify genes involved in the induction of UPRmt in response to a block in mitochondrial fusion, we performed a genome-wide RNAi screen in Caenorhabditis elegans mutants lacking the gene fzo-1, which encodes the ortholog of mammalian Mitofusin, and identified 299 suppressors and 86 enhancers. Approximately 90% of these 385 genes are conserved in humans, and one third of the conserved genes have been implicated in human disease. Furthermore, many have roles in developmental processes, which suggests that mitochondrial function and the response to stress are defined during development and maintained throughout life. Our dataset primarily contains mitochondrial enhancers and non-mitochondrial suppressors of UPRmt, indicating that the maintenance of mitochondrial homeostasis has evolved as a critical cellular function, which, when disrupted, can be compensated for by many different cellular processes. Analysis of the subsets 'non-mitochondrial enhancers' and 'mitochondrial suppressors' suggests that organellar contact sites, especially between the ER and mitochondria, are of importance for mitochondrial homeostasis. In addition, we identified several genes involved in IP3 signaling that modulate UPRmt in fzo-1 mutants and found a potential link between pre-mRNA splicing and UPRmt activation.

In silico candidate variant and gene identification using inbred mouse strains

Mice are the most widely used animal model to study genotype to phenotype relationships. Inbred mice are genetically identical, which eliminates genetic heterogeneity and makes them particularly useful for genetic studies. Many different strains have been bred over decades and a vast amount of phenotypic data has been generated. In addition, recently whole genome sequencing-based genome-wide genotype data for many widely used inbred strains has been released. Here, we present an approach for in silico fine-mapping that uses genotypic data of 37 inbred mouse strains together with phenotypic data provided by the user to propose candidate variants and genes for the phenotype under study. Public genome-wide genotype data covering more than 74 million variant sites is queried efficiently in real-time to provide those variants that are compatible with the observed phenotype differences between strains. Variants can be filtered by molecular consequences and by corresponding molecular impact. Candidate gene lists can be generated from variant lists on the fly. Fine-mapping together with annotation or filtering of results is provided in a Bioconductor package called MouseFM. In order to characterize candidate variant lists under various settings, MouseFM was applied to two expression data sets across 20 inbred mouse strains, one from neutrophils and one from CD4⁺ T cells. Fine-mapping was assessed for about 10,000 genes, respectively, and identified candidate variants and haplotypes for many expression quantitative trait loci (eQTLs) reported previously based on these data. For albinism, MouseFM reports only one variant allele of moderate or high molecular impact that only albino mice share: a missense variant in the Tyr gene, reported previously to be causal for this phenotype. Performing in silico fine-mapping for interfrontal bone formation in mice using four strains with and five strains without interfrontal bone results in 12 genes. Of these, three are related to skull shaping abnormality. Finally performing fine-mapping for dystrophic cardiac calcification by comparing 9 strains showing the phenotype with eight strains lacking it, we identify only one moderate impact variant in the known causal gene Abcc6. In summary, this illustrates the benefit of using MouseFM for candidate variant and gene identification.

Coping with brain amyloid: genetic heterogeneity and cognitive resilience to Alzheimer's pathophysiology

Although abnormal accumulation of amyloid in the brain is an early biomarker of Alzheimer's disease (AD), wide variation in cognitive trajectories during life can be seen in the setting of brain amyloidosis, ranging from maintenance of normal function to progression to dementia. It is widely presumed that cognitive resilience (i.e., coping) to amyloidosis may be influenced by environmental, lifestyle, and inherited factors, but relatively little in specifics is known about this architecture. Here, we leveraged multimodal longitudinal data from a large, population-based sample of older adults to discover genetic factors associated with differential cognitive resilience to brain amyloidosis determined by positron emission tomography (PET). Among amyloid-PET positive older adults, the AD risk allele APOE ɛ4 was associated with worse longitudinal memory trajectories as expected, and was thus covaried in the main analyses. Through a genome-wide association study (GWAS), we uncovered a novel association with cognitive resilience on chromosome 8 at the MTMR7/CNOT7/ZDHHC2/VPS37A locus (p = 4.66 × 10-8, β = 0.23), and demonstrated replication in an independent cohort. Post-hoc analyses confirmed this association as specific to the setting of elevated amyloid burden and not explained by differences in tau deposition or cerebrovascular disease. Complementary gene-based analyses and publically available functional data suggested that the causative variant at this locus may tag CNOT7 (CCR4-NOT Transcription Complex Subunit 7), a gene linked to synaptic plasticity and hippocampal-dependent learning and memory. Pathways related to cell adhesion and immune system activation displayed enrichment of association in the GWAS. Our findings, resulting from a unique study design, support the hypothesis that genetic heterogeneity is one of the factors that explains differential cognitive resilience to brain amyloidosis. Further characterization of the underlying biological mechanisms influencing cognitive resilience may facilitate improved prognostic counseling, therapeutic application, and trial enrollment in AD.

Integrated network analysis reveals new genes suggesting COVID-19 chronic effects and treatment

The COVID-19 disease led to an unprecedented health emergency, still ongoing worldwide. Given the lack of a vaccine or a clear therapeutic strategy to counteract the infection as well as its secondary effects, there is currently a pressing need to generate new insights into the SARS-CoV-2 induced host response. Biomedical data can help to investigate new aspects of the COVID-19 pathogenesis, but source heterogeneity represents a major drawback and limitation. In this work, we applied data integration methods to develop a Unified Knowledge Space (UKS) and used it to identify a new set of genes associated with SARS-CoV-2 host response, both in vitro and in vivo. Functional analysis of these genes reveals possible long-term systemic effects of the infection, such as vascular remodelling and fibrosis. Finally, we identified a set of potentially relevant drugs targeting proteins involved in multiple steps of the host response to the virus.

QTL Analysis of Adult Plant Resistance to Stripe Rust in a Winter Wheat Recombinant Inbred Population

Stripe rust, caused by the fungus Puccinia striiformis f. sp. tritici, is a worldwide disease of wheat that causes devastating crop losses. Resistant cultivars have been developed over the last 40 years that have significantly reduced the economic impact of the disease on growers, but in heavy infection years it is mostly controlled through the intensive application of fungicides. The Pacific Northwest of the United States has an ideal climate for stripe rust and has one of the most diverse race compositions in the country. This has resulted in many waves of epidemics that have overcome most of the resistance genes traditionally used in elite germplasm. The best way to prevent high yield losses, reduce production costs to growers, and reduce the heavy application of fungicides is to pyramid multiple stripe rust resistance genes into new cultivars. Using genotyping-by-sequencing, we identified 4662 high quality variant positions in a recombinant inbred line population of 196 individuals derived from a cross between Skiles, a highly resistant winter wheat cultivar, and Goetze, a moderately to highly susceptible winter wheat cultivar, both developed at Oregon State University. A subsequent genome wide association study identified two quantitative trait loci (QTL) on chromosomes 3B and 3D within the predicted locations of stripe rust resistance genes. Resistance QTL, when combined together, conferred high levels of stripe rust resistance above the level of Skiles in some locations, indicating that these QTL would be important additions to future breeding efforts of Pacific Northwest winter wheat cultivars.

RNA-seq of human T cells after hematopoietic stem cell transplantation identifies Linc00402 as a regulator of T cell alloimmunity

Mechanisms governing allogeneic T cell responses after solid organ and allogeneic hematopoietic stem cell transplantation (HSCT) are incompletely understood. To identify lncRNAs that regulate human donor T cells after clinical HSCT, we performed RNA sequencing on T cells from healthy individuals and donor T cells from three different groups of HSCT recipients that differed in their degree of major histocompatibility complex (MHC) mismatch. We found that lncRNA differential expression was greatest in T cells after MHC-mismatched HSCT relative to T cells after either MHC-matched or autologous HSCT. Differential expression was validated in an independent patient cohort and in mixed lymphocyte reactions using ex vivo healthy human T cells. We identified Linc00402, an uncharacterized lncRNA, among the lncRNAs differentially expressed between the mismatched unrelated and matched unrelated donor T cells. We found that Linc00402 was conserved and exhibited an 88-fold increase in human T cells relative to all other samples in the FANTOM5 database. Linc00402 was also increased in donor T cells from patients who underwent allogeneic cardiac transplantation and in murine T cells. Linc00402 was reduced in patients who subsequently developed acute graft-versus-host disease. Linc00402 enhanced the activity of ERK1 and ERK2, increased FOS nuclear accumulation, and augmented expression of interleukin-2 and Egr-1 after T cell receptor engagement. Functionally, Linc00402 augmented the T cell proliferative response to an allogeneic stimulus but not to a nominal ovalbumin peptide antigen or polyclonal anti-CD3/CD28 stimulus. Thus, our studies identified Linc00402 as a regulator of allogeneic T cell function.

Omics data integration identifies ELOVL7 and MMD gene regions as novel loci for adalimumab response in patients with Crohn's disease

Response to anti-TNF therapy is of pivotal importance in patients with Crohn’s disease (CD). Here we integrated our and previously reported PBMC derived transcriptomic and genomic data for identification of biomarkers for discrimination between responders and non-responders to anti-TNF therapy. CD patients, who were naïve with respect to the treatment with biologicals, were enrolled in the study. DNA and RNA were extracted from peripheral blood mononuclear cells. RNA-seq was performed using BGISEQ-500. Genotyping was performed using Infinium Global Screening Array. Association regressions were carried out with 12 week response to adalimumab as an outcome variable. RNA-seq analysis confirmed 7 out of 65 previously suggested genes involved in anti-TNF response. Subsequently, analysis of single nucleotide variants in regions of confirmed genes identified 5 variants near MMD and two in ELOVL7 intronic regions associated with treatment response to anti-TNF. Functional analysis has shown that rs1465352, rs4422035 and rs78620886 are listed at H3K9ac_Pro histone modification epigenetic mark. The present study confirmed MMD and ELOVL7 involvement in anti-TNF response and revealed that the regulation of MMD and ELOVL7 gene regions in ADA response may be a part of a complex interplay extending from genetic to epigenetic and to transcriptomic level.

The myosin regulatory light chain Myl5 localizes to mitotic spindle poles and is required for proper cell division

Myosins are ATP-dependent actin-based molecular motors critical for diverse cellular processes like intracellular trafficking, cell motility, and cell invasion. During cell division, myosin MYO10 is important for proper mitotic spindle assembly, the anchoring of the spindle to the cortex, and positioning of the spindle to the cell mid-plane. However, myosins are regulated by myosin regulatory light chains (RLCs), and whether RLCs are important for cell division has remained unexplored. Here, we have determined that the previously uncharacterized myosin RLC Myl5 associates with the mitotic spindle and is required for cell division. We show that Myl5 localizes to the leading edge and filopodia during interphase and to mitotic spindle poles and spindle microtubules during early mitosis. Importantly, depletion of Myl5 led to defects in mitotic spindle assembly, chromosome congression, and chromosome segregation and to a slower transition through mitosis. Furthermore, Myl5 bound to MYO10 in vitro and co-localized with MYO10 at the spindle poles. These results suggest that Myl5 is important for cell division and that it may be performing its function through MYO10.

Whole-Genome Sequence Data Suggest Environmental Adaptation of Ethiopian Sheep Populations

Great progress has been made over recent years in the identification of selection signatures in the genomes of livestock species. This work has primarily been carried out in commercial breeds for which the dominant selection pressures are associated with artificial selection. As agriculture and food security are likely to be strongly affected by climate change, a better understanding of environment-imposed selection on agricultural species is warranted. Ethiopia is an ideal setting to investigate environmental adaptation in livestock due to its wide variation in geo-climatic characteristics and the extensive genetic and phenotypic variation of its livestock. Here, we identified over three million single nucleotide variants across 12 Ethiopian sheep populations and applied landscape genomics approaches to investigate the association between these variants and environmental variables. Our results suggest that environmental adaptation for precipitation-related variables is stronger than that related to altitude or temperature, consistent with large-scale meta-analyses of selection pressure across species. The set of genes showing association with environmental variables was enriched for genes highly expressed in human blood and nerve tissues. There was also evidence of enrichment for genes associated with high-altitude adaptation although no strong association was identified with hypoxia-inducible-factor (HIF) genes. One of the strongest altitude-related signals was for a collagen gene, consistent with previous studies of high-altitude adaptation. Several altitude-associated genes also showed evidence of adaptation with temperature, suggesting a relationship between responses to these environmental factors. These results provide a foundation to investigate further the effects of climatic variables on small ruminant populations.

A Workflow for Selection of Single Nucleotide Polymorphic Markers for Studying of Genetics of Ischemic Stroke Outcomes

In this paper we propose a workflow for studying the genetic architecture of ischemic stroke outcomes. It develops further the candidate gene approach. The workflow is based on the animal model of brain ischemia, comparative genomics, human genomic variations, and algorithms of selection of tagging single nucleotide polymorphisms (tagSNPs) in genes which expression was changed after ischemic stroke. The workflow starts from a set of rat genes that changed their expression in response to brain ischemia and results in a set of tagSNPs, which represent other SNPs in the human genes analyzed and influenced on their expression as well.

Polymorphisms in canine immunoglobulin heavy chain gene cluster: a double-edged sword for diabetes mellitus in the dog

Insulin deficiency diabetes (IDD) in dogs is an endocrine disease similar to human type 1 diabetes. There are breeds more commonly affected, such as Yorkshire Terrier and Samoyed, suggesting an underlying genetic component. However, the genetic basis for canine diabetes mellitus (DM) is not fully established. We conducted both whole-genome scans for selection signatures and GWASs to compare the genomes of 136 dogs belonging to 29 breeds previously described at low or high risk for developing DM. Candidate variants were tested in dogs with a diagnosis of IDD and controls attending the Complutense Veterinary Teaching Hospital. The only genomic region under selection (CFA8:72 700 000-74 600 000; CanFam3.1) retrieved by our analyses is included in the immunoglobulin heavy chain gene cluster, which has already been related to human human type 1 diabetes susceptibility. This region contains two non-synonymous variants, rs852072969 and rs851728071, showing significant associations with high or low risk for IDD, respectively. The first variant, rs852072969, alters a protein poorly characterised in the dog. In contrast, rs851728071 was predicted to block the synthesis of an immunoglobulin variable (V) domain in breeds at low risk for DM. Although a large and diverse V gene repertoire is thought to offer a fitness advantage, we suggest that rs851728071 prevents the formation of an auto-reactive immunoglobulin V domain probably involved in the pathophysiology of IDD and, thus, decreases the risk for the disease. These results should be interpreted with caution until the functional roles of the proposed variants have been proved in larger studies.

Trans-ethnic meta-analysis identifies new loci associated with longitudinal blood pressure traits

Genome-wide association studies (GWAS) have identified thousands of genetic loci associated with cross-sectional blood pressure (BP) traits; however, GWAS based on longitudinal BP have been underexplored. We performed ethnic-specific and trans-ethnic GWAS meta-analysis using longitudinal and cross-sectional BP data of 33,720 individuals from five cohorts in the US and one in Brazil. In addition to identifying several known loci, we identified thirteen novel loci with nine based on longitudinal and four on cross-sectional BP traits. Most of the novel loci were ethnic- or study-specific, with the majority identified in African Americans (AA). Four of these discoveries showed additional evidence of association in independent datasets, including an intergenic variant (rs4060030, p = 7.3 × 10^–9) with reported regulatory function. We observed a high correlation between the meta-analysis results for baseline and longitudinal average BP (rho = 0.48). BP trajectory results were more correlated with those of average BP (rho = 0.35) than baseline BP(rho = 0.18). Heritability estimates trended higher for longitudinal traits than for cross-sectional traits, providing evidence for different genetic architectures. Furthermore, the longitudinal data identified up to 20% more BP known associations than did cross-sectional data. Our analyses of longitudinal BP data in diverse ethnic groups identified novel BP loci associated with BP trajectory, indicating a need for further longitudinal GWAS on BP and other age-related traits.

CONQUER: an interactive toolbox to understand functional consequences of GWAS hits

Numerous large genome-wide association studies have been performed to understand the influence of genetics on traits. Many identified risk loci are in non-coding and intergenic regions, which complicates understanding how genes and their downstream pathways are influenced. An integrative data approach is required to understand the mechanism and consequences of identified risk loci. Here, we developed the R-package CONQUER. Data for SNPs of interest are acquired from static- and dynamic repositories (build GRCh38/hg38), including GTExPortal, Epigenomics Project, 4D genome database and genome browsers. All visualizations are fully interactive so that the user can immediately access the underlying data. CONQUER is a user-friendly tool to perform an integrative approach on multiple SNPs where risk loci are not seen as individual risk factors but rather as a network of risk factors.

Mendelian randomization for studying the effects of perturbing drug targets

Drugs whose targets have genetic evidence to support efficacy and safety are more likely to be approved after clinical development. In this paper, we provide an overview of how natural sequence variation in the genes that encode drug targets can be used in Mendelian randomization analyses to offer insight into mechanism-based efficacy and adverse effects. Large databases of summary level genetic association data are increasingly available and can be leveraged to identify and validate variants that serve as proxies for drug target perturbation. As with all empirical research, Mendelian randomization has limitations including genetic confounding, its consideration of lifelong effects, and issues related to heterogeneity across different tissues and populations. When appropriately applied, Mendelian randomization provides a useful empirical framework for using population level data to improve the success rates of the drug development pipeline.

Describing variability in pig genes involved in coronavirus infections for a One Health perspective in conservation of animal genetic resources

Coronaviruses silently circulate in human and animal populations, causing mild to severe diseases. Therefore, livestock are important components of a “One Health” perspective aimed to control these viral infections. However, at present there is no example that considers pig genetic resources in this context. In this study, we investigated the variability of four genes (ACE2, ANPEP and DPP4 encoding for host receptors of the viral spike proteins and TMPRSS2 encoding for a host proteinase) in 23 European (19 autochthonous and three commercial breeds and one wild boar population) and two Asian Sus scrofa populations. A total of 2229 variants were identified in the four candidate genes: 26% of them were not previously described; 29 variants affected the protein sequence and might potentially interact with the infection mechanisms. The results coming from this work are a first step towards a “One Health” perspective that should consider conservation programs of pig genetic resources with twofold objectives: (i) genetic resources could be reservoirs of host gene variability useful to design selection programs to increase resistance to coronaviruses; (ii) the described variability in genes involved in coronavirus infections across many different pig populations might be part of a risk assessment including pig genetic resources.

Cryo-EM reveals structural breaks in a patient-derived amyloid fibril from systemic AL amyloidosis

Systemic AL amyloidosis is a debilitating and potentially fatal disease that arises from the misfolding and fibrillation of immunoglobulin light chains (LCs). The disease is patient-specific with essentially each patient possessing a unique LC sequence. In this study, we present two ex vivo fibril structures of a λ3 LC. The fibrils were extracted from the explanted heart of a patient (FOR005) and consist of 115-residue fibril proteins, mainly from the LC variable domain. The fibril structures imply that a 180° rotation around the disulfide bond and a major unfolding step are necessary for fibrils to form. The two fibril structures show highly similar fibril protein folds, differing in only a 12-residue segment. Remarkably, the two structures do not represent separate fibril morphologies, as they can co-exist at different z-axial positions within the same fibril. Our data imply the presence of structural breaks at the interface of the two structural forms.

Manually curated transcriptomics data collection for toxicogenomic assessment of engineered nanomaterials

Toxicogenomics (TGx) approaches are increasingly applied to gain insight into the possible toxicity mechanisms of engineered nanomaterials (ENMs). Omics data can be valuable to elucidate the mechanism of action of chemicals and to develop predictive models in toxicology. While vast amounts of transcriptomics data from ENM exposures have already been accumulated, a unified, easily accessible and reusable collection of transcriptomics data for ENMs is currently lacking. In an attempt to improve the FAIRness of already existing transcriptomics data for ENMs, we curated a collection of homogenized transcriptomics data from human, mouse and rat ENM exposures in vitro and in vivo including the physicochemical characteristics of the ENMs used in each study.

Mendelian randomization for studying the effects of perturbing drug targets

Drugs whose targets have genetic evidence to support efficacy and safety are more likely to be approved after clinical development. In this paper, we provide an overview of how natural sequence variation in the genes that encode drug targets can be used in Mendelian randomization analyses to offer insight into mechanism-based efficacy and adverse effects. Large databases of summary level genetic association data are increasingly available and can be leveraged to identify and validate variants that serve as proxies for drug target perturbation. As with all empirical research, Mendelian randomization has limitations including genetic confounding, its consideration of lifelong effects, and issues related to heterogeneity across different tissues and populations. When appropriately applied, Mendelian randomization provides a useful empirical framework for using population level data to improve the success rates of the drug development pipeline.

Chloroquine and Hydroxychloroquine Interact Differently with ACE2 Domains Reported to Bind with the Coronavirus Spike Protein: Mediation by ACE2 Polymorphism

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection inducing coronavirus disease 2019 (COVID-19) is still an ongoing challenge. To date, more than 95.4 million have been infected and more than two million deaths have been officially reported by the WHO. Angiotensin-converting enzyme (ACE) plays a key role in the disease pathogenesis. In this computational study, seventeen coding variants were found to be important for ACE2 binding with the coronavirus spike protein. The frequencies of these allele variants range from 3.88 × 10^-3 to 5.47 × 10^-6 for rs4646116 (K26R) and rs1238146879 (P426A), respectively. Chloroquine (CQ) and its metabolite hydroxychloroquine (HCQ) are mainly used to prevent and treat malaria and rheumatic diseases. They are also used in several countries to treat SARS-CoV-2 infection inducing COVID-19. Both CQ and HCQ were found to interact differently with the various ACE2 domains reported to bind with coronavirus spike protein. A molecular docking approach revealed that intermolecular interactions of both CQ and HCQ exhibited mediation by ACE2 polymorphism. Further explorations of the relationship and the interactions between ACE2 polymorphism and CQ/HCQ would certainly help to better understand the COVID-19 management strategies, particularly their use in the absence of specific vaccines or drugs.

Ensembl 2021

The Ensembl project (https://www.ensembl.org) annotates genomes and disseminates genomic data for vertebrate species. We create detailed and comprehensive annotation of gene structures, regulatory elements and variants, and enable comparative genomics by inferring the evolutionary history of genes and genomes. Our integrated genomic data are made available in a variety of ways, including genome browsers, search interfaces, specialist tools such as the Ensembl Variant Effect Predictor, download files and programmatic interfaces. Here, we present recent Ensembl developments including two new website portals. Ensembl Rapid Release (http://rapid.ensembl.org) is designed to provide core tools and services for genomes as soon as possible and has been deployed to support large biodiversity sequencing projects. Our SARS-CoV-2 genome browser (https://covid-19.ensembl.org) integrates our own annotation with publicly available genomic data from numerous sources to facilitate the use of genomics in the international scientific response to the COVID-19 pandemic. We also report on other updates to our annotation resources, tools and services. All Ensembl data and software are freely available without restriction.

Hypoxia-inducible miR-196a modulates glioblastoma cell proliferation and migration through complex regulation of NRAS

BACKGROUND

Glioblastoma (GBM) is the most common and aggressive malignant brain tumor in humans. Hypoxia has been correlated with the aggressive form of glial tumors, poor prognosis, recurrence and resistance to various therapies. MicroRNAs (miRNAs) have emerged as critical mediators of hypoxic responses and have shown great potential for cancer diagnostics and therapeutics. Here, we focus on the regulatory and functional characterization of miR-196a, a hypoxia-inducible miRNA, in GBM.

METHODS

Hypoxia/HIF regulation of miR-196a was assessed by RT-qPCR, promoter-luciferase and ChIP assays in GBM cell lines. miR-196a levels were analyzed in The Cancer Genome Atlas (TCGA)-GBM, Chinese Glioma Genome Atlas (CGGA) and Indian GBM patient cohorts. miR-target interactions were studied using RNA/protein quantification and 3'UTR luciferase assays. The effect of miR-196a overexpression/inhibition was assessed on cellular viability, migration and apoptosis under hypoxia and normoxia. Microarray-based gene expression profiling studies were performrd to study the effect of miR-196a on the GBM cellular transcriptome under hypoxia.

RESULTS

We identified miR-196a as a hypoxia-inducible and hypoxia-inducible factor (HIF)-regulated miRNA that plays an oncogenic role in GBM. miR-196a was found to be significantly up-regulated in TCGA-GBM, CGGA glioma as well as Indian GBM patient cohorts. miR-196a overexpression was found to induce cellular proliferation, migration, spheroid formation and colony formation and to inhibit apoptosis, while miR-196a inhibition using anti-miR-196a yielded opposite results, suggesting an oncogenic role of miR-196a in GBM. We further unveiled NRAS, AJAP1, TAOK1 and COL24A1 as direct targets of miR-196a. We also report a complex competitive regulation of oncogenic NRAS by miR-196a, miR-146a and let-7 in GBM. Analysis of microarray-based gene expression data obtained by miR-196a inhibition under hypoxia revealed a role of miR-196a in HIF, calcium adhesion, Wnt and cell adhesion pathways. Interestingly, miR-196a was found to positively regulate the expression of various genes involved in the induction or stabilization of HIFs and in maintenance of hypoxic conditions, thereby suggesting the existence of an indirect miR-196a/HIF positive feedback loop under hypoxia.

CONCLUSIONS

Overall, our work identifies a novel association between hypoxia/HIF signalling and miR-196a in GBM and suggests its therapeutic significance.

Suprabasin-A Review

Among the ~22,000 human genes, very few remain that have unknown functions. One such example is suprabasin (SBSN). Originally described as a component of the cornified envelope, the function of stratified epithelia-expressed SBSN is unknown. Both the lack of knowledge about the gene role under physiological conditions and the emerging link of SBSN to various human diseases, including cancer, attract research interest. The association of SBSN expression with poor prognosis of patients suffering from oesophageal carcinoma, glioblastoma multiforme, and myelodysplastic syndromes suggests that SBSN may play a role in human tumourigenesis. Three SBSN isoforms code for the secreted proteins with putative function as signalling molecules, yet with poorly described effects. In this first review about SBSN, we summarised the current knowledge accumulated since its original description, and we discuss the potential mechanisms and roles of SBSN in both physiology and pathology.

Application of Machine-Learning Methods to Recognize mitoBK Channels from Different Cell Types Based on the Experimental Patch-Clamp Results

(1) Background: In this work, we focus on the activity of large-conductance voltage- and Ca2+-activated potassium channels (BK) from the inner mitochondrial membrane (mitoBK). The characteristic electrophysiological features of the mitoBK channels are relatively high single-channel conductance (ca. 300 pS) and types of activating and deactivating stimuli. Nevertheless, depending on the isoformal composition of mitoBK channels in a given membrane patch and the type of auxiliary regulatory subunits (which can be co-assembled to the mitoBK channel protein) the characteristics of conformational dynamics of the channel protein can be altered. Consequently, the individual features of experimental series describing single-channel activity obtained by patch-clamp method can also vary. (2) Methods: Artificial intelligence approaches (deep learning) were used to classify the patch-clamp outputs of mitoBK activity from different cell types. (3) Results: Application of the K-nearest neighbors algorithm (KNN) and the autoencoder neural network allowed to perform the classification of the electrophysiological signals with a very good accuracy, which indicates that the conformational dynamics of the analyzed mitoBK channels from different cell types significantly differs. (4) Conclusion: We displayed the utility of machine-learning methodology in the research of ion channel gating, even in cases when the behavior of very similar microbiosystems is analyzed. A short excerpt from the patch-clamp recording can serve as a “fingerprint” used to recognize the mitoBK gating dynamics in the patches of membrane from different cell types.

Transcriptome-wide high-throughput mapping of protein-RNA occupancy profiles using POP-seq

Interaction between proteins and RNA is critical for post-transcriptional regulatory processes. Existing high throughput methods based on crosslinking of the protein–RNA complexes and poly-A pull down are reported to contribute to biases and are not readily amenable for identifying interaction sites on non poly-A RNAs. We present Protein Occupancy Profile-Sequencing (POP-seq), a phase separation based method in three versions, one of which does not require crosslinking, thus providing unbiased protein occupancy profiles on whole cell transcriptome without the requirement of poly-A pulldown. Our study demonstrates that ~ 68% of the total POP-seq peaks exhibited an overlap with publicly available protein–RNA interaction profiles of 97 RNA binding proteins (RBPs) in K562 cells. We show that POP-seq variants consistently capture protein–RNA interaction sites across a broad range of genes including on transcripts encoding for transcription factors (TFs), RNA-Binding Proteins (RBPs) and long non-coding RNAs (lncRNAs). POP-seq identified peaks exhibited a significant enrichment (p value < 2.2e−16) for GWAS SNPs, phenotypic, clinically relevant germline as well as somatic variants reported in cancer genomes, suggesting the prevalence of uncharacterized genomic variation in protein occupied sites on RNA. We demonstrate that the abundance of POP-seq peaks increases with an increase in expression of lncRNAs, suggesting that highly expressed lncRNA are likely to act as sponges for RBPs, contributing to the rewiring of protein–RNA interaction network in cancer cells. Overall, our data supports POP-seq as a robust and cost-effective method that could be applied to primary tissues for mapping global protein occupancies.

Machine learning suggests polygenic risk for cognitive dysfunction in amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a multi-system disease characterized primarily by progressive muscle weakness. Cognitive dysfunction is commonly observed in patients; however, factors influencing risk for cognitive dysfunction remain elusive. Using sparse canonical correlation analysis (sCCA), an unsupervised machine-learning technique, we observed that single nucleotide polymorphisms collectively associate with baseline cognitive performance in a large ALS patient cohort (N = 327) from the multicenter Clinical Research in ALS and Related Disorders for Therapeutic Development (CReATe) Consortium. We demonstrate that a polygenic risk score derived using sCCA relates to longitudinal cognitive decline in the same cohort and also to in vivo cortical thinning in the orbital frontal cortex, anterior cingulate cortex, lateral temporal cortex, premotor cortex, and hippocampus (N = 90) as well as post-mortem motor cortical neuronal loss (N = 87) in independent ALS cohorts from the University of Pennsylvania Integrated Neurodegenerative Disease Biobank. Our findings suggest that common genetic polymorphisms may exert a polygenic contribution to the risk of cortical disease vulnerability and cognitive dysfunction in ALS.

Rad18 mediates specific mutational signatures and shapes the genomic landscape of carcinogen-induced tumors in vivo

The E3 ubiquitin ligase Rad18 promotes a damage-tolerant and error-prone mode of DNA replication termed trans-lesion synthesis that is pathologically activated in cancer. However, the impact of vertebrate Rad18 on cancer genomes is not known. To determine how Rad18 affects mutagenesis in vivo, we have developed and implemented a novel computational pipeline to analyze genomes of carcinogen (7, 12-Dimethylbenz[a]anthracene, DMBA)-induced skin tumors from Rad18^+/+ and Rad18^{- / -} mice. We show that Rad18 mediates specific mutational signatures characterized by high levels of A(T)>T(A) single nucleotide variations (SNVs). In Rad18^{- /-} tumors, an alternative mutation pattern arises, which is characterized by increased numbers of deletions >4 bp. Comparison with annotated human mutational signatures shows that COSMIC signature 22 predominates in Rad18^+/+ tumors whereas Rad18^{- / -} tumors are characterized by increased contribution of COSMIC signature 3 (a hallmark of BRCA-mutant tumors). Analysis of The Cancer Genome Atlas shows that RAD18 expression is strongly associated with high SNV burdens, suggesting RAD18 also promotes mutagenesis in human cancers. Taken together, our results show Rad18 promotes mutagenesis in vivo, modulates DNA repair pathway choice in neoplastic cells, and mediates specific mutational signatures that are present in human tumors.

Bioinformatics: new tools and applications in life science and personalized medicine

While we have a basic understanding of the functioning of the gene when coding sequences of specific proteins, we feel the lack of information on the role that DNA has on specific diseases or functions of thousands of proteins that are produced. Bioinformatics combines the methods used in the collection, storage, identification, analysis, and correlation of this huge and complex information. All this work produces an "ocean" of information that can only be "sailed" with the help of computerized methods. The goal is to provide scientists with the right means to explain normal biological processes, dysfunctions of these processes which give rise to disease and approaches that allow the discovery of new medical cures. Recently, sequencing platforms, a large scale of genomes and transcriptomes, have created new challenges not only to the genomics but especially for bioinformatics. The intent of this article is to compile a list of tools and information resources used by scientists to treat information from the massive sequencing of recent platforms to new generations and the applications of this information in different areas of life sciences including medicine. KEY POINTS: • Biological data mining • Omic approaches • From genotype to phenotype.

Interpretation of psychiatric genome-wide association studies with multispecies heterogeneous functional genomic data integration

Genome-wide association studies and other discovery genetics methods provide a means to identify previously unknown biological mechanisms underlying behavioral disorders that may point to new therapeutic avenues, augment diagnostic tools, and yield a deeper understanding of the biology of psychiatric conditions. Recent advances in psychiatric genetics have been made possible through large-scale collaborative efforts. These studies have begun to unearth many novel genetic variants associated with psychiatric disorders and behavioral traits in human populations. Significant challenges remain in characterizing the resulting disease-associated genetic variants and prioritizing functional follow-up to make them useful for mechanistic understanding and development of therapeutics. Model organism research has generated extensive genomic data that can provide insight into the neurobiological mechanisms of variant action, but a cohesive effort must be made to establish which aspects of the biological modulation of behavioral traits are evolutionarily conserved across species. Scalable computing, new data integration strategies, and advanced analysis methods outlined in this review provide a framework to efficiently harness model organism data in support of clinically relevant psychiatric phenotypes.

Genomewide alteration of histone H3K4 methylation underlies genetic vulnerability to psychopathology

Dysregulated histone methylation has emerged as a recurring theme in multiple neuropsychiatric disorders. However, it is yet unclear whether the altered histone methylation is associated with aetiologic mechanisms or an outcome of disease manifestation. In this study, we examined the genomewide association studies datasets of three major psychiatric disorders, schizophrenia (SCZ), bipolar disorder (BIP), and major depression disorder (MDD), which represents a total of 231,783 cases and 425,444 controls, to clarify the relationship. Our gene-set enrichment analysis results identified statistically significant association of genes involved in three histone methylation biological processes with the three adult-onset psychiatric disorders, which is mainly driven by the histone H3K4 methylation pathway (GO: 0051568). Further analysis of histone H3K4 methylation pathway genes revealed a widespread role of the genes in brain function and disease; 29 (52%) and 41 genes (73.2%) were associated with at least one brain-related trait or brain disorder, respectively. Spatiotemporal gene expression analysis suggests that these pathway genes play a critical role during the prenatal period and are consistent regulators in the cerebral cortex throughout an individual's life. AUTS2, DNMT1 and TET2 are genes of particular interest due to their pervasive role in various aspects of brain function. Our findings support a critical aetiologic role of H3K4 methylation genes shared across SCZ, BIP and MDD, providing new direction for the development of epigenetically-focussed drugs targeting common causal factors of these devastating disorders.

GC and Repeats Profiling along Chromosomes-The Future of Fish Compositional Cytogenomics

The study of fish cytogenetics has been impeded by the inability to produce G-bands that could assign chromosomes to their homologous pairs. Thus, the majority of karyotypes published have been estimated based on morphological similarities of chromosomes. The reason why chromosome G-banding does not work in fish remains elusive. However, the recent increase in the number of fish genomes assembled to the chromosome level provides a way to analyse this issue. We have developed a Python tool to visualize and quantify GC percentage (GC%) of both repeats and unique DNA along chromosomes using a non-overlapping sliding window approach. Our tool profiles GC% and simultaneously plots the proportion of repeats (rep%) in a color scale (or vice versa). Hence, it is possible to assess the contribution of repeats to the total GC%. The main differences are the GC% of repeats homogenizing the overall GC% along fish chromosomes and a greater range of GC% scattered along fish chromosomes. This may explain the inability to produce G-banding in fish. We also show an occasional banding pattern along the chromosomes in some fish that probably cannot be detected with traditional qualitative cytogenetic methods.

Transcription-translation error: In-silico investigation of the structural and functional impact of deleterious single nucleotide polymorphisms in GULP1 gene

Nonsynonymous single nucleotide polymorphisms (nsSNPs) are one of the most common forms of mutations known to disrupt the product of translation thereby altering the protein structure-function relationship. GULP1 (PTB domain-containing engulfment adaptor protein 1) is an evolutionarily conserved adaptor protein that has been associated with glycated hemoglobin (HbA1c) in Genome-Wide Association Studies (GWAS). In order to understand the role of GULP1 in the etiology of diabetes, it is important to study some functional nsSNPs present within the GULP1 protein. We, therefore, used a SNPinformatics approach to retrieve, classify, and determine the stability effect of some nsSNPs. Y27C, G142D, A144T, and Y149C were jointly predicted by the pathogenic-classifying tools to be disease-causing, however, only G142D, A144T, and Y149C had their structural architecture perturbed as predicted by I-MUTANT and MuPro. Interestingly, G142D and Y149C occur at positions 142 and 149 of GULP1 which coincidentally are found within the binding site of GULP1. Protein-Protein interaction analysis also revealed that GULP1 interacted with 10 proteins such as Cell division cycle 5-like protein (CDC5L), ADP-ribosylation factor 6 (ARF6), Arf-GAP with coiled-coil (ACAP1), and Multiple epidermal growth factor-like domains protein 10 (MEGF10), etc. Taken together, rs1357922096, rs1264999716, and rs128246649 could be used as genetic biomarkers for the diagnosis of diabetes. However, being a computational study, these nsSNPs require experimental validation to explore their metabolic involvement in the pathogenesis of diseases.

The functions of CAP superfamily proteins in mammalian fertility and disease

BACKGROUND

Members of the cysteine-rich secretory proteins (CRISPS), antigen 5 (Ag5) and pathogenesis-related 1 (Pr-1) (CAP) superfamily of proteins are found across the bacterial, fungal, plant and animal kingdoms. Although many CAP superfamily proteins remain poorly characterized, over the past decade evidence has accumulated, which provides insights into the functional roles of these proteins in various processes, including fertilization, immune defence and subversion, pathogen virulence, venom toxicology and cancer biology.

OBJECTIVE AND RATIONALE

The aim of this article is to summarize the current state of knowledge on CAP superfamily proteins in mammalian fertility, organismal homeostasis and disease pathogenesis.

SEARCH METHODS

The scientific literature search was undertaken via PubMed database on all articles published prior to November 2019. Search terms were based on following keywords: 'CAP superfamily', 'CRISP', 'Cysteine-rich secretory proteins', 'Antigen 5', 'Pathogenesis-related 1', 'male fertility', 'CAP and CTL domain containing', 'CRISPLD1', 'CRISPLD2', 'bacterial SCP', 'ion channel regulator', 'CatSper', 'PI15', 'PI16', 'CLEC', 'PRY proteins', 'ASP proteins', 'spermatogenesis', 'epididymal maturation', 'capacitation' and 'snake CRISP'. In addition to that, reference lists of primary and review article were reviewed for additional relevant publications.

OUTCOMES

In this review, we discuss the breadth of knowledge on CAP superfamily proteins with regards to their protein structure, biological functions and emerging significance in reproduction, health and disease. We discuss the evolution of CAP superfamily proteins from their otherwise unembellished prokaryotic predecessors into the multi-domain and neofunctionalized members found in eukaryotic organisms today. At least in part because of the rapid evolution of these proteins, many inconsistencies in nomenclature exist within the literature. As such, and in part through the use of a maximum likelihood phylogenetic analysis of the vertebrate CRISP subfamily, we have attempted to clarify this confusion, thus allowing for a comparison of orthologous protein function between species. This framework also allows the prediction of functional relevance between species based on sequence and structural conservation.

WIDER IMPLICATIONS

This review generates a picture of critical roles for CAP proteins in ion channel regulation, sterol and lipid binding and protease inhibition, and as ligands involved in the induction of multiple cellular processes.

Present and Future Salmonid Cytogenetics

Salmonids are extremely important economically and scientifically; therefore, dynamic developments in their research have occurred and will continue occurring in the future. At the same time, their complex phylogeny and taxonomy are challenging for traditional approaches in research. Here, we first provide discoveries regarding the hitherto completely unknown cytogenetic characteristics of the Anatolian endemic flathead trout, Salmo platycephalus, and summarize the presently known, albeit highly complicated, situation in the genus Salmo. Secondly, by outlining future directions of salmonid cytogenomics, we have produced a prototypical virtual karyotype of Salmo trutta, the closest relative of S. platycephalus. This production is now possible thanks to the high-quality genome assembled to the chromosome level in S. trutta via soft-masking, including a direct labelling of repetitive sequences along the chromosome sequence. Repetitive sequences were crucial for traditional fish cytogenetics and hence should also be utilized in fish cytogenomics. As such virtual karyotypes become increasingly available in the very near future, it is necessary to integrate both present and future approaches to maximize their respective benefits. Finally, we show how the presumably repetitive sequences in salmonids can change the understanding of the overall relationship between genome size and G+C content, creating another outstanding question in salmonid cytogenomics waiting to be resolved.

Polymorphisms in TLR4 Gene Associated With Somatic Cell Score in Water Buffaloes (Bubalus bubalis)

Considering the importance of the diseases affecting the productive performance of animals in the dairy industry worldwide, it is necessary to implement tools that help to control and limit the occurrence of such diseases. As the increased somatic cell counts (SCC) are a direct expression of the inflammatory process, they are candidates to become the usual parameter for assessing udder health regarding milk quality and for monitoring mastitis incidences. Toll-Like Receptors are membrane proteins that play a key role in immunity, recognizing pathogens and, subsequently, activating immune responses. The present study was conducted to identify single nucleotide polymorphisms in the TLR4 gene of buffaloes and to analyze its associations with somatic cell counts. DNA samples of 120 Murrah buffaloes were used. The whole coding region of the TLR4 gene was amplified by polymerase chain reaction reactions and sequenced for polymorphism scanning. A total of 13 polymorphisms were identified for the sequenced regions of the TLR4, most of which are in the coding region. The association with the somatic cell score was highly significant (p < 0.001) for all identified polymorphisms of TLR4 gene (g.54621T>A, g.54429G>T, g.54407T>A, g.46616C>A, g.46613T>G, g.46612A>G, g.46611C>A, g.46609T>G, g.46541C>G, g.46526C>A, g.46516T>C, g.46376C>T, g.46372T>C). Therefore, it is suggested that the markers of the TLR4 gene can be used as molecular markers for mastitis resistance in buffaloes, due to their association with somatic cell counts.

In silico mutagenesis of human ACE2 with S protein and translational efficiency explain SARS-CoV-2 infectivity in different species

The coronavirus disease COVID-19 constitutes the most severe pandemic of the last decades having caused more than 1 million deaths worldwide. The SARS-CoV-2 virus recognizes the angiotensin converting enzyme 2 (ACE2) on the surface of human cells through its spike protein. It has been reported that the coronavirus can mildly infect cats, and ferrets, and perhaps dogs while not pigs, mice, chicken and ducks. Differences in viral infectivity among different species or individuals could be due to amino acid differences at key positions of the host proteins that interact with the virus, the immune response, expression levels of host proteins and translation efficiency of the viral proteins among other factors. Here, first we have addressed the importance that sequence variants of different animal species, human individuals and virus isolates have on the interaction between the RBD domain of the SARS-CoV-2 spike S protein and human angiotensin converting enzyme 2 (ACE2). Second, we have looked at viral translation efficiency by using the tRNA adaptation index. We find that integration of both interaction energy with ACE2 and translational efficiency explains animal infectivity. Humans are the top species in which SARS-CoV-2 is both efficiently translated as well as optimally interacting with ACE2. We have found some viral mutations that increase affinity for hACE and some hACE2 variants affecting ACE2 stability and virus binding. These variants suggest that different sensitivities to coronavirus infection in humans could arise in some cases from allelic variability affecting ACE2 stability and virus binding.

Characterization of reproductive gene diversity in the endangered Tasmanian devil

Interindividual variation at genes known to play a role in reproduction may impact reproductive fitness. The Tasmanian devil is an endangered Australian marsupial with low genetic diversity. Recent work has shown concerning declines in productivity in both wild and captive populations over time. Understanding whether functional diversity exists at reproductive genes in the Tasmanian devil is a key first step in identifying genes that may influence productivity. We characterized single nucleotide polymorphisms (SNPs) at 214 genes involved in reproduction in 37 Tasmanian devils. Twenty genes contained nonsynonymous substitutions, with genes involved in embryogenesis, fertilization and hormonal regulation of reproduction displaying greater numbers of nonsynonymous SNPs than synonymous SNPs. Two genes, ADAMTS9 and NANOG, showed putative signatures of balancing selection indicating that natural selection is maintaining diversity at these genes despite the species exhibiting low overall levels of genetic diversity. We will use this information in future to examine the interplay between reproductive gene variation and reproductive fitness in Tasmanian devil populations.

Determinants of vitamin D status in Kenyan calves

Vitamin D plays a critical role in calcium homeostasis and in the maintenance and development of skeletal health. Vitamin D status has increasingly been linked to non-skeletal health outcomes such as all-cause mortality, infectious diseases and reproductive outcomes in both humans and veterinary species. We have previously demonstrated a relationship between vitamin D status, assessed by the measurement of serum concentrations of the major vitamin D metabolite 25 hydroxyvitamin D (25(OH)D), and a wide range of non-skeletal health outcomes in companion and wild animals. The aims of this study were to define the host and environmental factors associated with vitamin D status in a cohort of 527 calves from Western Kenya which were part of the Infectious Disease of East African Livestock (IDEAL) cohort. A secondary aim was to explore the relationship between serum 25(OH)D concentrations measured in 7-day old calves and subsequent health outcomes over the following 12 months. A genome wide association study demonstrated that both dietary and endogenously produced vitamin D metabolites were under polygenic control in African calves. In addition, we found that neonatal vitamin D status was not predictive of the subsequent development of an infectious disease event or mortality over the 12 month follow up period.

nMAGMA: a network-enhanced method for inferring risk genes from GWAS summary statistics and its application to schizophrenia

MOTIVATION

Annotating genetic variants from summary statistics of genome-wide association studies (GWAS) is crucial for predicting risk genes of various disorders. The multimarker analysis of genomic annotation (MAGMA) is one of the most popular tools for this purpose, where MAGMA aggregates signals of single nucleotide polymorphisms (SNPs) to their nearby genes. In biology, SNPs may also affect genes that are far away in the genome, thus missed by MAGMA. Although different upgrades of MAGMA have been proposed to extend gene-wise variant annotations with more information (e.g. Hi-C or eQTL), the regulatory relationships among genes and the tissue specificity of signals have not been taken into account.

RESULTS

We propose a new approach, namely network-enhanced MAGMA (nMAGMA), for gene-wise annotation of variants from GWAS summary statistics. Compared with MAGMA and H-MAGMA, nMAGMA significantly extends the lists of genes that can be annotated to SNPs by integrating local signals, long-range regulation signals (i.e. interactions between distal DNA elements), and tissue-specific gene networks. When applied to schizophrenia (SCZ), nMAGMA is able to detect more risk genes (217% more than MAGMA and 57% more than H-MAGMA) that are involved in SCZ compared with MAGMA and H-MAGMA, and more of nMAGMA results can be validated with known SCZ risk genes. Some disease-related functions (e.g. the ATPase pathway in Cortex) are also uncovered in nMAGMA but not in MAGMA or H-MAGMA. Moreover, nMAGMA provides tissue-specific risk signals, which are useful for understanding disorders with multitissue origins.

The APOE gene cluster responds to air pollution factors in mice with coordinated expression of genes that differs by age in humans

Little is known of gene-environment interactions for Alzheimer's disease (AD) risk factors. Apolipoprotein E (APOE) and neighbors on chromosome 19q13.3 have variants associated with risks of AD, but with unknown mechanism. This study describes novel links among the APOE network, air pollution, and age-related diseases. Mice exposed to air pollution nano-sized particulate matter (nPM) had coordinate responses of Apoe-Apoc1-Tomm40 in the cerebral cortex. In humans, the AD vulnerable hippocampus and amygdala had stronger age decline in APOE cluster expression than the AD-resistant cerebellum and hypothalamus. Using consensus weighted gene co-expression network, we showed that APOE has a conserved co-expressed network in rodent and primate brains. SOX1, which has AD-associated single nucleotide polymorphisms, was among the co-expressed genes in the human hippocampus. Humans and mice shared 87% of potential binding sites for transcription factors in APOE cluster promoter, suggesting similar inducibility and a novel link among environment, APOE cluster, and risk of AD.

Prevalence of PKD1 gene mutation in cats in Turkey and pathogenesis of feline polycystic kidney disease

Polycystic kidney disease (PKD) is one of the most common hereditary diseases in cats, with high prevalence in Persian and Persian-related cats. PKD is caused mainly by an inherited autosomal dominant (AD) mutation, and animals may be asymptomatic for years. We screened 16 cats from various breeds exhibiting a renal abnormality by ultrasound examination and genotyped them for the c.10063C>A transversion on exon 29 of the polycystin-1 (PKD1) gene, by PCR-restriction fragment length polymorphism (PCR-RFLP). Among these cats, a Siamese nuclear family of 4 cats with ancestral hereditary renal failure were screened by whole-genome sequencing (WGS) to determine novel variations in genes associated with both AD and autosomal recessive PKD in humans. During the study period, one cat died as a result of renal failure and was forwarded for autopsy. Additionally, we screened 294 cats asymptomatic for renal disease (Angora, Van, Persian, Siamese, Scottish Fold, Exotic Shorthair, British Shorthair, and mixed breeds) to determine the prevalence of the mutation in cats in Turkey. Ten of the symptomatic and 2 of the asymptomatic cats carried the heterozygous C → A transversion, indicating a prevalence of 62.5% and 0.68%, respectively. In the WGS analysis of 4 cats in the Siamese nuclear family, novel variations were determined in the fibrocystin gene (PKHD1), which was not compatible with dominant inheritance of PKD.

Identification of synthetic chemosensitivity genes paired with BRAF for BRAF/MAPK inhibitors

Molecular-targeted approaches are important for personalised cancer treatment, which requires knowledge regarding drug target specificity. Here, we used the synthetic lethality concept to identify candidate gene pairs with synergistic effects on drug responses. A synergistic chemo-sensitivity response was identified if a drug had a significantly lower half-maximal inhibitory concentration (IC₅₀) in cell lines with a pair of mutated genes compared with those in other cell lines (wild-type or one mutated gene). Among significantly damaging mutations in the Genomics of Drug Sensitivity in Cancer database, we found 580 candidate synergistic chemo-sensitivity interaction sets for 456 genes and 54 commercial drugs. Clustering analyses according to drug/gene and drug/tissue interactions showed that BRAF/MAPK inhibitors clustered together; 11 partner genes for BRAF were identified. The combined effects of these partners on IC₅₀ values were significant for both drug-specific and drug-combined comparisons. Survival analysis using The Cancer Genome Atlas data showed that patients who had mutated gene pairs in synergistic interaction sets had longer overall survival compared with that in patients with other mutation profiles. Overall, this analysis demonstrated that synergistic drug-responsive gene pairs could be successfully used as predictive markers of drug sensitivity and patient survival, offering new targets for personalised medicine.

Methamphetamine Increases the Proportion of SIV-Infected Microglia/Macrophages, Alters Metabolic Pathways, and Elevates Cell Death Pathways: A Single-Cell Analysis

Both substance use disorder and HIV infection continue to affect many individuals. Both have untoward effects on the brain, and the two conditions often co-exist. In the brain, macrophages and microglia are infectable by HIV, and these cells are also targets for the effects of drugs of abuse, such as the psychostimulant methamphetamine. To determine the interaction of HIV and methamphetamine, we isolated microglia and brain macrophages from SIV-infected rhesus monkeys that were treated with or without methamphetamine. Cells were subjected to single-cell RNA sequencing and results were analyzed by statistical and bioinformatic analysis. In the animals treated with methamphetamine, a significantly increased proportion of the microglia and/or macrophages were infected by SIV. In addition, gene encoding functions in cell death pathways were increased, and the brain-derived neurotropic factor pathway was inhibited. The gene expression patterns in infected cells did not cluster separately from uninfected cells, but clusters comprised of microglia and/or macrophages from methamphetamine-treated animals differed in neuroinflammatory and metabolic pathways from those comprised of cells from untreated animals. Methamphetamine increases CNS infection by SIV and has adverse effects on both infected and uninfected microglia and brain macrophages, highlighting the dual and interacting harms of HIV infection and drug abuse on the brain.

Deep learning based DNA:RNA triplex forming potential prediction

BACKGROUND

Long non-coding RNAs (lncRNAs) can exert functions via forming triplex with DNA. The current methods in predicting the triplex formation mainly rely on mathematic statistic according to the base paring rules. However, these methods have two main limitations: (1) they identify a large number of triplex-forming lncRNAs, but the limited number of experimentally verified triplex-forming lncRNA indicates that maybe not all of them can form triplex in practice, and (2) their predictions only consider the theoretical relationship while lacking the features from the experimentally verified data.

RESULTS

In this work, we develop an integrated program named TriplexFPP (Triplex Forming Potential Prediction), which is the first machine learning model in DNA:RNA triplex prediction. TriplexFPP predicts the most likely triplex-forming lncRNAs and DNA sites based on the experimentally verified data, where the high-level features are learned by the convolutional neural networks. In the fivefold cross validation, the average values of Area Under the ROC curves and PRC curves for removed redundancy triplex-forming lncRNA dataset with threshold 0.8 are 0.9649 and 0.9996, and these two values for triplex DNA sites prediction are 0.8705 and 0.9671, respectively. Besides, we also briefly summarize the cis and trans targeting of triplexes lncRNAs.

CONCLUSIONS

The TriplexFPP is able to predict the most likely triplex-forming lncRNAs from all the lncRNAs with computationally defined triplex forming capacities and the potential of a DNA site to become a triplex. It may provide insights to the exploration of lncRNA functions.

DGCR8 microprocessor defect characterizes familial multinodular goiter with schwannomatosis

BACKGROUNDDICER1 is the only miRNA biogenesis component associated with an inherited tumor syndrome, featuring multinodular goiter (MNG) and rare pediatric-onset lesions. Other susceptibility genes for familial forms of MNG likely exist.METHODSWhole-exome sequencing of a kindred with early-onset MNG and schwannomatosis was followed by investigation of germline pathogenic variants that fully segregated with the disease. Genome-wide analyses were performed on 13 tissue samples from familial and nonfamilial DGCR8-E518K-positive tumors, including MNG, schwannomas, papillary thyroid cancers (PTCs), and Wilms tumors. miRNA profiles of 4 tissue types were compared, and sequencing of miRNA, pre-miRNA, and mRNA was performed in a subset of 9 schwannomas, 4 of which harbor DGCR8-E518K.RESULTSWe identified c.1552G>A;p.E518K in DGCR8, a microprocessor component located in 22q, in the kindred. The variant identified is a somatic hotspot in Wilms tumors and has been identified in 2 PTCs. Copy number loss of chromosome 22q, leading to loss of heterozygosity at the DGCR8 locus, was found in all 13 samples harboring c.1552G>A;p.E518K. miRNA profiling of PTCs, MNG, schwannomas, and Wilms tumors revealed a common profile among E518K hemizygous tumors. In vitro cleavage demonstrated improper processing of pre-miRNA by DGCR8-E518K. MicroRNA and RNA profiling show that this variant disrupts precursor microRNA production, impacting populations of canonical microRNAs and mirtrons.CONCLUSIONWe identified DGCR8 as the cause of an unreported autosomal dominant mendelian tumor susceptibility syndrome: familial multinodular goiter with schwannomatosis.FUNDINGCanadian Institutes of Health Research, Compute Canada, Alex's Lemonade Stand Foundation, the Mia Neri Foundation for Childhood Cancer, Cassa di Sovvenzioni e Risparmio fra il Personale della Banca d'Italia, and the KinderKrebsInitiative Buchholz/Holm-Seppensen.

A graph-based approach identifies dynamic H-bond communication networks in spike protein S of SARS-CoV-2

Corona virus spike protein S is a large homo-trimeric protein anchored in the membrane of the virion particle. Protein S binds to angiotensin-converting-enzyme 2, ACE2, of the host cell, followed by proteolysis of the spike protein, drastic protein conformational change with exposure of the fusion peptide of the virus, and entry of the virion into the host cell. The structural elements that govern conformational plasticity of the spike protein are largely unknown. Here, we present a methodology that relies upon graph and centrality analyses, augmented by bioinformatics, to identify and characterize large H-bond clusters in protein structures. We apply this methodology to protein S ectodomain and find that, in the closed conformation, the three protomers of protein S bring the same contribution to an extensive central network of H-bonds, and contribute symmetrically to a relatively large H-bond cluster at the receptor binding domain, and to a cluster near a protease cleavage site. Markedly different H-bonding at these three clusters in open and pre-fusion conformations suggest dynamic H-bond clusters could facilitate structural plasticity and selection of a protein S protomer for binding to the host receptor, and proteolytic cleavage. From analyses of spike protein sequences we identify patches of histidine and carboxylate groups that could be involved in transient proton binding.

The Y Chromosome: A Complex Locus for Genetic Analyses of Complex Human Traits

The Human Y chromosome (ChrY) has been demonstrated to be a powerful tool for phylogenetics, population genetics, genetic genealogy and forensics. However, the importance of ChrY genetic variation in relation to human complex traits is less clear. In this review, we summarise existing evidence about the inherent complexities of ChrY variation and their use in association studies of human complex traits. We present and discuss the specific particularities of ChrY genetic variation, including Y chromosomal haplogroups, that need to be considered in the design and interpretation of genetic epidemiological studies involving ChrY.

Prioritizing long range interactions in noncoding regions using GWAS and deletions perturbed TADs

Genome-wide association studies (GWAS) have contributed significantly to predisposing the disease etiology by associating single nucleotide polymorphisms (SNPs) with complex diseases. However, most GWAS-SNPs are in the noncoding regions that may affect distal genes via long range enhancer-promoter interactions. Thus, the common practice on GWAS discoveries cannot fully reveal the molecular mechanisms underpinning complex diseases. It is known that perturbations of topological associated domains (TADs) lead to long range interactions which underlie disease etiology. To identify the probable long range interactions in noncoding regions via GWAS and TADs perturbed by deletions, we integrated datasets from GWAS-SNPs, enhancers, TADs, and deletions. After ranking and clustering, we prioritized 201,132 high confident pairs of GWAS-SNPs and target genes. In this study, we performed a systematic inference on noncoding regions via GWAS-SNPs and deletion-perturbed TADs to boost GWAS discovery power. The high confident pairs of GWAS-SNPs and target genes (SE-Gs) provide the promising candidates to understand the molecular mechanisms underlying complex diseases with emphasis on the three-dimensional genome.