Unraveling tuberculosis patient cluster transmission chains: integrating WGS-based network with clinical and epidemiological insights.

Background: Tuberculosis remains a global health threat, and the World Health Organization reports a limited reduction in disease incidence rates, including both new and relapse cases. Therefore, studies targeting tuberculosis transmission chains and recurrent episodes are crucial for developing the most effective control measures. Herein, multiple tuberculosis clusters were retrospectively investigated by integrating patients' epidemiological and clinical information with median-joining networks recreated based on whole genome sequencing (WGS) data of Mycobacterium tuberculosis isolates. Methods: Epidemiologically linked tuberculosis patient clusters were identified during the source case investigation for pediatric tuberculosis patients. Only M. tuberculosis isolate DNA samples with previously determined spoligotypes identical within clusters were subjected to WGS and further median-joining network recreation. Relevant clinical and epidemiological data were obtained from patient medical records. Results: We investigated 18 clusters comprising 100 active tuberculosis patients 29 of whom were children at the time of diagnosis; nine patients experienced recurrent episodes. M. tuberculosis isolates of studied clusters belonged to Lineages 2 (sub-lineage 2.2.1) and 4 (sub-lineages 4.3.3, 4.1.2.1, 4.8, and 4.2.1), while sub-lineage 4.3.3 (LAM) was the most abundant. Isolates of six clusters were drug-resistant. Within clusters, the maximum genetic distance between closely related isolates was only 5-11 single nucleotide variants (SNVs). Recreated median-joining networks, integrated with patients' diagnoses, specimen collection dates, sputum smear microscopy, and epidemiological investigation results indicated transmission directions within clusters and long periods of latent infection. It also facilitated the identification of potential infection sources for pediatric patients and recurrent active tuberculosis episodes refuting the reactivation possibility despite the small genetic distance of ≤5 SNVs between isolates. However, unidentified active tuberculosis cases within the cluster, the variable mycobacterial mutation rate in dormant and active states, and low M. tuberculosis genetic variability inferred precise transmission chain delineation. In some cases, heterozygous SNVs with an allelic frequency of 10-73% proved valuable in identifying direct transmission events. Conclusion: The complex approach of integrating tuberculosis cluster WGS-data-based median-joining networks with relevant epidemiological and clinical data proved valuable in delineating epidemiologically linked patient transmission chains and deciphering causes of recurrent tuberculosis episodes within clusters.

Comparative Microbiome Analysis of Three Epidemiologically Important Tick Species in Latvia.

(1) Background: Amplicon-based 16S rRNA profiling is widely used to study whole communities of prokaryotes in many niches. Here, we comparatively examined the microbial composition of three tick species, Ixodes ricinus, Ixodes persulcatus and Dermacentor reticulatus, which were field-collected in Latvia. (2) Methods: Tick DNA samples were used for microbiome analysis targeting bacterial 16S rDNA using next-generation sequencing (NGS). (3) Results: The results showed significant differences in microbial species diversity and composition by tick species and life stage. A close similarity between microbiomes of I. ricinus and I. persulcatus ticks was observed, while the D. reticulatus microbiome composition appeared to be more distinct. Significant differences in alpha and beta microbial diversity were observed between Ixodes tick life stages and sexes, with lower taxa richness indexes obtained for female ticks. The Francisella genus was closely associated with D. reticulatus ticks, while endosymbionts Candidatus Midichlorii and Candidatus Lariskella were associated with I. ricinus and I. persulcatus females, respectively. In I. ricinus females, the endosymbiont load negatively correlated with the presence of the Rickettsia genus. (4) Conclusions: The results of this study revealed important associations between ticks and their microbial community and highlighted the microbiome features of three tick species in Latvia.

Tracing microbial communities associated with archaeological human samples in Latvia, 7-11th centuries AD.

In the grave environment, microorganisms are major ecological participants in the successional decomposition of vertebrates and could infiltrate the skeleton/skeletal material during taphonomic processes. The diversity of archaeological skeleton-associated microbial assemblages and the impact of various factors are poorly understood. This study aimed to evaluate the taxonomic microbial composition of archaeological human bone and teeth samples from the 7th to 11th centuries AD from two burial sites in Latvia. Samples were analysed by a shotgun metagenomics-based approach. The results showed a strong presence of the environmental DNA in the samples, and variability in microbial community structure between individual samples. Differences in microbial composition were observed between bone and tooth samples, as well as between different burial sites. Furthermore, the presence of endogenous ancient DNA (aDNA) in tooth samples was detected. Overall, compositions of microbial communities associated with archaeological human remains in Latvia dated 7-11th century AD were influenced by the sample type and burial location. These findings indicate that, while the content of historical DNA in archaeological samples is low, the comparison of archaeological skeleton-associated microbial assemblages across time and space, along with aDNA damage profile analysis, is important and could help to reveal putative ancient microorganisms.

Oral microbiome variations related to ageing: possible implications beyond oral health.

The global population is getting older due to a combination of longer life expectancy and declining birth rates. Growing evidence suggests that the oral microbiota composition and distribution may have a profound effect on how well we age. The purpose of this study was to investigate age-related oral microbiome variations of supragingival plaque and buccal mucosa samples in the general population in Latvia. Our results indicated significant difference between supragingival plaque bacterial profiles of three age groups (20-40; 40-60; 60 + years). Within supragingival plaque samples, age group 20-40 showed the highest bacterial diversity with a decline during the 40-60 age period and uprise again after the age of 60. Among other differences, the important oral commensal Neisseria had declined after the age of 40. Additionally, prevalence of two well-documented opportunistic pathogens Streptococcus anginosus and Gemella sanguinis gradually rose with age within our samples. Furthermore, supragingival plaque and buccal mucosa samples significantly differed in overall bacterial composition.

Advantages of analysing both pairwise SNV-distance and differing SNVs between Mycobacterium tuberculosis isolates for recurrent tuberculosis cause determination.

Endogenous reactivation and exogenous reinfection are two possible causes of recurrent tuberculosis (TB). However, in some cases, precise cause determination can be challenging. In this study, we used whole genome sequencing to determine pairwise SNV distances and detect differing SNVs in initial and subsequent isolates for recurrent TB cases when the first and second episodes were caused by Mycobacterium tuberculosis (Mtb) strains with an identical spoligotype pattern. In total, 104 Mtb isolates from 36 recurrent TB and 16 single TB episode patients were included in the study. Most isolate pairs belonged to the SIT1 (n=21), SIT42 (n=9), SIT53 (n=9), and SIT254 (n=7) spoligotypes, and in 27 cases, resistance to at least one anti-TB drug was found in either isolate. Drug susceptibility was more common in the recurrent TB patient cohort, and longitudinal single TB episode isolates were more prone to be drug-resistant (p=0.03), while the association between patient cohort and spoligotype was not statistically significant (p=0.07). The pairwise SNV-distance between the longitudinal single TB episode isolates was small (0-7 SNVs). Among the recurrent TB isolates, based on the high SNV-distance (38-273 SNVs), six reinfection cases (16.7%) were identified. This distance was small (<10 SNVs) in the remaining 30 isolate pairs. Further analysis of differing SNVs revealed that 22 (61.1%) cases could be classified as possible reactivation. Notably, despite the small distance of 2-7 SNVs, initial isolates of eight patients (22.2%) had several SNVs that were not found in the second isolates; therefore, these cases were classified as reinfection with a closely related Mtb strain. No statistically significant difference in the time interval between specimen collection in the reactivation and reinfection Mtb sample groups (p=0.13) or an association between recurrence cause and drug resistance status (p=0.62) or spoligotype (p=0.79) could be detected. The mycobacterial median mutation rate of longitudinal single TB episodes and possible reactivation isolate pairs (n=37) was 0.12 SNVs/genome/year (IQR 0-0.39), and in 18 cases (48.6%), it was equal to zero. No statistically significant differences in mutation rate were found between recurrent TB and longitudinal single TB episode isolates (p=0.087), drug-susceptible and resistant isolates (p=0.37) or isolates of Beijing and other genotype families (p=0.33). Furthermore, four cases of fluoroquinolone resistance development through the acquired SNVs in the gyrA gene were identified. To conclude, this study highlighted the complexity of recurrent episode cause determination and showed the usefulness of differing SNV identification in both Mtb isolates in such cases. Expected drug susceptibility was the only discriminative factor for recurrent TB episode-causing mycobacterial strains, while no differences between reactivation and reinfection sample groups could be identified.

Application of natural sciences methodology in archaeological study of Iron Age burials in Latvia: pilot study.

Natural sciences provide several modern methodologies that could be successfully applied in archaeological studies. In this pilot study, archaeological human remains from two Iron Age cemeteries (7th-twelfth centuries AD), Lejasbiteni and Cunkani-Drengeri, which are located in different regions of Latvia, were studied. We applied ancient DNA (aDNA) and tooth enamel peptide analysis to determine the biological sex of the individuals. In addition, aDNA analysis was used to perform mtDNA haplogroup analysis. In most cases, the results of aDNA analysis regarding the biological sex of individuals coincided with the gender assigned based on grave orientation and grave goods. The results of sex determination using peptide analysis in all four individuals for whom data were available matched the possible gender. Of the 17 samples that had sufficient DNA for sequencing, seven samples had enough reads to perform mtDNA haplogroup analysis. The H2a2a, I4a1, H2a2a1, and H16c mtDNA haplogroups were identified in the individuals from the Lejasbiteni cemetery, while the T2b and K1a + 150 mtDNA haplogroups were identified in the individuals from the Cunkani-Drengeri cemetery. Overall, the obtained results demonstrated the feasibility of applying aDNA and tooth enamel peptide analysis for biological sex determination within archaeological studies. The availability of human aDNA data will be highly useful for investigating the demographic history and social structures in Iron Age Latvia.

Implementation of a next-generation sequencing-based targeted approach for full-length CYP3A4 gene sequencing.

Aim: To evaluate the application of next-generation sequencing-based targeted protocol for full-length CYP3A4 gene sequencing analysis. Materials & methods: The developed sequencing protocol was applied to analyze human DNA samples (n = 7) obtained from tuberculosis patients admitted to the Riga East University Hospital, Center of Tuberculosis and Lung diseases. Results: The sequencing data quality was sufficient for the detection of already known genetic variants, as well as for identifying rare and novel variants dispersed throughout the CYP3A4 gene with a high degree of confidence. Conclusion: Developed protocol can be applied in subpopulation level association studies to determine whether specific genetic variants or variant combinations from multiple regions of the CYP3A4 gene are of clinical significance.

The Postmedieval Latvian Oral Microbiome in the Context of Modern Dental Calculus and Modern Dental Plaque Microbial Profiles.

Recent advantages in paleomicrobiology have provided an opportunity to investigate the composition of ancient microbial ecologies. Here, using metagenome analysis, we investigated the microbial profiles of historic dental calculus retrieved from archaeological human remains from postmedieval Latvia dated 16-17th century AD and examined the associations of oral taxa and microbial diversity with specific characteristics. We evaluated the preservation of human oral microbiome patterns in historic samples and compared the microbial composition of historic dental calculus, modern human dental plaque, modern human dental calculus samples and burial soil microbiota. Overall, the results showed that the majority of microbial DNA in historic dental calculus originated from the oral microbiome with little impact of the burial environment. Good preservation of ancient DNA in historical dental calculus samples has provided reliable insight into the composition of the oral microbiome of postmedieval Latvian individuals. The relative stability of the classifiable oral microbiome composition was observed. Significant differences between the microbiome profiles of dental calculus and dental plaque samples were identified, suggesting microbial adaptation to a specific human body environment.

Large-scale countrywide screening for tick-borne pathogens in field-collected ticks in Latvia during 2017-2019.

BACKGROUND: Tick-borne diseases are of substantial concern worldwide in both humans and animals. Several hard tick species are of medical and veterinary interest in Europe, and changes in the range of tick species can affect the spread of zoonotic pathogens. The aim of the present study was to map the current prevalence and distribution pattern of ticks and related tick-borne pathogens in Latvia, a Baltic state in northern Europe. METHODS: Nearly 4600 Ixodes ricinus, I. persulcatus and Dermacentor reticulatus tick samples were collected in all regions of Latvia during 2017-2019 and were screened by molecular methods to reveal the prevalence and distribution pattern of a wide spectrum of tick-borne pathogens. RESULTS: New localities of D. reticulatus occurrence were found in western and central Latvia, including the Riga region, indicating that the northern border of D. reticulatus in Europe has moved farther to the north. Among the analyzed ticks, 33.42% carried at least one tick-borne pathogen, and 5.55% of tick samples were positive for two or three pathogens. A higher overall prevalence of tick-borne pathogens was observed in I. ricinus (34.92%) and I. persulcatus (31.65%) than in D. reticulatus (24.2%). The molecular analysis revealed the presence of tick-borne encephalitis virus, Babesia spp., Borrelia spp., Anaplasma phagocytophilum and Rickettsia spp. Overall, 15 and 7 tick-borne pathogen species were detected in Ixodes spp. and D. reticulatus ticks, respectively. This is the first report of Borrelia miyamotoi in Latvian field-collected ticks. CONCLUSIONS: This large-scale countrywide study provides a snapshot of the current distribution patterns of Ixodes and Dermacentor ticks in Latvia and gives us a reliable overview of tick-borne pathogens in Latvian field-collected ticks.

Application of whole-genome sequencing in a case study of renal tuberculosis in a child.

BACKGROUND: Urogenital tuberculosis (TB) is rare in children and usually develops due to reactivation of the foci in the genitourinary tract after the latency period following initial infection. Urogenital TB in children has no pathognomonic clinical features that can result in overlooking or misdiagnosing this clinical entity. Here, we report important findings regarding the pathogenesis and transmission of TB by using genotyping and whole-genome sequencing (WGS) in a study of renal TB case in a child. CASE PRESENTATION: A 13-year-old boy was admitted to the hospital because of high fever, severe dry cough, flank pain and painful urination. Abdominal ultrasonography and CT revealed an 8 mm calculus in the kidney, and clinical findings were initially interpreted as nephrolithiasis. Nevertheless, due to the atypical clinical presentation of kidney stone disease, additional investigations for possible TB were performed. The QuantiFERON®-TB Gold Plus test was positive, and the Mantoux test resulted in 15 mm of induration, confirming infection with Mycobacterium tuberculosis (Mtb). Chest X-ray was normal. Chest CT revealed calcified intrathoracic lymph nodes. The urine sample tested positive for acid-fast bacilli, and Mtb cultures were obtained from urine and bronchial aspirate samples, resulting in a final diagnosis of intrathoracic lymph node and renal TB. Contact investigation revealed that the child's father was diagnosed with TB when the child was 1 year old. Genotyping and WGS analysis of Mtb isolates of the child and his father confirmed the epidemiological link and pointed to the latency of infection in the child. CONCLUSIONS: This case report confirmed the development of active TB from calcified lesions in adolescent after 12 years of exposure, demonstrated the absence of microevolutionary changes in the Mtb genome during the period of latency, and proved the importance of appropriate evaluation and management to prevent the progression of TB infection to active TB disease. The use of WGS provided the ultimate resolution for the detection of TB transmission and reactivation events.

Comparison of telomere length between population-specific mitochondrial haplogroups among different age groups in a Latvian population.

Population studies have demonstrated that telomere length (TL) displays great diversity among different populations. Previously described controversial findings associated longevity with specific mitochondrial DNA haplogroups (hgs) (e.g., J and U). These observations may be influenced by population diversity, geographic location, and/or specific historic background. The aims of this study were to identify a specific hg which correlates with aging in a Latvian populating and to evaluate the possible association of TL variability with specific mitochondrial hgs. The results show no significant correlation between TL, mitochondrial DNA hgs and longevity. A slight increase in frequency was observed among centenarians of hg H; however, these findings were not statistically significant. TL did not show any statically significant difference, only hg W had slightly longer telomeres among others. An insignificant increase in TL was observed in the 55-89 age group of hg W but in the <90 age group for hg J which also had the longest TL in the 20-45 age group. In conclusion this study indicates that specific mitochondrial DNA hgs do not have a significant, if any, influence on the variation of TL in Latvians.