Displaying publications 41 - 57 of 57 in total

Abstract:
Sort:
  1. Teo YY, Sim X, Ong RT, Tan AK, Chen J, Tantoso E, et al.
    Genome Res, 2009 Nov;19(11):2154-62.
    PMID: 19700652 DOI: 10.1101/gr.095000.109
    The Singapore Genome Variation Project (SGVP) provides a publicly available resource of 1.6 million single nucleotide polymorphisms (SNPs) genotyped in 268 individuals from the Chinese, Malay, and Indian population groups in Southeast Asia. This online database catalogs information and summaries on genotype and phased haplotype data, including allele frequencies, assessment of linkage disequilibrium (LD), and recombination rates in a format similar to the International HapMap Project. Here, we introduce this resource and describe the analysis of human genomic variation upon agglomerating data from the HapMap and the Human Genome Diversity Project, providing useful insights into the population structure of the three major population groups in Asia. In addition, this resource also surveyed across the genome for variation in regional patterns of LD between the HapMap and SGVP populations, and for signatures of positive natural selection using two well-established metrics: iHS and XP-EHH. The raw and processed genetic data, together with all population genetic summaries, are publicly available for download and browsing through a web browser modeled with the Generic Genome Browser.
    Matched MeSH terms: Genomics/methods
  2. Abdullah-Zawawi MR, Ahmad-Nizammuddin NF, Govender N, Harun S, Mohd-Assaad N, Mohamed-Hussein ZA
    Sci Rep, 2021 10 04;11(1):19678.
    PMID: 34608238 DOI: 10.1038/s41598-021-99206-y
    Transcription factors (TFs) form the major class of regulatory genes and play key roles in multiple plant stress responses. In most eukaryotic plants, transcription factor (TF) families (WRKY, MADS-box and MYB) activate unique cellular-level abiotic and biotic stress-responsive strategies, which are considered as key determinants for defense and developmental processes. Arabidopsis and rice are two important representative model systems for dicot and monocot plants, respectively. A comprehensive comparative study on 101 OsWRKY, 34 OsMADS box and 122 OsMYB genes (rice genome) and, 71 AtWRKY, 66 AtMADS box and 144 AtMYB genes (Arabidopsis genome) showed various relationships among TFs across species. The phylogenetic analysis clustered WRKY, MADS-box and MYB TF family members into 10, 7 and 14 clades, respectively. All clades in WRKY and MYB TF families and almost half of the total number of clades in the MADS-box TF family are shared between both species. Chromosomal and gene structure analysis showed that the Arabidopsis-rice orthologous TF gene pairs were unevenly localized within their chromosomes whilst the distribution of exon-intron gene structure and motif conservation indicated plausible functional similarity in both species. The abiotic and biotic stress-responsive cis-regulatory element type and distribution patterns in the promoter regions of Arabidopsis and rice WRKY, MADS-box and MYB orthologous gene pairs provide better knowledge on their role as conserved regulators in both species. Co-expression network analysis showed the correlation between WRKY, MADs-box and MYB genes in each independent rice and Arabidopsis network indicating their role in stress responsiveness and developmental processes.
    Matched MeSH terms: Genomics/methods*
  3. Yu G, Hatta A, Periyannan S, Lagudah E, Wulff BBH
    Methods Mol Biol, 2017;1659:207-213.
    PMID: 28856653 DOI: 10.1007/978-1-4939-7249-4_18
    DNA is widely used in plant genetic and molecular biology studies. In this chapter, we describe how to extract DNA from wheat tissues. The tissue samples are ground to disrupt the cell wall. Then cetyltrimethylammonium bromide (CTAB) or sodium dodecyl sulfate (SDS) is used to disrupt the cell and nuclear membranes to release the DNA into solution. A reducing agent, β-mercaptoethanol, is added to break the disulfide bonds between the cysteine residues and to help remove the tanins and polyphenols. A high concentration of salt is employed to remove polysaccharides. Ethylenediaminetetraacetic acid (EDTA) stops DNase activity by chelating the magnesium ions. The nucleic acid solution is extracted with chloroform-isoamyl alcohol (24:1) or 6 M ammonium acetate. The DNA in aqueous phase is precipated with ethanol or isopropanol, which makes DNA less hydrophilic in the presence of sodium ions (Na+).
    Matched MeSH terms: Genomics/methods
  4. Doni F, Suhaimi NSM, Mispan MS, Fathurrahman F, Marzuki BM, Kusmoro J, et al.
    Int J Mol Sci, 2022 Jan 10;23(2).
    PMID: 35054923 DOI: 10.3390/ijms23020737
    Rice, the main staple food for about half of the world's population, has had the growth of its production stagnate in the last two decades. One of the ways to further improve rice production is to enhance the associations between rice plants and the microbiome that exists around, on, and inside the plant. This article reviews recent developments in understanding how microorganisms exert positive influences on plant growth, production, and health, focusing particularly on rice. A variety of microbial species and taxa reside in the rhizosphere and the phyllosphere of plants and also have multiple roles as symbiotic endophytes while living within plant tissues and even cells. They alter the morphology of host plants, enhance their growth, health, and yield, and reduce their vulnerability to biotic and abiotic stresses. The findings of both agronomic and molecular analysis show ways in which microorganisms regulate the growth, physiological traits, and molecular signaling within rice plants. However, many significant scientific questions remain to be resolved. Advancements in high-throughput multi-omics technologies can be used to elucidate mechanisms involved in microbial-rice plant associations. Prospectively, the use of microbial inoculants and associated approaches offers some new, cost-effective, and more eco-friendly practices for increasing rice production.
    Matched MeSH terms: Genomics/methods
  5. Wu D, Shen E, Jiang B, Feng Y, Tang W, Lao S, et al.
    Nat Commun, 2022 02 03;13(1):689.
    PMID: 35115514 DOI: 10.1038/s41467-022-28359-9
    As one of the great survivors of the plant kingdom, barnyard grasses (Echinochloa spp.) are the most noxious and common weeds in paddy ecosystems. Meanwhile, at least two Echinochloa species have been domesticated and cultivated as millets. In order to better understand the genomic forces driving the evolution of Echinochloa species toward weed and crop characteristics, we assemble genomes of three Echinochloa species (allohexaploid E. crus-galli and E. colona, and allotetraploid E. oryzicola) and re-sequence 737 accessions of barnyard grasses and millets from 16 rice-producing countries. Phylogenomic and comparative genomic analyses reveal the complex and reticulate evolution in the speciation of Echinochloa polyploids and provide evidence of constrained disease-related gene copy numbers in Echinochloa. A population-level investigation uncovers deep population differentiation for local adaptation, multiple target-site herbicide resistance mutations of barnyard grasses, and limited domestication of barnyard millets. Our results provide genomic insights into the dual roles of Echinochloa species as weeds and crops as well as essential resources for studying plant polyploidization, adaptation, precision weed control and millet improvements.
    Matched MeSH terms: Genomics/methods*
  6. Kato T, Azegami J, Yokomori A, Dohra H, El Enshasy HA, Park EY
    BMC Genomics, 2020 Apr 23;21(1):319.
    PMID: 32326906 DOI: 10.1186/s12864-020-6709-7
    BACKGROUND: Ashbya gossypii naturally overproduces riboflavin and has been utilized for industrial riboflavin production. To improve riboflavin production, various approaches have been developed. In this study, to investigate the change in metabolism of a riboflavin-overproducing mutant, namely, the W122032 strain (MT strain) that was isolated by disparity mutagenesis, genomic analysis was carried out.

    RESULTS: In the genomic analysis, 33 homozygous and 1377 heterozygous mutations in the coding sequences of the genome of MT strain were detected. Among these heterozygous mutations, the proportion of mutated reads in each gene was different, ranging from 21 to 75%. These results suggest that the MT strain may contain multiple nuclei containing different mutations. We tried to isolate haploid spores from the MT strain to prove its ploidy, but this strain did not sporulate under the conditions tested. Heterozygous mutations detected in genes which are important for sporulation likely contribute to the sporulation deficiency of the MT strain. Homozygous and heterozygous mutations were found in genes encoding enzymes involved in amino acid metabolism, the TCA cycle, purine and pyrimidine nucleotide metabolism and the DNA mismatch repair system. One homozygous mutation in AgILV2 gene encoding acetohydroxyacid synthase, which is also a flavoprotein in mitochondria, was found. Gene ontology (GO) enrichment analysis showed heterozygous mutations in all 22 DNA helicase genes and genes involved in oxidation-reduction process.

    CONCLUSION: This study suggests that oxidative stress and the aging of cells were involved in the riboflavin over-production in A. gossypii riboflavin over-producing mutant and provides new insights into riboflavin production in A. gossypii and the usefulness of disparity mutagenesis for the creation of new types of mutants for metabolic engineering.

    Matched MeSH terms: Genomics/methods*
  7. Patel VL, Busch EL, Friebel TM, Cronin A, Leslie G, McGuffog L, et al.
    Cancer Res, 2020 Feb 01;80(3):624-638.
    PMID: 31723001 DOI: 10.1158/0008-5472.CAN-19-1840
    Pathogenic sequence variants (PSV) in BRCA1 or BRCA2 (BRCA1/2) are associated with increased risk and severity of prostate cancer. We evaluated whether PSVs in BRCA1/2 were associated with risk of overall prostate cancer or high grade (Gleason 8+) prostate cancer using an international sample of 65 BRCA1 and 171 BRCA2 male PSV carriers with prostate cancer, and 3,388 BRCA1 and 2,880 BRCA2 male PSV carriers without prostate cancer. PSVs in the 3' region of BRCA2 (c.7914+) were significantly associated with elevated risk of prostate cancer compared with reference bin c.1001-c.7913 [HR = 1.78; 95% confidence interval (CI), 1.25-2.52; P = 0.001], as well as elevated risk of Gleason 8+ prostate cancer (HR = 3.11; 95% CI, 1.63-5.95; P = 0.001). c.756-c.1000 was also associated with elevated prostate cancer risk (HR = 2.83; 95% CI, 1.71-4.68; P = 0.00004) and elevated risk of Gleason 8+ prostate cancer (HR = 4.95; 95% CI, 2.12-11.54; P = 0.0002). No genotype-phenotype associations were detected for PSVs in BRCA1. These results demonstrate that specific BRCA2 PSVs may be associated with elevated risk of developing aggressive prostate cancer. SIGNIFICANCE: Aggressive prostate cancer risk in BRCA2 mutation carriers may vary according to the specific BRCA2 mutation inherited by the at-risk individual.
    Matched MeSH terms: Genomics/methods*
  8. van Velzen R, Holmer R, Bu F, Rutten L, van Zeijl A, Liu W, et al.
    Proc Natl Acad Sci U S A, 2018 May 15;115(20):E4700-E4709.
    PMID: 29717040 DOI: 10.1073/pnas.1721395115
    Nodules harboring nitrogen-fixing rhizobia are a well-known trait of legumes, but nodules also occur in other plant lineages, with rhizobia or the actinomycete Frankia as microsymbiont. It is generally assumed that nodulation evolved independently multiple times. However, molecular-genetic support for this hypothesis is lacking, as the genetic changes underlying nodule evolution remain elusive. We conducted genetic and comparative genomics studies by using Parasponia species (Cannabaceae), the only nonlegumes that can establish nitrogen-fixing nodules with rhizobium. Intergeneric crosses between Parasponia andersonii and its nonnodulating relative Trema tomentosa demonstrated that nodule organogenesis, but not intracellular infection, is a dominant genetic trait. Comparative transcriptomics of P. andersonii and the legume Medicago truncatula revealed utilization of at least 290 orthologous symbiosis genes in nodules. Among these are key genes that, in legumes, are essential for nodulation, including NODULE INCEPTION (NIN) and RHIZOBIUM-DIRECTED POLAR GROWTH (RPG). Comparative analysis of genomes from three Parasponia species and related nonnodulating plant species show evidence of parallel loss in nonnodulating species of putative orthologs of NIN, RPG, and NOD FACTOR PERCEPTION Parallel loss of these symbiosis genes indicates that these nonnodulating lineages lost the potential to nodulate. Taken together, our results challenge the view that nodulation evolved in parallel and raises the possibility that nodulation originated ∼100 Mya in a common ancestor of all nodulating plant species, but was subsequently lost in many descendant lineages. This will have profound implications for translational approaches aimed at engineering nitrogen-fixing nodules in crop plants.
    Matched MeSH terms: Genomics/methods*
  9. Shi W, Louzada S, Grigorova M, Massaia A, Arciero E, Kibena L, et al.
    Hum Mol Genet, 2019 Aug 15;28(16):2785-2798.
    PMID: 31108506 DOI: 10.1093/hmg/ddz101
    Human RBMY1 genes are located in four variable-sized clusters on the Y chromosome, expressed in male germ cells and possibly associated with sperm motility. We have re-investigated the mutational background and evolutionary history of the RBMY1 copy number distribution in worldwide samples and its relevance to sperm parameters in an Estonian cohort of idiopathic male factor infertility subjects. We estimated approximate RBMY1 copy numbers in 1218 1000 Genomes Project phase 3 males from sequencing read-depth, then chose 14 for valid ation by multicolour fibre-FISH. These fibre-FISH samples provided accurate calibration standards for the entire panel and led to detailed insights into population variation and mutational mechanisms. RBMY1 copy number worldwide ranged from 3 to 13 with a mode of 8. The two larger proximal clusters were the most variable, and additional duplications, deletions and inversions were detected. Placing the copy number estimates onto the published Y-SNP-based phylogeny of the same samples suggested a minimum of 562 mutational changes, translating to a mutation rate of 2.20 × 10-3 (95% CI 1.94 × 10-3 to 2.48 × 10-3) per father-to-son Y-transmission, higher than many short tandem repeat (Y-STRs), and showed no evidence for selection for increased or decreased copy number, but possible copy number stabilizing selection. An analysis of RBMY1 copy numbers among 376 infertility subjects failed to replicate a previously reported association with sperm motility and showed no significant effect on sperm count and concentration, serum follicle stimulating hormone (FSH), luteinizing hormone (LH) and testosterone levels or testicular and semen volume. These results provide the first in-depth insights into the structural rearrangements underlying RBMY1 copy number variation across diverse human lineages.
    Matched MeSH terms: Genomics/methods
  10. Graham NS, Hammond JP, Lysenko A, Mayes S, O Lochlainn S, Blasco B, et al.
    Plant Cell, 2014 Jul;26(7):2818-30.
    PMID: 25082855 DOI: 10.1105/tpc.114.128603
    Although Ca transport in plants is highly complex, the overexpression of vacuolar Ca(2+) transporters in crops is a promising new technology to improve dietary Ca supplies through biofortification. Here, we sought to identify novel targets for increasing plant Ca accumulation using genetical and comparative genomics. Expression quantitative trait locus (eQTL) mapping to 1895 cis- and 8015 trans-loci were identified in shoots of an inbred mapping population of Brassica rapa (IMB211 × R500); 23 cis- and 948 trans-eQTLs responded specifically to altered Ca supply. eQTLs were screened for functional significance using a large database of shoot Ca concentration phenotypes of Arabidopsis thaliana. From 31 Arabidopsis gene identifiers tagged to robust shoot Ca concentration phenotypes, 21 mapped to 27 B. rapa eQTLs, including orthologs of the Ca(2+) transporters At-CAX1 and At-ACA8. Two of three independent missense mutants of BraA.cax1a, isolated previously by targeting induced local lesions in genomes, have allele-specific shoot Ca concentration phenotypes compared with their segregating wild types. BraA.CAX1a is a promising target for altering the Ca composition of Brassica, consistent with prior knowledge from Arabidopsis. We conclude that multiple-environment eQTL analysis of complex crop genomes combined with comparative genomics is a powerful technique for novel gene identification/prioritization.
    Matched MeSH terms: Genomics/methods*
  11. Barnard RT
    Expert Rev Vaccines, 2010 May;9(5):461-3.
    PMID: 20450319 DOI: 10.1586/erv.10.48
    The Recombinant Vaccines: Strategies for Candidate Discovery and Vaccine Delivery conference, organized by EuroSciCon, hosted a group of UK-based and international scientists from as far afield as Malaysia and Australia. Genomic analyses of pathogens and elucidation of mechanisms of pathogenesis has advanced candidate discovery and development of vaccines. Therefore, it was timely that this conference featured, in addition to detailed expositions of target selection and clinical trials, presentations on manufacturability, scale-up and delivery of vaccines. Ten talks were presented. This meeting report describes the key topics presented and the themes that emerged from this conference.
    Matched MeSH terms: Genomics/methods
  12. Masomian M, Rahman RN, Salleh AB, Basri M
    PLoS One, 2016;11(3):e0149851.
    PMID: 26934700 DOI: 10.1371/journal.pone.0149851
    Thermostable and organic solvent-tolerant enzymes have significant potential in a wide range of synthetic reactions in industry due to their inherent stability at high temperatures and their ability to endure harsh organic solvents. In this study, a novel gene encoding a true lipase was isolated by construction of a genomic DNA library of thermophilic Aneurinibacillus thermoaerophilus strain HZ into Escherichia coli plasmid vector. Sequence analysis revealed that HZ lipase had 62% identity to putative lipase from Bacillus pseudomycoides. The closely characterized lipases to the HZ lipase gene are from thermostable Bacillus and Geobacillus lipases belonging to the subfamily I.5 with ≤ 57% identity. The amino acid sequence analysis of HZ lipase determined a conserved pentapeptide containing the active serine, GHSMG and a Ca(2+)-binding motif, GCYGSD in the enzyme. Protein structure modeling showed that HZ lipase consisted of an α/β hydrolase fold and a lid domain. Protein sequence alignment, conserved regions analysis, clustal distance matrix and amino acid composition illustrated differences between HZ lipase and other thermostable lipases. Phylogenetic analysis revealed that this lipase represented a new subfamily of family I of bacterial true lipases, classified as family I.9. The HZ lipase was expressed under promoter Plac using IPTG and was characterized. The recombinant enzyme showed optimal activity at 65 °C and retained ≥ 97% activity after incubation at 50 °C for 1h. The HZ lipase was stable in various polar and non-polar organic solvents.
    Matched MeSH terms: Genomics/methods
  13. Chung FF, Tan PF, Raja VJ, Tan BS, Lim KH, Kam TS, et al.
    Sci Rep, 2017 02 15;7:42504.
    PMID: 28198434 DOI: 10.1038/srep42504
    Precursor mRNA (pre-mRNA) splicing is catalyzed by a large ribonucleoprotein complex known as the spliceosome. Numerous studies have indicated that aberrant splicing patterns or mutations in spliceosome components, including the splicing factor 3b subunit 1 (SF3B1), are associated with hallmark cancer phenotypes. This has led to the identification and development of small molecules with spliceosome-modulating activity as potential anticancer agents. Jerantinine A (JA) is a novel indole alkaloid which displays potent anti-proliferative activities against human cancer cell lines by inhibiting tubulin polymerization and inducing G2/M cell cycle arrest. Using a combined pooled-genome wide shRNA library screen and global proteomic profiling, we showed that JA targets the spliceosome by up-regulating SF3B1 and SF3B3 protein in breast cancer cells. Notably, JA induced significant tumor-specific cell death and a significant increase in unspliced pre-mRNAs. In contrast, depletion of endogenous SF3B1 abrogated the apoptotic effects, but not the G2/M cell cycle arrest induced by JA. Further analyses showed that JA stabilizes endogenous SF3B1 protein in breast cancer cells and induced dissociation of the protein from the nucleosome complex. Together, these results demonstrate that JA exerts its antitumor activity by targeting SF3B1 and SF3B3 in addition to its reported targeting of tubulin polymerization.
    Matched MeSH terms: Genomics/methods
  14. Kanniappan P, Ahmed SA, Rajasekaram G, Marimuthu C, Ch'ng ES, Lee LP, et al.
    J Cell Mol Med, 2017 10;21(10):2276-2283.
    PMID: 28756649 DOI: 10.1111/jcmm.13148
    Technological advances in RNA biology greatly improved transcriptome profiling during the last two decades. Besides the discovery of many small RNAs (sRNA) that are involved in the physiological and pathophysiological regulation of various cellular circuits, it becomes evident that the corresponding RNA genes might also serve as potential biomarkers to monitor the progression of disease and treatment. sRNA gene candidate npcTB_6715 was previously identified via experimental RNomic (unpublished data), and we report its application as potential biomarker for the detection of Mycobacterium tuberculosis (MTB) in patient samples. For proof of principle, we developed a multiplex PCR assay and report its validation with 500 clinical cultures, positive for Mycobacteria. The analysis revealed 98.9% sensitivity, 96.1% specificity, positive and negative predictive values of 98.6% and 96.8%, respectively. These results underscore the diagnostic value of the sRNA gene as diagnostic marker for the specific detection of MTB in clinical samples. Its successful application and the general ease of PCR-based detection compared to standard bacterial culture techniques might be the first step towards 'point-of-care' diagnostics of Mycobacteria. To the best of our knowledge, this is the first time for the design of diagnostic applications based on sRNA genes, in Mycobacteria.
    Matched MeSH terms: Genomics/methods
  15. Saunus JM, Quinn MC, Patch AM, Pearson JV, Bailey PJ, Nones K, et al.
    J Pathol, 2015 Nov;237(3):363-78.
    PMID: 26172396 DOI: 10.1002/path.4583
    Treatment options for patients with brain metastases (BMs) have limited efficacy and the mortality rate is virtually 100%. Targeted therapy is critically under-utilized, and our understanding of mechanisms underpinning metastatic outgrowth in the brain is limited. To address these deficiencies, we investigated the genomic and transcriptomic landscapes of 36 BMs from breast, lung, melanoma and oesophageal cancers, using DNA copy-number analysis and exome- and RNA-sequencing. The key findings were as follows. (a) Identification of novel candidates with possible roles in BM development, including the significantly mutated genes DSC2, ST7, PIK3R1 and SMC5, and the DNA repair, ERBB-HER signalling, axon guidance and protein kinase-A signalling pathways. (b) Mutational signature analysis was applied to successfully identify the primary cancer type for two BMs with unknown origins. (c) Actionable genomic alterations were identified in 31/36 BMs (86%); in one case we retrospectively identified ERBB2 amplification representing apparent HER2 status conversion, then confirmed progressive enrichment for HER2-positivity across four consecutive metastatic deposits by IHC and SISH, resulting in the deployment of HER2-targeted therapy for the patient. (d) In the ERBB/HER pathway, ERBB2 expression correlated with ERBB3 (r(2)  = 0.496; p < 0.0001) and HER3 and HER4 were frequently activated in an independent cohort of 167 archival BM from seven primary cancer types: 57.6% and 52.6% of cases were phospho-HER3(Y1222) or phospho-HER4(Y1162) membrane-positive, respectively. The HER3 ligands NRG1/2 were barely detectable by RNAseq, with NRG1 (8p12) genomic loss in 63.6% breast cancer-BMs, suggesting a microenvironmental source of ligand. In summary, this is the first study to characterize the genomic landscapes of BM. The data revealed novel candidates, potential clinical applications for genomic profiling of resectable BMs, and highlighted the possibility of therapeutically targeting HER3, which is broadly over-expressed and activated in BMs, independent of primary site and systemic therapy.
    Matched MeSH terms: Genomics/methods*
  16. Yap PSX, Ahmad Kamar A, Chong CW, Ngoi ST, Teh CSJ
    Microb Drug Resist, 2020 Mar;26(3):190-203.
    PMID: 31545116 DOI: 10.1089/mdr.2019.0199
    Background:
    Klebsiella pneumoniae is a major opportunistic pathogen frequently associated with nosocomial infections, and often poses a major threat to immunocompromised patients. In our previous study, two K. pneumoniae (K36 and B13), which displayed resistance to almost all major antibiotics, including colistin, were isolated. Both isolates were not associated with infection and isolated from the stools of two preterm neonates admitted to the neonatal intensive care unit (NICU) during their first week of life.
    Materials and Methods:
    In this study, whole genome sequencing was performed on these two clinical multidrug resistant K. pneumoniae. We aimed to determine the genetic factors that underline the antibiotic-resistance phenotypes of these isolates.
    Results:
    The strains harbored blaSHV-27, blaSHV-71, and oqxAB genes conferring resistance to cephalosporins, carbapenems, and fluoroquinolones, respectively, but not harboring any known plasmid-borne colistin resistance determinants such as mcr-1. However, genome analysis discovered interruption of mgrB gene by insertion sequences gaining insight into the development of colistin resistance.
    Conclusion:
    The observed finding that points to a scenario of potential gut-associated resistance genes to Gram negative (K. pneumoniae) host in the NICU environment warrants attention and further investigation.
    Matched MeSH terms: Genomics/methods
  17. Jelen MM, Chen Z, Kocjan BJ, Burt FJ, Chan PK, Chouhy D, et al.
    J Virol, 2014 Jul;88(13):7307-16.
    PMID: 24741079 DOI: 10.1128/JVI.00621-14
    Human papillomavirus type 6 (HPV6) is the major etiological agent of anogenital warts and laryngeal papillomas and has been included in both the quadrivalent and nonavalent prophylactic HPV vaccines. This study investigated the global genomic diversity of HPV6, using 724 isolates and 190 complete genomes from six continents, and the association of HPV6 genomic variants with geographical location, anatomical site of infection/disease, and gender. Initially, a 2,800-bp E5a-E5b-L1-LCR fragment was sequenced from 492/530 (92.8%) HPV6-positive samples collected for this study. Among them, 130 exhibited at least one single nucleotide polymorphism (SNP), indel, or amino acid change in the E5a-E5b-L1-LCR fragment and were sequenced in full. A global alignment and maximum likelihood tree of 190 complete HPV6 genomes (130 fully sequenced in this study and 60 obtained from sequence repositories) revealed two variant lineages, A and B, and five B sublineages: B1, B2, B3, B4, and B5. HPV6 (sub)lineage-specific SNPs and a 960-bp representative region for whole-genome-based phylogenetic clustering within the L2 open reading frame were identified. Multivariate logistic regression analysis revealed that lineage B predominated globally. Sublineage B3 was more common in Africa and North and South America, and lineage A was more common in Asia. Sublineages B1 and B3 were associated with anogenital infections, indicating a potential lesion-specific predilection of some HPV6 sublineages. Females had higher odds for infection with sublineage B3 than males. In conclusion, a global HPV6 phylogenetic analysis revealed the existence of two variant lineages and five sublineages, showing some degree of ethnogeographic, gender, and/or disease predilection in their distribution.

    IMPORTANCE: This study established the largest database of globally circulating HPV6 genomic variants and contributed a total of 130 new, complete HPV6 genome sequences to available sequence repositories. Two HPV6 variant lineages and five sublineages were identified and showed some degree of association with geographical location, anatomical site of infection/disease, and/or gender. We additionally identified several HPV6 lineage- and sublineage-specific SNPs to facilitate the identification of HPV6 variants and determined a representative region within the L2 gene that is suitable for HPV6 whole-genome-based phylogenetic analysis. This study complements and significantly expands the current knowledge of HPV6 genetic diversity and forms a comprehensive basis for future epidemiological, evolutionary, functional, pathogenicity, vaccination, and molecular assay development studies.

    Matched MeSH terms: Genomics/methods
Filters
Contact Us

Please provide feedback to Administrator (afdal@afpm.org.my)

External Links