Molecular alterations in retinoblastoma beyond RB1.

Retinoblastoma is the most common malignant ocular tumor in children. Although RB1 alterations are most frequently involved in the etiology of retinoblastoma, candidate driver events and somatic alterations leading to cell transformation, tumor onset and progression remain poorly understood. In this study, we identified novel genomic alterations in tumors with a panel of 160 genes. Sanger sequencing and Multiplex Ligation-dependent Probe Amplification (MLPA) were initially performed for identifying patients without apparent RB1 alterations in blood DNA. Subsequently, NGS analyses of 24 paired (blood/tumor) samples of these patients were carried out for identifying somatic mutations and copy number variation in RB1 and other 159 genes. One novel pathogenic RB1 mutation and seven novel VUS were identified as well as 90 novel pathogenic mutations in 61 other genes. Twenty-three genes appeared exclusively mutated in tumors without altered RB1 alleles and three frequently affected biological pathways while five other tumors did not show pathogenic RB1 alterations or SNV/indels in 159 other genes. Curiously, deletion of GATA2, AKT1, ARID1A, DNMT3A, MAP2K2, MEN1, MTOR, PTCH1 and SUFU (in homo- or heterozygosity) were exclusively found in these tumors when compared to those with any pathogenic alterations, probably indicating genes that might be essential for the development of retinoblastoma regardless of a functional RB1. Identification of genes associated with retinoblastoma will contribute to understanding presently unknown aspects of this malignancy, which might be essential for its initiation and progression, as well as providing valuable molecular markers.

Large ancestral effective population size explains the difficult phylogenetic placement of owl monkeys.

The phylogenetic position of owl monkeys, grouped in the genus Aotus, has been a controversial issue for understanding Neotropical primate evolution. Explanations of the difficult phylogenetic assignment of owl monkeys have been elusive, frequently relying on insufficient data (stochastic error) or scenarios of rapid speciation (adaptive radiation) events. Using a coalescent-based approach, we explored the population-level mechanisms likely explaining these topological discrepancies. We examined the topological variance of 2,192 orthologous genes shared between representatives of the three major Cebidae lineages and the outgroup. By employing a methodological framework that allows for reticulated tree topologies, our analysis explicitly tested for non-dichotomous evolutionary processes impacting the finding of the position of owl monkeys in the cebid phylogeny. Our findings indicated that Aotus is a sister lineage of the callitrichines. Most gene trees (>50%) failed to recover the species tree topology, although the distribution of gene trees mismatching the true species topology followed the standard expectation of the multispecies coalescent without reticulation. We showed that the large effective population size of the common ancestor of Aotus and callitrichines was the most likely factor responsible for generating phylogenetic uncertainty. On the other hand, fast speciation scenarios or introgression played minor roles. We propose that the difficult phylogenetic placement of Aotus is explained by population-level processes associated with the large ancestral effective size. These results shed light on the biogeography of the early cebid diversification in the Miocene, highlighting the relevance of evaluating phylogenetic relationships employing population-aware approaches.

Assessment of the gorilla gut virome in association with natural simian immunodeficiency virus infection.

BACKGROUND: Simian immunodeficiency viruses (SIVs) of chimpanzees and gorillas from Central Africa crossed the species barrier at least four times giving rise to human immunodeficiency virus type 1 (HIV-1) groups M, N, O and P. The paradigm of non-pathogenic lentiviral infections has been challenged by observations of naturally infected chimpanzees with SIVcpz associated with a negative impact on their life span and reproduction, CD4+ T-lymphocyte loss and lymphoid tissue destruction. With the advent and dissemination of new generation sequencing technologies, novel promising markers of immune deficiency have been explored in human and nonhuman primate species, showing changes in the microbiome (dysbiosis) that might be associated with pathogenic conditions. The aim of the present study was to identify and compare enteric viromes of SIVgor-infected and uninfected gorillas using noninvasive sampling and ultradeep sequencing, and to assess the association of virome composition with potential SIVgor pathogenesis in their natural hosts. RESULTS: We analyzed both RNA and DNA virus libraries of 23 fecal samples from 11 SIVgor-infected (two samples from one animal) and 11 uninfected western lowland gorillas from Campo-Ma'an National Park (CP), in southwestern Cameroon. Three bacteriophage families (Siphoviridae, Myoviridae and Podoviridae) represented 67.5 and 68% of the total annotated reads in SIVgor-infected and uninfected individuals, respectively. Conversely, mammalian viral families, such as Herpesviridae and Reoviridae, previously associated with gut- and several mammalian diseases were significantly more abundant (p < 0.003) in the SIVgor-infected group. In the present study, we analyzed, for the first time, the enteric virome of gorillas and their association with SIVgor status. This also provided the first evidence of association of specific mammalian viral families and SIVgor in a putative dysbiosis context. CONCLUSIONS: Our results suggested that viromes might be potentially used as markers of lentiviral disease progression in wild gorilla populations. The diverse mammalian viral families, herein described in SIVgor-infected gorillas, may play a pivotal role in a disease progression still unclear in these animals but already well characterized in pathogenic lentiviral infections in other organisms. Larger sample sets should be further explored to reduce intrinsic sampling variation.

Changes in selection intensity on the mitogenome of subterranean and fossorial rodents respective to aboveground species.

Several rodent lineages independently acquired the ability to dig complex networks of tunnels where fossorial and subterranean species spend part or their whole life, respectively. Their underground lifestyles imposed harsh physiological demands, presumably triggering strong selective pressures on genes involved in energy metabolism like those coding for mitochondrial proteins. Moreover, underground lifestyles must have increased inbreeding and susceptibility to population bottlenecks as well as restricted migration, leading to small effective population size (Ne) that, in turn, must have reduced the effectiveness of selection. These stringent environmental conditions and small Ne might be still operating as antagonist factors of selection efficacy in these rodents. In this report, we tested, in a phylogenetic framework, how the intensity of selection on protein-coding mitochondrial genes (mt-genes) fluctuated along the evolution of fossorial and subterranean rodents respective to aboveground lineages. Our findings showed significant selection relaxation in most mt-genes of subterranean hystricomorphs (African mole-rats, tuco-tucos, and coruro), while only in three mt-genes of fossorial hystricomorphs (degus, red vizcacha rat, and fossorial spiny rats) selection efficacy was strongly reduced, probably due to demographic constraints. Conversely, selection intensification was found to have occurred in three mt-genes in fossorial sciuromorphs (ground squirrels, chipmunks, marmot, and allies). Our findings indicated that evolution of mitogenomes in fossorial and, mainly, in subterranean rodents was a complex output of a balance between intense ecological and physiological pressures, together with demographic constraints leading to genetic drift that, in turn, might have resulted in relaxed selection in hystricomorphs.

Investigation of major genetic alterations in neuroblastoma.

Neuroblastoma (NB) is the most common extracranial solid tumor in childhood. This malignancy shows a wide spectrum of clinical outcome and its prognosis is conditioned by manifold biological and genetic factors. We investigated the tumor genetic profile and clinical data of 29 patients with NB by multiplex ligation-dependent probe amplification (MLPA) to assess therapeutic risk. In 18 of these tumors, MYCN status was assessed by fluorescence in situ hybridization (FISH). Copy number variation was also determined for confirming MLPA findings in two 6p loci. We found 2p, 7q and 17q gains, and 1p and 11q losses as the most frequent chromosome alterations in this cohort. FISH confirmed all cases of MYCN amplification detected by MLPA. In view of unexpected 6p imbalance, copy number variation of two 6p loci was assessed for validating MLPA findings. Based on clinical data and genetic profiles, patients were stratified in pretreatment risk groups according to international consensus. MLPA proved to be effective for detecting multiple genetic alterations in all chromosome regions as requested by the International Neuroblastoma Risk Group (INRG) for therapeutic stratification. Moreover, this technique proved to be cost effective, reliable, only requiring standard PCR equipment, and attractive for routine analysis. However, the observed 6p imbalances made PKHD1 and DCDC2 inadequate for control loci. This must be considered when designing commercial MLPA kits for NB. Finally, four patients showed a normal MLPA profile, suggesting that NB might have a more complex genetic pattern than the one assessed by presently available MLPA kits.

The human retinoblastoma susceptibility gene (RB1): an evolutionary story in primates.

The tumor suppressor gene RB1 (Human Retinoblastoma Susceptibility Gene) plays a prominent role in normal development, gene transcription, DNA replication, repair, and mitosis. Its complete biallelic dysfunction in retinoblasts is the main cause of retinoblastoma in the human. Although this gene has been evolutionary conserved, comparisons between the reference and human RB1 coding region with its counterparts in 19 non-human primates showed 359 sites where nucleotide replacements took place during the radiation of these species. These resulted in missense substitutions in 97 codons, 91 of which by amino acids with radically different physicochemical properties. Several in frame deletions and two insertions were also observed in the N-terminal region of the pRB protein where the highest number of amino acid substitutions and radical amino changes were found. Fifty-six codons were inferred to be under negative selection and five under positive selection. Differences in codon usage showed evident phylogenetic signals, with hominids generally presenting higher indices of codon bias than other catarrhines. The lineage leading to platyrrhines and, within platyrrhines, the lineage leading to Saimiri boliviensis showed a high rate of nucleotide substitutions and amino acids. Finally, several RB1 alterations associated to retinoblastoma in the human were present in several non-human primates without an apparent pathological effect.

Bonafide, type-specific human papillomavirus persistence among HIV-positive pregnant women: predictive value for cytological abnormalities, a longitudinal cohort study.

This study investigated the rate of human papillomavirus (HPV) persistence, associated risk factors, and predictors of cytological alteration outcomes in a cohort of human immunodeficiency virus-infected pregnant women over an 18-month period. HPV was typed through L1 gene sequencing in cervical smears collected during gestation and at 12 months after delivery. Outcomes were defined as nonpersistence (clearance of the HPV in the 2nd sample), re-infection (detection of different types of HPV in the 2 samples), and type-specific HPV persistence (the same HPV type found in both samples). An unfavourable cytological outcome was considered when the second exam showed progression to squamous intraepithelial lesion or high squamous intraepithelial lesion. Ninety patients were studied. HPV DNA persistence occurred in 50% of the cases composed of type-specific persistence (30%) or re-infection (20%). A low CD4+T-cell count at entry was a risk factor for type-specific, re-infection, or HPV DNA persistence. The odds ratio (OR) was almost three times higher in the type-specific group when compared with the re-infection group (OR = 2.8; 95% confidence interval: 0.43-22.79). Our findings show that bonafide (type-specific) HPV persistence is a stronger predictor for the development of cytological abnormalities, highlighting the need for HPV typing as opposed to HPV DNA testing in the clinical setting.

MECP2, a gene associated with Rett syndrome in humans, shows conserved coding regions, independent Alu insertions, and a novel transcript across primate evolution.

BACKGROUND: The methyl-CpG Binding Protein two gene (MECP2) encodes a multifunctional protein comprising two isoforms involved in nuclear organization and regulation of splicing and mRNA template activity. This gene is normally expressed in all tissues, with a higher expression level in the brain during neuronal maturation. Loss of MECP2 function is the primary cause of Rett syndrome (RTT) in humans, a dominant, X-linked disorder dramatically affecting neural and motor development. RESULTS: We investigated the molecular evolution of MECP2 in several primate taxa including 36 species in 16 genera of neotropical (platyrrhine) primates. The coding region of the MECP2_e2 isoform showed a high level of evolutionary conservation among humans and other primates, with amino acid substitutions in 14 codons and one in-frame insertion of a single serine codon, between codons 357 and 358, in Ateles paniscus. Most substitutions occurred in noncritical regions of MECP2 and the majority of the algorithms used for analyzing selection did not provide evidence of positive selection. Conversely, we found 48 sites under negative selection in different regions, 23 of which were consistently found by three different algorithms. Similar to an inverted Alu insert found previously in a lesser ape at a parallel location, one Alu insertion of approximately 300 bp in Cebus and Sapajus was found in intron 3. Phylogenetic reconstruction of the intron 3 data provided a topology that was coincident with the consensus arrangement of the primate taxa. RNAseq data in the neotropical primate Callimico goeldii revealed a novel transcript consisting of a noncontinuous region of the human-homologous intron 2 in this species; this transcript accounted for two putative polypeptides. CONCLUSIONS: Despite the remarkable evolutionary conservation of MECP2, one in-frame codon insertion was observed in A. paniscus, and one region of intron 3 was affected by a trans-specific Alu retrotransposition in two neotropical primate genera. Moreover, identification of novel MECP2 transcripts in Callimico suggests that part of a homologous human intronic region might be expressed, and that the potential open reading frame in this region might be a subject of interest in RTT patients who carry an apparently normal MECP2 sequence.

Family-based genome-wide association study in Patagonia confirms the association of the DMD locus and cleft lip and palate.

The etiology of cleft lip with or without cleft palate (CL±P) is complex and heterogeneous, and multiple genetic and environmental factors are involved. Some candidate genes reported to be associated with oral clefts are located on the X chromosome. At least three genes causing X-linked syndromes [midline 1 (MID1), oral-facial-digital syndrome 1 (OFD1), and dystrophin (DMD)] were previously found to be associated with isolated CL±P. We attempted to confirm the role of X-linked genes in the etiology of isolated CL±P in a South American population through a family-based genome-wide scan. We studied 27 affected children and their mothers, from 26 families, in a Patagonian population with a high prevalence of CL±P. We conducted an exploratory analysis of the X chromosome to identify candidate regions associated with CL±P. Four genomic segments were identified, two of which showed a statistically significant association with CL±P. One is an 11-kb region of Xp21.1 containing the DMD gene, and the other is an intergenic region (8.7 kb; Xp11.4). Our results are consistent with recent data on the involvement of the DMD gene in the etiology of CL±P. The MID1 and OFD1 genes were not included in the four potential CL±P-associated X-chromosome genomic segments.

Rare nasal cleft in a patient with holoprosencephaly due to a mutation in the ZIC2 gene.

BACKGROUND: Holoprosencephaly (HPE) is a spectrum of midline malformations of the prosencephalon generally reflected in a continuum of midline facial anomalies. Patients with mutation in the ZIC2 gene usually present a normal or mildly dysmorphic face associated with a severe brain malformation. Here we present a rare unilateral nasal cleft (Tessier cleft n. 1) with holoprosencephaly in a patient with a ZIC2 mutation. CASE: The male newborn presented with alobar HPE, microcephaly, ocular hypertelorism, upslanting palpebral fissures, a bulky nose with a left paramedian alar cleft. Mutational screening for HPE genes revealed the occurrence of a frameshift mutation in the ZIC2 gene. The mutation was inherited from the father who presented only mild ocular hypotelorism but had an affected child with HPE from his first marriage. CONCLUSION: The occurrence of oral clefts is common in patients with HPE, but unusual in patients with mutation in the ZIC2 gene. To our knowledge, clefts of the nasal alae have been reported only once or twice in patients with ZIC2 mutations. In documented patients from the literature, only 2% of individuals with described pathogenic mutations in the ZIC2 gene (3/171) presented facial clefts, one of them a nasal cleft, while common oral clefts were observed in 27% of individuals (7/26) described with nonpathogenic ZIC2 mutations or presenting a concomitant mutation in another HPE gene. When compared with the general population, nasal clefts are common in ZIC2 mutations and these mutations must be searched for in undiagnosed cases.

A molecular phylogeny of living primates.

Comparative genomic analyses of primates offer considerable potential to define and understand the processes that mold, shape, and transform the human genome. However, primate taxonomy is both complex and controversial, with marginal unifying consensus of the evolutionary hierarchy of extant primate species. Here we provide new genomic sequence (~8 Mb) from 186 primates representing 61 (~90%) of the described genera, and we include outgroup species from Dermoptera, Scandentia, and Lagomorpha. The resultant phylogeny is exceptionally robust and illuminates events in primate evolution from ancient to recent, clarifying numerous taxonomic controversies and providing new data on human evolution. Ongoing speciation, reticulate evolution, ancient relic lineages, unequal rates of evolution, and disparate distributions of insertions/deletions among the reconstructed primate lineages are uncovered. Our resolution of the primate phylogeny provides an essential evolutionary framework with far-reaching applications including: human selection and adaptation, global emergence of zoonotic diseases, mammalian comparative genomics, primate taxonomy, and conservation of endangered species.

Molecular analysis of holoprosencephaly in South America.

Holoprosencephaly (HPE) is a spectrum of brain and facial malformations primarily reflecting genetic factors, such as chromosomal abnormalities and gene mutations. Here, we present a clinical and molecular analysis of 195 probands with HPE or microforms; approximately 72% of the patients were derived from the Latin American Collaborative Study of Congenital Malformations (ECLAMC), and 82% of the patients were newborns. Alobar HPE was the predominant brain defect in almost all facial defect categories, except for patients without oral cleft and median or lateral oral clefts. Ethmocephaly, cebocephaly, and premaxillary agenesis were primarily observed among female patients. Premaxillary agenesis occurred in six of the nine diabetic mothers. Recurrence of HPE or microform was approximately 19%. The frequency of microdeletions, detected using Multiplex Ligation-dependant Probe Amplification (MLPA) was 17% in patients with a normal karyotype. Cytogenetics or QF-PCR analyses revealed chromosomal anomalies in 27% of the probands. Mutational analyses in genes SHH, ZIC2, SIX3 and TGIF were performed in 119 patients, revealing eight mutations in SHH, two mutations in SIX3 and two mutations in ZIC2. Thus, a detailed clinical description of new HPE cases with identified genetic anomalies might establish genotypic and phenotypic correlations and contribute to the development of additional strategies for the analysis of new cases.

HIV/HPV co-infection during pregnancy in southeastern Brazil: prevalence, HPV types, cytological abnormalities and risk factors.

OBJECTIVE: HIV(+) pregnant women are at a higher risk of HPV infection and development of cervical cancer. Our objectives were to assess the prevalence and HPV types in HIV(+) pregnant women and to identify risk factors for HPV infection and cytological abnormalities. METHODS: Cervicovaginal smears were collected during pregnancy from 140 women. Partial HPV L1 gene and the exon 4 of the human TP53 gene (containing codon 72) were PCR-amplified and sequenced. Amplified products indicating multiple HPV infection were further cloned and sequenced. The association of demographic, obstetric and HIV-related clinical variables with HPV infection and cervical lesions was tested by univariate analyses, and significant factors were subsequently tested by logistic regression multivariate analysis. RESULTS: HPV DNA tested positive for 118 patients and HPV types were identified in 104 samples. Twenty-eight different types were found, HPV-16 and HPV-58 being the most prevalent. High-risk types were present in 79.8% of samples and multiple infections in 16.3%. Abnormal cervical smears were found in 44 patients (31.4%). Absolute CD4(+) T-cell counts below 350 were associated with HPV infection. Younger age was associated with cervical abnormalities and higher CD4(+) T-cell count was an apparent protective factor. CONCLUSIONS: We found a high prevalence of HPV infection and high-risk types in this cohort. Our results highlighted the relevance of immune system integrity rather than TP53 variants for protecting this highly vulnerable population to HPV infection and carcinogenesis.

Identification of bacterial infection in neotropical primates.

Emerging infectious diseases usually arise from wild animal populations. In the present work, we performed a screening for bacterial infection in natural populations of New World primates. The blood cell bulk DNAs from 181 individuals of four Platyrrhini genera were PCR screened for eubacterial 16S rRNA genes. Bacteria were detected and identified in 13 distinct individuals of Alouatta belzebul, Alouatta caraya, and Cebus apella monkeys from geographically distant regions in the states of Mato Grosso and Pará, Brazil. Sequence analyses showed that these Platyrrhini bacteria are closely related not only to human pathogens Pseudomonas spp. but also to Pseudomonas simiae and sheep-Acari infecting Pseudomonas spp. The identified Pseudomonas possibly represents a group of bacteria circulating in natural monkey populations.

Influence of MDM2 and MDM4 on development and survival in hereditary retinoblastoma.

BACKGROUND: Retinoblastoma (RB) accounts for 3% of all childhood malignancies, with different incidences around the world. This malignancy results from loss-of-function of both RB1 alleles although other genes, like MDM2 and MDM4, have been proposed to be involved in tumor development. PROCEDURE: We genotyped rs2279744T>G and rs937283A>G in MDM2, and rs4252668T>C and rs116197192G>A in MDM4, in 104 unrelated RB patients and 104 controls. Sixty-month survival Kaplan-Meier curves and χ(2)-tests were performed for estimating the putative effect of MDM2 and MDM4 alleles on disease progression and survival of RB patients. RESULTS: MDM2 rs2279744G was significantly more frequent in controls, indicating an apparently protective effect on RB development. However, survival of patients who carried a constitutional RB1 mutation was significantly lower with rs2279744TG or GG than with rs2279744TT. Presence of rs2279744G and a constitutional RB1 mutation was sixfold more frequent in the 0-12 month age group than other age groups at onset of symptoms (P = 0.0401). MDM4 rs4252668C was present at a significantly higher frequency in controls while the frequency of MDM4 rs116197192G was significantly higher in RB patients, suggesting that this allele might increase the risk of developing RB. CONCLUSION: Our results indicate that MDM2 and MDM4 polymorphisms may influence development and/or survival in RB.

Constitutional and somatic methylation status of DMRH19 and KvDMR in Wilms tumor patients.

The most frequent epigenetic alterations in Wilms tumor (WT) occur at WT2, assigned to 11p15. WT2 consists of two domains: telomeric domain 1 (DMRH19) that contains the IGF2 gene and an imprinted maternally expressed transcript (H19) and centromeric domain 2 (KvDMR) that contains the genes KCNQ1, KCNQ1OT1 and CDKN1C. In this work, we used pyrosequencing and MS-MLPA to compare the methylation patterns of DMRH19/KvDMR in blood and tumor samples from 40 WT patients. Normal constitutional KvDMR methylation indicated that most of the epigenetic alterations in WT occur at DMRH19. Constitutional DMRH19 hypermethylation (HM DMRH19) was observed in two patients with Beckwith-Wiedemann syndrome. Pyrosequencing and MS-MLPA showed HM DMRH19 in 28/34 tumor samples: 16/34 with isolated HM DMRH19 and 12/34 with concomitant HM DMRH19 and KvDMR hypomethylation, indicating paternal uniparental disomy. With the exception of one blood sample, the MS-MLPA and pyrosequencing findings were concordant. Diffuse or focal anaplasia was present in five tumor samples and was associated with isolated somatic HM DMRH19 in four of them. Constitutional 11p15 methylation abnormalities were present in 5% of the samples and somatic abnormalities in the majority of tumors. Combined analysis of DMRH19/KvDMR by pyrosequencing and MS-MLPA is beneficial for characterizing epigenetic anomalies in WT, and MS-MLPA is useful and reliable for estimation of DNA methylation in a clinical setting.

Identification, classification and evolution of owl monkeys (Aotus, Illiger 1811).

BACKGROUND: Owl monkeys, belonging to the genus Aotus, have been extensively used as animal models in biomedical research but few reports have focused on the taxonomy and phylogeography of this genus. Moreover, the morphological similarity of several Aotus species has led to frequent misidentifications, mainly at the boundaries of their distribution. In this study, sequence data from five mitochondrial regions and the nuclear, Y-linked, SRY gene were used for species identification and phylogenetic reconstructions using well characterized specimens of Aotus nancymaae, A. vociferans, A. lemurinus, A. griseimembra, A. trivirgatus, A. nigriceps, A. azarae boliviensis and A. infulatus. RESULTS: The complete MT-CO1, MT-TS1, MT-TD, MT-CO2, MT-CYB regions were sequenced in 18 Aotus specimens. ML and Bayesian topologies of concatenated data and separate regions allowed for the proposition of a tentative Aotus phylogeny, indicating that Aotus diverged some 4.62 Million years before present (MYBP). Similar analyses with included GenBank specimens were useful for assessing species identification of deposited data. CONCLUSIONS: Alternative phylogenetic reconstructions, when compared with karyotypic and biogeographic data, led to the proposition of evolutionary scenarios questioning the conventional diversification of this genus in monophyletic groups with grey and red necks. Moreover, genetic distance estimates and haplotypic differences were useful for species validations.

Mitochondrial divergence between 2 populations of the hooded capuchin, Cebus (Sapajus) cay (Platyrrhini, Primates).

We analyzed the molecular divergence of 2 separate populations of Cebus apella paraguayanus, recently considered a junior synonym of Cebus cay, and estimated its time of separation from C. apella. Cytochrome b DNA from 23 C. cay from Brazil and 9 from Paraguay showed 24 haplotypes (20 and 4, respectively), accounting for 29 variable sites (19 transitions and 10 transversions), with 40.0%, 26.7%, and 33.0% replacements at first, second, and third codon positions, respectively. Genetic distance between haplotypes averaged 0.5%, with 1.1% between C. cay populations. Phylogenetic reconstructions and median joining separated C. cay from Brazil and Paraguay. Neighbor joining showed C. cay and C. apella as sister groups, although C. cay and C. apella collapsed in maximum parsimony and maximum likelihood topologies. Analysis of molecular variance showed the highest variance component between C. cay populations, and mismatch distribution indicated that this species suffered a recent demographic expansion. Divergence time estimates suggested that the 2 populations of C. cay split in the Pleistocene, a period of repeated glaciation events leading to drastic changes in the vegetation composition of different biomes.

The characterization of two novel neotropical primate papillomaviruses supports the ancient within-species diversity model.

Papillomaviruses (PVs) are non-enveloped icosahedral viruses with a circular double-stranded DNA genome of ∼8,000 base pairs (bp). More than 200 different PV types have been identified to date in humans, which are distributed in five genera, with several strains associated with cancer development. Although widely distributed in vertebrates, Neotropical Primates (NP) PV infection was described for the first time only in 2016. Currently, four complete genomes of NP PVs have been characterized, three from Saimiri sciureus (SscPV1 to SscPV3) and one from Alouatta guariba (AgPV1). In this work, we describe two novel PV strains infecting Callithrix penicillata (provisionally named CpenPV1 and CpenPV2), using anal swab samples from animals residing at the Brasilia Primatology Center and next generation sequencing. The genomes of CpenPV1 (7,288 bp; 41.5% guanine-cytosine content - GC) and CpenPV2 (7,250 bp; 40.7% GC) contain the characteristic open reading frames (ORFs) for the early (E6, E7, E1, E2, and E4) and late (L2 and L1) PV genes. The L1 ORFs, commonly used for phylogenetic identification, share 76 per cent similarity with each other and differ 32 per cent from any other known PV, indicating that these new strains meet the criteria for defining novel species. PV genes phylogenetic variance was analyzed and different degrees of saturation revealed similar levels of topological heterogeneity, ruling out saturation as primary etiological factor for this phenomenon. Interestingly, the two CpenPV strains form a monophyletic clade within the Gammapapillomavirus genus (provisionally named gammapapillomavirus 32). Unlike for other NP PV strains, which grouped into a new sister genus of Alphapapillomavirus, this is the first report of NP PV strains grouping into a genus previously considered to exclusively comprise Old World Primates (OWP) PVs, including human PVs. These findings confirm the existence of a common ancestor for Gammapapillomavirus already infecting primates before the split of OWP and NP at ∼40 million years ago. Finally, our findings are consistent with an ancient within-species diversity model and emphasize the importance of increasing sampling to help understanding the PV-primate codivergence dynamics and pathogenic potential.