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.

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.

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.

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.

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.

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.

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.

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.

Aging Effects of Caenorhabditis elegans Ryanodine Receptor Variants Corresponding to Human Myopathic Mutations.

Delaying the decline in skeletal muscle function will be critical to better maintenance of an active lifestyle in old age. The skeletal muscle ryanodine receptor, the major intracellular membrane channel through which calcium ions pass to elicit muscle contraction, is central to calcium ion balance and is hypothesized to be a significant factor for age-related decline in muscle function. The nematode Caenorhabditis elegans is a key model system for the study of human aging, and strains were generated with modified C. elegans ryanodine receptors corresponding to human myopathic variants linked with malignant hyperthermia and related conditions. The altered response of these strains to pharmacological agents reflected results of human diagnostic tests for individuals with these pathogenic variants. Involvement of nerve cells in the C. elegans responses may relate to rare medical symptoms concerning the central nervous system that have been associated with ryanodine receptor variants. These single amino acid modifications in C. elegans also conferred a reduction in lifespan and an accelerated decline in muscle integrity with age, supporting the significance of ryanodine receptor function for human aging.

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.

The heritability of abortion in pedigree Charollais flocks.

Foetal death, or abortion at term, in sheep is of major significance to the livestock industry, accounting for more than £24million lost per annum. We have investigated whether there is a genetic component to abortion within two flocks of pedigree Charollais sheep, one followed from 1989 to 2006, the other from 1992 to 2006. Abortion occurred at a rate of 5.74-8.78% per annum against a total mortality rate of 14-24%. By model covariate analysis we have shown that 15.5% aborting ewes went on to have one or more abortions and that this risk increased with parity (p=0.006). Heritability estimates were approximately 0.08 as calculated by SOLAR, pedigreemm and ASReml3, with sire and dam components of 0.046 and 0.048, respectively. Where the lamb was aborted, heritability estimates were highly variable according to the method employed, 0.046-0.378, with sex of the lamb being a significant covariate. This variability indicated one or more underlying, significant factors that were not measured in these analyses, potentially including infectious agents that may be involved. Nevertheless, the ASReml3 estimate (0.179) resolved to 0.074 variance attributable to the sire and 0.092 attributable to the dam, which, while not significant, was suggestive that genetic variants passed by the dam to the lamb may be of more weight than that from the sire in determining whether a lamb will abort.

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.

Genetic diversity and transmission characteristics of Beijing family strains of Mycobacterium tuberculosis in Peru.

Beijing family strains of Mycobacterium tuberculosis have attracted worldwide attention because of their wide geographical distribution and global emergence. Peru, which has a historical relationship with East Asia, is considered to be a hotspot for Beijing family strains in South America. We aimed to unveil the genetic diversity and transmission characteristics of the Beijing strains in Peru. A total of 200 Beijing family strains were identified from 2140 M. tuberculosis isolates obtained in Lima, Peru, between December 2008 and January 2010. Of them, 198 strains were classified into sublineages, on the basis of 10 sets of single nucleotide polymorphisms (SNPs). They were also subjected to variable number tandem-repeat (VNTR) typing using an international standard set of 15 loci (15-MIRU-VNTR) plus 9 additional loci optimized for Beijing strains. An additional 70 Beijing family strains, isolated between 1999 and 2006 in Lima, were also analyzed in order to make a longitudinal comparison. The Beijing family was the third largest spoligotyping clade in Peru. Its population structure, by SNP typing, was characterized by a high frequency of Sequence Type 10 (ST10), which belongs to a modern subfamily of Beijing strains (178/198, 89.9%). Twelve strains belonged to the ancient subfamily (ST3 [n=3], ST25 [n=1], ST19 [n=8]). Overall, the polymorphic information content for each of the 24 loci values was low. The 24 loci VNTR showed a high clustering rate (80.3%) and a high recent transmission index (RTI(n-1)=0.707). These strongly suggest the active and on-going transmission of Beijing family strains in the survey area. Notably, 1 VNTR genotype was found to account for 43.9% of the strains. Comparisons with data from East Asia suggested the genotype emerged as a uniquely endemic clone in Peru. A longitudinal comparison revealed the genotype was present in Lima by 1999.

CCL3L1 copy number and susceptibility to malaria.

Copy number variation can contribute to the variation observed in susceptibility to complex diseases. Here we present the first study to investigate copy number variation of the chemokine gene CCL3L1 with susceptibility to malaria. We present a family-based genetic analysis of a Tanzanian population (n=922), using parasite load, mean number of clinical infections of malaria and haemoglobin levels as phenotypes. Copy number of CCL3L1 was measured using the paralogue ratio test (PRT) and the dataset exhibited copy numbers ranging between 1 and 10 copies per diploid genome (pdg). Association between copy number and phenotypes was assessed. Furthermore, we were able to identify copy number haplotypes in some families, using microsatellites within the copy variable region, for transmission disequilibrium testing. We identified a high level of copy number haplotype diversity and find some evidence for an association of low CCL3L1 copy number with protection from anaemia.

Inflammatory response of porcine epithelial IPEC J2 cells to enterotoxigenic E. coli infection is modulated by zinc supplementation.

Enterotoxigenic Escherichia coli (ETEC) is a major cause of diarrhoea in pigs and humans. The duration and severity of diarrhoea can be controlled using zinc supplementation, typically pharmacological levels of zinc oxide in pigs. In this study, IPEC J2 cells were used as an in vitro model of intestinal ETEC infection, with separate and simultaneous zinc treatment. Genomic analysis identified increased expression of a variety of innate immune response genes (NF-κB targets) in response to ETEC exposure, and several stress response genes in response to zinc exposure, provided as ZnO. Expression of genes involved in the innate immune response was reduced when cells were simultaneously exposed to ZnO, and it is suggested that ZnO treatment inhibits the induction of NF-κB in response to pathogens, possibly through up-regulated heat shock proteins. A similar response in vivo with consequent down-regulation in the inflammatory response would reduce further pathogen invasion, maintain normal gut function and maintain growth.