Coding Region Research Articles

Draft genome sequence of multidrug resistant Proteus mirabilis strain MAD23

Genome sequencing became one of the important tools for diagnosis of microorganisms. Proteus mirabilis is an enteric Gramnegative bacterium causes gastrointestinal and urinary tract infections in human and animals. Here, the draft sequencing of Proteus mirabilis (MAD23) strain was performed. The strain was collected from Chinese duck cecum suffering from diarrhea. Bacteria was identified by biochemical tests and characteristic swarming pattern on the plate. Antimicrobial susceptibility test was performed by VITEC-2 MS system. Also, draft genome sequencing was conducted using the Illumina platform NovaSeq sequencer. The genome is about 3,729,695 bp long, GC content is of 38.9%., number of contigs is 123, largest contig is 137547pb, smallest contig is 513, total of protein-coding sequences (CDSs) is 3410, and number of RNA genes is 72. Genome Pairwise comparisons of Proteus mirabilis MAD23 with other vs. type strain was detected based on is DDH, GC content, δ- value, genome size and number of proteins. Neighbour-Joining phylogenetic trees show relationship between Proteus mirabilis MAD23 and the close strains of Proteus species based on 16S rRNA sequences. In addition, analysis of the sequence reveals presence of many antibiotic resistance genes as dfrA1, sul2, dfrA32, aadA2b, aph(3’)-Ia, aadA1, aph(6)-Id, aph(3’’)-Ib, ere(A), aadA1, aadA2b, tet(J), tet(C), floR, and Cat genes.These genes encoding the resistance to different antibiotics including Aminoglycosides, Macrolides, Aminocyclitol, Tetracycline and Amphenicol. This work is considered as an important study can help in more understanding of WGS in Proteus spp. and recognize the genes that related to antibiotic resistance in this pathogen.

MHC-I-presented non-canonical antigens expand the cancer immunotherapy targets in acute myeloid leukemia

Identification of tumor neoantigens is indispensable for the development of cancer immunotherapies. However, we are still lacking knowledge about the potential neoantigens derived from sequences outside protein-coding regions. Here, we comprehensively characterized the immunopeptidome landscape by integrating multi-omics data in acute myeloid leukemia (AML). Both canonical and non-canonical MHC-associated peptides (MAPs) in AML were identified. We found that the quality and characteristics of ncMAPs are comparable or superior to cMAPs, suggesting ncMAPs are indispensable sources for tumor neoantigens. We further proposed a computational framework to prioritize the neoantigens by integrating additional transcriptome and immunopeptidome in normal tissues. Notably, 6 of prioritized 13 neoantigens were derived from ncMAPs. The expressions of corresponding source genes are highly related to infiltrations of immune cells. Finally, a risk model was developed, which exhibited good performance for clinical prognosis in AML. Our findings expand potential cancer immunotherapy targets and provide in-depth insights into AML treatment, laying a new foundation for precision therapies in AML.

Complete genome sequence of Bacillus subtilis phage fHSPT3.

Here, we describe a Bacillus subtilis bacteriophage isolated from sewage water. The phage fHSPT3 was isolated against the host Bacillus subtilis 168 and has a genome size of 150,187 bp with 221 protein-coding sequences. The fHSPT3 belongs to the genus Siophivirus, lacks virulence or antibiotic resistance genes, and shows a virulent life cycle predicted through PhageScope and PhageAI.

The first complete mitochondrial genome of Grossulariaceae: Molecular features, structure recombination, and genetic evolution

BackgroundMitochondria play crucial roles in the growth, development, and adaptation of plants. Blackcurrant (Ribes nigrum L.) stands out as a significant berry species due to its rich nutritional profile, medicinal properties, and health benefits. Despite its importance, the mitochondrial genome of blackcurrant remains unassembled.ResultsThis study presents the first assembly of the mitochondrial genome of R. nigrum in the Grossulariaceae family. The genome spans 450,227 base pairs (bp) and encompasses 39 protein-coding genes (PCGs), 19 transfer RNAs (tRNAs), and three ribosomal RNAs (rRNAs). Protein-coding regions constitute 8.88% of the entire genome. Additionally, we identified 180 simple sequence repeats, 12 tandem repeats, and 432 pairs of dispersed repeats. Notably, the dispersed sequence R1 (cotig3, 1,129 bp) mediated genome recombination, resulting in the formation of two major conformations, namely master and double circles. Furthermore, we identified 731 C-to-U RNA editing sites within the PCGs. Among these, cox1-2, nad1-2, and nad4L-2 were associated with the creation of start codons, whereas atp6-718 and rps10-391 were linked to termination codons. We also detected fourteen plastome fragments within the mitogenome, constituting 1.11% of the total length. Phylogenetic analysis suggests that R. nigrum might have undergone multiple genomic reorganization and/or gene transfer events, resulting in the loss of two PCGs (rps2 and rps11) during its evolutionary history.ConclusionsThis investigation unveils the molecular characteristics of the R. nigrum mitogenome, shedding light on its evolutionary trajectory and phylogenetic implications. Furthermore, it serves as a valuable reference for evolutionary research and germplasm identification within the genus.

Features of the CD1 gene family in rodents and the uniqueness of the immune system of naked mole-rat

Cluster of Differentiation 1 (CD1) proteins are widely expressed throughout jawed vertebrates and present lipid antigens to specific CD1-restricted T lymphocytes. CD1 molecules play an important role in immune defense with the presence or absence of particular CD1 proteins frequently associated with the functional characteristics of the immune system. Here, we show the evolution of CD1 proteins in the Rodentia family and the diversity among its members. Based on the analysis of CD1 protein-coding regions in rodent genomes and the reconstruction of protein structures, we found that Heterocephalus glaber represents a unique member of the suborder Hystricomorpha with significant changes in protein sequences and structures of the CD1 family. Multiple lines of evidence point to the absence of CD1d and CD1e and probably a dysfunctional CD1b protein in Heterocephalus glaber. In addition, the impact of CD1d loss on the CD1d/Natural killer T (NKT) cell axis in the naked mole-rat and its potential implications for immune system function are discussed in detail.

Explorer: efficient DNA coding by De Bruijn graph toward arbitrary local and global biochemical constraints.

With the exponential growth of digital data, there is a pressing need for innovative storage media and techniques. DNA molecules, due to their stability, storage capacity, and density, offer a promising solution for information storage. However, DNA storage also faces numerous challenges, such as complex biochemical constraints and encoding efficiency. This paper presents Explorer, a high-efficiency DNA coding algorithm based on the De Bruijn graph, which leverages its capability to characterize local sequences. Explorer enables coding under various biochemical constraints, such as homopolymers, GC content, and undesired motifs. This paper also introduces Codeformer, a fast decoding algorithm based on the transformer architecture, to further enhance decoding efficiency. Numerical experiments indicate that, compared with other advanced algorithms, Explorer not only achieves stable encoding and decoding under various biochemical constraints but also increases the encoding efficiency and bit rate by ¿10%. Additionally, Codeformer demonstrates the ability to efficiently decode large quantities of DNA sequences. Under different parameter settings, its decoding efficiency exceeds that of traditional algorithms by more than two-fold. When Codeformer is combined with Reed-Solomon code, its decoding accuracy exceeds 99%, making it a good choice for high-speed decoding applications. These advancements are expected to contribute to the development of DNA-based data storage systems and the broader exploration of DNA as a novel information storage medium.

Parkinson's Disease Gene Screening in Familial Cases from Central and South America.

Parkinson's disease (PD) is the second most common neurodegenerative disease following Alzheimer's disease. Nearly 30 causative genes have been identified for PD and related disorders. However, most of these genes were identified in European-derived families, and little is known about their role in Latin American populations. Our goal was to assess the spectrum and frequency of pathogenic variants in known PD genes in familial PD patients from Latin America. We selected 335 PD patients with a family history of PD from the Latin American Research Consortium on the Genetics of PD. We capture-sequenced the coding regions of 26 genes related to neurodegenerative parkinsonism. Of the 335 PD patients, 324 had sufficient sequencing coverage to be analyzed. We identified pathogenic variants in 41 individuals (12.7%) in FBXO7, GCH1, LRRK2, PARK7, PINK1, PLA2G6, PRKN, SNCA, and TARDBP, GBA1 risk variants in 25 individuals (7.7%), and variants of uncertain significance in another 24 individuals (7.4%) in ATP13A2, ATP1A3, DNAJC13, DNAJC6, GBA1, LRKK2, PINK1, VPS13C, and VPS35. Of the 70 unique variants identified, 19 were more frequent in Latin Americans than in any other population. This is the first screening of known PD genes in a large cohort of patients with familial PD from Latin America. There were substantial differences in the spectrum of variants observed in comparison to previous findings from PD families of European origin. Our data provide further evidence that differences exist between the genetic architecture of PD in Latinos and European-derived populations. © 2024 The Author(s). Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Viral hijacking of hnRNPH1 unveils a G-quadruplex-driven mechanism of stress control.

Viral genomes are enriched with G-quadruplexes (G4s), non-canonical structures formed in DNA or RNA upon assembly of four guanine stretches into stacked quartets. Because of their critical roles, G4s are potential antiviral targets, yet their function remains largely unknown. Here, we characterize the formation and functions of a conserved G4 within the polymerase coding region of orthoflaviviruses of the Flaviviridae family. Using yellow fever virus, we determine that this G4 promotes viral replication and suppresses host stress responses via interactions with hnRNPH1, a host nuclear protein involved in RNA processing. G4 binding to hnRNPH1 causes its cytoplasmic retention with subsequent impacts on G4-containing tRNA fragments (tiRNAs) involved in stress-mediated reductions in translation. As a result, these host stress responses and associated antiviral effects are impaired. These data reveal that the interplay between hnRNPH1 and both host and viral G4 targets controls the integrated stress response and viral replication.

Neuropathologically directed profiling of PRNP somatic and germline variants in sporadic human prion disease.

Creutzfeldt-Jakob Disease (CJD), the most common human prion disease, is associated with pathologic misfolding of the prion protein (PrP), encoded by the PRNP gene. Of human prion disease cases,<1% were transmitted by misfolded PrP, ~ 15% are inherited, and ~ 85% are sporadic (sCJD). While familial cases are inherited through germline mutations in PRNP, the cause of sCJD is unknown. Somatic mutations have been hypothesized as a cause of sCJD, and recent studies have revealed that somatic mutations accumulate in neurons during aging. To investigate the hypothesis that somatic mutations in PRNP may underlie sCJD, we performed deep DNA sequencing of PRNP in 205 sCJD cases and 170 age-matched non-disease controls. We included 5 cases of Heidenhain variant sporadic CJD (H-sCJD), where visual symptomatology and neuropathology implicate localized initiation of prion formation, and examined multiple regions across the brain including in the affected occipital cortex. We employed Multiple Independent Primer PCR Sequencing (MIPP-Seq) with a median depth of > 5000× across the PRNP coding region and analyzed for variants using MosaicHunter. An allele mixing experiment showed positive detection of variants in bulk DNA at a variant allele fraction (VAF) as low as 0.2%. We observed multiple polymorphic germline variants among individuals in our cohort. However, we did not identify bona fide somatic variants in sCJD, including across multiple affected regions in H-sCJD, nor in control individuals. Beyond our stringent variant-identification pipeline, we also analyzed VAFs from raw sequencing data, and observed no evidence of prion disease enrichment for the known germline pathogenic variants P102L, D178N, and E200K. The lack of PRNP pathogenic somatic mutations in H-sCJD or the broader cohort of sCJD suggests that clonal somatic mutations may not play a major role in sporadic prion disease. With H-sCJD representing a localized presentation of neurodegeneration, this serves as a test of the potential role of clonal somatic mutations in genes known to cause familial neurodegeneration.

Identification of a novel neurovirulence factor encoded by the cryptic orphan gene UL31.6 of herpes simplex virus 1.

Although the herpes simplex virus type 1 (HSV-1) genome was thought to contain approximately 80 different protein coding sequences (CDSs), recent multi-omics analyses reported HSV-1 encodes more than 200 potential CDSs. However, few of the newly identified CDSs were confirmed to be expressed at the peptide or protein level in HSV-1-infected cells. Furthermore, the impact of the proteins they encode on HSV-1 infection is largely unknown. This study focused on a newly identified CDS, UL31.6. Re-analyzation of our previous chemical proteomics data verified that UL31.6 was expressed at the peptide level in HSV-1-infected cells. Antisera raised against a viral protein encoded by UL31.6 (pUL31.6) reacted with a protein with an approximate molecular mass of 37 kDa in lysates of Vero cells infected with each of three HSV-1 strains. pUL31.6 was efficiently dissociated from virions in high-salt solution. A UL31.6-null mutation had a minimal effect on HSV-1 gene expression, replication, cell-to-cell spread, and morphogenesis in Vero cells; in contrast, it significantly reduced HSV-1 cell-to-cell spread in three neural cells but not in four non-neural cells including Vero cells. The UL31.6-null mutation also significantly reduced the mortality and viral replication in the brains of mice after intracranial infection, but had minimal effects on pathogenic manifestations in and around the eyes, and viral replication detected in the tear films of mice after ocular infection. These results indicated that pUL31.6 was a tegument protein and specifically acted as a neurovirulence factor by potentially promoting viral transmission between neuronal cells in the central nervous system.IMPORTANCERecent multi-omics analyses reported the herpes simplex virus type 1 (HSV-1) genome encodes an additional number of potential coding sequences (CDSs). However, the expressions of these CDSs at the peptide or protein levels and the biological effects of these CDSs on HSV-1 infection remain largely unknown. This study annotated a cryptic orphan CDS, termed UL31.6, an HSV-1 gene that encodes a tegument protein with an approximate molecular mass of 37 kDa, which specifically acts as a neurovirulence factor. Our study indicates that HSV-1 proteins important for viral pathogenesis remain to be identified and a comprehensive understanding of the pathogenesis of HSV-1 will require not only the identification of cryptic orphan CDSs using emerging technologies but also step-by-step and in-depth analyses of each of the cryptic orphan CDSs.

Global Update on Measles Molecular Epidemiology.

Molecular surveillance of circulating measles variants serves as a line of evidence for the absence of endemic circulation and provides a means to track chains of transmission. Molecular surveillance for measles (genotyping) is based on the sequence of 450 nucleotides at the end of the nucleoprotein coding region (N450) of the measles genome. Genotyping was established in 1998 and, with over 50,000 sequence submissions to the Measles Nucleotide Surveillance database, has proven to be an effective resource for countries attempting to trace pathways of transmission. This review summarizes the tools used for the molecular surveillance of measles and describes the challenge posed by the decreased number of circulating measles genotypes. The Global Measles and Rubella Laboratory Network addressed this challenge through the development of new tools such as named strains and distinct sequence identifiers that analyze the diversity within the currently circulating genotypes. The advantages and limitations of these approaches are discussed, together with the need to generate additional sequence data including whole genome sequences to ensure the continued utility of strain surveillance for measles.

Complete genome sequence of Escherichia coli O157:H7 phage Φ241.

We report the genome sequence of phage Φ241 infecting Escherichia coli O157:H7. Phage Φ241 was isolated from an industrial cucumber fermentation at high acidity (pH 3.7) and high salinity (5% NaCl). The phage genome consists of a 157,291 bp circular double-stranded DNA with 203 coding regions and 44.96% GC content.

African Swine Fever Virus Protein-Protein Interaction Prediction.

The African swine fever virus (ASFV) is an often deadly disease in swine and poses a threat to swine livestock and swine producers. With its complex genome containing more than 150 coding regions, developing effective vaccines for this virus remains a challenge due to a lack of basic knowledge about viral protein function and protein-protein interactions between viral proteins and between viral and host proteins. In this work, we identified ASFV-ASFV protein-protein interactions (PPIs) using artificial intelligence-powered protein structure prediction tools. We benchmarked our PPI identification workflow on the Vaccinia virus, a widely studied nucleocytoplasmic large DNA virus, and found that it could identify gold-standard PPIs that have been validated in vitro in a genome-wide computational screening. We applied this workflow to more than 18,000 pairwise combinations of ASFV proteins and were able to identify seventeen novel PPIs, many of which have corroborating experimental or bioinformatic evidence for their protein-protein interactions, further validating their relevance. Two protein-protein interactions, I267L and I8L, I267L__I8L, and B175L and DP79L, B175L__DP79L, are novel PPIs involving viral proteins known to modulate host immune response.

'Drifting' Buchnera genomes track the microevolutionary trajectories of their aphid hosts.

Evolution of Buchnera-aphid host symbioses is often studied among species at macroevolutionary scales. Investigations within species offer a different perspective about how eco-evolutionary processes shape patterns of genetic variation at microevolutionary scales. Our study leverages new and publicly available whole-genome sequencing data to study Buchnera-aphid host evolution in Myzus persicae, the peach potato aphid, a globally invasive and polyphagous pest. Across 43 different asexual, clonally reproducing isofemale strains, we examined patterns of genomic covariation between Buchnera and their aphid host and considered the distribution of mutations in protein-coding regions of the Buchnera genome. We found Buchnera polymorphisms within aphid strains, suggesting the presence of genetically different Buchnera strains within the same clonal lineage. Genetic distance between pairs of Buchnera samples was positively correlated to genetic distance between their aphid hosts, indicating shared evolutionary histories. However, there was no segregation of genetic variation for both M. persicae and Buchnera with plant host (Brassicaceae and non-tobacco Solanaceae) and no associations between genetic and geographic distance at global or regional spatial scales. Abundance patterns of non-synonymous mutations were similar to synonymous mutations in the Buchnera genome, and both mutation classes had similar site frequency spectra. We hypothesize that a predominance of neutral processes results in the Buchnera of M. persicae to simply 'drift' with the evolutionary trajectory of their aphid hosts. Our study presents a unique microevolutionary characterization of Buchnera-aphid host genomic covariation across multiple aphid clones. This provides a new perspective on the eco-evolutionary processes generating and maintaining polymorphisms in a major pest aphid species and its obligate primary endosymbiont.

Diverse origins of near-identical antifreeze proteins in unrelated fish lineages provide insights into evolutionary mechanisms of new gene birth and protein sequence convergence.

Determining the origins of novel genes and the genetic mechanisms underlying the emergence of new functions is challenging yet crucial for understanding evolutionary innovations. The convergently evolved fish antifreeze proteins provide excellent opportunities to investigate evolutionary origins and pathways of new genes. Particularly notable is the near-identical type I antifreeze proteins (AFPI) in four phylogenetically divergent fish taxa. This study tested the hypothesis of protein sequence convergence beyond functional convergence in three unrelated AFPI-bearing fish lineages, revealing different paths by which a similar protein arose from diverse genomic resources. Comprehensive comparative analyses of de novo sequenced genome of the winter flounder and grubby sculpin, available high-quality genome of the cunner and 14 other relevant species found that the near-identical AFPI originated from a distinct genetic precursor in each lineage. Each independently evolved a coding region for the novel ice-binding protein while retaining sequence identity in the regulatory regions with their respective ancestor. The deduced evolutionary processes and molecular mechanisms are consistent with the Innovation-Amplification-Divergence (IAD) model applicable to AFPI formation in all three lineages, a new Duplication-Degeneration-Divergence (DDD) model we propose for the sculpin lineage, and a DDD model with gene fission for the cunner lineage. This investigation illustrates the multiple ways by which a novel functional gene with sequence convergence at the protein level could evolve across divergent species, advancing our understanding of the mechanistic intricacies in new gene formation.

Comparative Analysis of Plastome Sequences of Seven Tulipa L. (Liliaceae Juss.) Species from Section Kolpakowskianae Raamsd. Ex Zonn and Veldk.

Tulipa L. is a genus of significant economic, environmental, and cultural importance in several parts of the world. The exact number of species in the genus remains uncertain due to inherent taxonomic challenges. We utilized next-generation sequencing technology to sequence and assemble the plastid genomes of seven Tulipa species collected in Kazakhstan and conducted a comparative analysis. The total number of annotated genes was 136 in all seven studied Tulipa species, 114 of which were unique, including 80 protein-coding, 30 tRNA, and 4 rRNA genes. Nine regions (petD, ndhH, ycf2-ycf3, ndhA, rpl16, clpP, ndhD-ndhF, rpoC2, and ycf1) demonstrated significant nucleotide variability, suggesting their potential as molecular markers. A total of 1388 SSRs were identified in the seven Tulipa plastomes, with mononucleotide repeats being the most abundant (60.09%), followed by dinucleotide (34.44%), tetranucleotide (3.90%), trinucleotide (1.08%), pentanucleotide (0.22%), and hexanucleotide (0.29%). The Ka/Ks values of the protein-coding genes ranged from 0 to 3.9286, with the majority showing values <1. Phylogenetic analysis based on a complete plastid genome and protein-coding gene sequences divided the species into three major clades corresponding to their subgenera. The results obtained in this study may contribute to understanding the phylogenetic relationships and molecular taxonomy of Tulipa species.

Expanding the genetics and phenotypes of ocular congenital cranial dysinnervation disorders

PurposeTo identify genetic etiologies and genotype/phenotype associations for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs). MethodsWe coupled phenotyping with exome or genome sequencing of 467 probands (550 affected and 1108 total individuals) with genetically unsolved oCCDDs, integrating analyses of pedigrees, human and animal model phenotypes, and de novo variants to identify rare candidate single-nucleotide variants, insertion/deletions, and structural variants disrupting protein-coding regions. Prioritized variants were classified for pathogenicity and evaluated for genotype/phenotype correlations. ResultsAnalyses elucidated phenotypic subgroups, identified pathogenic/likely pathogenic variant(s) in 43 of 467 probands (9.2%), and prioritized variants of uncertain significance in 70 of 467 additional probands (15.0%). These included known and novel variants in established oCCDD genes, genes associated with syndromes that sometimes include oCCDDs (eg, MYH10, KIF21B, TGFBR2, and TUBB6), genes that fit the syndromic component of the phenotype but had no prior oCCDD association (eg, CDK13, TGFB2), genes with no reported association with oCCDDs or the syndromic phenotypes (eg, TUBA4A, KIF5C, CTNNA1, KLB, and FGF21), and genes associated with oCCDD phenocopies that had resulted in misdiagnoses. ConclusionThis study suggests that unsolved oCCDDs are clinically and genetically heterogeneous disorders often overlapping other Mendelian conditions and nominates many candidates for future replication and functional studies.

Evolution of the regulation of developmental gene expression in blind Mexican cavefish.

Developmental evolution and diversification of morphology can arise through changes in the regulation of gene expression or protein-coding sequence. To unravel mechanisms underlying early developmental evolution in cavefish of the species Astyanax mexicanus, we compared transcriptomes of surface-dwelling and blind cave-adapted morphs at the end of gastrulation. Twenty percent of the transcriptome was differentially expressed. Allelic expression ratios in cave X surface hybrids showed that cis-regulatory changes are the quasi-exclusive contributors to inter-morph variations in gene expression. Among a list of 108 genes with change at the cis-regulatory level, we explored the control of expression of rx3, which is a master eye gene. We discovered that cellular rx3 levels are cis-regulated in a cell-autonomous manner, whereas rx3 domain size depends on non-autonomous Wnt and Bmp signalling. These results highlight how uncoupled mechanisms and regulatory modules control developmental gene expression and shape morphological changes. Finally, a transcriptome-wide search for fixed coding mutations and differential exon use suggested that variations in coding sequence have a minor contribution. Thus, during early embryogenesis, changes in gene expression regulation are the main drivers of cavefish developmental evolution.

Sequencing of Kaposi's Sarcoma Herpesvirus (KSHV) genomes from persons of diverse ethnicities and provenances with KSHV-associated diseases demonstrate multiple infections, novel polymorphisms, and low intra-host variance.

Recently published near full-length KSHV genomes from a Cameroon Kaposi sarcoma case-control study showed strong evidence of viral recombination and mixed infections, but no sequence variations associated with disease. Using the same methodology, an additional 102 KSHV genomes from 76 individuals with KSHV-associated diseases have been sequenced. Diagnoses comprise all KSHV-associated diseases (KAD): Kaposi sarcoma (KS), primary effusion lymphoma (PEL), KSHV-associated large cell lymphoma (KSHV-LCL), a type of multicentric Castleman disease (KSHV-MCD), and KSHV inflammatory cytokine syndrome (KICS). Participants originated from 22 different countries, providing the opportunity to obtain new near full-length sequences of a wide diversity of KSHV genomes. These include near full-length sequence of genomes with KSHV K1 subtypes A, B, C, and F as well as subtype E, for which no full sequence was previously available. High levels of recombination were observed. Fourteen individuals (18%) showed evidence of infection with multiple KSHV variants (from two to four unique genomes). Twenty-six comparisons of sequences, obtained from various sampling sites including PBMC, tissue biopsies, oral fluids, and effusions in the same participants, identified near complete genome conservation between different biological compartments. Polymorphisms were identified in coding and non-coding regions, including indels in the K3 and K15 genes and sequence inversions here reported for the first time. One such polymorphism in KSHV ORF46, specific to the KSHV K1 subtype E2, encoded a mutation in the leucine loop extension of the uracil DNA glycosylase that results in alteration of biochemical functions of this protein. This confirms that KSHV sequence variations can have functional consequences warranting further investigation. This study represents the largest and most diverse analysis of KSHV genome sequences to date among individuals with KAD and provides important new information on global KSHV genomics.

Approach for Phased Sequence-Based Genotyping of the Critical Pharmacogene Dihydropyrimidine Dehydrogenase (DPYD).

Pre-treatment genotyping of four well-characterized toxicity risk-variants in the dihydropyrimidine dehydrogenase gene (DPYD) has been widely implemented in Europe to prevent serious adverse effects in cancer patients treated with fluoropyrimidines. Current genotyping practices are largely limited to selected commonly studied variants and are unable to determine phasing when more than one variant allele is detected. Recent evidence indicates that common DPYD variants modulate the functional impact of deleterious variants in a phase-dependent manner, where a cis- or a trans-configuration translates into different toxicity risks and dosing recommendations. DPYD is a large gene with 23 exons spanning nearly a mega-base of DNA, making it a challenging candidate for full-gene sequencing in the diagnostic setting. Herein, we present a time- and cost-efficient long-read sequencing approach for capturing the complete coding region of DPYD. We demonstrate that this method can reliably produce phased genotypes, overcoming a major limitation with current methods. This method was validated using 21 subjects, including two cancer patients, each of whom carried multiple DPYD variants. Genotype assignments showed complete concordance with conventional approaches. Furthermore, we demonstrate that the method is robust to technical challenges inherent in long-range sequencing of PCR products, including reference alignment bias and PCR chimerism.