An Association between OXPHOS-Related Gene Expression and Malignant Hyperthermia Susceptibility in Human Skeletal Muscle Biopsies.

Malignant hyperthermia (MH) is a pharmacogenetic condition of skeletal muscle that manifests in hypermetabolic responses upon exposure to volatile anaesthetics. This condition is caused primarily by pathogenic variants in the calcium-release channel RYR1, which disrupts calcium signalling in skeletal muscle. However, our understanding of MH genetics is incomplete, with no variant identified in a significant number of cases and considerable phenotype diversity. In this study, we applied a transcriptomic approach to investigate the genome-wide gene expression in MH-susceptible cases using muscle biopsies taken for diagnostic testing. Baseline comparisons between muscle from MH-susceptible individuals (MHS, n = 8) and non-susceptible controls (MHN, n = 4) identified 822 differentially expressed genes (203 upregulated and 619 downregulated) with significant enrichment in genes associated with oxidative phosphorylation (OXPHOS) and fatty acid metabolism. Investigations of 10 OXPHOS target genes in a larger cohort (MHN: n = 36; MHS: n = 36) validated the reduced expression of ATP5MD and COQ6 in MHS samples, but the remaining 8 selected were not statistically significant. Further analysis also identified evidence of a sex-linked effect in SDHB and UQCC3 expression, and a difference in ATP5MD expression across individuals with MH sub-phenotypes (trigger from in vitro halothane exposure only, MHSh (n = 4); trigger to both in vitro halothane and caffeine exposure, MHShc (n = 4)). Our data support a link between MH-susceptibility and dysregulated gene expression associated with mitochondrial bioenergetics, which we speculate plays a role in the phenotypic variability observed within MH.

Comparison of Transcriptomic Changes in Survivors of Exertional Heat Illness with Malignant Hyperthermia Susceptible Patients.

Exertional heat illness (EHI) is an occupational health hazard for athletes and military personnel-characterised by the inability to thermoregulate during exercise. The ability to thermoregulate can be studied using a standardised heat tolerance test (HTT) developed by The Institute of Naval Medicine. In this study, we investigated whole blood gene expression (at baseline, 2 h post-HTT and 24 h post-HTT) in male subjects with either a history of EHI or known susceptibility to malignant hyperthermia (MHS): a pharmacogenetic condition with similar clinical phenotype. Compared to healthy controls at baseline, 291 genes were differentially expressed in the EHI cohort, with functional enrichment in inflammatory response genes (up to a four-fold increase). In contrast, the MHS cohort featured 1019 differentially expressed genes with significant down-regulation of genes associated with oxidative phosphorylation (OXPHOS). A number of differentially expressed genes in the inflammation and OXPHOS pathways overlapped between the EHI and MHS subjects, indicating a common underlying pathophysiology. Transcriptome profiles between subjects who passed and failed the HTT (based on whether they achieved a plateau in core temperature or not, respectively) were not discernable at baseline, and HTT was shown to elevate inflammatory response gene expression across all clinical phenotypes.

Bioenergetic defects in muscle fibers of RYR1 mutant knock-in mice associated with malignant hyperthermia.

Mutations in the skeletal muscle ryanodine receptor gene (RYR1) can cause susceptibility to malignant hyperthermia (MH), a potentially lethal genetic condition triggered by volatile anesthetics. MH is associated with hypermetabolism, which has directed research interest into oxidative phosphorylation and muscle bioenergetics. The most common cause of MH in the United Kingdom is the c.7300G>A RYR1 variant, which is present in ∼16% of MH families. Our study focuses on the MH susceptible G2435R-RYR1 knock-in mouse model, which is the murine equivalent of the human c.7300G>A genotype. Using a combination of transcriptomics, protein expression, and functional analysis, we investigated adult muscle fiber bioenergetics in this mouse model. RNA-Seq data showed reduced expression of genes associated with mitochondria and fatty acid oxidation in RYR1 mutants when compared with WT controls. Mitochondrial function was assessed by measuring oxygen consumption rates in permeabilized muscle fibers. Comparisons between WT and homozygous G2435R-RYR1 mitochondria showed a significant increase in complex I-facilitated oxidative phosphorylation in mutant muscle. Furthermore, we observed a gene-dose-specific increase in reactive oxygen species production in G2435R-RYR1 muscle fibers. Collectively, these findings provide evidence of metabolic defects in G2435R-RYR1 knock-in mouse muscle under basal conditions. Differences in metabolic profile could be the result of differential gene expression in metabolic pathways, in conjunction with mitochondrial damage accumulated from chronic exposure to increased oxidative stress.

The Malaria-Protective Human Glycophorin Structural Variant DUP4 Shows Somatic Mosaicism and Association with Hemoglobin Levels.

Glycophorin A and glycophorin B are red blood cell surface proteins and are both receptors for the parasite Plasmodium falciparum, which is the principal cause of malaria in sub-Saharan Africa. DUP4 is a complex structural genomic variant that carries extra copies of a glycophorin A-glycophorin B fusion gene and has a dramatic effect on malaria risk by reducing the risk of severe malaria by up to 40%. Using fiber-FISH and Illumina sequencing, we validate the structural arrangement of the glycophorin locus in the DUP4 variant and reveal somatic variation in copy number of the glycophorin B-glycophorin A fusion gene. By developing a simple, specific, PCR-based assay for DUP4, we show that the DUP4 variant reaches a frequency of 13% in the population of a malaria-endemic village in south-eastern Tanzania. We genotype a substantial proportion of that village and demonstrate an association of DUP4 genotype with hemoglobin levels, a phenotype related to malaria, using a family-based association test. Taken together, we show that DUP4 is a complex structural variant that may be susceptible to somatic variation and show that DUP4 is associated with a malarial-related phenotype in a longitudinally followed population.

Investigating the genetic susceptibility to exertional heat illness.

BACKGROUND: We aimed to identify rare (minor allele frequency ≤1%), potentially pathogenic non-synonymous variants in a well-characterised cohort with a clinical history of exertional heat illness (EHI) or exertional rhabdomyolysis (ER). The genetic link between malignant hyperthermia (MH) and EHI was investigated due to their phenotypic overlap. METHODS: The coding regions of 38 genes relating to skeletal muscle calcium homeostasis or exercise intolerance were sequenced in 64 patients (mostly military personnel) with a history of EHI, or ER and who were phenotyped using skeletal muscle in vitro contracture tests. We assessed the pathogenicity of variants using prevalence data, in silico analysis, phenotype and segregation evidence and by review of the literature. RESULTS: We found 51 non-polymorphic, potentially pathogenic variants in 20 genes in 38 patients. Our data indicate that RYR1 p.T3711M (previously shown to be likely pathogenic for MH susceptibility) and RYR1 p.I3253T are likely pathogenic for EHI. PYGM p.A193S was found in 3 patients with EHI, which is significantly greater than the control prevalence (p=0.000025). We report the second case of EHI in which a missense variant at CACNA1S p.R498 has been found. Combinations of rare variants in the same or different genes are implicated in EHI. CONCLUSION: We confirm a role of RYR1 in the heritability of EHI as well as ER but highlight the likely genetic heterogeneity of these complex conditions. We propose defects, or combinations of defects, in skeletal muscle calcium homeostasis, oxidative metabolism and membrane excitability are associated with EHI.

Permeabilised skeletal muscle reveals mitochondrial deficiency in malignant hyperthermia-susceptible individuals.

BACKGROUND: Individuals genetically susceptible to malignant hyperthermia (MH) exhibit hypermetabolic reactions when exposed to volatile anaesthetics. Mitochondrial dysfunction has previously been associated with the MH-susceptible (MHS) phenotype in animal models, but evidence of this in human MH is limited. METHODS: We used high resolution respirometry to compare oxygen consumption rates (oxygen flux) between permeabilised human MHS and MH-negative (MHN) skeletal muscle fibres with or without prior exposure to halothane. A substrate-uncoupler-inhibitor titration protocol was used to measure the following components of the electron transport chain under conditions of oxidative phosphorylation (OXPHOS) or after uncoupling the electron transport system (ETS): complex I (CI), complex II (CII), CI+CII and, as a measure of mitochondrial mass, complex IV (CIV). RESULTS: Baseline comparisons without halothane exposure showed significantly increased mitochondrial mass (CIV, P=0.021) but lower flux control ratios in CI+CII(OXPHOS) and CII(ETS) of MHS mitochondria compared with MHN (P=0.033 and 0.005, respectively) showing that human MHS mitochondria have a functional deficiency. Exposure to halothane triggered a hypermetabolic response in MHS mitochondria, significantly increasing mass-specific oxygen flux in CI(OXPHOS), CI+CII(OXPHOS), CI+CII(ETS), and CII(ETS) (P=0.001-0.012), while the rates in MHN samples were unaltered by halothane exposure. CONCLUSIONS: We present evidence of mitochondrial dysfunction in human MHS skeletal muscle both at baseline and after halothane exposure.

Mosaic maternal ancestry in the Great Lakes region of East Africa.

The Great Lakes lie within a region of East Africa with very high human genetic diversity, home of many ethno-linguistic groups usually assumed to be the product of a small number of major dispersals. However, our knowledge of these dispersals relies primarily on the inferences of historical, linguistics and oral traditions, with attempts to match up the archaeological evidence where possible. This is an obvious area to which archaeogenetics can contribute, yet Uganda, at the heart of these developments, has not been studied for mitochondrial DNA (mtDNA) variation. Here, we compare mtDNA lineages at this putative genetic crossroads across 409 representatives of the major language groups: Bantu speakers and Eastern and Western Nilotic speakers. We show that Uganda harbours one of the highest mtDNA diversities within and between linguistic groups, with the various groups significantly differentiated from each other. Despite an inferred linguistic origin in South Sudan, the data from the two Nilotic-speaking groups point to a much more complex history, involving not only possible dispersals from Sudan and the Horn but also large-scale assimilation of autochthonous lineages within East Africa and even Uganda itself. The Eastern Nilotic group also carries signals characteristic of West-Central Africa, primarily due to Bantu influence, whereas a much stronger signal in the Western Nilotic group suggests direct West-Central African ancestry. Bantu speakers share lineages with both Nilotic groups, and also harbour East African lineages not found in Western Nilotic speakers, likely due to assimilating indigenous populations since arriving in the region ~3000 years ago.

Next-generation Sequencing of RYR1 and CACNA1S in Malignant Hyperthermia and Exertional Heat Illness.

BACKGROUND: Variants in RYR1 are associated with the majority of cases of malignant hyperthermia (MH), a form of heat illness pharmacogenetically triggered by general anesthetics, and they have also been associated with exertional heat illness (EHI). CACNA1S has also been implicated in MH. The authors applied a targeted next-generation sequencing approach to identify variants in RYR1 and CACNA1S in a cohort of unrelated patients diagnosed with MH susceptibility. They also provide the first comprehensive report of sequencing of these two genes in a cohort of survivors of EHI. METHODS: DNA extracted from blood was genotyped using a "long" polymerase chain reaction technique, with sequencing on the Illumina GAII or MiSeq platforms (Illumina Inc., USA). Variants were assessed for pathogenicity using bioinformatic approaches. For further follow-up, DNA from additional family members and up to 211 MH normal and 556 MH-susceptible unrelated individuals was tested. RESULTS: In 29 MH patients, the authors identified three pathogenic and four novel RYR1 variants, with a further five RYR1 variants previously reported in association with MH. Three novel RYR1 variants were found in the EHI cohort (n = 28) along with two more previously reported in association with MH. Two other variants were reported previously associated with centronuclear myopathy. The authors found one and three rare variants of unknown significance in CACNA1S in the MH and EHI cohorts, respectively. CONCLUSIONS: Targeted next-generation sequencing proved efficient at identifying diagnostically useful and potentially implicated variants in RYR1 and CACNA1S in MH and EHI.

CCL3L1 copy number, CCR5 genotype and susceptibility to tuberculosis.

BACKGROUND: Tuberculosis is a major infectious disease and functional studies have provided evidence that both the chemokine MIP-1α and its receptor CCR5 play a role in susceptibility to TB. Thus by measuring copy number variation of CCL3L1, one of the genes that encode MIP-1α, and genotyping a functional promoter polymorphism -2459A > G in CCR5 (rs1799987) we investigate the influence of MIP-1α and CCR5, independently and combined, in susceptibility to clinically active TB in three populations, a Peruvian population (n = 1132), a !Xhosa population (n = 605) and a South African Coloured population (n = 221). The three populations include patients with clinically diagnosed pulmonary TB, as well as other, less prevalent forms of extrapulmonary TB. METHODS AND RESULTS: Copy number of CCL3L1 was measured using the paralogue ratio test and exhibited ranges between 0-6 copies per diploid genome (pdg) in Peru, between 0-12 pdg in !Xhosa samples and between 0-10 pdg in South African Coloured samples. The CCR5 promoter polymorphism was observed to differ significantly in allele frequency between populations (*A; Peru f = 0.67, !Xhosa f = 0.38, Coloured f = 0.48). CONCLUSIONS: The case-control association studies performed however find, surprisingly, no evidence for an influence of variation in genes coding for MIP-1α or CCR5 individually or together in susceptibility to clinically active TB in these populations.

Variants in ASPH cause exertional heat illness and are associated with malignant hyperthermia susceptibility.

Exertional heat illness (EHI) and malignant hyperthermia (MH) are life threatening conditions associated with muscle breakdown in the setting of triggering factors including volatile anesthetics, exercise, and high environmental temperature. To identify new genetic variants that predispose to EHI and/or MH, we performed genomic sequencing on a cohort with EHI/MH and/or abnormal caffeine-halothane contracture test. In five individuals, we identified rare, pathogenic heterozygous variants in ASPH, a gene encoding junctin, a regulator of excitation-contraction coupling. We validated the pathogenicity of these variants using orthogonal pre-clinical models, CRISPR-edited C2C12 myotubes and transgenic zebrafish. In total, we demonstrate that ASPH variants represent a new cause of EHI and MH susceptibility.

Karyomapping: a universal method for genome wide analysis of genetic disease based on mapping crossovers between parental haplotypes.

The use of genome wide single nucleotide polymorphism (SNP) arrays for high resolution molecular cytogenetic analysis using a combination of quantitative and genotype analysis is well established. This study demonstrates that by Mendelian analysis of the SNP genotypes of the parents and a sibling or other appropriate family member to establish phase, it is possible to identify informative loci for each of the four parental haplotypes across each chromosome and map the inheritance of these haplotypes and the position of any crossovers in the proband. The resulting 'karyomap', unlike a karyotype, identifies the parental and grandparental origin of each chromosome and chromosome segment and is unique for every individual being defined by the independent segregation of parental chromosomes and the pattern of non-recombinant and recombinant chromosomes. Karyomapping, therefore, enables both genome wide linkage based analysis of inheritance and detection of chromosome imbalance where either both haplotypes from one parent are present (trisomy) or neither are present (monosomy/deletion). The study also demonstrates that karyomapping is possible at the single cell level following whole genome amplification and, without any prior patient or disease specific test development, provides a universal linkage based methodology for preimplantation genetic diagnosis readily available worldwide.

One year into the pandemic: Short-term evolution of SARS-CoV-2 and emergence of new lineages.

The COVID-19 pandemic was officially declared on March 11th, 2020. Since the very beginning, the spread of the virus has been tracked nearly in real-time by worldwide genome sequencing efforts. As of March 2021, more than 830,000 SARS-CoV-2 genomes have been uploaded in GISAID and this wealth of data allowed researchers to study the evolution of SARS-CoV-2 during this first pandemic year. In parallel, nomenclatures systems, often with poor consistency among each other, have been developed to designate emerging viral lineages. Despite general fears that the virus might mutate to become more virulent or transmissible, SARS-CoV-2 genetic diversity has remained relatively low during the first ~ 8 months of sustained human-to-human transmission. At the end of 2020/beginning of 2021, though, some alarming events started to raise concerns of possible changes in the evolutionary trajectory of the virus. Specifically, three new viral variants associated with extensive transmission have been described as variants of concern (VOC). These variants were first reported in the UK (B.1.1.7), South Africa (B.1.351) and Brazil (P.1). Their designation as VOCs was determined by an increase of local cases and by the high number of amino acid substitutions harboured by these lineages. This latter feature is reminiscent of viral sequences isolated from immunocompromised patients with long-term infection, suggesting a possible causal link. Here we review the events that led to the identification of these lineages, as well as emerging data concerning their possible implications for viral phenotypes, reinfection risk, vaccine efficiency and epidemic potential. Most of the available evidence is, to date, provisional, but still represents a starting point to uncover the potential threat posed by the VOCs. We also stress that genomic surveillance must be strengthened, especially in the wake of the vaccination campaigns.

Single Amino Acid Changes in the Ryanodine Receptor in the Human Population Have Effects In Vivo on Caenorhabditis elegans Neuro-Muscular Function.

The ryanodine receptor mediates intracellular calcium ion release with excitation of nerve and muscle cells. Ryanodine receptor missense variants cause a number of myopathologies, such as malignant hyperthermia, and have been linked with various neuropathologies, including Alzheimer's disease. We characterized the consequences of ryanodine receptor variants in vivo. Eight Caenorhabditis elegans strains, with ryanodine receptor modifications equivalent to human myopathic RYR1 variants, were generated by genome editing. In humans, these variants are rare and confer sensitivity to the inhalational anaesthetic halothane when heterozygous. Increased sensitivity to halothane was found in both homozygous and heterozygous C. elegans. Close analysis revealed distinct subtle locomotion defects, due to the different single amino acid residue changes, even in the absence of the external triggering agent. Distinct pre- and postsynaptic consequences of the variants were characterized through the responses to cholinergic pharmacological agents. The range of phenotypes reflects the complexity of the regulatory inputs to the ryanodine receptor and the criticality of the calcium ion channel opening properties, in different cell types and with age. Ryanodine receptors with these single amino acid residue changes still function as calcium ion channels, but with altered properties which are likely to have subtle consequences for human carriers of such variants. The long-term consequences of subtly altered calcium ion signalling could be cumulative and may be focussed in the smaller nerve cells rather than the more robust muscle cells. It was important to assess phenotypes in vivo to properly appreciate consequences for a whole organism.

SARS-CoV-2 and COVID-19: A genetic, epidemiological, and evolutionary perspective.

In less than five months, COVID-19 has spread from a small focus in Wuhan, China, to more than 5 million people in almost every country in the world, dominating the concern of most governments and public health systems. The social and political distresses caused by this epidemic will certainly impact our world for a long time to come. Here, we synthesize lessons from a range of scientific perspectives rooted in epidemiology, virology, genetics, ecology and evolutionary biology so as to provide perspective on how this pandemic started, how it is developing, and how best we can stop it.

The role of CACNA1S in predisposition to malignant hyperthermia.

BACKGROUND: Malignant hyperthermia (MH) is an inherited pharmacogenetic disorder of skeletal muscle, characterised by an elevated calcium release from the skeletal muscle sarcoplasmic reticulum. The dihydropyridine receptor (DHPR) plays an essential role in excitation-contraction coupling and calcium homeostasis in skeletal muscle. This study focuses on the gene CACNA1S which encodes the alpha1 subunit of the DHPR, in order to establish whether CACNA1S plays a major role in MH susceptibility in the UK. METHODS: We investigate the CACNA1S locus in detail in 50 independent MH patients, the largest study to date, to identify novel variants that may predispose to disease and also to characterise the haplotype structure across CACNA1S. RESULTS: We present CACNA1S cDNA sequencing data from 50 MH patients in whom RYR1 mutations have been excluded, and subsequent mutation screening analysis. Furthermore we present haplotype analysis of unphased CACNA1S SNPs to (1) assess CACNA1S haplotype frequency differences between susceptible MH cases and a European control group and (2) analyse population-based association via clustering of CACNA1S haplotypes based on disease risk. CONCLUSION: The study identified a single potentially pathogenic change in CACNA1S (p.Arg174Trp), and highlights that the haplotype structure across CACNA1S is diverse, with a high degree of variability.

Genetics of susceptibility to malaria related phenotypes.

Previous studies have established a genetic component for susceptibility to malaria. Here we use a pedigree based approach, and transmission disequilibrium testing (TDT), to identify immune response genes that influence susceptibility to Plasmodium falciparum malarial phenotypes (parasite density and frequency of clinical episodes) in a Tanzanian population. Evidence for association was observed between markers in the TNF gene cluster and both the malarial phenotypes. There was weaker evidence for associations between HLA-DRB1*04, HLA-DRB1*10, and loci in the TCRBV region with parasite density. There was no evidence for association with polymorphisms in the IL10 promoter, IL1 gene cluster, or from the IL4/IL13 region.

Pathogen evolution and disease emergence in carnivores.

Emerging infectious diseases constitute some of the most pressing problems for both human and domestic animal health, and biodiversity conservation. Currently it is not clear whether the removal of past constraints on geographical distribution and transmission possibilities for pathogens alone are sufficient to give rise to novel host-pathogen combinations, or whether pathogen evolution is also generally required for establishment in novel hosts. Canine distemper virus (CDV) is a morbillivirus that is prevalent in the world dog population and poses an important conservation threat to a diverse range of carnivores. We performed an extensive phylogenetic and molecular evolution analysis on complete sequences of all CDV genes to assess the role of selection and recombination in shaping viral genetic diversity and driving the emergence of CDV in non-dog hosts. We tested the specific hypothesis that molecular adaptation at known receptor-binding sites of the haemagglutinin gene is associated with independent instances of the spread of CDV to novel non-dog hosts in the wild. This hypothesis was upheld, providing compelling evidence that repeated evolution at known functional sites (in this case residues 530 and 549 of the haemagglutinin molecule) is associated with multiple independent occurrences of disease emergence in a range of novel host species.

Immunogenetic control of antibody responsiveness in a malaria endemic area.

This study builds upon the established genetic control of antimalarial immune responses and prior association studies by using a family-based approach, transmission disequilibrium testing, to identify immune response genes that influence antibody responses to Plasmodium falciparum infection in an endemic Tanzanian population. Candidate polymorphisms are within the interleukin-1 (IL-1) gene cluster, the IL-10 promoter, Major histocompatibility complex class II and III, the 5q31-q33 region, and the T-Cell Receptor beta variable region. There was a significant association between the IL1RN alleles and total IgE. Weak evidence for association was present between polymorphisms in the IL10 promoter region and both anti-P falciparum IgE and IgG4 antibodies.