Evolution and diversity of nucleotide and dinucleotide composition in poxviruses.

Poxviruses (family Poxviridae) have long dsDNA genomes and infect a wide range of hosts, including insects, birds, reptiles and mammals. These viruses have substantial incidence, prevalence and disease burden in humans and in other animals. Nucleotide and dinucleotide composition, mostly CpG and TpA, have been largely studied in viral genomes because of their evolutionary and functional implications. We analysed here the nucleotide and dinucleotide composition, as well as codon usage bias, of a set of representative poxvirus genomes, with a very diverse host spectrum. After correcting for overall nucleotide composition, entomopoxviruses displayed low overall GC content, no enrichment in TpA and large variation in CpG enrichment, while chordopoxviruses showed large variation in nucleotide composition, no obvious depletion in CpG and a weak trend for TpA depletion in GC-rich genomes. Overall, intergenome variation in dinucleotide composition in poxviruses is largely accounted for by variation in overall genomic GC levels. Nonetheless, using vaccinia virus as a model, we found that genes expressed at the earliest times in infection are more CpG-depleted than genes expressed at later stages. This observation has parallels in betahepesviruses (also large dsDNA viruses) and suggests an antiviral role for the innate immune system (e.g. via the zinc-finger antiviral protein ZAP) in the early phases of poxvirus infection. We also analysed codon usage bias in poxviruses and we observed that it is mostly determined by genomic GC content, and that stratification after host taxonomy does not contribute to explaining codon usage bias diversity. By analysis of within-species diversity, we show that genomic GC content is the result of mutational biases. Poxvirus genomes that encode a DNA ligase are significantly AT-richer than those that do not, suggesting that DNA repair systems shape mutation biases. Our data shed light on the evolution of poxviruses and inform strategies for their genetic manipulation for therapeutic purposes.

Does exposure to different menstrual products affect the vaginal environment?

The vaginal ecosystem is a key component of women's health. It also represents an ideal system for ecologists to investigate the consequence of perturbations on species diversity and emerging properties between organizational levels. Here, we study how exposure to different types of menstrual products is linked to microbial, immunological, demographic, and behavioural measurements in a cohort of young adult women who reported using more often tampons (n = 107) or menstrual cups (n = 31). We first found that cup users were older and smoked less than tampon users. When analysing health indicators, we detected potential associations between cups use reporting and fungal genital infection. A multivariate analysis confirmed that in our cohort, reporting using cups over tampons was associated with the higher odds ratio to report a fungal genital infection diagnosis by a medical doctor within the last 3 months. We did not detect significant differences between groups in terms of their bacterial vaginal microbiota composition and found marginal differences in the level of expression of 20 cytokines. However, a multivariate analysis of these biological data identified some level of clustering based on the menstrual product type preferred (cups or tampons). These results suggest that exposure to different types of menstrual products could influence menstrual health. Larger studies and studies with a more powered setting are needed to assess the robustness of these associations and identify causal mechanisms.

Subfunctionalisation of paralogous genes and evolution of differential codon usage preferences: The showcase of polypyrimidine tract binding proteins.

Gene paralogs are copies of an ancestral gene that appear after gene or full genome duplication. When two sister gene copies are maintained in the genome, redundancy may release certain evolutionary pressures, allowing one of them to access novel functions. Here, we focused our study on gene paralogs on the evolutionary history of the three polypyrimidine tract binding protein genes (PTBP) and their concurrent evolution of differential codon usage preferences (CUPrefs) in vertebrate species. PTBP1-3 show high identity at the amino acid level (up to 80%) but display strongly different nucleotide composition, divergent CUPrefs and, in humans and in many other vertebrates, distinct tissue-specific expression levels. Our phylogenetic inference results show that the duplication events leading to the three extant PTBP1-3 lineages predate the basal diversification within vertebrates, and genomic context analysis illustrates that local synteny has been well preserved over time for the three paralogs. We identify a distinct evolutionary pattern towards GC3-enriching substitutions in PTBP1, concurrent with enrichment in frequently used codons and with a tissue-wide expression. In contrast, PTBP2s are enriched in AT-ending, rare codons, and display tissue-restricted expression. As a result of this substitution trend, CUPrefs sharply differ between mammalian PTBP1s and the rest of PTBPs. Genomic context analysis suggests that GC3-rich nucleotide composition in PTBP1s is driven by local substitution processes, while the evidence in this direction is thinner for PTBP2-3. An actual lack of co-variation between the observed GC composition of PTBP2-3 and that of the surrounding non-coding genomic environment would raise an interrogation on the origin of CUPrefs, warranting further research on a putative tissue-specific translational selection. Finally, we communicate an intriguing trend for the use of the UUG-Leu codon, which matches the trends of AT-ending codons. Our results are compatible with a scenario in which a combination of directional mutation-selection processes would have differentially shaped CUPrefs of PTBPs in vertebrates: the observed GC-enrichment of PTBP1 in placental mammals may be linked to genomic location and to the strong and broad tissue-expression, while AT-enrichment of PTBP2 and PTBP3 would be associated with rare CUPrefs and thus, possibly to specialized spatio-temporal expression. Our interpretation is coherent with a gene subfunctionalisation process by differential expression regulation associated with the evolution of specific CUPrefs.

Human DNA decays faster with time than viral dsDNA: an analysis on HPV16 using pathology archive samples spanning 85 years.

BACKGROUND: Quality of the nucleic acids extracted from Formalin Fixed Paraffin Embedded (FFPE) samples largely depends on pre-analytic, fixation and storage conditions. We assessed the differential sensitivity of viral and human double stranded DNA (dsDNA) to degradation with storage time. METHODS: We randomly selected forty-four HPV16-positive invasive cervical cancer (ICC) FFPE samples collected between 1930 and 1935 and between 2000 and 2004. We evaluated through qPCR the amplification within the same sample of two targets of the HPV16 L1 gene (69 bp, 134 bp) compared with two targets of the human tubulin-ß gene (65 bp, 149 bp). RESULTS: Both viral and human, short and long targets were amplified from all samples stored for 15 years. In samples archived for 85 years, we observed a significant decrease in the ability to amplify longer targets and this difference was larger in human than in viral DNA: longer fragments were nine times (CI 95% 2.6-35.2) less likely to be recovered from human DNA compared with 1.6 times (CI 95% 1.1-2.2) for viral DNA. CONCLUSIONS: We conclude that human and viral DNA show a differential decay kinetics in FFPE samples. The faster degradation of human DNA should be considered when assessing viral DNA prevalence in long stored samples, as HPV DNA detection remains a key biomarker of viral-associated transformation.

Genomic and phylogenetic characterization of ChPV2, a novel goat PV closely related to the Xi-PV1 species infecting bovines.

BACKGROUND: Papillomaviruses (PVs) infecting artiodactyls are very diverse, and only second in number to PVs infecting primates. PVs associated to lesions in economically important ruminant species have been isolated from cattle and sheep. METHODS: Potential PV DNA from teat lesions of a Damascus goat was isolated, cloned and sequenced. The PV genome was analyzed using bioinformatics approaches to detect open reading frames and to predict potential features of encoded proteins as well as putative regulatory elements. Sequence comparison and phylogenetic analyses using the concatenated E1E2L2L1 nucleotide and amino acid alignments was used to reveal the relationship of the new PV to the known PV diversity and its closest relevants. RESULTS: We isolated and characterized the full-genome of novel Capra hircus papillomavirus. We identified the E6, E7, E1, E2, L2, L1 open reading frames with protein coding potential and putative active elements in the ChPV2 proteins and putative regulatory genome elements. Sequence similarities of L1 and phylogenetic analyses using concatenated E1E2L2L1 nucleotide and amino acid alignments suggest the classification as a new PV type designated ChPV2 with a phylogenetic position within the XiPV genus, basal to the XiPV1 species. ChPV2 is not closely related to ChPV1, the other known goat PV isolated from healthy skin, although both of them belong confidently into a clade composed of PVs infecting cervids and bovids. Interestingly, ChPV2 contains an E6 open reading frame whereas all closely related PVs do not CONCLUSION: ChPV2 is a novel goat PV closely related to the Xi-PV1 species infecting bovines. Phylogenetic relationships and genome architecture of ChPV2 and closely related PV types suggest at least two independent E6 losses within the XiPV clade.

Distinct geographic clustering of oncogenic human papillomaviruses multiple infections in cervical cancers: Results from a worldwide cross-sectional study.

Coinfections by multiple Human Papillomaviruses (HPVs) are observed in approximately 6-8% of invasive cervical cancer (ICC) cases worldwide. But neither the presence of persistent HPVs coinfections nor their etiological role in the development of ICC is well understood. Cervical HPVs coinfections have been observed randomly, mostly in women with preneoplastic lesions, and only few studies have globally analyzed ICC cases. Here we explored the HPVs multiple infection patterns in a large worldwide sample of cross-sectional ICC cases. Paraffin-embedded ICC biopsy samples were tested using stringent HPV genotyping. Logistic regression models were used to identify the most likely pairwise HPV types in multiple infections. Multivariate analysis was applied to detect significant HPV coinfection patterns beyond pairwise HPVs comparison. Among 8780 HPV DNA-positive ICC cases worldwide, 6.7% (N = 587) contained multiple HPVs. Pairwise analysis revealed that HPV16|74, HPV31|33, HPV31|44, HPV33|44 and HPV45|70 pairs were significantly more frequently found together in multiple infections compared to any other HPV type combination, which supports the occasional role of Alpha-10 LR-HPVs in cervical cancers. In contrast, HPV16|31, HPV16|45, HPV16|51 and HPV18|HPV45 pairs were significantly less frequently found together than with any other HPV pair combination. Multivariate analysis sustained the results and revealed for the first time a distinct coinfection pattern in African ICCs stemming from the clustering of oncogenic HPV51/35/18/52 coinfections in African women. We suggest that the differential geographic HPVs coinfections clustering observed might be compatible with a specific modulation of the natural history/oncogenic potential of particular HPVs multiple infections and warrant monitoring for post-vaccinated.

Stably expressed APOBEC3H forms a barrier for cross-species transmission of simian immunodeficiency virus of chimpanzee to humans.

APOBEC3s (A3s) are potent restriction factors of human immunodeficiency virus type 1/simian immunodeficiency viruses (HIV-1/SIV), and can repress cross-species transmissions of lentiviruses. HIV-1 originated from a zoonotic infection of SIV of chimpanzee (SIVcpz) to humans. However, the impact of human A3s on the replication of SIVcpz remains unclear. By using novel SIVcpz reporter viruses, we identified that human APOBEC3B (A3B) and APOBEC3H (A3H) haplotype II strongly reduced the infectivity of SIVcpz, because both of them are resistant to SIVcpz Vifs. We further demonstrated that human A3H inhibited SIVcpz by deaminase dependent as well independent mechanisms. In addition, other stably expressed human A3H haplotypes and splice variants showed strong antiviral activity against SIVcpz. Moreover, most SIV and HIV lineage Vif proteins could degrade chimpanzee A3H, but no Vifs from SIVcpz and SIV of gorilla (SIVgor) lineages antagonized human A3H haplotype II. Expression of human A3H hapII in human T cells efficiently blocked the spreading replication of SIVcpz. The spreading replication of SIVcpz was also restricted by stable A3H in human PBMCs. Thus, we speculate that stably expressed human A3H protects humans against the cross-species transmission of SIVcpz and that SIVcpz spillover to humans may have started in individuals that harbor haplotypes of unstable A3H proteins.

Human papillomavirus DNA detected in fingertip, oral and bathroom samples from unvaccinated adolescent girls in Tanzania.

OBJECTIVE: Human papillomavirus (HPV) DNA has been detected in vaginal samples from adolescent girls who report no previous sex and, in high-income settings, from fingertips, raising the possibility of non-sexual transmission. No such studies originate from East Africa which bears among the highest cervical cancer incidence and HPV prevalence worldwide. HPV-related oral cancer incidence is increasing, but oral HPV prevalence data from East Africa are limited. We aimed to describe the HPV DNA prevalence in genital and non-genital sites and in the bathroom of unvaccinated adolescent girls, and examine genotype concordance between sites. METHODS: We nested a cross-sectional study of HPV in genital and extragenital sites within a cohort study of vaginal HPV acquisition. Unvaccinated girls age 16-18 years in Tanzania, who reported ever having had sex, were consented, enrolled and tested for the presence of HPV DNA in vaginal samples collected using self-administered swabs, oral samples collected using an oral rinse, and on fingertips and bathroom surfaces collected using a cytobrush. RESULTS: Overall, 65 girls were enrolled and 23 (35%, 95% CI 23% to 47%) had detectable vaginal HPV. Adequate (ß-globin positive) samples were collected from 36 girls' fingertips and HPV was detected in 7 (19%, 95% CI 6% to 33%). 63 girls provided adequate oral samples, 4 (6%, 95% CI 0% to 13%) of which had HPV DNA detected. In bathroom samples from 58 girls, 4 (7%, 95% CI 0% to 14%) had detectable HPV DNA. Of the 11 girls with extragenital HPV, six had the same genotype in >1 site. CONCLUSION: We found a high prevalence of HPV in non-genital sites in adolescent girls and in their bathrooms, in this region with a high cervical cancer incidence. Concordance of genotypes between sites supports the possibility of autoinoculation.

Transmission between Archaic and Modern Human Ancestors during the Evolution of the Oncogenic Human Papillomavirus 16.

Every human suffers through life a number of papillomaviruses (PVs) infections, most of them asymptomatic. A notable exception are persistent infections by Human papillomavirus 16 (HPV16), the most oncogenic infectious agent for humans and responsible for most infection-driven anogenital cancers. Oncogenic potential is not homogeneous among HPV16 lineages, and genetic variation within HPV16 exhibits some geographic structure. However, an in-depth analysis of the HPV16 evolutionary history was still wanting. We have analyzed extant HPV16 diversity and compared the evolutionary and phylogeographical patterns of humans and of HPV16. We show that codivergence with modern humans explains at most 30% of the present viral geographical distribution. The most explanatory scenario suggests that ancestral HPV16 already infected ancestral human populations and that viral lineages co-diverged with the hosts in parallel with the split between archaic Neanderthal-Denisovans and ancestral modern human populations, generating the ancestral HPV16A and HPV16BCD viral lineages, respectively. We propose that after out-of-Africa migration of modern human ancestors, sexual transmission between human populations introduced HPV16A into modern human ancestor populations. We hypothesize that differential coevolution of HPV16 lineages with different but closely related ancestral human populations and subsequent host-switch events in parallel with introgression of archaic alleles into the genomes of modern human ancestors may be largely responsible for the present-day differential prevalence and association with cancers for HPV16 variants.

Evolutionary analysis of Old World arenaviruses reveals a major adaptive contribution of the viral polymerase.

The Old World (OW) arenavirus complex includes several species of rodent-borne viruses, some of which (i.e., Lassa virus, LASV and Lymphocytic choriomeningitis virus, LCMV) cause human diseases. Most LCMV and LASV infections are caused by rodent-to-human transmissions. Thus, viral evolution is largely determined by events that occur in the wildlife reservoirs. We used a set of human- and rodent-derived viral sequences to investigate the evolutionary history underlying OW arenavirus speciation, as well as the more recent selective events that accompanied LASV spread in West Africa. We show that the viral RNA polymerase (L protein) was a major positive selection target in OW arenaviruses and during LASV out-of-Nigeria migration. No evidence of selection was observed for the glycoprotein, whereas positive selection acted on the nucleoprotein (NP) during LCMV speciation. Positively selected sites in L and NP are surrounded by highly conserved residues, and the bulk of the viral genome evolves under purifying selection. Several positively selected sites are likely to modulate viral replication/transcription. In both L and NP, structural features (solvent exposed surface area) are important determinants of site-wise evolutionary rate variation. By incorporating several rodent-derived sequences, we also performed an analysis of OW arenavirus codon adaptation to the human host. Results do not support a previously hypothesized role of codon adaptation in disease severity for non-Nigerian strains. In conclusion, L and NP represent the major selection targets and possible determinants of disease presentation; these results suggest that field surveys and experimental studies should primarily focus on these proteins.

Assessing parallel gene histories in viral genomes.

BACKGROUND: The increasing abundance of sequence data has exacerbated a long known problem: gene trees and species trees for the same terminal taxa are often incongruent. Indeed, genes within a genome have not all followed the same evolutionary path due to events such as incomplete lineage sorting, horizontal gene transfer, gene duplication and deletion, or recombination. Considering conflicts between gene trees as an obstacle, numerous methods have been developed to deal with these incongruences and to reconstruct consensus evolutionary histories of species despite the heterogeneity in the history of their genes. However, inconsistencies can also be seen as a source of information about the specific evolutionary processes that have shaped genomes. RESULTS: The goal of the approach here proposed is to exploit this conflicting information: we have compiled eleven variables describing phylogenetic relationships and evolutionary pressures and submitted them to dimensionality reduction techniques to identify genes with similar evolutionary histories. To illustrate the applicability of the method, we have chosen two viral datasets, namely papillomaviruses and Turnip mosaic virus (TuMV) isolates, largely dissimilar in genome, evolutionary distance and biology. Our method pinpoints viral genes with common evolutionary patterns. In the case of papillomaviruses, gene clusters match well our knowledge on viral biology and life cycle, illustrating the potential of our approach. For the less known TuMV, our results trigger new hypotheses about viral evolution and gene interaction. CONCLUSIONS: The approach here presented allows turning phylogenetic inconsistencies into evolutionary information, detecting gene assemblies with similar histories, and could be a powerful tool for comparative pathogenomics.

Determinants of FIV and HIV Vif sensitivity of feline APOBEC3 restriction factors.

BACKGROUND: Feline immunodeficiency virus (FIV) is a global pathogen of Felidae species and a model system for Human immunodeficiency virus (HIV)-induced AIDS. In felids such as the domestic cat (Felis catus), APOBEC3 (A3) genes encode for single-domain A3Z2s, A3Z3 and double-domain A3Z2Z3 anti-viral cytidine deaminases. The feline A3Z2Z3 is expressed following read-through transcription and alternative splicing, introducing a previously untranslated exon in frame, encoding a domain insertion called linker. Only A3Z3 and A3Z2Z3 inhibit Vif-deficient FIV. Feline A3s also are restriction factors for HIV and Simian immunodeficiency viruses (SIV). Surprisingly, HIV-2/SIV Vifs can counteract feline A3Z2Z3. RESULTS: To identify residues in feline A3s that Vifs need for interaction and degradation, chimeric human-feline A3s were tested. Here we describe the molecular direct interaction of feline A3s with Vif proteins from cat FIV and present the first structural A3 model locating these interaction regions. In the Z3 domain we have identified residues involved in binding of FIV Vif, and their mutation blocked Vif-induced A3Z3 degradation. We further identified additional essential residues for FIV Vif interaction in the A3Z2 domain, allowing the generation of FIV Vif resistant A3Z2Z3. Mutated feline A3s also showed resistance to the Vif of a lion-specific FIV, indicating an evolutionary conserved Vif-A3 binding. Comparative modelling of feline A3Z2Z3 suggests that the residues interacting with FIV Vif have, unlike Vif-interacting residues in human A3s, a unique location at the domain interface of Z2 and Z3 and that the linker forms a homeobox-like domain protruding of the Z2Z3 core. HIV-2/SIV Vifs efficiently degrade feline A3Z2Z3 by possible targeting the linker stretch connecting both Z-domains. CONCLUSIONS: Here we identified in feline A3s residues important for binding of FIV Vif and a unique protein domain insertion (linker). To understand Vif evolution, a structural model of the feline A3 was developed. Our results show that HIV Vif binds human A3s differently than FIV Vif feline A3s. The linker insertion is suggested to form a homeo-box domain, which is unique to A3s of cats and related species, and not found in human and mouse A3s. Together, these findings indicate a specific and different A3 evolution in cats and human.

Human papillomavirus DNA prevalence and type distribution in anal carcinomas worldwide.

Knowledge about human papillomaviruses (HPV) types involved in anal cancers in some world regions is scanty. Here, we describe the HPV DNA prevalence and type distribution in a series of invasive anal cancers and anal intraepithelial neoplasias (AIN) grades 2/3 from 24 countries. We analyzed 43 AIN 2/3 cases and 496 anal cancers diagnosed from 1986 to 2011. After histopathological evaluation of formalin-fixed paraffin-embedded samples, HPV DNA detection and genotyping was performed using SPF-10/DEIA/LiPA25 system (version 1). A subset of 116 cancers was further tested for p16(INK4a) expression, a cellular surrogate marker for HPV-associated transformation. Prevalence ratios were estimated using multivariate Poisson regression with robust variance in the anal cancer data set. HPV DNA was detected in 88.3% of anal cancers (95% confidence interval [CI]: 85.1-91.0%) and in 95.3% of AIN 2/3 (95% CI: 84.2-99.4%). Among cancers, the highest prevalence was observed in warty-basaloid subtype of squamous cell carcinomas, in younger patients and in North American geographical region. There were no statistically significant differences in prevalence by gender. HPV16 was the most frequent HPV type detected in both cancers (80.7%) and AIN 2/3 lesions (75.4%). HPV18 was the second most common type in invasive cancers (3.6%). p16(INK4a) overexpression was found in 95% of HPV DNA-positive anal cancers. In view of the results of HPV DNA and high proportion of p16(INK4a) overexpression, infection by HPV is most likely to be a necessary cause for anal cancers in both men and women. The large contribution of HPV16 reinforces the potential impact of HPV vaccines in the prevention of these lesions.

Time trends of human papillomavirus types in invasive cervical cancer, from 1940 to 2007.

Contribution over time of human papillomavirus (HPV) types in human cancers has been poorly documented. Such data is fundamental to measure current HPV vaccines impact in the years to come. We estimated the HPV type-specific distribution in a large international series of invasive cervical cancer (ICC) over 70 years prior to vaccination. Paraffin embedded ICC cases diagnosed between 1940 and 2007 were retrieved from eleven countries in Central-South America, Asia and Europe. Included countries reported to have low-medium cervical cancer screening uptake. Information on age at and year of diagnosis was collected from medical records. After histological confirmation, HPV DNA detection was performed by SPF-10/DEIA/LiPA25 (version1). Logistic regression models were used for estimating the adjusted relative contributions (RC) of HPV16 and of HPV18 over time. Among 4,771 HPV DNA positive ICC cases, HPV16 and HPV18 were the two most common HPVs in all the decades with no statistically significant variations of their adjusted-RC from 1940-59 to 2000-07 (HPV16-from 61.5 to 62.1%, and HPV18-from 6.9 to 7.2%). As well, the RC of other HPV types did not varied over time. In the stratified analysis by histology, HPV16 adjusted-RC significantly increased across decades in adenocarcinomas. Regarding age, cases associated to either HPV16, 18 or 45 were younger than those with other HPV types in all the evaluated decades. The observed stability on the HPV type distribution predicts a high and stable impact of HPV vaccination in reducing the cervical cancer burden in future vaccinated generations.

Two Human papillomavirus 11 complete genomes recovered from inverted sinonasal papillomas in humans.

We communicate here two complete Human papillomavirus 11 (HPV11) genomes recovered from one transitional and from one squamous inverted sinonasal papilloma, a rare proliferative disease in humans. Both genomes belong to the HPV11_A2 sublineage.

Methylation of human papillomavirus Type 16 CpG sites at E2-binding site 1 (E2BS1), E2BS2, and the Sp1-binding site in cervical cancer samples as determined by high-resolution melting analysis-PCR.

High-risk (HR) human papillomavirus (HPV)-associated carcinogenesis is driven mainly by the overexpression of E7 and E6 oncoproteins following viral DNA integration and the concomitant loss of the E2 open reading frame (ORF). However, the integration of HR-HPV DNA is not systematically observed in cervical cancers. The E2 protein acts as a transcription factor that governs viral oncogene expression. The methylation of CpGs in the E2-binding sites (E2BSs) in the viral long control region abrogates E2 binding, thus impairing the E2-mediated regulation of E7/E6 transcription. Here, high-resolution melting (HRM)-PCR was developed to quantitatively analyze the methylation statuses of E2BS1, E2BS2, and the specificity protein 1 (Sp1)-binding site in 119 HPV16-positive cervical smears. This is a rapid assay that is suitable for the analysis of cervical samples. The proportion of cancer samples with methylated E2BS1, E2BS2, and Sp1-binding site CpGs was 47%, whereas the vast majority of samples diagnosed as being within normal limits, low-grade squamous intraepithelial lesions (LSIL), or high-grade squamous intraepithelial lesions (HSIL) harbored unmethylated CpGs. Methylation levels varied widely, since some cancer samples harbored up to 60% of methylated HPV16 genomes. A pyrosequencing approach was used as a confirmation test and highlighted that quantitative measurement of methylation can be achieved by HRM-PCR. Its prognostic value deserves to be investigated alone or in association with other biomarkers. The reliability of this single-tube assay offers great opportunities for the investigation of HPV16 methylation in other HPV-related cancers, such as head and neck cancers, which are a major public health burden.

Transcriptomic, proteomic, and functional consequences of codon usage bias in human cells during heterologous gene expression.

Differences in codon frequency between genomes, genes, or positions along a gene, modulate transcription and translation efficiency, leading to phenotypic and functional differences. Here, we present a multiscale analysis of the effects of synonymous codon recoding during heterologous gene expression in human cells, quantifying the phenotypic consequences of codon usage bias at different molecular and cellular levels, with an emphasis on translation elongation. Six synonymous versions of an antibiotic resistance gene were generated, fused to a fluorescent reporter, and independently expressed in HEK293 cells. Multiscale phenotype was analyzed by means of quantitative transcriptome and proteome assessment, as proxies for gene expression; cellular fluorescence, as a proxy for single-cell level expression; and real-time cell proliferation in absence or presence of antibiotic, as a proxy for the cell fitness. We show that differences in codon usage bias strongly impact the molecular and cellular phenotype: (i) they result in large differences in mRNA levels and protein levels, leading to differences of over 15 times in translation efficiency; (ii) they introduce unpredicted splicing events; (iii) they lead to reproducible phenotypic heterogeneity; and (iv) they lead to a trade-off between the benefit of antibiotic resistance and the burden of heterologous expression. In human cells in culture, codon usage bias modulates gene expression by modifying mRNA availability and suitability for translation, leading to differences in protein levels and eventually eliciting functional phenotypic changes.

Digital Papillary Adenocarcinoma in Nonacral Skin: Clinicopathologic and Genetic Characterization of 5 Cases.

Digital papillary adenocarcinoma (DPA) is a rare sweat gland neoplasm that has exceptionally been reported outside acral locations. Recently, human papillomavirus 42 was identified as the main oncogenic driver of DPA. Herein, we report 5 tumors arising in extra-acral locations predominantly in the female anogenital skin. Four patients were female and 1 patient was male. The mean age at the diagnosis time was 65 years (range: 55 to 82 y). Tumors were located on the vulva (n=3), perianal area (n=1), and forearm (n=1). Histologically, all tumors were lobular and mainly solid and composed of sheets of cells with rare focal papillae and frequent glandular structures in a "back-to-back" pattern and lined by atypical basophilic cells. Immunohistochemistry showed diffuse positivity for SOX10. Epithelial membrane antigen and carcinoembryonic antigen highlighted the luminal cells and staining for p63 and p40 revealed a consistent and continuous myoepithelial component around glandular structures. Follow-up was available in 3 cases (mean duration: 12 mo [range: 8 to 16 mo]). One patient developed local recurrence and 1 experienced regional lymph node metastases. HPV Capture Next-generation sequencing revealed the presence of the HPV42 genome in all samples. Viral reads distributions were compatible in the 5 cases with an episomal nature of the viral genome, with a recurrent deletion in the E1 and/or E2 open reading frames. In conclusion, this study demonstrates that digital DPA may rarely present in nonacral locations mainly in the female anogenital area, usually with a more solid pattern as compared with those cases presenting on the digits and it is also associated with HPV42.

Sweat Gland Tumors Arising on Acral Sites: A Molecular Survey.

Recurrent oncogenic drivers have been identified in a variety of sweat gland tumors. Recently, integration of human papillomavirus type 42 (HPV42) has been reported in digital papillary adenocarcinoma (DPA). The main objectives of the present study were (i) to provide an overview of the prevalence of previously identified oncogenic drivers in acral sweat gland tumors and (ii) to genetically characterize tumors in which no recurrent genetic alteration has been identified yet. Cases of acral sweat gland tumors were identified from the database of the French network CARADERM. After histologic review, the presence of previously identified genetic alterations was investigated in the entire cohort (n=79) using a combination of immunohistochemistry and targeted DNA and RNA sequencing. Tumor entities with no recurrent genetic alterations were submitted to whole-transcriptome sequencing. CRTC1::MAML2 fusion was identified in cases of hidradenoma and hidradenocarcinoma (n=9/12 and n=9/12). A p.V600E mutation of BRAF was observed in all cases of tubular adenoma (n=4). YAP1:MAML2 and YAP1::NUTM1 fusions were observed in poroid tumors (n=15/25). ETV6::NTRK3 and TRPS1::PLAG1 fusion transcripts were identified in secretory carcinoma (n=1/1) and cutaneous mixed tumors (n=3/4), respectively. The HPV42 genome was detected in most cases of DPA (n=10/11) and in 1 adnexal adenocarcinoma not otherwise specified. Finally, whole-transcriptome analysis revealed BRD3::NUTM1 or NSD3::NUTM1 fusions in 2 cases of NUT adnexal carcinoma and NCOA4::RET and CCDC6::RET fusion transcripts in 2 cystadenoma/hidrocystoma-like tumors. Our study confirms distinctive cytogenetic abnormalities in a wide number of acral adnexal neoplasms and supports the use of molecular analysis as a valuable aid in the diagnosis of these rare and often difficult to diagnose group of neoplasms.

An ancient history of gene duplications, fusions and losses in the evolution of APOBEC3 mutators in mammals.

BACKGROUND: The APOBEC3 (A3) genes play a key role in innate antiviral defense in mammals by introducing directed mutations in the DNA. The human genome encodes for seven A3 genes, with multiple splice alternatives. Different A3 proteins display different substrate specificity, but the very basic question on how discerning self from non-self still remains unresolved. Further, the expression of A3 activity/ies shapes the way both viral and host genomes evolve. RESULTS: We present here a detailed temporal analysis of the origin and expansion of the A3 repertoire in mammals. Our data support an evolutionary scenario where the genome of the mammalian ancestor encoded for at least one ancestral A3 gene, and where the genome of the ancestor of placental mammals (and possibly of the ancestor of all mammals) already encoded for an A3Z1-A3Z2-A3Z3 arrangement. Duplication events of the A3 genes have occurred independently in different lineages: humans, cats and horses. In all of them, gene duplication has resulted in changes in enzyme activity and/or substrate specificity, in a paradigmatic example of convergent adaptive evolution at the genomic level. Finally, our results show that evolutionary rates for the three A3Z1, A3Z2 and A3Z3 motifs have significantly decreased in the last 100 Mya. The analysis constitutes a textbook example of the evolution of a gene locus by duplication and sub/neofunctionalization in the context of virus-host arms race. CONCLUSIONS: Our results provide a time framework for identifying ancestral and derived genomic arrangements in the APOBEC loci, and to date the expansion of this gene family for different lineages through time, as a response to changes in viral/retroviral/retrotransposon pressure.