Q fever is a zoonotic disease caused by the highly infectious Gram-negative coccobacillus, Coxiella burnetii (C. burnetii). The Q fever vaccine Q-VAX® is characterised by high reactogenicity, requiring individuals to be pre-screened for prior exposure before vaccination. To date it remains unclear whether vaccine side effects in pre-exposed individuals are associated with pre-existing adaptive immune responses to C. burnetii or are also a function of innate responses to Q-VAX®. In the current study, we measured innate and adaptive cytokine responses to C. burnetii and compared these among individuals with different pre-exposure status. Three groups were included: n=98 Dutch blood bank donors with unknown exposure status, n=95 Dutch village inhabitants with known natural exposure status to C. burnetii during the Dutch Q fever outbreak of 2007-2010, and n=96 Australian students receiving Q-VAX® vaccination in 2021. Whole blood cytokine responses following ex vivo stimulation with heat-killed C. burnetii were assessed for IFNγ, IL-2, IL-6, IL-10, TNFα, IL-1β, IP-10, MIP-1α and IL-8. Serological data were collected for all three cohorts, as well as data on skin test and self-reported vaccine side effects and clinical symptoms during past infection. IFNγ, IP-10 and IL-2 responses were strongly elevated in individuals with prior C. burnetii antigen exposure, whether through infection or vaccination, while IL-1β, IL-6 and TNFα responses were slightly increased in naturally exposed individuals only. High dimensional analysis of the cytokine data identified four clusters of individuals with distinct cytokine response signatures. The cluster with the highest levels of adaptive cytokines and antibodies comprised solely individuals with prior exposure to C. burnetii, while another cluster was characterized by high innate cytokine production and an absence of C. burnetii-induced IP-10 production paired with high baseline IP-10 levels. Prior exposure status was partially associated with these signatures, but could not be clearly assigned to a single cytokine response signature. Overall, Q-VAX® vaccination and natural C. burnetii infection were associated with comparable cytokine response signatures, largely driven by adaptive cytokine responses. Neither individual innate and adaptive cytokine responses nor response signatures were associated retrospectively with clinical symptoms during infection or prospectively with side effects post-vaccination.

Hereditary haemorrhagic telangiectasia (HHT) can result in challenging anaemia and thrombosis phenotypes. Clinical presentations of HHT vary for relatives with identical casual mutations, suggesting other factors may modify severity. To examine objectively, we developed unsupervised machine learning algorithms to test whether haematological data at presentation could be categorised into sub‐groupings and fitted to known biological factors. With ethical approval, we examined 10 complete blood count (CBC) variables, four iron index variables, four coagulation variables and eight iron/coagulation indices combined from 336 genotyped HHT patients (40% male, 60% female, 86.5% not using iron supplementation) at a single centre. T‐SNE unsupervised, dimension reduction, machine learning algorithms assigned each high‐dimensional datapoint to a location in a two‐dimensional plane. k‐Means clustering algorithms grouped into profiles, enabling visualisation and inter‐profile comparisons of patients’ clinical and genetic features. The unsupervised machine learning algorithms using t‐SNE and k‐Means identified two distinct CBC profiles, two iron profiles, four clotting profiles and three combined profiles. Validating the methodology, profiles for CBC or iron indices fitted expected patterns for haemorrhage. Distinct coagulation profiles displayed no association with age, sex, C‐reactive protein, pulmonary arteriovenous malformations (AVMs), ENG/ACVRL1 genotype or epistaxis severity. The most distinct profiles were from t‐SNE/k‐Means analyses of combined iron‐coagulation indices and mapped to three risk states – for venous thromboembolism in HHT; for ischaemic stroke attributed to paradoxical emboli through pulmonary AVMs in HHT; and for cerebral abscess attributed to odontogenic bacteremias in immunocompetent HHT patients with right‐to‐left shunting through pulmonary AVMs. In conclusion, unsupervised machine learning algorithms categorise HHT haematological indices into distinct, clinically relevant profiles which are independent of age, sex or HHT genotype. Further evaluation may inform prophylaxis and management for HHT patients’ haemorrhagic and thrombotic phenotypes.