Menstrual abnormalities after COVID-19 vaccination in the Netherlands: A description of spontaneous and longitudinal patient-reported data.

AIMS: During the COVID-19 vaccination campaigns, the number of reports of menstrual abnormalities increased rapidly. Here, we describe the nature and potential risk factors associated with menstrual abnormalities based on spontaneously reporting data as well as data from a prospective cohort event monitoring (CEM) study as these are poorly studied. METHODS: Reports of menstrual abnormalities received by the Netherlands Pharmacovigilance Centre Lareb in the spontaneous reporting system between February 2021 and April 2022 were summarized. In addition, logistic regression analysis was performed on the reported menstrual abnormalities in the CEM study to assess the association between person characteristics, prior SARS-CoV-2 infection and use of hormonal contraceptives and the occurrence of menstrual abnormalities after vaccination. RESULTS: We analysed over 24 000 spontaneous reports of menstrual abnormalities and over 500 episodes (among 16 929 included women) of menstrual abnormalities in the CEM study. The CEM study showed an incidence of 41.4 per 1000 women aged ≤54 years. Amenorrhoea/oligomenorrhoea and heavy menstrual bleeding collectively accounted for about half of all abnormalities reported. Significant associations were observed for the age group 25-34 years (odds ratio 2.18; 95% confidence interval 1.45-3.41) and the Pfizer vaccine (odds ratio 3.04; 95% confidence interval 2.36-3.93). No association was observed for body mass index and presence of most comorbidities assessed. CONCLUSION: The cohort study showed a high incidence of menstrual disorders among women aged ≤54 years, and this observation was supported by the analysis of spontaneous reports. This suggests that a relation between COVID-19 vaccination and menstrual abnormalities is plausible and should be further investigated.

Integrative transnational analysis to dissect tuberculosis transmission events along the migratory route from Africa to Europe.

BACKGROUND: Growing international migration has increased the complexity of tuberculosis transmission patterns. Italy's decision to close its borders in 2018 made of Spain the new European porte entrée for migration from the Horn of Africa (HA). In one of the first rescues of migrants from this region at the end of 2018, tuberculosis was diagnosed in eight subjects, mainly unaccompanied minors. METHODS: Mycobacterium tuberculosis isolates from these recently arrived migrants were analysed by Mycobacterial Interspersed Repetitive-Unit/Variable-Number of Tandem Repeat (MIRU-VNTR) and subsequent whole genome sequencing (WGS) analysis. Data were compared with those from collections from other European countries receiving migrants from the HA and a strain-specific PCR was applied for a fast searching of common strains. Infections in a cellular model were performed to assess strain virulence. RESULTS: MIRU-VNTR analysis allowed identifying an epidemiological cluster involving three of the eight cases from Somalia (0 single-nucleotide polymorphisms between isolates, HA cluster). Following detailed interviews revealed that two of these cases had shared the same migratory route in most of the trip and had spent a long time at a detention camp in Libya. To confirm potential en route transmission for the three cases, we searched the same strain in collections from other European countries receiving migrants from the HA. MIRU-VNTR, WGS and a strain-specific PCR for the HA strain were applied. The same strain was identified in 12 cases from Eritrea diagnosed soon after their arrival in 2018 to the Netherlands, Belgium and Italy. Intracellular replication rate of the strain did not reveal abnormal virulence. CONCLUSIONS: Our study suggests a potential en route transmission of a pan-susceptible strain, which caused at least 15 tuberculosis cases in Somalian and Eritrean migrants diagnosed in four different European countries.

Towards standardisation: comparison of five whole genome sequencing (WGS) analysis pipelines for detection of epidemiologically linked tuberculosis cases.

BackgroundWhole genome sequencing (WGS) is a reliable tool for studying tuberculosis (TB) transmission. WGS data are usually processed by custom-built analysis pipelines with little standardisation between them.AimTo compare the impact of variability of several WGS analysis pipelines used internationally to detect epidemiologically linked TB cases.MethodsFrom the Netherlands, 535 Mycobacterium tuberculosis complex (MTBC) strains from 2016 were included. Epidemiological information obtained from municipal health services was available for all mycobacterial interspersed repeat unit-variable number of tandem repeat (MIRU-VNTR) clustered cases. WGS data was analysed using five different pipelines: one core genome multilocus sequence typing (cgMLST) approach and four single nucleotide polymorphism (SNP)-based pipelines developed in Oxford, United Kingdom; Borstel, Germany; Bilthoven, the Netherlands and Copenhagen, Denmark. WGS clusters were defined using a maximum pairwise distance of 12 SNPs/alleles.ResultsThe cgMLST approach and Oxford pipeline clustered all epidemiologically linked cases, however, in the other three SNP-based pipelines one epidemiological link was missed due to insufficient coverage. In general, the genetic distances varied between pipelines, reflecting different clustering rates: the cgMLST approach clustered 92 cases, followed by 84, 83, 83 and 82 cases in the SNP-based pipelines from Copenhagen, Oxford, Borstel and Bilthoven respectively.ConclusionConcordance in ruling out epidemiological links was high between pipelines, which is an important step in the international validation of WGS data analysis. To increase accuracy in identifying TB transmission clusters, standardisation of crucial WGS criteria and creation of a reference database of representative MTBC sequences would be advisable.

WGS more accurately predicts susceptibility of Mycobacterium tuberculosis to first-line drugs than phenotypic testing.

BACKGROUND: Drug-susceptibility testing (DST) of Mycobacterium tuberculosis complex (MTBC) isolates by the Mycobacteria Growth Indicator Tube (MGIT) approach is the most widely applied reference standard. However, the use of WGS is increasing in many developed countries to detect resistance and predict susceptibility. We investigated the reliability of WGS in predicting drug susceptibility, and analysed the discrepancies between WGS and MGIT against the first-line drugs rifampicin, isoniazid, ethambutol and pyrazinamide. METHODS: DST by MGIT and WGS was performed on MTBC isolates received in 2016/2017. Nine genes and/or their promotor regions were investigated for resistance-associated mutations: rpoB, katG, fabG1, ahpC, inhA, embA, embB, pncA and rpsA. Isolates that were discrepant in their MGIT/WGS results and a control group with concordant results were retested in the MGIT, at the critical concentration and a lower concentration, and incubated for up to 45 days after the control tube became positive in the MGIT. RESULTS: In total, 1136 isolates were included, of which 1121 were routine MTBC isolates from the Netherlands. The negative predictive value of WGS was ≥99.3% for all four first-line antibiotics. The majority of discrepancies for isoniazid and ethambutol were explained by growth at the lower concentrations, and for rifampicin by prolonged incubation in the MGIT, both indicating low-level resistance. CONCLUSIONS: Applying WGS in a country like the Netherlands, with a low TB incidence and low prevalence of resistance, can reduce the need for phenotypic DST for ∼90% of isolates and accurately detect mutations associated with low-level resistance, often missed in conventional DST.

SNP-IT Tool for Identifying Subspecies and Associated Lineages of Mycobacterium tuberculosis Complex.

The clinical phenotype of zoonotic tuberculosis and its contribution to the global burden of disease are poorly understood and probably underestimated. This shortcoming is partly because of the inability of currently available laboratory and in silico tools to accurately identify all subspecies of the Mycobacterium tuberculosis complex (MTBC). We present SNPs to Identify TB (SNP-IT), a single-nucleotide polymorphism-based tool to identify all members of MTBC, including animal clades. By applying SNP-IT to a collection of clinical genomes from a UK reference laboratory, we detected an unexpectedly high number of M. orygis isolates. M. orygis is seen at a similar rate to M. bovis, yet M. orygis cases have not been previously described in the United Kingdom. From an international perspective, it is possible that M. orygis is an underestimated zoonosis. Accurate identification will enable study of the clinical phenotype, host range, and transmission mechanisms of all subspecies of MTBC in greater detail.

Prevalence and Characterization of Heterogeneous Variable-Number Tandem-Repeat Clusters Comprising Drug-Susceptible and/or Variable Resistant Mycobacterium tuberculosis Complex Isolates in the Netherlands from 2004 to 2016.

The variable-number tandem-repeat (VNTR) typing method is used to study tuberculosis (TB) transmission. Clustering of Mycobacterium tuberculosis isolates with identical VNTR patterns is assumed to reflect recent transmission. Hence, clusters are thought to be homogeneous regarding antibiotic resistance. In practice, however, heterogeneous clusters are also identified. This study investigates the prevalence and characteristics of heterogeneous VNTR clusters and assesses whether isolates in these clusters remain clustered when subjected to whole-genome sequencing (WGS). In the period from 2004 to 2016, 9,072 isolates were included. Demographic and epidemiological linkage data were obtained from the Netherlands Tuberculosis Register. VNTR clusters were defined as homogeneous when isolates shared identical resistance profiles or as heterogeneous if both susceptible and (variable) resistant isolates were found. Multivariate logistic regression analysis was performed to identify factors associated with heterogeneous clustering. Isolates from 2016 were subjected to WGS, and a genetic distance of 12 single nucleotide polymorphisms (SNPs) was used as the cutoff for WGS clustering. In total, 4,661/9,072 (51%) isolates were clustered into 985 different VNTR clusters, of which 217 (22%) were heterogeneous. Patient characteristics associated with heterogeneous clustering were non-Dutch ethnicity (odds ratio [OR], 1.46 [95% confidence interval {CI}, 1.22 to 1.75]), asylum seeker (OR, 1.51 [95% CI, 1.24 to 1.85]), extrapulmonary TB (OR, 1.26 [95% CI, 1.09 to 1.46]), previous TB diagnosis (OR, 1.38 [95% CI, 1.04 to 1.82]), and not being a contact of a TB patient (OR, 1.35 [95% CI, 1.08 to 1.69]). With WGS, 34% of heterogeneous and 78% of homogeneous isolates from 2016 remained clustered. Heterogeneous VNTR clusters are common but seem to be explained by a substantial degree of false clustering by VNTR typing compared to WGS.

Correction: Epidemiological links between tuberculosis cases identified twice as efficiently by whole genome sequencing than conventional molecular typing: A population-based study.

[This corrects the article DOI: 10.1371/journal.pone.0195413.].

Epidemiological links between tuberculosis cases identified twice as efficiently by whole genome sequencing than conventional molecular typing: A population-based study.

BACKGROUND: Patients with Mycobacterium tuberculosis isolates sharing identical DNA fingerprint patterns can be epidemiologically linked. However, municipal health services in the Netherlands are able to confirm an epidemiological link in only around 23% of the patients with isolates clustered by the conventional variable number of tandem repeat (VNTR) genotyping. This research aims to investigate whether whole genome sequencing (WGS) is a more reliable predictor of epidemiological links between tuberculosis patients than VNTR genotyping. METHODS: VNTR genotyping and WGS were performed in parallel on all Mycobacterium tuberculosis complex isolates received at the Netherlands National Institute for Public Health and the Environment in 2016. Isolates were clustered by VNTR when they shared identical 24-loci VNTR patterns; isolates were assigned to a WGS cluster when the pair-wise genetic distance was ≤ 12 single nucleotide polymorphisms (SNPs). Cluster investigation was performed by municipal health services on all isolates clustered by VNTR in 2016. The proportion of epidemiological links identified among patients clustered by either method was calculated. RESULTS: In total, 535 isolates were genotyped, of which 25% (134/535) were clustered by VNTR and 14% (76/535) by WGS; the concordance between both typing methods was 86%. The proportion of epidemiological links among WGS clustered cases (57%) was twice as common than among VNTR clustered cases (31%). CONCLUSION: When WGS was applied, the number of clustered isolates was halved, while all epidemiologically linked cases remained clustered. WGS is therefore a more reliable tool to predict epidemiological links between tuberculosis cases than VNTR genotyping and will allow more efficient transmission tracing, as epidemiological investigations based on false clustering can be avoided.

Occurrence and Nature of Double Alleles in Variable-Number Tandem-Repeat Patterns of More than 8,000 Mycobacterium tuberculosis Complex Isolates in The Netherlands.

Since 2004, variable-number tandem-repeat (VNTR) typing of Mycobacterium tuberculosis complex isolates has been applied on a structural basis in The Netherlands to study the epidemiology of tuberculosis (TB). Although this technique is faster and technically less demanding than the previously used restriction fragment length polymorphism (RFLP) typing, reproducibility remains a concern. In the period from 2004 to 2015, 8,532 isolates were subjected to VNTR typing in The Netherlands, with 186 (2.2%) of these exhibiting double alleles at one locus. Double alleles were most common in loci 4052 and 2163b. The variables significantly associated with double alleles were urban living (odds ratio [OR], 1.503; 95% confidence interval [CI], 1.084 to 2.084; P = 0.014) and pulmonary TB (OR, 1.703; 95% CI, 1.216 to 2.386; P = 0.002). Single-colony cultures of double-allele strains were produced and revealed single-allele profiles; a maximum of five single nucleotide polymorphisms (SNPs) was observed between the single- and double-allele isolates from the same patient when whole-genome sequencing (WGS) was applied. This indicates the presence of two bacterial populations with slightly different VNTR profiles in the parental population, related to genetic drift. This observation is confirmed by the fact that secondary cases from TB source cases with double-allele isolates sometimes display only one of the two alleles present in the source case. Double alleles occur at a frequency of 2.2% in VNTR patterns in The Netherlands. They are caused by biological variation rather than by technical aberrations and can be transmitted either as single- or double-allele variants.

A Predominant Variable-Number Tandem-Repeat Cluster of Mycobacterium tuberculosis Isolates among Asylum Seekers in the Netherlands and Denmark, Deciphered by Whole-Genome Sequencing.

In many countries, Mycobacterium tuberculosis isolates are routinely subjected to variable-number tandem-repeat (VNTR) typing to investigate M. tuberculosis transmission. Unexpectedly, cross-border clusters were identified among African refugees in the Netherlands and Denmark, although transmission in those countries was unlikely. Whole-genome sequencing (WGS) was applied to analyze transmission in depth and to assess the precision of VNTR typing. WGS was applied to 40 M. tuberculosis isolates from refugees in the Netherlands and Denmark (most of whom were from the Horn of Africa) that shared the exact same VNTR profile. Cluster investigations were undertaken to identify in-country epidemiological links. Combining WGS results for the isolates (all members of the central Asian strain [CAS]/Delhi genotype), from both European countries, an average genetic distance of 80 single-nucleotide polymorphisms (SNPs) (maximum, 153 SNPs) was observed. The few pairs of isolates with confirmed epidemiological links, except for one pair, had a maximum distance of 12 SNPs. WGS divided this refugee cluster into several subclusters of patients from the same country of origin. Although the M. tuberculosis cases, mainly originating from African countries, shared the exact same VNTR profile, most were clearly distinguished by WGS. The average genetic distance in this specific VNTR cluster was 2 times greater than that in other VNTR clusters. Thus, identical VNTR profiles did not represent recent direct M. tuberculosis transmission for this group of patients. It appears that either these strains from Africa are extremely conserved genetically or there is ongoing transmission of this genotype among refugees on their long migration routes from Africa to Europe.

Whole genome sequencing as the ultimate tool to diagnose tuberculosis.

In the past two decades, DNA techniques have been increasingly used in the laboratory diagnosis of tuberculosis (TB). The (sub) species of the Mycobacterium tuberculosis complex are usually identified using reverse line blot techniques. The resistance is predicted by the detection of mutations in genes associated with resistance. Nevertheless, all cases are still subjected to cumbersome phenotypic resistance testing. The production of a strain-characteristic DNA fingerprint, to investigate the epidemiology of TB, is done by the 24-locus variable number tandem repeat (VNTR) typing. However, most of the molecular techniques in the diagnosis of TB can eventually be replaced by whole genome sequencing (WGS). Many international TB reference laboratories are currently working on the introduction of WGS; however, standardization in the international context is lacking. The European Centre for Infectious Disease Prevention and Control in Stockholm, Sweden organizes a yearly round of quality control on VNTR typing and in 2015 for the first time also WGS. In this first proficiency study, only three out of eight international TB laboratories produced WGS results in line with those of the reference laboratory. The whole process of DNA isolation, purification, quantification, sequencing, and analysis/interpretation of data is still under development. In this presentation, many aspects will be covered that influence the quality and interpretation of WGS results. The turn-around-time, analysis, and utility of WGS will be discussed. Moreover, the experiences in the use of WGS in the molecular epidemiology of TB in The Netherlands are detailed. It can be concluded that many difficulties still have to be conquered. The state of the art is that bacteria still have to be cultured to have sufficient quality and quantity of DNA for succesful WGS. The quality of sequencing has improved significantly over the past 7years, and the detection of mutations has, therefore, become more reliable. The resistance mutations detected in WGS are in line with the ones visualized in reverse line blot techniques. The turnover in the genome of M. tuberculosis is very low, ∼0.3-0.5 mutations per genome per year. However, there is a wide variation in the occurrence of mutations per strain and genotype. Still, the resolution of WGS in epidemiological typing is higher than that in VNTR typing; previously suggested epidemiological links by VNTR typing are sometimes refuted on the basis of WGS. Although WGS offers the highest resolution in typing, in a country like The Netherlands, there are many strains with a limited genetic distance up to 100 mutations, without an apparent epidemiological link between the respective cases. These lookalikes are presumably even more prevalent in settings where predominant genotypes of M. tuberculosis are circulating. In summary, WGS seems to yield a more reliable prediction of resistance by the (lack of) detection of mutations in all 25 genes ever associated with resistance. This may within a short while prevent the need for many phenotypic resistance tests. Although more robust algorithms need to be developed, the recognition of the (sub) species in the M. tuberculosis complex seems possible. The first detailed studies on the population structure of M. tuberculosis strains in The Netherlands provide more resolution in typing but also an interesting observation that a part of the strains are genetically so conserved that they are separated by less than 100 mutations. This demands a more extended and accurate validation and understanding of the utility of WGS in the epidemiology of TB.

Persistent high antibody titres against Coxiella burnetii after acute Q fever not explained by continued exposure to the source of infection: a case-control study.

BACKGROUND: From 2007 to 2010, (the southern part of) the Netherlands experienced a large Q fever epidemic, with more than 4,000 reported symptomatic cases. Approximately 1 - 5% of the acute Q fever patients develop chronic Q fever. A high IgG antibody titre against phase I of Coxiella burnetii during follow-up is considered a marker of chronic Q fever. However, there is uncertainty about the significance and cause of persistence of high IgG phase I antibody titres in patients that do not have any additional manifestations of chronic Q fever. We studied whether continued or repeated exposure to the source of infection could explain elevated IgG phase I antibody levels. METHODS: A case-control study was performed to analyze predictors for possible chronic Q fever. Possible chronic Q fever cases (n = 53) are patients with phase I IgG antibody titre ≥1:1,024 at any point in the 9 - 18 months after acute Q fever diagnosis, with a negative PCR test result for C. burnetii DNA and without other disease manifestations. Controls (n = 110) are acute Q fever patients that did not develop chronic Q fever, and who consistently had phase I IgG antibody titre <1:1,024 during the 9 - 18 months follow-up. Binary logistic regression was performed to analyze the effect of living close to an infected farm on the high antibody titres. A longitudinal analysis described the serological profiles of cases and controls. RESULTS: Proximity to infected farms and contact with animal placental material were not associated with an increased risk for possible chronic Q fever. Possible chronic Q fever patients have high IgG phase II as well as IgG phase I antibody titres, even after 48 months of follow-up. CONCLUSION: We were unable to explain the cause of persistent high IgG phase I titres among possible chronic Q fever patients by being continuously exposed to the source of infection.