Genetic analysis of six communities of Mbyá-Guaraní inhabiting northeastern Argentina by means of nuclear and mitochondrial polymorphic markers

Autosomal STRs, Y-chromosome markers, and mitochondrial DNA sequences were investigated in six Mbyá-Guaraní villages (Fortín M'Bororé, Yryapu, Tabay, Kaaguy Poty, Jejy, and Yaboti), all of them settled within the province of Misiones, northeastern Argentina. One hundred twenty-one unrelated individuals were analyzed. The study involved typing fifteen autosomal STRs, nine Y-chromosome STRs, and four biallele loci in the nonrecombinant region of the Y chromosome, sequencing the mtDNA of hypervariable regions I and II, and detecting the 9-bp ins/del in region V of mtDNA. All autosomal STRs were in Hardy-Weinberg equilibrium. The four major native American mtDNA haplogroups were represented in the sample. Haplogroups A2 and D1 exhibited the highest frequencies (40.5% and 36.0%, respectively), and haplogroups B2 and C1 appeared to be less frequent (17.5% and 6.0%, respectively). The native American haplogroup Q1a3a was observed in a relevant proportion (88.8%). In addition, a nine-STR Y-chromosome haplo-type (DYS19*13, DYS389I*14, DYS389II*31, DYS390*24, DYS391*11, DYS392*14, DYS393*11, DYS385A*14, DYS385B*16) exhibited a frequency of more than 36%. Our results indicate that the analyzed Argentinean Guaraní individuals are genetically more closely related to Guaraní from Brazil [genetic distance (Δμ)(2) = 0.48] than to other related tribes that are geographically closer. Statistical approaches based on autosomal data do not support the hypothesis of genetic drift previously proposed; however, this apparent discrepancy might be due to the lack of sensitivity of the autosomal markers used here.     

Divergent Human Y-Chromosome Microsatellite Evolution Rates

In this work, we analyze several characteristics influencing the low variability of the microsatellite DYS19 in the major founder Amerindian Y chromosome lineage containing the point mutation DYS199-T. Variation of DYS19 was compared with that of five other Y-linked tetranucleotide repeat loci (DYS389A, DYS389B, DYS390, DYS391, and DYS393) in the DYS199-T lineage. All the other microsatellites showed significantly higher levels of variability than DYS19 as measured by gene diversity and repeat number variance. Moreover, we had previously shown that DYS19 had high diversity in Brazilians and in several other populations worldwide. Thus, the slow DYS19 evolution in the DYS199-T lineage seems to be both locus and allele specific. To understand the slow DYS19 evolutionary rate, the microsatellite loci were compared according to their mapping on the Y chromosome and also on the basis of structural aspects such as the base composition of the repeat motif and flanking regions and the degree of perfection and size (repeat number) of the variable blocks. The only observed difference that might be related to the low DYS19 variability is its small average number of repeats, a value expected to be closer to the founder DYS19 allele in the DYS199-T lineage. These data were also compared to other derived Y lineages. The Tat-C lineage displayed a lower DYS19 variability correlated to a small average repeat number, while in the DYS234-G lineage, a high DYS19 variability was found associated to a larger average repeat number. This approach reveals that evolution of Y microsatellites in lineages defined by slowly evolving markers, such as point mutations, can be greatly influenced by the size (number of repeats of the variable block) of the founder allele in each microsatellite locus. Thus lineage-dating methods using microsatellite variation should be practiced with great care. Received: 7 November 1998 / Accepted: 9 April 1999     

Genetic diversity of Y-chromosome microsatellites in the Fujian Han and the Sichuan Han populations of China

Y-chromosome short tandem repeats (STRs) are potentially useful for forensic, anthropological and evolutionary studies. In this study we chose the loci DYS 19, DYS 388, DYS 389 I, DYS 389 II, DYS 390, DYS 391, DYS 392, DYS 393, DYS 425 and DYS 426. Blood samples were taken from 46 unrelated male individuals from Fujian Han and 43 unrelated males from Sichuan Han in China. DNA was extracted by conventional chelex extraction procedure. PCR was carried out in two multiplex reactions. Fragment analysis was conducted on an ABI PRISM 310 Genetic Analyzer. Allele frequency distributions and discrimination indices were calculated, and the two populations were tested for genetic differences by means of analysis of molecular variance (AMOVA). Here we obtained 75 Y-STR haplotypes and the haplotype diversity for the complete haplotype was 0.9884 in Fujian Han and 0.9967 in Sichuan Han. A larger genetic difference became apparent between the two populations that belong to the Sino-Tibetan speaking populations.     

Genetic diversity of Y-specific STRs in chimpanzees (Pan troglodytes)

Using the primers described for humans, sequences for 11 Y-specific microsatellites (DYS434, DYS435, DYS436, DYS437, DYS438, DYS439, GATA A10, A7.1, A7.2, C4, and H4 [Gusm?o et al., in press]), previously described in 10 male chimpanzees (Pan troglodytes), were confirmed in nine additional male chimpanzees. Sequences for nine additional microsatellites (DYS19, DYS385I and II, DYS389I and II, DYS390, DYS391, DYS392, and DYS393) were determined in all 19 male chimpanzees; homology to human Y-Short Tandem Repeat (STRs) was confirmed by sequencing. Good amplification results were not obtained for DYS19 and DYS385I/II. Two amplicons were obtained for DYS389I/II, but in contrast to humans, the larger fragment was not Y-specific. Moreover, no polymorphism was observed for DYS434, DYS435, or GATA A10. Consequently, these eight STRs were eliminated from further analyses, and haplotype and allele frequencies were estimated for the remaining 12 STRs. A high haplotype diversity value was found (1.000 +/- 0.017), demonstrating the usefulness and informative power of these Y-STRs for future studies on chimpanzee population genetics.     

Characterization of a Small Family (CAIII) of Microsatellite-Containing Sequences with X-Y Homology

Four X-linked loci showing homology with a previously described Y-linked polymorphic locus (DYS413) were identified and characterized. By fluorescent in situ hybridization (FISH), somatic cell hybrids, and YAC screening, the X-linked members of this small family of sequences (CAIII) all map in Xp22, while the Y members map in Yq11. These loci contribute to the overall similarity of the two genomic regions. All of the CAIII loci contain an internal microsatellite of the (CA)_n type. The microsatellites display extensive length polymorphism in two of the X-linked members as well as in the Y members. In addition, common sequence variants are found in the portions flanking the microsatellites in two of the X-linked members. Our results indicate that, during the evolution of this family, length variation on the Y chromosome was accumulated at a rate not slower than that on the X chromosome. Finally, these sequences represent a model system with which to analyze human populations for similar X- and Y-linked polymorphisms. Received: 29 July 1996 / Accepted: 15 January 1997     

Y-chromosome polymorphisms and the origins of the European gene pool

Gradients of allele frequencies have long been considered the main genetic characteristic of the European population, but mitochondrial DNA diversity seems to be distributed differently. One Alu insertion (YAP), five tetranucleotide (DYS19, DYS389B, DYS390, DYS391 and DYS393) and one trinucleotide (DYS392) microsatellite loci of the Y chromosome were analysed for geographical patterns in 59 European populations. Spatial correlograms showed clines for most markers, which paralleled the gradients previously observed for two RFLP polymorphisms. Effective separation times between populations were estimated from genetic distances at microsatellite loci. Even after correcting for the possible effects of continuous local gene flow, the most distant Indo-European-speaking populations seem to have separated no more than 7000 years ago. The clinal patterns and the estimated, recent separation times between populations jointly suggest that Y-chromosome diversity in Europe largely reflects a directional demic expansion, which is unlikely to have occurred before the Neolithic period.     

A Concertedly Evolving Region in Chironomus, Unique Within the Telomere

Chromosome terminal, complex repeats in the dipteran Chironomus pallidivittatus show rapid concerted evolution during which there is remarkably efficient homogenization of the repeat units within and between chromosome ends. It has been shown previously that gene conversion is likely to be an important component during these changes. The sequence evolution could be a result of different processes—exchanges between repeats in the tandem array as well as information transfer between units in different chromosomes—and is therefore difficult to analyze in detail. In this study the concerted evolution of a region present only once per chromosome, at the junction between the telomeric complex repeats and the subtelomeric DNA was therefore investigated in the two sibling species C. pallidivittatus and C. tentans. Material from individual microdissected chromosome ends was used, as well as clones from bulk genomic DNA. On the telomeric side of the border pronounced species-specific sequence differences were observed, the patterns being similar for clones of different origin within each species. Mutations had been transmitted efficiently between chromosomes also when adjoining, more distally localized DNA showed great differences in sequence, suggesting that gene conversion had taken place. The evolving telomeric region bordered proximally to subtelomeric DNA with high evolutionary constancy. More proximally localized, subtelomeric DNA evolved more rapidly and showed heterogeneity between species and chromosomes. Received: 24 September 1997 / Accepted: 24 November 1997     

Male ancestry structure and interethnic admixture in African-descent communities from the Amazon as revealed by Y-chromosome Strs

Some genetic markers on both the Y chromosome and mtDNA are highly polymorphic and population‐specific in humans, representing useful tools for reconstructing the past history of populations with poor historical records. Such lack of information is usually true in the case of recent African‐descent populations of the New World founded by fugitive slaves throughout the slavery period in the Americas, particularly in Brazil, where those communities are known as quilombos. Aiming to recover male‐derived ethnic structure of nine quilombos from the Brazilian Amazon, a total of 300 individuals, belonging to Mazagão Velho (N = 24), Curiaú (N = 48), Mazagão (N = 36), Trombetas (N = 20), Itacoã (N = 22), Saracura (N = 46), Marajó (N = 58), Pitimandeua (N = 26), and Pontal (N = 20), were investigated for nine Y‐STRs (DYS393, DYS19, DYS390, DYS389 I, DYS389 II, DYS392, DYS391, DYS385 I/II). From the 169 distinct haplotypes obtained, 120 were singletons. The results suggest the West African coast as the main origin of slaves brought to Brazil (54% of male contribution); the European contribution was high (41%), while the Amerindian's was low (5%). Those results contrast with previous mtDNA data that showed high Amerindian female contribution (46.6%) in African‐descent populations. AMOVA suggests that the genetic differentiation among the quilombos is mainly influenced by admixture with European. However, when restricting AMOVA to African‐specific haplotypes, low differentiation was detected, suggesting great genetic homogeneity of the African founding populations and/or a later homogenization by intense slave trade inside Brazil. Am J Phys Anthropol, 2011. © 2010 Wiley‐Liss, Inc.     

A study of Y-chromosome microsatellite variation in sub-Saharan Africa: a comparison between F(ST) and R(ST) genetic distances

Seven Y-chromosome microsatellite loci (DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, and DYS393) were analyzed in three populations from sub-Saharan Africa: the Bamileke and Ewondo populations from Cameroon and the Hutu from Rwanda. Complete typing was obtained for 112 individuals, and a total of 53 different haplotypes was observed. The single-locus gene diversity, averaged across populations, ranges from 0.100 for the DYS392 locus to 0.610 for the DYS389I locus. The haplotype diversity ranges from 0.832 (Ewondo) to 0.965 (Hutu), with an intermediate value of 0.918 in the Bamileke. The diversity among Bamileke, Ewondo, Hutu, and other sub-Saharan populations selected from the literature was analyzed using both a classical (F(ST)) and a stepwise-based (R(ST)) genetic distance method. The pattern of interpopulational diversity based on F(ST) was congruent with anthropological knowledge, while that based on R(ST) revealed unexpected and unconvincing population affinities. From a practical point of view, our study indicates that Y-chromosome microsatellite data may provide useful information for analyses of interpopulational diversity among sub-Saharan populations if an adequate number of loci and individuals along with an appropriate genetic distance method are used. On a theoretical ground, we propose that the lesser performance of R(ST) compared to F(ST) could be explained by the important role played by genetic drift in shaping the relationships among examined populations.     

Population data for 17 short tandem repeat loci on Y chromosome in northern Croatia

Human Y-short tandem repeats (STRs) are tandem repeat arrays of two to seven base pair units on non-recombining region (NRY) of the human Y chromosome. Studies on Y-STR are interesting in both population genetics and forensics. The aim of this study was to investigate the population genetic properties of 17 STR loci on Y chromosome in the northern Croatia region. We carried out a statistical analysis of the data from previously performed genetic analysis collected during routine forensic work by the Forensic Science Centre “Ivan Vu?eti?”. A total of 220 unrelated healthy men from northern Croatia were selected for the purpose of this study. Genomic DNA was extracted using Chelex procedure from FTA^® cards. Y-chromosomal STRs were determined using the AmpFISTR Yfiler PCR amplification kit. The haplotype frequencies were determined by direct counting and analyzed using Arlequin 3.1 and analysis of molecular variance calculated with the Y chromosome haplotype reference database online analysis tool. A total of 210 haplotypes were identified, 200 of which were unique. Total haplotype diversity was 0.995. Locus diversity varied from 0.331 for DYS392 to 0.783 for DYS385 locus. Allele frequencies diversity was 0.662. Discrimination capacity was 95.7%. The use of European minimal haplotype set indicated the most resemblance of this population to the Croatian capital of Zagreb, with modest resemblance to Bosnia and Herzegovina, Serbia and Hungary. This article provides the first overview of the Y chromosome STR variability in northern Croatia, thus providing the referent point for any future forensic and genetic epidemiology efforts in this region.     

A short tandem repeat-based phylogeny for the human Y chromosome

Human Y-chromosomal short tandem repeat (STR) data provide a potential model system for the understanding of autosomal STR mutations in humans and other species. Yet, the reconstruction of STR evolution is rarely attempted, because of the absence of an appropriate methodology. We here develop and validate a phylogenetic-network approach. We have typed 256 Y chromosomes of indigenous descent from Africa, Asia, Europe, Australia, and highland Papua New Guinea, for the STR loci DYS19, DXYS156Y, DYS389, DYS390, DYS392, and DYS393, as well as for five ancient biallelic mutation events: two poly (A) length variants associated with the YAP insertion, two independent SRY-1532 mutations, and the 92R7 mutation. We have used our previously published pedigree data from 11,000 paternity-tested autosomal STR-allele transfers to produce a two-class weighting system for the Y-STR loci that is based on locus lengths and motif lengths. Reduced-median-network analysis yields a phylogeny that is independently supported by the five biallelic mutations, with an error of 6%. We find the earliest branch in our African San (Bushmen) sample. Assuming an age of 20,000 years for the Native American DYS199 T mutation, we estimate a mutation rate of 2.6x10-4 mutations/20 years for slowly mutating Y STRs, approximately 10-fold slower than the published average pedigree rate.     

DYS19 and DYS199 loci in a Chilean population of mixed ancestry

The current Chilean population originated from admixture between aboriginal populations (Amerindians) and Spanish conquerors of European origin. Consequently, the unions that gave rise to the Chilean population were chiefly between Spanish males and aboriginal females, and not the converse. To test the hypothesis that the Y chromosome of the Chilean population is mainly of Spanish origin, while the other chromosomes are from mixed (European and aboriginal) origin, we studied the DYS19 and DYS199 loci in two samples. One sample was obtained from a high socioeconomic stratum, while a second sample was from a low stratum. We studied male blood donors (N = 187) from Santiago, the capital of the country. Subjects were typed for the autosomal ABO and Rh (locus D) blood groups, and for the Y-linked DYS19 and the DYS199 loci, reported as Y-chromosome haplotypes. The aboriginal admixture was estimated for each genetic marker. The percentage of aboriginal admixture was 38.17% for the ABO system and 31.28% for the Rh system in the low socioeconomic stratum and 19.22% and 22.5%, respectively, in the high stratum. Y-chromosome haplotype frequencies constructed from the DYS19 and DYS199 loci demonstrated that the main haplotypes were DYS19*14/DYS199 C, as is often the case with many European populations, and DYS19*13/DYS199 C. The aboriginal admixture from Y-haplotype frequencies was estimated to be 15.83% in the low socioeconomic stratum and 6.91% in the high stratum. These values are lower than the values found using autosomal genetic markers, and are consistent with the historical background of the population studied. This study highlights the population genetic consequences of the asymmetric pattern of genome admixture between two ancestral populations (European and Amerindian).     

The distribution of human genetic diversity: a comparison of mitochondrial, autosomal, and Y-chromosome data

We report a comparison of worldwide genetic variation among 255 individuals by using autosomal, mitochondrial, and Y-chromosome polymorphisms. Variation is assessed by use of 30 autosomal restriction-site polymorphisms (RSPs), 60 autosomal short-tandem-repeat polymorphisms (STRPs), 13 Alu-insertion polymorphisms and one LINE-1 element, 611 bp of mitochondrial control-region sequence, and 10 Y-chromosome polymorphisms. Analysis of these data reveals substantial congruity among this diverse array of genetic systems. With the exception of the autosomal RSPs, in which an ascertainment bias exists, all systems show greater gene diversity in Africans than in either Europeans or Asians. Africans also have the largest total number of alleles, as well as the largest number of unique alleles, for most systems. GST values are 11%-18% for the autosomal systems and are two to three times higher for the mtDNA sequence and Y-chromosome RSPs. This difference is expected because of the lower effective population size of mtDNA and Y chromosomes. A lower value is seen for Y-chromosome STRs, reflecting a relative lack of continental population structure, as a result of rapid mutation and genetic drift. Africa has higher GST values than does either Europe or Asia for all systems except the Y-chromosome STRs and Alus. All systems except the Y-chromosome STRs show less variation between populations within continents than between continents. These results are reassuring in their consistency and offer broad support for an African origin of modern human populations.     

Genetic variation of the Y chromosome in Chibcha-speaking Amerindians of Costa Rica and Panama

Genetic variation of the Y chromosome in five Chibchan tribes (Bribri, Cabecar, Guaymi, Huetar, and Teribe) of Costa Rica and Panama was analyzed using six microsatellite loci (DYS19, DYS389A, DYS389B, DYS390, DYS391, and DYS393), the Y-chromosome-specific alphoid system (alphah), the Y-chromosome Alu polymorphism (YAP), and a specific pre-Columbian transition (C-->T) (M3 marker) in the DYS 199 locus that defines the Q-M3 haplogroup. Thirty-nine haplotypes were found, resulting in a haplotype diversity of 0.937. The Huetar were the most diverse tribe, probably because of their high levels of interethnic admixture. A candidate founder Y-chromosome haplotype was identified (15.1% of Chibchan chromosomes), with the following constitution: YAP-, DYS199*T, alphah-II, DYS19*13, DYS389A*17, DYS389B*10, DYS390*24, DYS391*10, and DYS393*13. This haplotype is the same as the one described previously as one of the most frequent founder paternal lineages in native American populations. Analysis of molecular variance indicated that the between-population variation was smaller than the within-population variation, and the comparison with mtDNA restriction data showed no evidence of differential structuring between maternally and paternally inherited genes in the Chibchan populations. The mismatch-distribution approach indicated estimated coalescence times of the Y chromosomes of the Q-M3 haplogroup of 3,113 and 13,243 years before present; for the mtDNA-restriction haplotypes the estimated coalescence time was between 7,452 and 9,834 years before present. These results are compatible with the suggested time for the origin of the Chibchan group based on archeological, linguistic, and genetic evidence.     

Partial Sequence of a Sponge Mitochondrial Genome Reveals Sequence Similarity to Cnidaria in Cytochrome Oxidase Subunit II and the Large Ribosomal RNA Subunit

A 2550-bp portion of the mitochondrial genome of a Demosponge, genus Tetilla, was amplified from whole genomic DNA extract and sequenced. The sequence was found to code for the 3′ end of the 16S rRNA gene, cytochrome c oxidase subunit II, a lysine tRNA, ATPase subunit 8, and a 5′ portion of ATPase subunit 6. The Porifera cluster distinctly within the eumetazoan radiation, as a sister group to the Cnidaria. Also, the mitochondrial genetic code of this sponge is likely identical to that found in the Cnidaria. Both the full COII DNA and protein sequences and a portion of the 16S rRNA gene were found to possess a striking similarity to published Cnidarian mtDNA sequences, allying the Porifera more closely to the Cnidaria than to any other metazoan phylum. The gene arrangement, COII—tRNA^Lys—ATP8—ATP6, is observed in many Eumetazoan phyla and is apparently ancestral in the metazoa. Received: 24 November 1997 / Accepted: 14 September 1998     

中国13个民族7个Y-STR基因座遗传关系的研究

选择具有高度遗传多态性与稳定性的7个Y-STR位点分析宁夏回族与国内17个群体之间的群体遗传关系。采用Y-STR复合扩增、ABI全自动测序仪测序方法,结合基因扫描和自动分型技术获得150名宁夏回族个体DYS19、DYS389Ⅰ、DYS389Ⅱ、DYS390、DYS391、DYS392和DYS393七个基因座的基因频率,同时还收集国内17个群体这7个基因座的基因频率数据,计算他们间的遗传距离,进行聚类分析。结果显示:宁夏回族与北京汉族、云南汉族、安徽汉族、天津汉族、湖南汉族先聚为一类,然后与福建汉族相聚,最后与其他少数民族相聚,揭示回族与汉族在遗传结构上具有很大的相似性,为回族起源和发展过程中融入了大量汉族血缘的论点提供了遗传学证据。     

Early population differentiation in extinct aborigines from Tierra del Fuego-Patagonia: ancient mtDNA sequences and Y-chromosome STR characterization

Ancient mtDNA was successfully recovered from 24 skeletal samples of a total of 60 ancient individuals from Patagonia-Tierra del Fuego, dated to 100-400 years BP, for which consistent amplifications and two-strand sequences were obtained. Y-chromosome STRs (DYS434, DYS437, DYS439, DYS393, DYS391, DYS390, DYS19, DYS389I, DYS389II, and DYS388) and the biallelic system DYS199 were also amplified, Y-STR alleles could be characterized in nine cases, with an average of 4.1 loci per sample correctly typed. In two samples of the same ethnic group (Aonikenk), an identical and complete eight-loci haplotype was recovered. The DYS199 biallelic system was used as a control of contamination by modern DNA and, along with DYS19, as a marker of American origin. The analysis of both mtDNA and Y-STRs revealed DNA from Amerindian ancestry. The observed polymorphisms are consistent with the hypothesis that the ancient Fuegians are close to populations from south-central Chile and Argentina, but their high nucleotide diversity and the frequency of single lineages strongly support early genetic differentiation of the Fuegians through combined processes of population bottleneck, isolation, and/or migration, followed by strong genetic drift. This suggests an early genetic diversification of the Fuegians right after their arrival at the southernmost extreme of South America.     

Genomic organization of the extracellular coding region of the human FGFR4 and FLT4 genes: Evolution of the genes encoding receptor tyrosine kinases with immunoglobulin-like domains

Receptor tyrosine kinases with five, seven, and three Ig-like domains in their extracellular region are grouped in subclasses IIIa, IIIb, and IIIc, respectively. Here, we describe the genomic organization of the extracellular coding region of the human FGFR4 (IIIc) and FLT4 (IIIb) genes and compare it to that of the human FGFR1(IIIc), KIT, and FMS (IIIa). The results show that while genes belonging to the same subclass have an identical exon/intron structure in their extracellular coding region—as they do in their intracellular coding region—genes of related subclasses only have a similar exon/intron structure. These results strongly support the hypothesis that the genes of the three subclasses evolved from a common ancestor by duplications involving entire genes, already in pieces. Hypotheses on the origin of introns and on the difference in the number of extracellular Ig-like domains in the three gene subclasses are discussed. Received: 19 August 1996 / Accepted: 2 January 1997     

Accumulation of Transposable Elements in the Heterochromatin and on the Y Chromosome of Drosophila simulans and Drosophila melanogaster

The elements of the transposon families G, copia, mdg 1, 412, and gypsy that are located in the heterochromatin and on the Y chromosome have been identified by the Southern blotting technique in Drosophila simulans and D. melanogaster populations. Within species, the abundance of such elements differs between transposon families. Between species, the abundance in the heterochromatin and on the Y chromosome of the elements of the same family can differ greatly suggesting that differences within a species are unrelated to structural features of elements. By shedding some new light on the mechanism of accumulation of transposable elements in the heterochromatin, these data appear relevant to the understanding of the long-term interaction between transposable elements and the host genome. Received: 8 August 1997 / Accepted: 11 December 1997     

Elucidating the Population Histories and Transmission Dynamics of Papillomaviruses Using Phylogenetic Trees

Using gene genealogies constructed from gene sequence data, we show that both the mucosal and cutaneous papillomaviruses (PV)—supergroups A and B—appear to have been transmitted through susceptible populations faster than exponentially. The data and methods involved (1) examining the PV database for phylogenetic signal in an L1 open reading frame (ORF) fragment and an E1 ORF segment, (2) demonstrating that the same two fragments have evolved in a way consistent with a molecular clock, and (3) applying methods of phylogenetic tree analysis that test different scenarios for the dynamics of viral transmission within populations. The results indicate increases in PV populations of both supergroups A and B in the recent past. This form of the increases, which fit a null model of population growth with an exponent increasing in time, is compatible with the fact that human populations have grown at a faster than exponential rate, thus increasing the numbers of susceptible hosts for HPVs. There are, however, indications that the population of supergroup A has now stopped increasing in size. Received: 4 June 1996 / Accepted: 12 August 1996