Evidence supports a causal association between allele-specific vitamin D receptor binding and multiple sclerosis among Europeans.

Although evidence exists for a causal association between 25-hydroxyvitamin D (25(OH)D) serum levels, and multiple sclerosis (MS), the role of variation in vitamin D receptor (VDR) binding in MS is unknown. Here, we leveraged previously identified variants associated with allele imbalance in VDR binding (VDR-binding variant; VDR-BV) in ChIP-exo data from calcitriol-stimulated lymphoblastoid cell lines and 25(OH)D serum levels from genome-wide association studies to construct genetic instrumental variables (GIVs). GIVs are composed of one or more genetic variants that serve as proxies for exposures of interest. Here, GIVs for both VDR-BVs and 25(OH)D were used in a two-sample Mendelian Randomization study to investigate the relationship between VDR binding at a locus, 25(OH)D serum levels, and MS risk. Data for 13,598 MS cases and 38,887 controls of European ancestry from Kaiser Permanente Northern California, Swedish MS studies, and the UK Biobank were included. We estimated the association between each VDR-BV GIV and MS. Significant interaction between a VDR-BV GIV and a GIV for serum 25OH(D) was evidence for a causal association between VDR-BVs and MS unbiased by pleiotropy. We observed evidence for associations between two VDR-BVs (rs2881514, rs2531804) and MS after correction for multiple tests. There was evidence of interaction between rs2881514 and a 25(OH)D GIV, providing evidence of a causal association between rs2881514 and MS. This study is the first to demonstrate evidence that variation in VDR binding at a locus contributes to MS risk. Our results are relevant to other autoimmune diseases in which vitamin D plays a role.

Comparison of T-cell receptor diversity of people with myalgic encephalomyelitis versus controls.

OBJECTIVE: Myalgic Encephalomyelitis (ME; sometimes referred to as Chronic Fatigue Syndrome) is a chronic disease without laboratory test, detailed aetiological understanding or effective therapy. Its symptoms are diverse, but it is distinguished from other fatiguing illnesses by the experience of post-exertional malaise, the worsening of symptoms even after minor physical or mental exertion. Its frequent onset after infection suggests autoimmune involvement or that it arises from abnormal T-cell activation. RESULTS: To test this hypothesis, we sequenced the genomic loci of α/δ, ß and γ T-cell receptors (TCR) from 40 human blood samples from each of four groups: severely affected people with ME; mildly or moderately affected people with ME; people diagnosed with Multiple Sclerosis, as disease controls; and, healthy controls. Seeking to automatically classify these individuals' samples by their TCR repertoires, we applied P-SVM, a machine learning method. However, despite working well on a simulated data set, this approach did not allow statistically significant partitioning of samples into the four subgroups. Our findings do not support the hypothesis that blood samples from people with ME frequently contain altered T-cell receptor diversity.

The histone chaperone SPT2 regulates chromatin structure and function in Metazoa.

Histone chaperones control nucleosome density and chromatin structure. In yeast, the H3-H4 chaperone Spt2 controls histone deposition at active genes but its roles in metazoan chromatin structure and organismal physiology are not known. Here we identify the Caenorhabditis elegans ortholog of SPT2 (CeSPT-2) and show that its ability to bind histones H3-H4 is important for germline development and transgenerational epigenetic gene silencing, and that spt-2 null mutants display signatures of a global stress response. Genome-wide profiling showed that CeSPT-2 binds to a range of highly expressed genes, and we find that spt-2 mutants have increased chromatin accessibility at a subset of these loci. We also show that SPT2 influences chromatin structure and controls the levels of soluble and chromatin-bound H3.3 in human cells. Our work reveals roles for SPT2 in controlling chromatin structure and function in Metazoa.

Typing myalgic encephalomyelitis by infection at onset: A DecodeME study.

Background: People with myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) experience core symptoms of post-exertional malaise, unrefreshing sleep, and cognitive impairment. Despite numbering 0.2-0.4% of the population, no laboratory test is available for their diagnosis, no effective therapy exists for their treatment, and no scientific breakthrough regarding pathogenesis has been made. It remains unknown, despite decades of small-scale studies, whether individuals experience different types of ME/CFS separated by onset-type, sex or age. Methods: DecodeME is a large population-based study of ME/CFS that recruited 17,074 participants in the first 3 months following full launch. Detailed questionnaire responses from UK-based participants who all reported being diagnosed with ME/CFS by a health professional provided an unparalleled opportunity to investigate, using logistic regression, whether ME/CFS severity or onset type is significantly associated with sex, age, illness duration, comorbid conditions or symptoms. Results: The well-established sex-bias among ME/CFS patients is evident in the initial DecodeME cohort: 83.5% of participants were females. What was not known previously was that females tend to have more comorbidities than males. Moreover, being female, being older and being over 10 years from ME/CFS onset are significantly associated with greater severity. Five different ME/CFS onset types were examined in the self-reported data: those with ME/CFS onset (i) after glandular fever (infectious mononucleosis); (ii) after COVID-19 infection; (iii) after other infections; (iv) without an infection at onset; and, (v) where the occurrence of an infection at or preceding onset is not known. Among other findings, ME/CFS onset with unknown infection status was significantly associated with active fibromyalgia. Conclusions: DecodeME participants differ in symptoms, comorbid conditions and/or illness severity when stratified by their sex-at-birth and/or infection around the time of ME/CFS onset.

Myalgic Encephalomyelitis / Chronic Fatigue Syndrome (ME/CFS) is a chronic disease that affects an estimated 250,000 people in the UK. Its defining symptom is post-exertional malaise, an excessive delayed worsening of symptoms following even minor physical or mental exertion. For those with it, ME/CFS means disability and poor quality of life. DecodeME is a research study which is looking for DNA differences between people with ME/CFS and people without any health problems. People with ME/CFS who take part in DecodeME complete a questionnaire that assesses their symptoms and whether they will then be invited to donate a DNA sample. This paper analyses the answers to this questionnaire; we will publish results of the DNA analysis separately. So far, more than 17 thousand people with ME/CFS have completed the DecodeME questionnaire. Their answers help us to address the question: "Are there different types of ME/CFS linked to different causes and how severe it becomes?" Results show that people with ME/CFS do not form a single group reporting similar symptoms and additional medical conditions. Instead, participants who had an infection at the start of their ME/CFS reported a different pattern of symptoms and conditions compared to those without an infection. It is well known that most people with ME/CFS are females. What was not clear previously was that females tend to have more additional health conditions. Also, being female, being older and being over 10 years from ME/CFS onset all make it more likely that someone is more severely affected by their ME/CFS. These findings could indicate that by studying people with different ME/CFS onset-types separately ­ rather than analysing all people with ME/CFS together ­ it will be easier to understand what is going wrong.

GWAS and meta-analysis identifies 49 genetic variants underlying critical COVID-19.

Critical illness in COVID-19 is an extreme and clinically homogeneous disease phenotype that we have previously shown1 to be highly efficient for discovery of genetic associations2. Despite the advanced stage of illness at presentation, we have shown that host genetics in patients who are critically ill with COVID-19 can identify immunomodulatory therapies with strong beneficial effects in this group3. Here we analyse 24,202 cases of COVID-19 with critical illness comprising a combination of microarray genotype and whole-genome sequencing data from cases of critical illness in the international GenOMICC (11,440 cases) study, combined with other studies recruiting hospitalized patients with a strong focus on severe and critical disease: ISARIC4C (676 cases) and the SCOURGE consortium (5,934 cases). To put these results in the context of existing work, we conduct a meta-analysis of the new GenOMICC genome-wide association study (GWAS) results with previously published data. We find 49 genome-wide significant associations, of which 16 have not been reported previously. To investigate the therapeutic implications of these findings, we infer the structural consequences of protein-coding variants, and combine our GWAS results with gene expression data using a monocyte transcriptome-wide association study (TWAS) model, as well as gene and protein expression using Mendelian randomization. We identify potentially druggable targets in multiple systems, including inflammatory signalling (JAK1), monocyte-macrophage activation and endothelial permeability (PDE4A), immunometabolism (SLC2A5 and AK5), and host factors required for viral entry and replication (TMPRSS2 and RAB2A).

Human genetics seen through an evolutionary lens.

Most DNA bases crucial for species perpetuation are marked by a dearth of sequence change among species related over long evolutionary time. Recently, Christmas et al.1 and Sullivan et al.2 cast light on human DNA and its variants through comparison with 239 other mammalian species' genomes.

An ependymal cell census identifies heterogeneous and ongoing cell maturation in the adult mouse spinal cord that changes dynamically on injury.

The adult spinal cord stem cell potential resides within the ependymal cell population and declines with age. Ependymal cells are, however, heterogeneous, and the biological diversity this represents and how it changes with age remain unknown. Here, we present a single-cell transcriptomic census of spinal cord ependymal cells from adult and aged mice, identifying not only all known ependymal cell subtypes but also immature as well as mature cell states. By comparing transcriptomes of spinal cord and brain ependymal cells, which lack stem cell abilities, we identify immature cells as potential spinal cord stem cells. Following spinal cord injury, these cells re-enter the cell cycle, which is accompanied by a short-lived reversal of ependymal cell maturation. We further analyze ependymal cells in the human spinal cord and identify widespread cell maturation and altered cell identities. This in-depth characterization of spinal cord ependymal cells provides insight into their biology and informs strategies for spinal cord repair.

DecodeME: community recruitment for a large genetics study of myalgic encephalomyelitis / chronic fatigue syndrome.

BACKGROUND: Myalgic encephalomyelitis / chronic fatigue syndrome (ME/CFS) is a common, long-term condition characterised by post-exertional malaise, often with fatigue that is not significantly relieved by rest. ME/CFS has no confirmed diagnostic test or effective treatment and we lack knowledge of its causes. Identification of genes and cellular processes whose disruption adds to ME/CFS risk is a necessary first step towards development of effective therapy. METHODS: Here we describe DecodeME, an ongoing study co-produced by people with lived experience of ME/CFS and scientists. Together we designed the study and obtained funding and are now recruiting up to 25,000 people in the UK with a clinical diagnosis of ME/CFS. Those eligible for the study are at least 16 years old, pass international study criteria, and lack any alternative diagnoses that can result in chronic fatigue. These will include 5,000 people whose ME/CFS diagnosis was a consequence of SARS-CoV-2 infection. Questionnaires are completed online or on paper. Participants' saliva DNA samples are acquired by post, which improves participation by more severely-affected individuals. Digital marketing and social media approaches resulted in 29,000 people with ME/CFS in the UK pre-registering their interest in participating. We will perform a genome-wide association study, comparing participants' genotypes with those from UK Biobank as controls. This should generate hypotheses regarding the genes, mechanisms and cell types contributing to ME/CFS disease aetiology. DISCUSSION: The DecodeME study has been reviewed and given a favourable opinion by the North West - Liverpool Central Research Ethics Committee (21/NW/0169). Relevant documents will be available online ( www.decodeme.org.uk ). Genetic data will be disseminated as associated variants and genomic intervals, and as summary statistics. Results will be reported on the DecodeME website and via open access publications.

Genome-Wide Analysis of Human Long Noncoding RNAs: A Provocative Review.

Do long noncoding RNAs (lncRNAs) contribute little or substantively to human biology? To address how lncRNA loci and their transcripts, structures, interactions, and functions contribute to human traits and disease, we adopt a genome-wide perspective. We intend to provoke alternative interpretation of questionable evidence and thorough inquiry into unsubstantiated claims. We discuss pitfalls of lncRNA experimental and computational methods as well as opposing interpretations of their results. The majority of evidence, we argue, indicates that most lncRNA transcript models reflect transcriptional noise or provide minor regulatory roles, leaving relatively few human lncRNAs that contribute centrally to human development, physiology, or behavior. These important few tend to be spliced and better conserved but lack a simple syntax relating sequence to structure and mechanism, and so resist simple categorization. This genome-wide view should help investigators prioritize individual lncRNAs based on their likely contribution to human biology.

Mapping the developing human cardiac endothelium at single-cell resolution identifies MECOM as a regulator of arteriovenous gene expression.

AIMS: Coronary vasculature formation is a critical event during cardiac development, essential for heart function throughout perinatal and adult life. However, current understanding of coronary vascular development has largely been derived from transgenic mouse models. The aim of this study was to characterize the transcriptome of the human foetal cardiac endothelium using single-cell RNA sequencing (scRNA-seq) to provide critical new insights into the cellular heterogeneity and transcriptional dynamics that underpin endothelial specification within the vasculature of the developing heart. METHODS AND RESULTS: We acquired scRNA-seq data of over 10 000 foetal cardiac endothelial cells (ECs), revealing divergent EC subtypes including endocardial, capillary, venous, arterial, and lymphatic populations. Gene regulatory network analyses predicted roles for SMAD1 and MECOM in determining the identity of capillary and arterial populations, respectively. Trajectory inference analysis suggested an endocardial contribution to the coronary vasculature and subsequent arterialization of capillary endothelium accompanied by increasing MECOM expression. Comparative analysis of equivalent data from murine cardiac development demonstrated that transcriptional signatures defining endothelial subpopulations are largely conserved between human and mouse. Comprehensive characterization of the transcriptional response to MECOM knockdown in human embryonic stem cell-derived EC (hESC-EC) demonstrated an increase in the expression of non-arterial markers, including those enriched in venous EC. CONCLUSIONS: scRNA-seq of the human foetal cardiac endothelium identified distinct EC populations. A predicted endocardial contribution to the developing coronary vasculature was identified, as well as subsequent arterial specification of capillary EC. Loss of MECOM in hESC-EC increased expression of non-arterial markers, suggesting a role in maintaining arterial EC identity.

Refining the domain architecture model of the replication origin firing factor Treslin/TICRR.

Faithful genome duplication requires appropriately controlled replication origin firing. The metazoan origin firing regulation hub Treslin/TICRR and its yeast orthologue Sld3 share the Sld3-Treslin domain and the adjacent TopBP1/Dpb11 interaction domain. We report a revised domain architecture model of Treslin/TICRR. Protein sequence analyses uncovered a conserved Ku70-homologous ß-barrel fold in the Treslin/TICRR middle domain (M domain) and in Sld3. Thus, the Sld3-homologous Treslin/TICRR core comprises its three central domains, M domain, Sld3-Treslin domain, and TopBP1/Dpb11 interaction domain, flanked by non-conserved terminal domains, the CIT (conserved in Treslins) and the C terminus. The CIT includes a von Willebrand factor type A domain. Unexpectedly, MTBP, Treslin/TICRR, and Ku70/80 share the same N-terminal domain architecture, von Willebrand factor type A and Ku70-like ß-barrels, suggesting a common ancestry. Binding experiments using mutants and the Sld3-Sld7 dimer structure suggest that the Treslin/Sld3 and MTBP/Sld7 ß-barrels engage in homotypic interactions, reminiscent of Ku70-Ku80 dimerization. Cells expressing Treslin/TICRR domain mutants indicate that all Sld3-core domains and the non-conserved terminal domains fulfil important functions during origin firing in human cells. Thus, metazoa-specific and widely conserved molecular processes cooperate during metazoan origin firing.

Citrullination Was Introduced into Animals by Horizontal Gene Transfer from Cyanobacteria.

Protein posttranslational modifications add great sophistication to biological systems. Citrullination, a key regulatory mechanism in human physiology and pathophysiology, is enigmatic from an evolutionary perspective. Although the citrullinating enzymes peptidylarginine deiminases (PADIs) are ubiquitous across vertebrates, they are absent from yeast, worms, and flies. Based on this distribution PADIs were proposed to have been horizontally transferred, but this has been contested. Here, we map the evolutionary trajectory of PADIs into the animal lineage. We present strong phylogenetic support for a clade encompassing animal and cyanobacterial PADIs that excludes fungal and other bacterial homologs. The animal and cyanobacterial PADI proteins share functionally relevant primary and tertiary synapomorphic sequences that are distinct from a second PADI type present in fungi and actinobacteria. Molecular clock calculations and sequence divergence analyses using the fossil record estimate the last common ancestor of the cyanobacterial and animal PADIs to be less than 1 billion years old. Additionally, under an assumption of vertical descent, PADI sequence change during this evolutionary time frame is anachronistically low, even when compared with products of likely endosymbiont gene transfer, mitochondrial proteins, and some of the most highly conserved sequences in life. The consilience of evidence indicates that PADIs were introduced from cyanobacteria into animals by horizontal gene transfer (HGT). The ancestral cyanobacterial PADI is enzymatically active and can citrullinate eukaryotic proteins, suggesting that the PADI HGT event introduced a new catalytic capability into the regulatory repertoire of animals. This study reveals the unusual evolution of a pleiotropic protein modification.

OAF: a new member of the BRICHOS family.

Summary: The 10 known BRICHOS domain-containing proteins in humans have been linked to an unusually long list of pathologies, including cancer, obesity and two amyloid-like diseases. BRICHOS domains themselves have been described as intramolecular chaperones that act to prevent amyloid-like aggregation of their proteins' mature polypeptides. Using structural comparison of coevolution-based AlphaFold models and sequence conservation, we identified the Out at First (OAF) protein as a new member of the BRICHOS family in humans. OAF is an experimentally uncharacterized protein that has been proposed as a candidate biomarker for clinical management of coronavirus disease 2019 infections. Our analysis revealed how structural comparison of AlphaFold models can discover remote homology relationships and lead to a better understanding of BRICHOS domain molecular mechanism. Supplementary information: Supplementary data are available at Bioinformatics Advances online.

PERCC1, a new member of the Yap/TAZ/FAM181 transcriptional co-regulator family.

Motivation: Disrupted PERCC1 gene expression causes an intractable congenital diarrhoea in infants. However, this gene's molecular mechanism is unknown and no homologous proteins have been reported. Results: Our detailed evolutionary analysis of PERCC1 sequence reveals it to be a previously unappreciated member of the YAP/TAZ/FAM181 family of homologous transcriptional regulators. Like YAP and TAZ, PERCC1 likely interacts with DNA via binding to TEA/ATTS domain transcription factors (TEADs) using its conserved interface-2 and -3 sequences. We compare the expression patterns of PERCC1 with those of YAP, TAZ, TEADs. Our report provides the identification and first in-depth bioinformatic analysis of a YAP/TAZ homologue, and a likely new regulator of the YAP/TAZ-TEAD transcriptional complex. Availability and implementation: The data underlying this article are available in UniProt Database. Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Extending the Horizon of Homology Detection with Coevolution-based Structure Prediction.

Traditional sequence analysis algorithms fail to identify distant homologies when they lie beyond a detection horizon. In this review, we discuss how co-evolution-based contact and distance prediction methods are pushing back this homology detection horizon, thereby yielding new functional insights and experimentally testable hypotheses. Based on correlated substitutions, these methods divine three-dimensional constraints among amino acids in protein sequences that were previously devoid of all annotated domains and repeats. The new algorithms discern hidden structure in an otherwise featureless sequence landscape. Their revelatory impact promises to be as profound as the use, by archaeologists, of ground-penetrating radar to discern long-hidden, subterranean structures. As examples of this, we describe how triplicated structures reflecting longin domains in MON1A-like proteins, or UVR-like repeats in DISC1, emerge from their predicted contact and distance maps. These methods also help to resolve structures that do not conform to a "beads-on-a-string" model of protein domains. In one such example, we describe CFAP298 whose ubiquitin-like domain was previously challenging to perceive owing to a large sequence insertion within it. More generally, the new algorithms permit an easier appreciation of domain families and folds whose evolution involved structural insertion or rearrangement. As we exemplify with α1-antitrypsin, coevolution-based predicted contacts may also yield insights into protein dynamics and conformational change. This new combination of structure prediction (using innovative co-evolution based methods) and homology inference (using more traditional sequence analysis approaches) shows great promise for bringing into view a sea of evolutionary relationships that had hitherto lain far beyond the horizon of homology detection.

Ageing compromises mouse thymus function and remodels epithelial cell differentiation.

Ageing is characterised by cellular senescence, leading to imbalanced tissue maintenance, cell death and compromised organ function. This is first observed in the thymus, the primary lymphoid organ that generates and selects T cells. However, the molecular and cellular mechanisms underpinning these ageing processes remain unclear. Here, we show that mouse ageing leads to less efficient T cell selection, decreased self-antigen representation and increased T cell receptor repertoire diversity. Using a combination of single-cell RNA-seq and lineage-tracing, we find that progenitor cells are the principal targets of ageing, whereas the function of individual mature thymic epithelial cells is compromised only modestly. Specifically, an early-life precursor cell population, retained in the mouse cortex postnatally, is virtually extinguished at puberty. Concomitantly, a medullary precursor cell quiesces, thereby impairing maintenance of the medullary epithelium. Thus, ageing disrupts thymic progenitor differentiation and impairs the core immunological functions of the thymus.

Genetic risk factors of ME/CFS: a critical review.

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a complex multisystem illness that lacks effective therapy and a biomedical understanding of its causes. Despite a prevalence of ∼0.2-0.4% and its high public health burden, and evidence that it has a heritable component, ME/CFS has not yet benefited from the advances in technology and analytical tools that have improved our understanding of many other complex diseases. Here we critically review existing evidence that genetic factors alter ME/CFS risk before concluding that most ME/CFS candidate gene associations are not replicated by the larger CFS cohort within the UK Biobank. Multiple genome-wide association studies of this cohort also have not yielded consistently significant associations. Ahead of upcoming larger genome-wide association studies, we discuss how these could generate new lines of enquiry into the DNA variants, genes and cell types that are causally involved in ME/CFS disease.