Phylogeny of ladybirds (Coleoptera: Coccinellidae): are the subfamilies monophyletic?

The Coccinellidae (ladybirds) is a highly speciose family of the Coleoptera. Ladybirds are well known because of their use as biocontrol agents, and are the subject of many ecological studies. However, little is known about phylogenetic relationships of the Coccinellidae, and a precise evolutionary framework is needed for the family. This paper provides the first phylogenetic reconstruction of the relationships within the Coccinellidae based on analysis of five genes: the 18S and 28S rRNA nuclear genes and the mitochondrial 12S, 16S rRNA and cytochrome oxidase subunit I (COI) genes. The phylogenetic relationships of 67 terminal taxa, representative of all the subfamilies of the Coccinellidae (61 species, 37 genera), and relevant outgroups, were reconstructed using multiple approaches, including Bayesian inference with partitioning strategies. The recovered phylogenies are congruent and show that the Coccinellinae is monophyletic but the Coccidulinae, Epilachninae, Scymninae and Chilocorinae are paraphyletic. The tribe Chilocorini is identified as the sister-group of the Coccinellinae for the first time.

Origin and number of founders in an introduced insular primate: estimation from nuclear genetic data.

Cynomolgus macaques (Macaca fascicularis) were introduced on the island of Mauritius between 400 and 500 years ago and underwent a strong population expansion after a probable initial founding event. However, in practice, little is known of the geographical origin of the individuals that colonized the island, on how many individuals were introduced, and of whether the following demographic expansion erased any signal of this putative bottleneck. In this study, we asked whether the current nuclear genome of the Mauritius population retained a signature that would allow us to answer these questions. Altogether, 21 polymorphic autosomal and sex-linked microsatellites were surveyed from 81 unrelated Mauritius individuals and 173 individuals from putative geographical sources in Southeast Asia: Java, the Philippines islands and the Indochinese peninsula. We found that (i) the Mauritius population was closer to different populations depending on the markers we used, which suggests a possible mixed origin with Java playing most probably a major role; and (ii) the level of diversity was lower than the other populations but there was no clear and consistent bottleneck signal using either summary statistics or full-likelihood methods. However, summary statistics strongly suggest that Mauritius is not at mutation-drift equilibrium and favours an expansion rather than a bottleneck. This suggests that on a short time scale, population decline followed by growth can be difficult to deduce from genetic data based on mutation-drift theory. We then used a simple Bayesian rejection algorithm to estimate the number of founders under different demographic models (exponential, logistic and logistic with lag) and pure genetic drift. This new method uses current population size estimates and expected heterozygosity of Mauritius and source population(s). Our results indicate that a simple exponential growth is unlikely and that, under the logistic models, the population may have expanded from an initial effective number of individuals of 10-15. The data are also consistent with a logistic growth with different lag values, indicating that we cannot exclude past population fluctuation.

Population studies on an endemic troglobitic beetle: geographical patterns of genetic variation, gene flow and genetic structure compared with morphometric data.

Highly specialized obligatory cave beetles endemic to the French Pyrenees offer an opportunity to investigate the relative importance of environmental conditions and ecological characteristics on the organization of genetic variability, to describe the genetic structure of populations, and to assess the extent of gene flow between local populations in relation to geologic structure. Twenty-three geographically close populations of the beetle Speonomus hydrophilus occurring both in caves (reduced fluctuations in many abiotic parameters) and under the deepest layer of soil in mountains (more exposed to climatic variations) were studied. Significant genetic differentiation at 17 allozyme loci was found among populations in close proximity, as well as among those from distant parts of range. On a larger scale, genetic differences among populations appear to result from low dispersal rates between populations. The spatial patterning observed suggests that allozyme frequencies are not responding to environmentally controlled selection. Substantial genetic divergence (F(ST) = 0.112) occurred throughout the range, with important variation in levels of genetic variability (H: 0.065-0.184) among populations. A significant level of substructuring has occurred among the populations with four major geographic areas of similarity indicated. The substructuring of the species into regions suggests an influence of paleoclimatic gradient and paleoenvironment on the population's genetic structure. Also, founder effect and reduced gene flow appear to have influenced populations in the southeastern portion of the range.

Analysis of HLA-A/B recombinant families with new polymorphic markers.

The MHC in humans is a much studied region of the genome, the genes of which display a high rate of polymorphism and a high rate of linkage disequilibrium. Four families in which intra-class-I recombination has occurred have been analyzed with six polymorphic markers between HLA-A and -B in order to determine the full haplotypes of the whole families and to localize the points of crossover. The previously proposed order of the markers was confirmed by recombination mapping. In one family, the crossover was shown to have occurred in the 20-kb stretch of DNA bounded by the two markers (P3B and P5) between which no evidence of linkage disequilibrium was found, a region which constitutes of only about 1% of the distance between HLA-A and HLA-B. Although supportive of the suggestion of a hot spot of recombination in this region, based on the apparent lack of linkage disequilibrium between the markers P3B and P5, more such families need to be tested in order to confirm or refute this hypothesis.

Tumor necrosis factor microsatellites in four European populations.

The human genome contains a large number of interspersed simple repeat sequences that vary in length among individuals and can therefore serve as highly informative polymorphic markers. Several such variable sites (microsatellites) have been described within the TNF genes within the MHC. In this study, individuals from four Caucasian populations have been typed for three TNF-associated microsatellites in order to define their haplotypes. Of the 208 possible haplotypes, eight exist at a high frequency in all populations and account for approximately 60% of the haplotypes studied, but with marked variations in their frequencies among populations. A few population/sample-specific haplotypes have been identified. The ability of alleles to define haplotypes uniquely varies not only among the loci, but also among the alleles: some alleles displaying complete gametic association (linkage disequilibrium) and others displaying very little.

New highly polymorphic microsatellite marker in linkage disequilibrium with HLA-B.

The difficulty of molecular typing of the HLA class I genes and the relevance of the genes of this region to disease susceptibility and transplantation have provided an impetus to develop useful typing markers. We have characterized by polymerase chain reaction analysis a new highly informative CA repeat localized approximately 25-kb centromeric to the gene HLA-B and 10-kb telomeric to the gene MICA. Twelve alleles defined by length were found in a sample of French Basques, with the PIC being 0.82. A detailed haplotype analysis was performed to investigate the association between this microsatellite and two others markers of the region (HLA-B gene and TNF region microsatellite). The 10 haplotypes with the highest estimated frequencies show evidence of a gametic association or linkage disequilibrium. A very strong association between the expressed HLA-B polymorphism and microsatellite alleles was also revealed in this sample and confirmed in the workshop cells lines of the Fourth Asia-Oceania Histocompatibility Workshop. This marker can be used in the fine mapping of this region and the association with some alleles of HLA-B may allow the replacement of HLA-B typing at least in a preliminary study. Moreover, these studies support the hypothesis of a high mutability for large alleles in microsatellite loci.

Microsatellite, restriction fragment-length polymorphism, and sequence-specific oligonucleotide typing of the tumor necrosis factor region. Comparisons of the 4AOHW cell panel.

The location of the TNF and other genes in the central MHC and their possible relevance to disease susceptibility provided an impetus to develop useful typing markers. The 4AOHW undertook to assess the various markers available, including DNA sequence-based systems. A panel of well-characterized lymphoblastoid cell lines were typed by Nco I RLFP analysis, SSO typing, and TNF microsatellite typing. RFLP and SSO typing were relatively reproducible as judged by the blind replicates. The two techniques provided the same results with only one exception, and it would be reasonable to prefer SSO typing because of its advantages in terms of cost and time. Microsatellite typing was much more discriminating but, as expected, less robust in that some discrepancies were apparent. As a result of the workshop and subsequent testing, alleles and haplotypes were allocated to most cells within the 4AOHW panel, including 10W cells typed in previous studies. While there was evidence that microsatellites may be relatively stable, they have the potential to identify recent mutations within ancestral haplotypes.

Factors shaping genetic variation in the MHC of natural non-human primate populations.

Across a large distribution range, population-specific factors as well as pathogen-mediated selection may shape species genetic diversity in the major histocompatibility complex (MHC). We have studied genetic diversity and population differentiation in the MHC region of the Southeast Asian cynomolgus macaque (Macaca fascicularis fascicularis), a species with large and discontinuous range, in order to investigate the role of demography vs selection. Genetic variation was assessed at seven MHC microsatellites on 272 individuals from five populations (Indochina, Java, Borneo, Philippines, and Mauritius). A high genetic diversity was observed in all populations and the Philippines but also the Mauritius populations were the most genetically differentiated. The strength and extent of linkage disequilibrium (LD) (up to 4 Mb) varies across populations mainly because of demographic factors. In Indochina, the complete lack of LD could be the signature of ancient hybridization between cynomolgus and rhesus macaques in the Indochinese peninsula. With the additional support of seven autosomal microsatellites, tests for outlier loci based on intrapopulation diversity and interpopulation differentiation (using F-statistic) allowed to dissociate demographic from selective histories: (i) demographic history may itself explain levels of MHC variability in the Mauritius populations and (ii) positive selection could be responsible for the Philippines population differentiation, especially in the MHC class II region. Among various pathogens, Plasmodium knowlesi and Plasmodium coatneyi are two likely candidates to explain the higher frequency of some MHC haplotypes. Indeed, literature describes low parasitemia in the Philippines individuals, contrasting with fatal infections provoked by these parasites in other cynomolgus macaque populations.

Unbiased interpretation of haplotypes at duplicated microsatellites.

The Y-chromosome is a powerful tool for population geneticists to study human evolutionary history. Haploid and largely non-recombining, it should contain a simple record of past mutational events. However, this apparent simplicity is compromised by Y-linked duplicons, which make up approximately 35% of this chromosome; 25% of these duplicons are large inverted repeats (palindromes). For microsatellites lying in these palindromes, two loci cannot be easily distinguished due to PCR co-amplification, and this order misspecification of alleles generates an additional variance component. Due to this ambiguity, population geneticists have traditionally used an arbitrary method to assign the alleles (shorter allele to locus 1, larger allele to locus 2). Here, we simulate these posterior estimate distributions under three different novel allele assignment priors and compare this with the original method. We use a sample of 33 human populations, typed for duplicated microsatellites lying within palindrome P8, to illustrate our approach. We show that both intra- and inter-population statistics can be dramatically affected by order misspecification. Surprisingly, matrices of pairwise F-statistics or distance estimates appear far less sensitive to order misspecification and remain relatively unchanged under the priors considered, suggesting that these microsatellites can be considered as useful markers for population genetic studies using an appropriate data treatment. Duplicated microsatellites represent an attractive source of information to investigate the extensive structural polymorphism observed among human Y chromosomes, as well as processes of intra-chromosomal gene conversion acting between duplicons.

Multiplexed microsatellites for rapid identification and characterization of individuals and populations of Cercopithecidae.

Cross-amplification of 15 human microsatellites was performed successfully in cynomolgus (Macaca fascicularis) and rhesus (M. mulatta) macaques and 11 other Cercopithecidae species of biomedical and conservation relevance. To allow for quick, efficient, and high-throughput genotyping to assess intra- and interspecific genetic variation, we performed three multiplex sets, each comprised of five markers from different parts of the genome (i.e., autosomes, the MHC region, and the X-chromosome). These multiplex sets are likely to reveal allelic divergence between taxa, which could be used for their discrimination. Population studies on three regional populations of M. fascicularis and one of M. mulatta revealed that most of the loci, with the exception of one monomorphic locus, displayed polymorphisms (the expected heterozygosities were 0.48-0.91 for M. fascicularis, and 0.61-0.93 for M. mulatta), which makes them useful for population genetics. For the multiplex set M1, including the nonlinked autosomal markers, low probabilities of identity were observed: P(ID) values ranged from 8 x 10(-7) to 3 x 10(-5). This multiplex set is reliable for forensic applications, such as individual identification, parentage testing, and kinship analysis, in wild and captive populations.

Trans-speciation maintenance in the MHC region of a polymorphism which includes a polymorphic dinucleotide locus, and the de novo arisal of a polymorphic tetranucleotide microsatellite.

Alleles and the surrounding regions of DQCAR, a dinucleotide repeat tightly linked to HLA-DQB1, were sequenced in a range of primate species including man. Three polymorphic regions can usefully be defined in the description of these sequences: the dinucleotide GT repeat itself, the anonymous region 5' of this repeat, and a variable CTGT repeat in the 3' region. The 5' sequence displayed six alleles in the individuals studied. One of these alleles was invariably associated with substitutions in the GT repeat and absence of the CTGT repeat, the others with pure, polymorphic GT repeats and variation in the numbers of CTGT repeats. Haplotypes can be classified by the allele in the 5' region. Those carrying allele 1 were only found in man, those with allele 2 in man, chimpanzee and gorilla. The third haplotype (indicated by the presence of allele 3) was found in chimpanzee, gorilla and orang-utan, the fourth in chimpanzee and gibbon, the fifth in baboon, guenon and mangabey and the sixth in guenon and macaque. The alleles in the 5' region, but from different species, are thus often more similar than alleles from the same species, a phenomenon already shown for some HLA genes. This suggests that major histocompatibility sequences and surrounding sequences shared a correlated evolutionary history. The new polymorphic tetranucleotide microsatellite (CTGT, 3rd region) has possibly arisen de novo from the pre-existing dinucleotide GT. This study provides information not only on the molecular evolution of this particular microsatellite but also of the trans-speciation maintenance of polymorphism of its surrounding sequences.

Microsatellite allelic homoplasy due to variable flanking sequences.

Microsatellite DNA sequences have become the dominant source of nuclear genetic markers for most applications. It is important to investigate the basis of variation between alleles and to know if current assumptions about the mechanisms of microsatellite mutation (that is to say, variations involving simple changes in the number of repeat) are correct. We have characterized, by DNA sequencing, the human alleles of a new highly informative (CA)n repeat localized approximately 20 kb centromeric to the HLA-B gene. Although 12 alleles were identified based on conventional length criteria, sequencing of the alleles demonstrated that differences between alleles were found to be more complex than previously assumed: A high degree of microsatellite variability is due to variation in the region immediately flanking the repeat. These data indicate that the mutational process which generates polymorphism in this region has involved not only simple changes in the number of dinucleotide CA repeats but also perturbations in the nonrepeated 5' and 3' flanking sequences. Three families of alleles (not visible from the overall length of the alleles), with presumably separate evolutionary histories, exist and can yield to homoplasy of size. Effectively, we can observe alleles of the same size with different internal structures which are separated by a significant amount of variation. Although allelic homoplasy for non-interrupted microsatellite loci has been suggested between different species, it has not been unequivocally demonstrated within species. A strong association is noted between alleles defined at the sequence level and HLA-B alleles. The observation of several families of alleles at the population level provides information about the evolutionary history and mutation processes of microsatellites and may have implications for the use of these markers in phylogenetic, linkage disequilibrium studies, and gene mapping.

Evolution of an Alu DNA element of type Sx in the lineage of primates and the origin of an associated tetranucleotide microsatellite.

A 394-bp DNA fragment, which in human is on chromosome 6 near the MOG (myelin oligodendrocyte glycoprotein) gene and encompasses an Alu element and an associated tetranucleotide microsatellite, was sequenced from a large range of primate species to follow its evolutionary divergence and to understand the origin of the microsatellite. This Alu element is found at the same orthologous position in all primates sequenced, but the tetranucleotide repeat is present only in Catarrhini between the 3'-oligo(dA) of the Alu element and the 3' flanking direct repeat. Little intraspecific variation was found. Sequence identity values for this orthologous primate Alu averaged 90% (82-99%) with transitions comprising between 70% and 100% of the observed nucleotide substitutions. Although the insertion of the Alu element predates the separation of these species, the original sequence of the site of integration can still be identified. This identification of the direct repeats suggests an active role of the oligo(dA) of the Alu element in the origin of the tetranucleotide repeats. The microsatellite probably appeared after the insertion of the Alu element, early in the lineage leading to the common ancestor of the hominoids and the Old World monkeys.

Polymorphism, monomorphism, and sequences in conserved microsatellites in primate species.

Dimeric short tandem repeats are a source of highly polymorphic markers in the mammalian genome. Genetic variation at these hypervariable loci is extensively used for linkage analysis, for the identification of individuals, and may be useful for interpopulation and interspecies studies. In this paper, we analyze the variability and the sequences of a segment including three microsatellites, first described in man, in several species of primates (chimpanzee, orangutan, gibbon, and macaque) using the heterologous primers (man primers). This region is located on the human chromosome 6p, near the tumor necrosis factor genes, in the major histocompatibility complex. The fact that these primers work in all species studied indicates that they are conserved throughout the different lineages of the two superfamilies, the Hominoidea and the Cercopithecidea, represented by the macaques. However, the intervening sequence displays intraspecific and interspecific variability. The sites of base substitutions and the insertion/deletion events are not evenly distributed within this region. The data suggest that it is necessary to have a minimal number of repeats to increase the rate of mutation sufficiently to allow the development of polymorphism. In some species, the microsatellites present single base variations which reduce the number of contiguous repeats, thus apparently slowing the rate of additional slippage events. Species with such variations or a low number of repeats are monomorphic. These microsatellite sequences are informative in the comparison of closely related species and reflect the phylogeny of the Old World monkeys, apes, and man.

Haplotypes without children: PCR applied to close loci on individual human sperm.

Indirect evidence of interindividual variation of the recombination fraction is considerable, but as yet direct evidence is lacking. Such interindividual variation could explain the widely observed gametic association in the major histocompatibility complex (MHC) region. A possible approach to this is the application of the polymerase chain reaction to closely linked markers on individual sperm. The feasibility of this approach is investigated here. Three loci in the human MHC (chromosome 6) were used in this study. Two are polymorphisms identifiable by oligonucleotide dot blotting, and one is a microsatellite polymorphism. The genetic distances between these loci are 0.5 and 1 centi-Morgan. We were able to double-type most cases for the two markers HLA-DPB and HLA-DRB in single haploid cells. The rate of double typings is comparable with that expected from the hypothesis of independent failures at two loci. However, this rate of failure remains too great to allow an analysis of segregation because allele specific failure cannot be ruled out. All double and triple typings (implicating the microsatellite) were perfectly correlated with each other; thus no recombinations were identified. With such close markers recombination would have been surprising and may indicate an inherent accuracy of positive haplotyping.

Conservation and evolution of microsatellite loci in primate taxa.

Microsatellites are promising genetic markers for the study of demographic structure and phylogenetic history in populations. However, little information exists on the molecular nature of the repeats and their flanking sequences of a same microsatellite in a large range of species. In this study, we report polymorphism and consensus sequences of eight microsatellite loci using human primers in 20 primate species. The results show size polymorphism in almost all species and microsatellites. These loci are therefore useful markers for population genetic studies between populations of the same species. Insertion/deletion events are frequent in the flanking regions, the majority concerning several contiguous bases. This is in contrast with the more usual single base pair events in non-coding regions. The ranges of allele lengths in non-human primates often show no overlap with that of human, usually due to the deletion/insertion events in the flanking sequences, producing smaller allele lengths rather than smaller numbers of repeats. The use of length of PCR product will bias the inter-species interpretation reducing the number of observable alleles and treating as the same allele very divergent molecular sequences. Caution should be used when employing microsatellites in cross-species comparisons in which the species under study are separated by significant amounts of evolutionary time: in such cases allele comparison cannot be based on lengths alone.

A precise meiotic map in the class I region of the human major histocompatibility complex.

Human families in which recombinant meiotic event(s) are known to have occurred are powerful tools with which to analyze more precisely the structures of defined genomic regions, especially unstable areas. Such families allow the determination of the haplotypes of each member and, taking into account the recombinant event, it is possible to localize very precisely the point of crossover. Using families in which crossovers between the genes HLA-A and -B have occurred, we have constructed a meiotic map localizing the meiotic breakpoint events with respect to both anonymous markers and the principal genes of the region. Such mapping, which depends on the direct analysis of genomic DNA, is essential for fine structural analysis and is a powerful means of verification of the order and the localization of markers: physical mapping alone, using yeast artificial chromosomes, presents some uncertainties due to the numerous chimeras and inversions that can be produced. The establishment of this map will allow us to determine efficiently the precise location for new markers already localized to the map region. Three microsatellites (D6S265, D6S276, and D6S306), localized in the HLA region by linkage analysis, have been precisely located with respect to the points of recombination in the class I region. The sites of meiotic recombination in the MHC class I region seem to be not randomly distributed but in the majority of cases occurred between HLA-C and the microsatellite D6S265. This study also shows two cases of abnormal segregation of alleles. We discuss how these mutations correspond to a spontaneous mutation event at the somatic or germinal level.

A fine-scale comparison of the human and chimpanzee genomes: linkage, linkage disequilibrium and sequence analysis.

We have performed a fine-scale comparative study of the human and chimpanzee genomes, using linkage, linkage disequilibrium and sequence analyses on microsatellite loci spanning a region of approximately 30 cM on human chromosome 4p. Our results extend the findings of previous studies that indicated virtually complete conservation between the human and chimpanzee genomes at the chromosomal and sub-chromosomal level and support the hypothesis, derived from previous analyses of mitochondrial DNA, that chimpanzee populations are more diverse than human ones. By sequencing several human and chimpanzee alleles of two microsatellites we showed that base substitutions that diminish the length of perfect repeats (but do not change allele sizes) are probably responsible for the low heterozygosity of these loci in chimpanzees; our results suggest that the evolutionary history of microsatellites should not be inferred from comparisons of mean allele lengths between populations or species.

Testing the reliability of noninvasive genetic sampling by comparing analyses of blood and fecal samples in Barbary macaques (Macaca sylvanus).

Genetic studies of wild animal populations are often hindered by difficulties in obtaining blood samples. Recent advances in molecular biology have allowed the use of noninvasive samples as sources of DNA (e.g., hair or feces), but such samples may provide low-quality DNA and prevent the determination of true genotypes in subsequent DNA analysis. We present a preliminary study aimed at assessing the reliability of using fecal samples for genotyping in Barbary macaques (Macaca sylvanus). The test was performed on samples of blood and feces from 11 captive animals, using three dinucleotide microsatellites. The CTAB DNA extraction method was found to be the most relevant for Barbary macaque feces, yielding successful amplification at all loci for 70% of PCRs. All the fecal samples tested gave correct genotypes at least once for each locus when referenced against blood-derived genotypes. An average of 18.3% of PCRs displayed spurious genotypes (false homozygous or false allele). The minimum theoretical probability required to obtain a 100% accurate genotype is 0.74, based on the criterion that a correct genotype is assessed only if it was observed at least twice. The observed probability of obtaining a correct genotype from three PCRs, based on our genotyping results, was greater (0.81 on average) than the minimum threshold. In conclusion, our comparison of blood and fecal samples showed that fecal sampling is a reliable tool for the further study of wild Barbary macaque populations.

Microsatellites in the HLA region: on overview.

Microsatellites are repeats of a DNA base motif (1-6 bp, mostly CA repeats) up to 100 times; they are distributed regularly all over the genome. Many of them are polymorphic and their high polymorphism is based upon a variable number of repeats. They are widely used for genetic mapping, linkage analysis, population genetics, evolutionary studies and in forensic medicine. Such markers have also been described in the HLA region since 1991, and a growing interest in their potential applications is being expressed. The aims of this review are: 1) to outline the presently available information from literature and molecular databases concerning 53 microsatellites in the HLA region (localization, type of repeat, number of alleles, heterozygosity, primers used for amplification); 2) to address the question of technical pitfalls when using such markers; 3) to discuss specific features such as their mutation rate (10 (-3) to 10 (-6), which is higher than that reported for HLA genes, and their linkage disequilibrium with HLA alleles; 4) to present an integrated map of microsatellites and genes of this region; and 5) to provide a synopsis of their different applications in HLA-related fields (disease studies, population genetics, recombination point studies, HLA region mapping, transplantation) along with perspectives for future use. Although some HLA region microsatellites have already been applied to the analysis of more than 10 diseases, it is now evident that their use in population genetics and the determination of genomic compatibility in bone marrow transplantation represent growing areas of application.