Pol Gene Research Articles

A preliminary study of gene expression changes in Koalas Infected with Koala Retrovirus (KoRV) and identification of potential biomarkers for KoRV pathogenesis

BackgroundKoala retrovirus (KoRV), a major pathogen of koalas, exists in both endogenous (KoRV-A) and exogenous forms (KoRV-A to I and K to M) and causes multiple disease phenotypes, including carcinomas and immunosuppression. However, the direct association between the different KoRV subtypes and carcinogenesis remains unknown. Differentially expressed gene (DEG) analysis of peripheral blood mononuclear cells (PBMCs) of koalas carrying both endogenous (KoRV-A) and exogenous (KoRV-A, B, and C) subtypes was performed using a high-throughput RNA-seq approach. PBMCs were obtained from three healthy koalas: one infected with endogenous (KoRV-A; Group I) and two infected with exogenous (KoRV-B and/or KoRV-C; Group II) subtypes. Additionally, spleen samples (n = 6) from six KoRV-infected deceased koalas (K1- K6) and blood samples (n = 1) from a live koala (K7) were collected and examined to validate the findings.ResultsAll koalas were positive for the endogenous KoRV-A subtype, and eight koalas were positive for KoRV-B and/or KoRV-C. Transcription of KoRV gag, pol, and env genes was detected in all koalas. Upregulation of cytokine and immunosuppressive genes was observed in koalas infected with KoRV-B or KoRV-B and -C subtypes, compared to koalas infected with only KoRV-A. We found 550 DEG signatures with significant (absolute p < 0.05, and absolute log2 Fold Change (FC) > 1.5) dysregulation, out of which 77.6% and 22.4% DEGs were upregulated (log2FC > 1.5) and downregulated (log2FC < − 1.5), and downregulated (log2 FC < − 1), respectively. We identified 17 unique hub genes (82.3% upregulated and 17.7% down-regulated), with KIF23, CCNB2, POLR3F, and RSL24D1 detected as the potential hub genes modified with KoRV infection. Real-time RT-qPCR was performed on seven koalas to ascertain the expression levels of four potential hub genes, which were subsequently normalized to actin copies. Notably, all seven koalas exhibited distinct expression signatures for the hub genes, especially, KIF23 and CCNB2 show the highest expression in healthy koala PBMC, and POLR3F shows the highest expression in koala with lymphoma (K1).ConclusionThus, it can be concluded that multiple KoRV subtypes affect disease progression in koalas and that the predicted hub genes could be promising prognostic biomarkers for pathogenesis.

Prevalence of major INSTI and HIV-1 drug resistance mutations in pre- and antiretroviral-treated patients in Indonesia

Indonesia has one of the highest HIV infection rates in Southeast Asia. The use of dolutegravir, an integrase strand transfer inhibitor (INSTI), as a first-line treatment underscores the need for detailed data on INSTI drug resistance mutations (DRMs). Currently, there is a lack of comprehensive data on DRMs INSTI and other HIV drug resistance in Indonesian patients, both pre- and post-treatment. The aim of this study was to identify the subtypes and drug resistance mutations of the protease, reverse transcriptase, and integrase genes in both treatment-naive and ARV-treated patients in Bandung, West Java, Indonesia. A cross-sectional study was conducted involving HIV-positive patients at Hasan Sadikin Hospital, Bandung, Indonesia, from September 2022 to January 2023. The patients were categorized into two groups: ARV-treated and pre-treatment patients. Peripheral blood mononuclear cells (PBMCs) were processed for DNA extraction, followed by amplification and sequencing of the pol gene to detect mutations and subtypes. The study found that the predominant subtype was CRF01_AE, accounting for 85.4% and 69% of pre-treatment and treated patients, respectively, followed by recombinant forms such as A1/CRF01_AE, CRF01_AE/CRF02_AG, subtype B, and other subtypes. Among ARV-treated/INSTI-naive patients, major INSTI DRMs R263K and Y143H were identified, while pre-treatment patients exhibited accessory integrase DRMs. The most common DRMs detected were non-nucleoside reverse transcriptase inhibitor (NNRTI) DRMs, with prevalences of 14.6% and 7% in pre-treatment and ARV-treated patients, respectively. In conclusion, CRF01_AE emerged as the predominant subtype in both pre-treatment and ARV-treated patients in Bandung, underscoring the necessity for ongoing surveillance of integrase DRMs, particularly given the presence of major INSTI DRMs in patients undergoing INSTI treatment.

Molecular detection and phylogenetic analysis of ovine herpesvirus-2 in subclinical infections of cattle and sheep.

Ovine herpesvirus-2 (OvHV-2) is the causative agent of malignant catarrhal fever (MCF), a serious and often fatal disease that affects cattle and other ruminants. This study aimed to investigate the molecular epidemiology and genetic diversity of OvHV-2 strains circulating in sheep and cattle populations in the Jammu and Kashmir region of India. Screening of 150 sheep and 57 cattle blood samples revealed the presence of the OvHV-2 polymerase (pol) gene in 8.6% of sheep, 10% of apparently healthy cattle, and 29.7% of cattle exhibiting MCF-like symptoms. The full-length glycoprotein B (gB) gene (2800bp) and an 875bp internal fragment were successfully amplified, cloned, and sequenced from pol-positive samples. Comparative sequence analysis of the deduced gB amino acid sequences identified seven substitutions at positions 278, 341, 390, 440, 468, 539, and 566 compared to reference strains. Phylogenetic analysis based on the gB nucleotide sequences clustered the OvHV-2 strains from this study within the Indian clade, distinct from strains reported in the UK and US. These findings provide insights into the genetic diversity of OvHV-2 strains circulating in Jammu and Kashmir, with the identified mutations potentially influencing virus-host interactions. Further investigations into the functional implications of these mutations are warranted to understand their role in viral pathogenesis and tropism.

Enhanced Transcription of Human Endogenous Retroviruses and TRIM28 Downregulation in Patients with Inflammatory Bowel Disease.

Inflammatory bowel disease (IBD) includes patients affected by Crohn's disease or ulcerative colitis. IBD is thought to be a chronic immune-mediated disease, but its origin remains elusive, and this limits new therapeutic approaches. Human endogenous retroviruses (HERVs) originate from ancestral infections and represent 8% of the human genome. HERVs are no longer infectious, but some retroviral sequences can be activated, and their aberrant expressions have been implicated in inflammatory and autoimmune disorders. HERV transcription is regulated by TRIM28 and SETDB1, which are also directly involved in epigenetic processes and modulation of the immune response. Using a PCR real-time Taqman amplification assay, we assessed, for the first time, the transcription levels of pol genes of HERV-H, -K, and -W families of env genes of syncytin 1 (SYN1), SYN2, and HERV-W, as well as of TRIM28 and SETDB1 in the whole blood of 48 patients with Crohn's disease (CD), 20 with ulcerative colitis (UC), and in healthy controls (HC) of comparable age. The transcriptional levels of HERV-H-pol (p = 0.0003) and HERV-K-pol (p = 0.001) were significantly higher in IBD patients compared with HC, with no differences between patients with CD and UC. No significant differences were found for the remaining HERVs between IBD patients and HC. The transcript levels of TRIM28 were significantly downregulated in IBD patients (p < 0.001), without differences between CD and UC, while the SETDB1 levels were preserved. The enhanced transcription of HERV-H-pol and HERV-K-pol, as well as the impaired activation of TRIM28, were not influenced by clinical disease activity and type of treatment. The overexpression of HERVs and impaired transcription of TRIM28 in patients affected by CD or UC suggest that they might be the main actors in the pathophysiology of IBD, opening the way to innovative targeted interventions.

Comparison of short-read and long-read next-generation sequencing technologies for determining HIV-1 drug resistance.

Accurate HIV-1 genome sequencing is necessary to identify drug resistance mutations (DRMs) in people with HIV-1 (PWH). Next-generation-sequencing (NGS) allows the detection of minor variants and is now available in many laboratories. Our study aimed to compare two NGS approaches, a "short read" sequencing protocol using DeepChek® Whole Genome HIV-1 Assay on Illumina, and a "long read" sequencing protocol of HIV-1 pol and env single-molecule real-time sequencing (SMRT) on Pacific Biosciences (PacBio). We analyzed 16 plasma samples and 13 cellular samples from PWH. HIV-1 whole genome was amplified into five amplicons using DeepChek® Whole Genome HIV-1 Assay and sequenced on an iSeq. 100. In parallel, HIV-1 pol and env genes were separately amplified and sequenced using PacBio SMRT system with the circular consensus sequencing mode on a Sequel IIe. Concordance rates for determining DRMs with both approaches varied depending on the HIV-1 region, with higher concordance in the integrase region compared to the reverse transcriptase and protease regions. DeepChek® Whole Genome HIV-1 Assay exhibited better sensitivity in HIV-1 RNA sequencing of plasmas with lower viral loads. In cell HIV-1 DNA sequencing, the DeepChek® Whole Genome HIV-1 Assay performed better in pol and env sequencing but detected more APOBEC-induced DRMs, which can represent defective proviruses. Our findings indicate that both DeepChek® Whole Genome HIV-1 Assay and PacBio SMRT sequencing exhibit good performance for subtype determination, detection, and quantification of DRMs of the HIV-1 genome. However, some discrepancies were found in cellular samples, highlighting the challenges of interpreting HIV-1 DNA DRMs.

Characterizing HIV-1 transmission by genetic cluster analysis among newly diagnosed patients in the China-Myanmar border region from 2020 to 2023

ABSTRACT Cluster analysis of HIV sequence can provide insights into viral transmission patterns in border regions. This study aims to illuminate the HIV-1 subtype distribution and transmission dynamics among newly diagnosed individuals in Dehong prefecture, a region along the China-Myanmar border. Among 948 participants with pol gene sequences, 36 HIV-1 subtypes were identified, with URFs (18.8%, 178/948) being the dominant strain, followed by CRF01_AE (18.5%, 175/948) and CRF07_BC (10.9%, 103/948). Additionally, 287 sequences (30.3%, 287/948) were grouped into 91 clusters, 31 of which contained both Chinese and Burmese individuals. Multivariable logistic regression indicated that men who have sex with men (MSM), CD4 + cell count of 200∼499 cells/μl, and 500 cells/μl and above, and CRF01_AE were risk factors for entering the network. Through Chord diagram, we found frequent transmission relationships among heterosexual China male group, especially those over 35 years of age. Additionally, the correlation between heterosexual Myanmar female group and heterosexual China male group among cross-risk groups deserved to be emphasized. Furthermore, the network exhibited a growing trend over time, with the largest active transmission cluster identified in Ruili county. In conclusion, the HIV-1 subtype landscape in Dehong has become increasingly complex, and the region has faced risks of transmission from both domestic and international sources. Targeted intervention strategies should be implemented for MSM, heterosexual Chinese middle-aged and elderly men, and heterosexual Burmese young adults to mitigate these risks. These findings provided evidence-based insights for local government to formulate coordinated transnational intervention approaches.

Combining molecular transmission network analysis and spatial epidemiology to reveal HIV-1 transmission pattern among the older people in Nanjing, China

BackgroundIn China, the problem of HIV infection among the older people has become increasingly prominent. This study aimed to analyze the pattern and influencing factors of HIV transmission based on a genomic and spatial epidemiological analysis among this population.MethodsA total of 432 older people who were aged ≥ 50 years, newly diagnosed with HIV-1 between January 2018 and December 2021 and without a history of ART were enrolled. HIV-1 pol gene sequence was obtained by viral RNA extraction and nested PCR. The molecular transmission network was constructed using HIV-TRACE and the spatial distribution analyses were performed in ArcGIS. The multivariate logistic regression analysis was performed to analyze the factors associated with clustering.ResultsA total of 382 sequences were successfully sequenced, of which CRF07_BC (52.3%), CRF01_AE (32.5%), and CRF08_BC (6.8%) were the main HIV-1 strains. A total of 176 sequences entered the molecular network, with a clustering rate of 46.1%. Impressively, the clustering rate among older people infected through commercial heterosexual contact was as high as 61.7% and three female sex workers (FSWs) were observed in the network. The individuals who were aged ≥ 60 years and transmitted the virus by commercial heterosexual contact had a higher clustering rate, while those who were retirees or engaged other occupations and with higher education degree were less likely to cluster. There was a positive spatial correlation of clustering rate (Global Moran I = 0.206, P < 0.001) at the town level and the highly aggregated regions were mainly distributed in rural area. We determined three large clusters which mainly spread in the intra-region of certain towns in rural areas. Notably, 54.5% of cases in large clusters were transmitted through commercial heterosexual contact.ConclusionsOur joint analysis of molecular and spatial epidemiology effectively revealed the spatial aggregation of HIV transmission and highlighted that towns of high aggregation were mainly located in rural area. Also, we found vital role of commercial heterosexual contact in HIV transmission among older people. Therefore, health resources should be directed towards highly aggregated rural areas and prevention strategy should take critical persons as entry points.

Molecular cluster, transmission characteristics, origin and dynamics analysis of HIV-1 CRF59_01B in China: A molecular epidemiology study

PurposeThis study investigated for the HIV-1 CRF59_01B epidemic's spatiotemporal dynamics and its transmission networks in China. MethodsBetween 2007 and 2020, a total of 250 partial pol gene sequences of HIV-1 CRF59_01B were collected from four regions (10 Chinese provinces). Phylogenetic tree construction and cluster identification were then performed. The Bayesian skyline and birth–death susceptible–infected–removed models were employed for the phylodynamic analyses of subtypes and large clusters, respectively. Phylogenetic analyses and trait diffusion of these sequences were performed using Bayesian phylogenetic methods (beast-classic package). Distance-based molecular network analyses were performed to identify putative relationships. ResultsUsing a genetic distance threshold of 1.3 %, We identified 45 clusters that included 62.40 % (156/250) of the sequences. Three clusters (6.67 %, 3/45) had 10 or more sequences, and were considered "large clusters". Six clusters (13.33 %) included sequences from different regions (Southeast, Northeast, Southeast, and Central China). Thirteen clusters (28.89 %) included sequences of men who had sex with men only, three clusters (6.67 %) included sequences of heterosexuals only, and 12 clusters (26.67 %) included sequences of both groups. The substitution rate of CRF59_01B was 1.91 × 10–3 substitutions per site per year [95 % highest posterior density (HPD) interval: 1.39 × 10−3–2.49 × 10−3)], the time to the most recent common ancestor of CRF59_01B was to be 1992.83 (95 % HPD: 1977.97–2002.81). A Bayesian skyline plot revealed that the effective population size of CRF59_01B increased from 2000 to 2015 and remained stable after 2015. The large clusters showed continuous growth from 2013 to 2020. Phylogeographic analysis showed that CRF59_01B B most likely originated in Southeast China, with a posterior probability of 97.44 %, and then spread to other regions. ConclusionsOur study revealed the temporal and geographical origins of HIV-1 CRF59_01B as well as the process of transmission among various regions and risk groups in China, which can help develop targeted HIV prevention strategies.

Molecular genetic monitoring of HIV-1 variants circulating in St. Petersburg

The aim of the study: to assess the genetic diversity of HIV-1 variants circulating in St. Petersburg.Materials and methods. The study included 289 patients with virological ART failure in 2022 in the St. Petersburg AIDS Center. Fragments of the pol gene encoding integrase, reverse transcriptase and protease were analyzed by polymerase chain reaction and Sanger sequencing. Phylogenetic tree created by the Neighbor-joining method with 1000 repeats of nucleotide sequences, bootstrap values &gt;70. To assess the circulation of HIV genovariants in dynamics in St. Petersburg, sequences and clinical and laboratory parameters obtained from 544 patients since 2018. The total sample included 833 samples (289 were collected in 2022), compared with the results of a study from 1104 HIV-infected patients in 2006–2011.Results and discussion. Monovariants of HIV dominated in the examined patients (95.1%, 275 people), subtype A6 — 88.2% (255 people); subtype B — 5.9% (17 people); C — 0.3% (1 person); G — 0.7% (2 people), the proportion of recombinant forms — 4.9% (14 people). In the sample collection, the proportion of new HIV cases of non-A subtype was 13.3% in 2006–2011, and 11.1% in 2018–2022. A significant increase in the incidence of recombinant forms of HIV-1 was revealed from 1.6% to 3.5% (c2=6.111; p=0.014). In the group (2018–2022), the recombinant form of CRF63_02A6 was more common (15/29 people).Conclusion. Phylogenetic analyses makes it possible to determine HIV subtypes, but also to establish the potential geographical origin of the virus, to identify transmission clusters taking into account the socio-demographic indicators of HIV-infected patients. Molecular epidemiological monitoring can be used to develop and implement programs to counter the spread of HIV among the population. The dominant genetic variant of HIV circulating in St. Petersburg is sub-subtype A6, as it was 10 years ago. The proportion of new cases of infection with non-A subtype of HIV remains stable, with a downward trend. The increase in the frequency of detection of recombinant forms of HIV-1 is probably related to the migration processes of the population.

Molecular transmission network analysis reveals the challenge of HIV-1 in ageing patients in China: elderly people play a crucial role in the transmission of subtypes and high pretreatment drug resistance in developed Eastern China, 2019–2023

BackgroundThe number and proportion of HIV/AIDS patients among older people are continuously and rapidly increasing in China. We conducted a detailed molecular epidemiological analysis of HIV-1 epidemic strains in a developed city in eastern China and found that elderly people play a crucial role in the transmission of subtypes and high pretreatment drug resistance (PDR).MethodsA total of 1048 samples were obtained from 1129 (92.8%) newly confirmed HIV-1-positive and treatment-naive patients between 2019 and 2023. The 1316 bp target fragment of the pol gene was amplified by reverse transcription polymerase chain reaction (RT‒PCR) and nested PCR, and Maximum-likelihood (ML) phylogenetic trees and molecular transmission network were constructed to analyse the subtypes and transmission clusters. Molecular transmission network was visualized using Cytoscape with the distance threshold of 0.0075. PDR-associated mutations were determined according to the Stanford University HIV Drug Resistance Database.ResultsA total of 933 pol sequences (89.0%, 933/1048) were successfully obtained, and twelve HIV-1 subtypes were detected. CRF07_BC was the predominant subtype, accounting for 48.1% (449/933) of sequences, followed by CRF01_AE (29.4%, 274/933). A total of 398 individuals (42.7%, 398/933) formed 89 clusters in the network. Multivariable logistic regression analysis revealed that age, nationality, subtype, and PDR were the most significant factors associated with clustering in the transmission network. The prevalence of PDR was 14.6% (136/933).PDR associated with non-nucleoside reverse transcriptase inhibitors (10.0%, 93/933) was much more common than that associated with nucleoside reverse transcriptase inhibitors (1.8%, 17/933) and protease inhibitors (3.2%, 30/933) (χ2 = 77.961, p < 0.001). The most frequent NNRTI mutations were K103N/S/KN/NS (52.2%, 71/136), which led to the highest proportion of high-level resistance to nevirapine and efavirenz (52.2%).ConclusionsOur study revealed the important influence of elderly people on CRF07_BC transmission and the high prevalence of PDR. The clustering of drug-resistant cases was significant, which suggested the potential for localized widespread transmission of drug-resistant strains. HIV screening and the determination of PDR are recommended for older patients to improve early detection and reduce treatment failure and second-generation transmission.

Transmitted Antiretroviral Drug Resistance to Integrase Strand Transfer Inhibitors Class in São Paulo Metropolitan Area, Brazil.

A newer integrase strand transfer inhibitor (INSTI) cabotegravir was recently approved for both therapy and prophylaxis and can play an essential role in the fight against AIDS. It shares similar resistance profile to dolutegravir, the cornerstone of Brazilian antiretroviral (ARV) treatment, with about 600 thousand people living with HIV in Brazil currently on regimens that contain this INSTI. Health services in the São Paulo metropolitan area are responsible for a large proportion of ARV dispensation in the country. Estimating transmitted drug resistance mutation (TDRM) in the area before cabotegravir introduction may provide a useful baseline information. Partial HIV-1 pol gene was sequenced (Sanger) from 192 newly diagnosed individuals from São Paulo and nearby cities (2020 to March 2023) at integrase, with 85 also at protease/reverse transcriptase regions. Retrotranscribed plasma RNA, amplified with nested PCR, was edited (Recall or Sequencher) and analyzed at Rega and Stanford db. Surveillance drug resistance mutations (SDRM) to INSTI class was detected in three cases (1.6%; 95% CI: 0.5%-5%), two E138K and one R263K, with 7.8% (95% CI: 5%-13%) with resistance mutations (major or accessory). SDRM for nucleoside reverse transcriptase inhibitors, non-nucleoside reverse transcriptase inhibitors, and PI classes were identified in 7 (8.2% CI: 95% 4%-16%) cases. Subtype B predominated (69%), followed by subtype C (16%), now the second most prevalent infection in this area. Among 131 patients treated for over 6 months, 92% were virally suppressed below 200 copies/mL, with low TCD4 counts independently associated to failure. SDRM to INSTI class is rare in the area. Intermediate rates of transmitted resistance to other ARV classes are comparable to previous estimates. Viral suppression rates may depend on TCD4 counts, another negative impact of late diagnosis in care that deserves more attention.

Analysis of virus gene subtypes and drug resistance monitoring results of newly reported HIV/AIDS population in Anhui Province from 2020 to 2023

Objective: To investigate the genetic subtypes and drug resistance monitoring of newly reported human immunodeficiency virus (HIV) infection/AIDS virus in Anhui Province from 2020 to 2023. Methods: An observational design study was used to collect blood samples from patients diagnosed with HIV/AIDS in the AIDS Prevention and Control Department of Anhui Provincial Center for Disease Control and Prevention from January 2020 to December 2023.The HIV-1 pol gene was amplified by reverse transcription-nested PCR, and the genetic subtypes were identified by phylogenetic tree analysis using MEGA 7.0 software. The mutation sites of drug resistance were analyzed by the online software tool of Stanford University's HIV Drug resistance database. The influencing factors of drug resistance before treatment were analyzed by multivariate logistic analysis. Results: A total of 335 plasma samples were collected, and 332 HIV-1 pol gene sequences were obtained successfully. The main gene subtypes were CRF01-AE, accounting for 35.55% (118/332), followed by CRF07-BC, B and B+C types [29.22% (97/332), 11.74% (39/332), 9.93% (33/332)]. The total drug resistance rate before treatment was 30.12%(32/100), and the drug resistance rate of protease inhibitor (PIs) in HIV-1 was 6.33% (21/332). The drug resistance rate of nucleoside reverse transcriptase inhibitors (NRTI) before treatment was 6.33% (21/332). The drug resistance rate of non-nucleoside reverse transcriptase inhibitors (NNRTI) before treatment was 17.47% (58/332).The comparison of drug resistance rate of different drug types showed statistical significance (χ2=30.435, P<0.05).Among the 100 cases of drug resistance, the main mutation point of HIV-1 protease inhibitor was Q58E (21.00%), and the main mutation point of nucleoside reverse transcriptase inhibitor was M184V/I (6.00%). Non-nucleoside reverse transcriptase inhibitor resistance mutation points mainly K103N (22.00%).There were statistically significant differences in the starting time of antiviral therapy, the number of CD4+T cells at baseline and the drug resistance rate of gene subtypes (the chi-square values are respectively 24.152, 32.516, 11.652, P<0.05).Multivariate logistic analysis showed that the baseline CD4+T cell count was <200/μl, subtype B, subtype B+C, CRF01-AE subtype, CRF55-01B subtype and 01-BC subtype was the influential factor of drug resistance before treatment (the chi-square values are respectively 4.577, 8.202, 4.416, 5.206, 7.603 and 4.804, P<0.05). Conclusion: The newly reported HIV/AIDS population in Anhui Province from 2020 to 2023 has a variety of viral gene subtypes, and NNRTIs are the main types of drug resistance gene mutations before treatment. Attention should be paid to the number of baseline CD4+T cells, the duration of antiviral treatment, and the distribution of gene subtypes to reduce the drug resistance of HIV/AIDS patients before treatment.

Exploring Codon Usage Patterns and Influencing Factors in Ranavirus DNA Polymerase Genes.

Ranaviruses, members of the genus Ranavirus within the family Iridoviridae, have become a significant concern for amphibian populations globally, along with other cold-blooded vertebrates, due to their emergence as a significant threat. We employed bioinformatics tools to examine the codon usage patterns in 61 DNA pol genes from Ranavirus, Lymphocystivirus, Megalocytivirus, and two unclassified ranaviruses, as no prior studies had been conducted on this topic. The results showed a slight or low level of codon usage bias (CUB) in the DNA pol genes of Ranavirus. Relative synonymous codon usage (RSCU) analysis indicated that the predominant codons favored in Ranavirus DNA pol genes terminate with C or G. Correlation analysis examining nucleotide content, third codon position, effective number of codons (ENC), correspondence analysis (COA), Aroma values, and GRAVY values indicated that the CUB across DNA pol genes could be influenced by both mutation pressure and natural selection. The neutrality plot indicated that natural selection is the primary factor driving codon usage. Furthermore, the analysis of the codon adaptation index (CAI) illustrated the robust adaptability of Ranavirus DNA pol genes to their hosts. Analysis of the relative codon deoptimization index (RCDI) suggested that Ranavirus DNA pol genes underwent greater selection pressure from their hosts. These findings will aid in comprehending the factors influencing the evolution and adaptation of Ranavirus to its hosts.

Complete genome sequence analysis of Reticuloendotheliosis virus integrated in nonhomologous Avipoxvirus

Integration of nucleic acid sequences of Reticuloendotheliosis virus (REV) in Avipoxvirus（APV） has become commonplace. In this study, 4 strains of suspected Fowlpox virus (FPV) and 1 strain of suspected Pigeonpox virus (PPV) collected in Taiyuan, Shanxi Province were cultured in chicken embryos, and the 4b core protein gene was amplified by PCR, and the identity and genome similarity were determined by sequence analysis. The sequences between the end of ORF201 and the beginning of ORF203 of FPV and PPV were then amplified, sequenced, and subjected to sequence comparison to determine genome similarity. The results showed that the isolates were 4 strains of FPV and 1 strain of PPV. The 4 isolated strains of FPV belong to type A1 virus, with 100 % identity to each other and to the FWPV-09-Jilin strain isolated in Jilin, China, and the lowest identity to the type B2 virus TNPV5/NZL/2009, which is only 74 %. PPV belongs to type A2 virus, and its identity with local strain of fowlpox virus was 90.1 %, with the highest identity of 100 % with PPLH and ROPI/W370/ON/2012 and ow_2017_3 strains, which also belong to type A2 pigeonpox virus, and the lowest identity of 73.7 % with TNPV5/NZL/2009, a type B2 virus. The complete genome of REV sequences integrated into FPV and PPV were amplified, and 5 REV nucleic acid sequences were obtained after sequencing and concatenation, with lengths ranging from 7942 to 8005 bp. The identity analysis results indicate that it has high identity with isolates from Northeast China, Guangdong, and Guangxi regions in China. Based on its gp90 protein gene, the REV integrated into the poxvirus belong to type III, with the highest identity of 99.9% with strains such as APC-566 and CY1111, and the lowest identity with REV-Anhui1, at 95.4 %. The length of the pol gene varies among different strains of REV, and its encoded amino acid changes significantly after position 675, with deletions and alterations. This study indicates that all fowlpox viruses isolated in Taiyuan, Shanxi Province have integrated the entire REV gene sequence, with high identity between them. At the same time, it indicates that the pigeonpox virus isolate has also integrated the entire REV gene sequence, and has the highest identity with the integrated REV gene sequence in fowlpox virus.

Characteristics of Subtype and Molecular Transmission Networks among Newly Diagnosed HIV-1 Infections in Patients Residing in Taiyuan City, Shanxi Province, China, from 2021 to 2023.

The HIV-1 pandemic, spanning four decades, presents a significant challenge to global public health. This study aimed to understand the molecular transmission characteristics of newly reported HIV infections in Taiyuan, Shanxi Province, China, to analyze the characteristics of subtypes and the risk factors of the transmission network, providing a scientific basis for precise prevention and intervention measures. A total of 720 samples were collected from newly diagnosed HIV-1 patients residing in Taiyuan between 2021 and 2023. Sequencing of partial genes of the HIV-1 pol gene resulted in multiple sequence acquisitions and was conducted to analyze their subtypes and molecular transmission networks. Out of the samples, 584 pol sequences were obtained, revealing 17 HIV-1 subtypes, with CRF07_BC (48.29%), CRF01_AE (31.34%), and CRF79_0107 (7.19%) being the dominant subtypes. Using a genetic distance threshold of 1.5%, 49 molecular transmission clusters were generated from the 313 pol gene sequences. Univariate analysis showed significant differences in the HIV transmission molecular network in terms of HIV subtype and household registration (p < 0.05). Multivariate logistic regression analysis showed that CRF79_0107 subtype and its migrants were associated with higher proportions of sequences in the HIV transmission network. These findings provide a scientific foundation for the development of localized HIV-specific intervention strategies.

High level of genomic divergence in orf-I p12 and hbz genes of HTLV-1 subtype-C in Central Australia

BackgroundHuman T cell lymphotropic virus type 1 (HTLV-1) infection remains a largely neglected public health problem, particularly in resource-poor areas with high burden of communicable and non-communicable diseases, such as some remote populations in Central Australia where an estimated 37% of adults are infected with HTLV-1. Most of our understanding of HTLV-1 infection comes from studies of the globally spread subtype-A (HTLV-1a), with few molecular studies reported with the Austral-Melanesian subtype-C (HTLV-1c) predominant in the Indo-Pacific and Oceania regions.ResultsUsing a primer walking strategy and direct sequencing, we constructed HTLV-1c genomic consensus sequences from 22 First Nations participants living with HTLV-1c in Central Australia. Phylogenetic and pairwise analysis of this subtype-C proviral gDNA showed higher levels of genomic divergence in comparison to previously published HTLV-1a genomes. While the overall genomic homology between subtypes was 92.5%, the lowest nucleotide and amino acid sequence identity occurred near the 3′ end of the proviral genome coding regulatory genes, especially overlapping hbz (85.37%, 77.46%, respectively) and orf-I product p12 (82.00%, 70.30%, respectively). Strikingly, the HTLV-1c genomic consensus sequences uniformly showed a defective translation start codon for the immune regulatory proteins p12/p8 encoded by the HTLV-1A orf-I. Deletions in the proviral genome were detected in many subjects, particularly in the structural gag, pol and env genes. Similarly, using a droplet digital PCR assay measuring the copies of gag and tax per reference host genome, we quantitatively confirmed that provirus retains the tax gene region at higher levels than gag.ConclusionsOur genomic analysis of HTLV-1c in Central Australia in conjunction with earlier Melanesian HTLV-1c sequences, elucidate substantial differences with respect to the globally spread HTLV-1a. Future studies should address the impact these genomic differences have on infection and the regionally distinctive frequency of associated pulmonary disease. Understanding the host and virus subtype factors which contribute to the differential morbidity observed, is crucial for the development of much needed therapeutics and vaccine strategies against this highly endemic infection in remote First Nations communities in Central Australia.

AUTO-TUNE: selecting the distance threshold for inferring HIV transmission clusters.

Molecular surveillance of viral pathogens and inference of transmission networks from genomic data play an increasingly important role in public health efforts, especially for HIV-1. For many methods, the genetic distance threshold used to connect sequences in the transmission network is a key parameter informing the properties of inferred networks. Using a distance threshold that is too high can result in a network with many spurious links, making it difficult to interpret. Conversely, a distance threshold that is too low can result in a network with too few links, which may not capture key insights into clusters of public health concern. Published research using the HIV-TRACE software package frequently uses the default threshold of 0.015 substitutions/site for HIV pol gene sequences, but in many cases, investigators heuristically select other threshold parameters to better capture the underlying dynamics of the epidemic they are studying. Here, we present a general heuristic scoring approach for tuning a distance threshold adaptively, which seeks to prevent the formation of giant clusters. We prioritize the ratio of the sizes of the largest and the second largest cluster, maximizing the number of clusters present in the network. We apply our scoring heuristic to outbreaks with different characteristics, such as regional or temporal variability, and demonstrate the utility of using the scoring mechanism's suggested distance threshold to identify clusters exhibiting risk factors that would have otherwise been more difficult to identify. For example, while we found that a 0.015 substitutions/site distance threshold is typical for US-like epidemics, recent outbreaks like the CRF07_BC subtype among men who have sex with men (MSM) in China have been found to have a lower optimal threshold of 0.005 to better capture the transition from injected drug use (IDU) to MSM as the primary risk factor. Alternatively, in communities surrounding Lake Victoria in Uganda, where there has been sustained heterosexual transmission for many years, we found that a larger distance threshold is necessary to capture a more risk factor-diverse population with sparse sampling over a longer period of time. Such identification may allow for more informed intervention action by respective public health officials.

The origin and spread of HIV-1 CRF59_01B epidemic in China: A molecular network and phylogeographic analysis.

Human immunodeficiency virus type 1 CRF59_01B, identified in China in 2013, has been detected nationwide, exhibiting notably high prevalence in Guangzhou and its vicinity. This study aimed to unravel its origin and migration. A data set was established, incorporating all available CRF59_01B pol gene sequences and their metadata from Guangzhou and the public database. Bayesian phylogeographic analysis demonstrated that CRF59_01B originated in Shenzhen, the neighboring city of Guangzhou, around 1998 with posterior probability of 0.937. Molecular network analysis detected 1131 transmission links and showed a remarkably high clustering rate (78.9%). Substantial inter-city transmissions (26.5%, 300/1131) were observed between Shenzhen and Guangzhou while inter-region transmissions linked Guangzhou with South (46) and Southwest (64) China. The centre of Guangzhou was the hub of CRF59_01B transmission, including the inflow from Shenzhen (3.57 events/year) and outflow to the outskirts of Guangzhou (>2 events/year). The large-scale analysis revealed significant migration from Shenzhen to Guangzhou (5.08 events/year) and North China (0.59 events/year), and spread from Guangzhou to Central (0.47 events/year), East (0.42 events/year), South (0.76 events/year), Southwest China (0.76 events/year) and Shenzhen (1.89 events/year). Shenzhen and Guangzhou served as the origin and the hub of CRF59_01B circulation, emphasizing inter-city cooperation and data sharing to confine its nationwide diffusion.

Human Immunodeficiency Virus Type-1 Genetic Diversity and Drugs Resistance Mutations among People Living with HIV in Karachi, Pakistan.

The human immunodeficiency virus type-1 epidemic in Pakistan has significantly increased over the last two decades. In Karachi, Pakistan, there is a lack of updated information on the complexity of HIV-1 genetic diversity and the burden of drug resistance mutations (DRMs) that can contribute to ART failure and poor treatment outcomes. This study aimed to determine HIV-1 genetic diversity and identify drug-resistance mutations among people living with HIV in Karachi. A total of 364 HIV-positive individuals, with a median age of 36 years, were enrolled in the study. The HIV-1 partial pol gene was successfully sequenced from 268 individuals. The sequences were used to generate phylogenetic trees to determine clade diversity and also to assess the burden of DRMs. Based on the partial pol sequences, 13 distinct HIV-1 subtypes and recombinant forms were identified. Subtype A1 was the most common clade (40%), followed by CRF02_AG (33.2%). Acquired DRMs were found in 30.6% of the ART-experienced patients, of whom 70.7%, 20.7%, and 8.5% were associated with resistance to NNRTIs, NRTIs, and PIs, respectively. Transmitted DRMs were found in 5.6% of the ART-naïve patients, of whom 93% were associated with resistance against NNRTIs and 7% to PIs. The high prevalence of DRMs in ART-experienced patients poses significant challenges to the long-term benefits and sustainability of the ART program. This study emphasizes the importance of continuous HIV molecular epidemiology and drug resistance surveillance to support evidence-based HIV prevention, precise ART, and targeted AIDS care.

Molecular transmission network analysis of newly diagnosed HIV-1 infections in Nanjing from 2019 to 2021

ObjectiveThe objective of this study was to conduct a comprehensive analysis of the molecular transmission networks and transmitted drug resistance (TDR) patterns among individuals newly diagnosed with HIV-1 in Nanjing.MethodsPlasma samples were collected from newly diagnosed HIV patients in Nanjing between 2019 and 2021. The HIV pol gene was amplified, and the resulting sequences were utilized for determining TDR, identifying viral subtypes, and constructing molecular transmission network. Logistic regression analyses were employed to investigate the epidemiological characteristics associated with molecular transmission clusters.ResultsA total of 1161 HIV pol sequences were successfully extracted from newly diagnosed individuals, each accompanied by reliable epidemiologic information. The analysis revealed the presence of multiple HIV-1 subtypes, with CRF 07_BC (40.57%) and CRF01_AE (38.42%) being the most prevalent. Additionally, six other subtypes and unique recombinant forms (URFs) were identified. The prevalence of TDR among the newly diagnosed cases was 7.84% during the study period. Employing a genetic distance threshold of 1.50%, the construction of the molecular transmission network resulted in the identification of 137 clusters, encompassing 613 nodes, which accounted for approximately 52.80% of the cases. Multivariate analysis indicated that individuals within these clusters were more likely to be aged ≥ 60, unemployed, baseline CD4 cell count ≥ 200 cells/mm3, and infected with the CRF119_0107 (P < 0.05). Furthermore, the analysis of larger clusters revealed that individuals aged ≥ 60, peasants, those without TDR, and individuals infected with the CRF119_0107 were more likely to be part of these clusters.ConclusionsThis study revealed the high risk of local HIV transmission and high TDR prevalence in Nanjing, especially the rapid spread of CRF119_0107. It is crucial to implement targeted interventions for the molecular transmission clusters identified in this study to effectively control the HIV epidemic.