Output list
Preprint
Posted to a preprint site 28/05/2025
medRxiv
Co-trimoxazole is a leading global cause of severe cutaneous adverse drug reactions (SCAR) including Stevens-Johnson syndrome/toxic epidermal necrolysis (SJS/TEN) and drug reaction with eosinophilia and systemic symptoms (DRESS). Co-trimoxazole-induced SCAR are associated with HLA class I alleles including HLA-B*13:01 and HLA-B*38:02 in Southeast Asian (SEA) populations. However, the global generalizability of these associations is unknown but critical for population-appropriate risk stratification and diagnosis.
To determine HLA risk factors associated with co-trimoxazole-induced SJS/TEN and DRESS in populations from the United States (US) and South Africa (SA).
We performed high-resolution HLA typing on dermatologist-adjudicated co-trimoxazole-induced SCAR patients in the US (n=63) and SA (n=26) compared to population controls. Peptide binding and docking analyses were performed using MHCcluster2.0 and CB-Dock2.
In a multiple logistic regression model, HLA-B*44:03 (Pc<0.001, OR: 4.08), HLA-B*38:01 (Pc<0.001, OR: 5.66), and HLA-C*04:01 (Pc=0.003, OR: 2.50) were independently associated with co-trimoxazole-induced SJS/TEN in the US. HLA-B*44:03 was also associated with co-trimoxazole-induced DRESS in SA (Pc=0.019, OR: 10.69). Distinct HLA-B variants with shared peptide binding specificities (SPBS) and HLA-C*04:01 identified 94% and 78% of co-trimoxazole-induced SJS/TEN and DRESS in the US, respectively. The SEA risk allele HLA-B*13:01, with SPBS to HLA-B*44:03, was identified in just 1/63 US SCAR patients.
HLA alleles with SPBS to SEA-related risk alleles including HLA-B*44:03 (SPBS with HLA-B*13:01) and HLA-B*38:01 (SPBS with HLA-B*38:02) but also HLA-C*04:01 predisposed to co-trimoxazole-induced SCAR in the US and SA. These findings provide biological plausibility and strategies for global risk prediction and diagnosis of co-trimoxazole-induced SCAR.
HLA alleles including HLA-B*13:01 and HLA-B*38:02 are risk factors for co-trimoxazole-induced SCAR in Asian populations. However, the generalizability of these associations to other global populations is unknown but critical for population-appropriate risk stratification and diagnosis.
HLA alleles with shared peptide binding specificities (SPBS) to Asian-related risk alleles including HLA-B*44:03 (SPBS with HLA-B*13:01) and HLA-B*38:01 (SPBS with HLA-B*38:02) but also HLA-C*04:01 predisposed to co-trimoxazole-induced SCAR in the US and South Africa.
HLA alleles previously associated with co-trimoxazole-induced SCAR do not identify risk across populations. However, HLA alleles with SPBS provide biological plausibility and strategies for global and population-appropriate clinical risk stratification and diagnosis of cotrimoxazole-induced SCAR.
Journal article
HIV-1 adapts to HLA class II-associated selection pressure exerted by CD4 + and CD8 + T cells
Published 2025
Science advances, 11, 7, eadr4238
Developing an effective HIV-1 vaccine is a global health priority, but HIV-1 mutational escape from T cells poses a challenge. While escape from human leukocyte antigen class I (HLA-I)–restricted CD8+ T cells is well characterized, less is known about HLA-II–restricted T cell escape. We used computational methods to identify 149 sites across the HIV-1 clade B genome under HLA-II–associated selection. Functional assays, including activation-induced intracellular cytokine staining and enzyme-linked immunospot for interferon-γ, revealed diverse mechanisms of HIV-1 adaptation to HLA-II–associated immune pressure, ranging from loss to sustained antigen recognition. T cell receptor and RNA sequencing demonstrated variable clonotype overlap of T cell clones to recognize adapted versus non-adapted peptides, with cells targeting adapted peptides exhibiting a dysfunctional transcriptomic state. Moreover, incorporating HLA-II–associated adaptation strengthened the correlation between Gag-specific viral adaptation and poor disease outcomes. Last, we mapped viral regions prone to HLA-II–associated adaptation and found that these adaptations can increase in frequency within populations.
Journal article
Published 2023
Journal of autoimmunity, 142, 103150
OBJECTIVES: Inclusion body myositis (IBM) is a progressive inflammatory-degenerative muscle disease of older individuals, with some patients producing anti-cytosolic 5'-nucleotidase 1A (NT5C1A, aka cN1A) antibodies. Human Leukocyte Antigens (HLA) is the highest genetic risk factor for developing IBM. In this study, we aimed to further define the contribution of HLA alleles to IBM and the production of anti-cN1A antibodies.
METHODS: We HLA haplotyped a Western Australian cohort of 113 Caucasian IBM patients and 112 ethnically matched controls using Illumina next-generation sequencing. Allele frequency analysis and amino acid alignments were performed using the Genentech/MiDAS bioinformatics package. Allele frequencies were compared using Fisher's exact test. Age at onset analysis was performed using the ggstatsplot package. All analysis was carried out in RStudio version 1.4.1717.
RESULTS: Our findings validated the independent association of HLA-DRB1*03:01:01 with IBM and attributed the risk to an arginine residue in position 74 within the DRβ1 protein. Conversely, DRB4*01:01:01 and DQA1*01:02:01 were found to have protective effects; the carriers of DRB1*03:01:01 that did not possess these alleles had a fourteenfold increased risk of developing IBM over the general Caucasian population. Furthermore, patients with the abovementioned genotype developed symptoms on average five years earlier than patients without. We did not find any HLA associations with anti-cN1A antibody production.
CONCLUSIONS: High-resolution HLA sequencing more precisely characterised the alleles associated with IBM and defined a haplotype linked to earlier disease onset. Identification of the critical amino acid residue by advanced biostatistical analysis of immunogenetics data offers mechanistic insights and future directions into uncovering IBM aetiopathogenesis.
Conference presentation
A SQSTM1 polymorphism confers risk for sporadic inclusion body myositis disease expression
Date presented 06/2022
4th Global Conference on Myositis (GCOM), 06/06/2022–09/06/2022, Prague, Czech Republic
Background. Inclusions of Sequestosome1 (SQSTM1/p62) within muscle fibers are a pathological hallmark of sporadic inclusion body myositis (sIBM) with p62 overexpression reported in patients. Mounting evidence suggests a role for p62 expression and/or variation in sIBM pathology, due to the presence of rare and potentially pathogenic missense variants (A117V, G194R, K238E, P392L). Consequently, we hypothesized that genetic modifiers of SQSTM1 may present a critical missing link for sIBM pathology and contribute to disease expres¬sion. Short structural variants (SSVs) are a class of genetic variation that can be difficult to characterized due to their highly repetitive and complex nature. Evidence that SSVs play an important role in complex diseases such as Alzhei¬mer’s Disease, Amyotrophic lateral sclerosis, Spinocerebellar Ataxia type 2, and Huntington’s disease is now confirmed and further investigations of this type of genetic variation is necessary to uncover missing heritability in complex dis¬eases. We and others have previously reported an SSV within SQSTM1 that is associated with altered expression of p62. The SSV rs60327661 is a CAAA in¬sertion/deletion within intron 5 of SQSTM1, which also confers risk for familial Amyotrophic lateral sclerosis. Due to the role of the SSV in Amyotrophic lateral sclerosis and altered p62 expression, we hypothesized the SSV rs60327661 may have disease-modifying effects in a longitudinal cohort of sIBM patients.
Methods. DNA samples from 218 sIBM patients and 242 healthy controls were received from The Institute for Immunology and Infectious Diseases, Murdoch, Western Australia, and the NINDS Repository, Coriell Institute for Medical Research, New Jersey. Genomic DNA samples were systematically assessed through polymerase chain reaction, capillary separation, and Sanger sequencing for rs60327661 allele genotyping.
Results: In the present study, when controlling for self-declared ancestry, car¬riage of the D/D genotype is associated with sIBM disease expression (p<0.05). Both the case and control groups did not violate Hardy-Weinberg equilibrium (p=0.99, p=0.98; respectively. Intriguingly, patients who were CN1A seroposi¬tive were more likely carry the D allele (n=18) when compared to patients with-out a D allele (n=3; p<0.047). Patients classified as fast progressors (n=2) carried only the D/D genotype.
Conclusion. In this study, we present the first report of an association between the SQSTM1 insertion/deletion and sIBM disease expression. We provide evi¬dence that the investigation of genetic variants outside of the HLA region is war¬ranted, and that such investigations are likely to uncover critical information for sIBM. We present the SSV rs60327661 as a novel disease modifying variant for sIBM which is functionally linked to p62 by altering protein expression. Our data adds to growing evidence that examination of SSVs may uncover novel genetic risk markers, and consequently further understanding of the pathogenic mechanisms at play.
Journal article
Published 2022
Scientific Reports, 12, 1, Art. 14739
Neurofilament heavy (NEFH) is one of the critical proteins required for the formation of the neuronal cytoskeleton and polymorphisms in NEFH are reported as a rare cause of sporadic ALS (sALS). In the current study, a candidate tetranucleotide (TTTA) repeat variant in NEFH was selected using an in-silico short structural variant (SSV) evaluation algorithm and investigated in two cohorts of North American sALS patients, both separately and combined (Duke cohort n = 138, Coriell cohort n = 333; combined cohort n = 471), compared to a group of healthy controls from the Coriell Institute biobank (n = 496). Stratification according to site of disease onset revealed that the 9 TTTA allele was associated with reduced disease risk, specifically confined to spinal-onset sALS patients in the Duke cohort (p = 0.001). Furthermore, carriage of the 10 TTTA allele was associated with a 2.7 year later age of disease onset in the larger combined sALS cohort (p = 0.02). These results suggest that the 9 and 10 TTTA motif length may have a protective advantage for potentially lowering the risk of sALS and delaying the age of disease onset, however, these results need to be replicated in larger multicenter and multi-ethnic cohorts.
Journal article
Published 2021
Pathogens and Immunity, 6, 2, 27 - 49
Background: Genetic variations across the SARS-CoV-2 genome may influence transmissibility of the virus and the host’s anti-viral immune response, in turn affecting the frequency of variants over time. In this study, we examined the adjacent amino acid polymorphisms in the nucleocapsid (R203K/G204R) of SARS-CoV-2 that arose on the background of the spike D614G change and describe how strains harboring these changes became dominant circulating strains globally. Methods: Deep-sequencing data of SARS-CoV-2 from public databases and from clinical samples were analyzed to identify and map genetic variants and sub-genomic RNA transcripts across the genome. Results: Sequence analysis suggests that the 3 adjacent nucleotide changes that result in the K203/R204 variant have arisen by homologous recombination from the core sequence of the leader transcription-regulating sequence (TRS) rather than by stepwise mutation. The resulting sequence changes generate a novel sub-genomic RNA transcript for the C-terminal dimerization domain of nucleocapsid. Deep-sequencing data from 981 clinical samples confirmed the presence of the novel TRS-CS-dimerization domain RNA in individuals with the K203/R204 variant. Quantification of sub-genomic RNA indicates that viruses with the K203/R204 variant may also have increased expression of sub-genomic RNA from other open reading frames. Conclusions: The finding that homologous recombination from the TRS may have occurred since the introduction of SARS-CoV-2 in humans, resulting in both coding changes and novel sub-genomic RNA transcripts, suggests this as a mechanism for diversification and adaptation within its new host.
Other
Published 2021
bioRxiv: The preprint server for Biology
Background Genetic variations across the SARS-CoV-2 genome may influence transmissibility of the virus and the host’s anti-viral immune response, in turn affecting the frequency of variants over-time. In this study, we examined the adjacent amino acid polymorphisms in the nucleocapsid (R203K/G204R) of SARS-CoV-2 that arose on the background of the spike D614G change and describe how strains harboring these changes became dominant circulating strains globally. Methods Deep sequencing data of SARS-CoV-2 from public databases and from clinical samples were analyzed to identify and map genetic variants and sub-genomic RNA transcripts across the genome. Results Sequence analysis suggests that the three adjacent nucleotide changes that result in the K203/R204 variant have arisen by homologous recombination from the core sequence (CS) of the leader transcription-regulating sequence (TRS) rather than by stepwise mutation. The resulting sequence changes generate a novel sub-genomic RNA transcript for the C-terminal dimerization domain of nucleocapsid. Deep sequencing data from 981 clinical samples confirmed the presence of the novel TRS-CS-dimerization domain RNA in individuals with the K203/R204 variant. Quantification of sub-genomic RNA indicates that viruses with the K203/R204 variant may also have increased expression of sub-genomic RNA from other open reading frames. Conclusions The finding that homologous recombination from the TRS may have occurred since the introduction of SARS-CoV-2 in humans resulting in both coding changes and novel sub-genomic RNA transcripts suggests this as a mechanism for diversification and adaptation within its new host.
Journal article
Published 2021
Thrombosis and Haemostasis, 121, 11, 1483 - 1496
Background: High estradiol (E2) levels are linked to an increased risk of venous thromboembolism; however, the underlying molecular mechanism(s) remain poorly understood. We previously identified an E2-responsive microRNA (miR), miR-494–3p, that downregulates protein S expression, and posited additional coagulation factors, such as tissue factor, may be regulated in a similar manner via miRs. Objectives: To evaluate the coagulation capacity of cohorts with high physiological E2, and to further characterize novel E2-responsive miR and miR regulation on tissue factor in E2-related hypercoagulability. Methods: Ceveron Alpha thrombin generation assay (TGA) was used to assess plasma coagulation profile of three cohorts. The effect of physiological levels of E2, 10 nM, on miR expression in HuH-7 cells was compared using NanoString nCounter and validated with independent assays. The effect of tissue factor-interacting miR was confirmed by dual-luciferase reporter assays, immunoblotting, flow cytometry, biochemistry assays, and TGA. Results: Plasma samples from pregnant women and women on the contraceptive pill were confirmed to be hypercoagulable (compared with sex-matched controls). At equivalent and high physiological levels of E2, miR-365a-3p displayed concordant E2 downregulation in two independent miR quantification platforms, and tissue factor protein was upregulated by E2 treatment. Direct interaction between miR-365a-3p and F3-3′UTR was confirmed and overexpression of miR-365a-3p led to a decrease of (1) tissue factor mRNA transcripts, (2) protein levels, (3) activity, and (4) tissue factor-initiated thrombin generation. Conclusion: miR-365a-3p is a novel tissue factor regulator. High E2 concentrations induce a hypercoagulable state via a miR network specific for coagulation factors.
Journal article
Published 2021
Tuberculosis, 126, Art. 102023
Background To determine the association of human leukocyte antigen (HLA) alleles as correlates of risk for and protection against tuberculin skin test (TST) positivity and active TB disease amongst HIV-infected adults. Methods Genomic DNA was extracted from 754 HIV-infected adults whole-blood. HLA-A, -B, -C and -DRB1 loci were genotyped by next generation sequencing methods. HLA alleles were analysed by the presence/absence of TST immune conversion and active TB disease and further stratified by exposure to a household TB contact, CD4 T-cell count and, for active TB disease, TST-positivity. Results HLA-A*29:11 and - B*45:01/07 were associated with TST-positivity, while HLA-A*24:02, -A*29:02 and -B*15:16 with TST-negativity. In participants with a household TB contact, HLA-A*66:01, -A*68:02 and -B*49:01 were associated with TST-negativity. For TB disease, HLA-B*41:01, -C*06:02, -DRB1*04:01 and -DRB1*15:01 were associated with susceptibility, while HLA-B*07:02 and -DRB1*11:01 were protective, even for CD4 T-cell count <350 cells/mm3. For initial TST-positivity and subsequent TB disease, HLA-A*01:01 and -DRB1*11:01 conveyed protection including those with CD4 T-cell count <350 cells/mm3. Conclusion Several HLA alleles are noted as correlates of TB infection, risk and natural protection in HIV-infected individuals. HLA associations may enable risk stratification of those with HIV infection. Protective alleles may assist in future TB vaccine development.
Journal article
Novel STMN2 variant linked to amyotrophic lateral sclerosis risk and clinical phenotype
Published 2021
Frontiers in Aging Neuroscience, 13, 658226
Objective: There is a critical need to establish genetic markers that explain the complex phenotypes and pathogenicity of ALS. This study identified a polymorphism in the Stathmin-2 gene and investigated its association with sporadic ALS (sALS) disease risk, age-of onset and survival duration.
Methods: The candidate CA repeat was systematically analyzed using PCR, Sanger sequencing and high throughput capillary separation for genotyping. Stathmin-2 expression was investigated using RT-PCR in patient olfactory neurosphere-derived (ONS) cells and RNA sequencing in laser-captured spinal motor neurons.
Results: In a case-control analysis of a combined North American sALS cohort (n = 321) and population control group (n = 332), long/long CA genotypes were significantly associated with disease risk (p = 0.042), and most strongly when one allele was a 24 CA repeat (p = 0.0023). In addition, longer CA allele length was associated with earlier age-of-onset (p = 0.039), and shorter survival duration in bulbar-onset cases (p = 0.006). In an Australian longitudinal sALS cohort (n = 67), ALS functional rating scale scores were significantly lower in carriers of the long/long genotype (p = 0.034). Stathmin-2 mRNA expression was reduced in sporadic patient ONS cells. Additionally, sALS patients and controls exhibited variable expression of Stathmin-2 mRNA according to CA genotype in laser-captured spinal motor neurons.
Conclusions: We report a novel non-coding CA repeat in Stathmin-2 which is associated with sALS disease risk and has disease modifying effects. The potential value of this variant as a disease marker and tool for cohort enrichment in clinical trials warrants further investigation.