Phylogeny of the macaques (Cercopithecidae: Macaca) based on Alu elements

Genus Macaca (Cercopithecidae: Papionini) is one of the most successful primate radiations. Despite previous studies on morphology and mitochondrial DNA analysis, a number of issues regarding the details of macaque evolution remain unsolved. Alu elements are a class of non-autonomous retroposons belonging to short interspersed elements that are specific to the primate lineage. Because retroposon insertions show very little homoplasy, and because the ancestral state (absence of the SINE) is known, Alu elements are useful genetic markers and have been utilized for analyzing primate phylogenentic relationships and human population genetic relationships. Using PCR display methodology, 298 new Alu insertions have been identified from ten species of macaques. Together with 60 loci reported previously, a total of 358 loci are used to infer the phylogenetic relationships of genus Macaca. With regard to earlier unresolved issues on the macaque evolution, the topology of our tree suggests that: 1) genus Macaca contains four monophyletic species groups; 2) within the Asian macaques, the silenus group diverged first, and members of the sinica and fascicularis groups share a common ancestor; 3) Macaca arctoides are classified in the sinica group. Our results provide a robust molecular phylogeny for genus Macaca with stronger statistical support than previous studies. The present study also illustrates that SINE-based approaches are a powerful tool in primate phylogenetic studies and can be used to successfully resolve evolutionary relationships between taxa at scales from the ordinal level to closely related species within one genus.     

Evolution of assortative mating following selective introgression of pigmentation genes between two Drosophila species

Adaptive introgression is ubiquitous in animals, but experimental support for its role in driving speciation remains scarce. In the absence of conscious selection, admixed laboratory strains of Drosophila asymmetrically and progressively lose alleles from one parental species and reproductive isolation against the predominant parent ceases after 10 generations. Here, we selectively introgressed during 1 year light pigmentation genes of D. santomea into the genome of its dark sibling D. yakuba, and vice versa. We found that the pace of phenotypic change differed between the species and the sexes and identified through genome sequencing common as well as distinct introgressed loci in each species. Mating assays showed that assortative mating between introgressed flies and both parental species persisted even after 4 years (~60 generations) from the end of the selection. Those results indicate that selective introgression of as low as 0.5% of the genome can beget morphologically distinct and reproductively isolated strains, two prerequisites for the delimitation of new species. Our findings hence represent a significant step toward understanding the genome‐wide dynamics of speciation‐through‐introgression.     

G-banding in the barbary macaque (Macaca sylvanus, Linnaeus, 1758) of Morocco

G-banding was performed on the lymphocytes of six barbary macaques (Macaca sylvanus). Four chromosomes appear to have polymorphic variants: number 1, number 4, number 16, and the “marker” chromosome, number 9. These polymorphic forms differe from the G-banding reported for Macaca mulatta, Macaca fascicularis, and Macaca arctoides. It may be that evolutionary relationships among various macaque species can be ascertained by using such chromosomal polymorphisms.     

Ancient genome‐wide admixture extends beyond the current hybrid zone between Macaca fascicularis and M. mulatta

Macaca fascicularis and Macaca mulatta are two of the most commonly used laboratory macaques, yet their genetic differences at a genome‐wide level remain unclear. We analysed the multilocus DNA sequence data of 54 autosomal loci obtained from M. fascicularis samples from three different geographic origins and M. mulatta samples of Burmese origin. M. fascicularis shows high nucleotide diversity, four to five times higher than humans, and a strong geographic population structure between Indonesian‐Malaysian and Philippine macaques. The pattern of divergence and polymorphism between M. fascicularis and M. mulatta shows a footprint of genetic exchange not only within their current hybrid zone but also across a wider range for more than 1 million years. However, genetic admixture may not be a random event in the genome. Whereas randomly selected genic and intergenic regions have the same evolutionary dynamics between the species, some cytochrome oxidase P450 (CYP) genes (major chemical metabolizing genes and potential target genes for local adaptation) have a significantly larger species divergence than other genes. By surveying CYP3A5 gene sequences of more than a hundred macaques, we identified three nonsynonymous single nucleotide polymorphisms that were highly differentiated between the macaques. The mosaic pattern of species divergence in the genomes may be a consequence of genetic differentiation under ecological adaptation and may be a salient feature in the genomes of nascent species under parapatry.     

Male external genitalia and systematic relationships of the Japanese macaque (Macaca fuscata Blyth, 1875)

Penial morphology was studied in six preserved specimens of the Japanese macaque,Macaca fuscata. Structure of the glans penis and baculum inM. fuscata is essentially similar to that in the rhesus macaque,M. mulatta. With age the skin of the penis inM. fuscata apparently becomes progressively more darkly and extensively pigmented, and the baculum becomes progressively longer. One immature specimen examined, with penial skin unpigmented and baculum cartilaginous, appears to be abnormally retarded in penial development. Penial morphology suggests that the allopatric speciesM. fuscata, M. mulatta, M. cyclopis, andM. fascicularis constitute a closely interrelated group. Contrary to a previous report, there is no evidence that penial morphology inM. fuscata is transitional to that inM. radiata andM. sinica. InM. fuscata, M. mulatta, andM. fascicularis, latitudinal variation of body size and relative tail length conforms toBergmann's rule andAllen's rule; external proportions in three specimens ofM. cyclopis only partly conform to this pattern.     

Genetic polymorphism of alpha-1-antitrypsin in macaque species

Thirty-two phenotypes of alpha-1-antitrypsin (al-AT) controlled by 18 codominant alleles of PI_mac locus were identified by isoelectric focusing in 1.121 macaques of nine species. In terms of al-AT polymorphism macaque species vary markedly from nearly monomorphic (Macaca mulatta, M. arctoides) to highly polymorphic (M. nemestrina, M. fascicularis). Only 4 of 18 PI_mac alleles are common for two or more macaque species, the rest of the alleles are species-specific.     

Species boundaries in the messy middle—A genome‐scale validation of species delimitation in a recently diverged lineage of coastal fog desert lichen fungi

Species delimitation among closely related species is challenging because traditional phenotype‐based approaches, for example, using morphology, ecological, or chemical characteristics, may not coincide with natural groupings. With the advent of high‐throughput sequencing, it has become increasingly cost‐effective to acquire genome‐scale data which can resolve previously ambiguous species boundaries. As the availability of genome‐scale data has increased, numerous species delimitation analyses, such as BPP and SNAPP+Bayes factor delimitation (BFD*), have been developed to delimit species boundaries. However, even empirical molecular species delimitation approaches can be biased by confounding evolutionary factors, for example, hybridization/introgression and incomplete lineage sorting, and computational limitations. Here, we investigate species boundaries and the potential for micro‐endemism in a lineage of lichen‐forming fungi, Niebla Rundel & Bowler, in the family Ramalinaceae by analyzing single‐locus and genome‐scale data consisting of (a) single‐locus species delimitation analysis using ASAP, (b) maximum likelihood‐based phylogenetic tree inference, (c) genome‐scale species delimitation models, e.g., BPP and SNAPP+BFD, and (d) species validation using the genealogical divergence index (gdi). We specifically use these methods to cross‐validate results between genome‐scale and single‐locus datasets, differently sampled subsets of genomic data and to control for population‐level genetic divergence. Our species delimitation models tend to support more speciose groupings that were inconsistent with traditional taxonomy, supporting a hypothesis of micro‐endemism, which may include morphologically cryptic species. However, the models did not converge on robust, consistent species delimitations. While the results of our analysis are somewhat ambiguous in terms of species boundaries, they provide a valuable perspective on how to use these empirical species delimitation methods in a nonmodel system. This study thus highlights the challenges inherent in delimiting species, particularly in groups such as Niebla, with complex, relatively recent phylogeographic histories.     

Genetic differentiation and connectivity of morphological types of the broadcast‐spawning coral Galaxea fascicularis in the Nansei Islands,Japan

Population connectivity resulting from larval dispersal is essential for the maintenance or recovery of populations in marine ecosystems, including coral reefs. Studies of species diversity and genetic connectivity within species are essential for the conservation of corals and coral reef ecosystems. We analyzed mitochondrial DNA sequence types and microsatellite genotypes of the broadcast‐spawning coral, Galaxea fascicularis, from four regions in the subtropical Nansei Islands in the northwestern Pacific Ocean. Two types (soft and hard types) of nematocyst morphology are known in G. fascicularis and are significantly correlated with the length of a mitochondrial DNA noncoding sequence (soft type: mt‐L; hard type: mt‐S type). Using microsatellites, significant genetic differentiation was detected between the mitochondrial DNA sequence types in all regions. We also found a third genetic cluster (mt‐L+), and this unexpected type may be a cryptic species of Galaxea. High clonal diversity was detected in both mt‐L and mt‐S types. Significant genetic differentiation, which was found among regions within a given type (F _ST = 0.009–0.024, all Ps ≤ 0.005 in mt‐L; 0.009–0.032, all Ps ≤ 0.01 in mt‐S), may result from the shorter larval development than in other broadcast‐spawning corals, such as the genus Acropora. Nevertheless, intraspecific genetic diversity and connectivity have been maintained, and with both sexual and asexual reproduction, this species appears to have a potential for the recovery of populations after disturbance.     

Characteristics of the mitochondrial genome of Rana omeimontis and related species in Ranidae: Gene rearrangements and phylogenetic relationships

The Omei wood frog (Rana omeimontis), endemic to central China, belongs to the family Ranidae. In this study, we achieved detail knowledge about the mitogenome of the species. The length of the genome is 20,120 bp, including 13 protein‐coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a noncoding control region. Similar to other amphibians, we found that only nine genes (ND6 and eight tRNA genes) are encoded on the light strand (L) and other genes on the heavy strand (H). Totally, The base composition of the mitochondrial genome included 27.29% A, 28.85% T, 28.87% C, and 15.00% G, respectively. The control regions among the Rana species were found to exhibit rich genetic variability and A + T content. R. omeimontis was clustered together with R. chaochiaoensis in phylogenetic tree. Compared to R. amurensis and R. kunyuensi, it was more closely related to R. chaochiaoensis, and a new way of gene rearrangement (ND6‐trnE‐Cytb‐D‐loop‐trnL2 (CUN)‐ND5‐D‐loop) was also found in the mitogenome of R. amurensis and R. kunyuensi. Our results about the mitochondrial genome of R. omeimontis will contribute to the future studies on phylogenetic relationship and the taxonomic status of Rana and related Ranidae species.     

The coevolution of male and female genitalia in a mammal: A quantitative genetic insight

Male genitalia are among the most phenotypically diverse morphological traits, and sexual selection is widely accepted as being responsible for their evolutionary divergence. Studies of house mice suggest that the shape of the baculum (penis bone) affects male reproductive fitness and experimentally imposed postmating sexual selection has been shown to drive divergence in baculum shape across generations. Much less is known of the morphology of female genitalia and its coevolution with male genitalia. In light of this, we used a paternal half-sibling design to explore patterns of additive genetic variation and covariation underlying baculum shape and female vaginal tract size in house mice (Mus musculus domesticus). We applied a landmark-based morphometrics approach to measure baculum size and shape in males and the length of the vaginal tract and width of the cervix in females. Our results reveal significant additive genetic variation in house mouse baculum morphology and cervix width, as well as evidence for genetic covariation between male and female genital measures. Our data thereby provide novel insight into the potential for the coevolutionary divergence of male and female genital traits in a mammal.     

Complete Taiwanese Macaque (Macaca cyclopis) Mitochondrial Genome: Reference-Assisted de novo Assembly with Multiple k-mer Strategy

The Taiwanese (Formosan) macaque (Macaca cyclopis) is the only nonhuman primate endemic to Taiwan. This primate species is valuable for evolutionary studies and as subjects in medical research. However, only partial fragments of the mitochondrial genome (mitogenome) of this primate species have been sequenced, not mentioning its nuclear genome. We employed next-generation sequencing to generate 2 x 90 bp paired-end reads, followed by reference-assisted de novo assembly with multiple k-mer strategy to characterize the M. cyclopis mitogenome. We compared the assembled mitogenome with that of other macaque species for phylogenetic analysis. Our results show that, the M. cyclopis mitogenome consists of 16,563 nucleotides encoding for 13 protein-coding genes, 2 ribosomal RNAs and 22 transfer RNAs. Phylogenetic analysis indicates that M. cyclopis is most closely related to M. mulatta lasiota (Chinese rhesus macaque), supporting the notion of Asia-continental origin of M. cyclopis proposed in previous studies based on partial mitochondrial sequences. Our work presents a novel approach for assembling a mitogenome that utilizes the capabilities of de novo genome assembly with assistance of a reference genome. The availability of the complete Taiwanese macaque mitogenome will facilitate the study of primate evolution and the characterization of genetic variations for the potential usage of this species as a non-human primate model for medical research.     

Genetic markers in Andean Puya species (Bromeliaceae) with implications on plastome evolution and phylogeny

The Andean plant endemic Puya is a striking example of recent and rapid diversification from central Chile to the northern Andes, tracking mountain uplift. This study generated 12 complete plastomes representing nine Puya species and compared them to five published plastomes for their features, genomic evolution, and phylogeny. The total size of the Puya plastomes ranged from 159,542 to 159,839 bp with 37.3%–37.4% GC content. The Puya plastomes were highly conserved in organization and structure with a typical quadripartite genome structure. Each of the 17 consensus plastomes harbored 133 genes, including 87 protein‐coding genes, 38 tRNA (transfer RNA) genes, and eight rRNA (ribosomal RNA) genes; we found 69–78 tandem repeats, 45–60 SSRs (simple sequence repeats), and 8–22 repeat structures among 13 species. Four protein‐coding genes were identified under positive site‐specific selection in Puya. The complete plastomes and hypervariable regions collectively provided pronounced species discrimination in Puya and a practical tool for future phylogenetic studies. The reconstructed phylogeny and estimated divergence time for the lineage suggest that the diversification of Puya is related to Andean orogeny and Pleistocene climatic oscillations. This study provides plastome resources for species delimitation and novel phylogenetic and biogeographic studies.     

Copulatory behavior in relation to anatomical characteristics of three macaques

Anatomy of copulatory organs and patterns of copulatory behavior were studied in three species of Macaca. The copulatory behavior of the rhesus monkey (M. mulatta) and crab-eating macaque (M. fascicularis) were similar. Mounting time was longer in stump-tailed macaque (M. arctoides) (48 seconds) than in the rhesus monkey or crab-eating macaque (8.5 seconds). Species differences were correlated with anatomical differences in the female and male reproductive organs. In the stump-tailed macaque the male has a relatively long penis, and the female has an elaborate vestibular colliculus in the vestibule. Both rhesus and crab-eating macaques have a short penis, and the female has no vestibular colliculus. Under scanning electron microscope, cone-like projections (3–4 mm) were observed on the surface of the glans and corpus penis of M. arctoides. In M. fascicularis shorter projections (0.05–0.1 mm) were observed on the glans penis but not the corpus penis.     

Do Heliconius butterfly species exchange mimicry alleles?

Hybridization has the potential to transfer beneficial alleles across species boundaries, and there are a growing number of examples in which this has apparently occurred. Recent studies suggest that Heliconius butterflies have transferred wing pattern mimicry alleles between species via hybridization, but ancestral polymorphism could also produce a signature of shared ancestry around mimicry genes. To distinguish between these alternative hypotheses, we measured DNA sequence divergence around putatively introgressed mimicry loci and compared this with the rest of the genome. Our results reveal that putatively introgressed regions show strongly reduced sequence divergence between co-mimetic species, suggesting that their divergence times are younger than the rest of the genome. This is consistent with introgression and not ancestral variation. We further show that this signature of introgression occurs at sites throughout the genome, not just around mimicry genes.     

Genome sequencing and comparison of two nonhuman primate animal models, the cynomolgus and Chinese rhesus macaques

The nonhuman primates most commonly used in medical research are from the genus Macaca. To better understand the genetic differences between these animal models, we present high-quality draft genome sequences from two macaque species, the cynomolgus/crab-eating macaque and the Chinese rhesus macaque. Comparison with the previously sequenced Indian rhesus macaque reveals that all three macaques maintain abundant genetic heterogeneity, including millions of single-nucleotide substitutions and many insertions, deletions and gross chromosomal rearrangements. By assessing genetic regions with reduced variability, we identify genes in each macaque species that may have experienced positive selection. Genetic divergence patterns suggest that the cynomolgus macaque genome has been shaped by introgression after hybridization with the Chinese rhesus macaque. Macaque genes display a high degree of sequence similarity with human disease gene orthologs and drug targets. However, we identify several putatively dysfunctional genetic differences between the three macaque species, which may explain functional differences between them previously observed in clinical studies.     

Baculum length and copulatory behavior in primates

The length of the baculum (os penis) was measured in 74 adult males representing 46 primate species. These data, and a review of previously published measurements, indicate that variation in baculum length among primates is related to taxonomic and behavioral differences. Thus, many New World monkeys have shorter bacula, relative to body weight, than Old World monkeys. The baculum is shorter in colobine monkeys than in cercopithecines. Among the great apes, reduction of the baculum is more pronounced in Pan and Gorilla than in Pongo. Very long bacula are found in some nocturnal prosimians (eg, Lorisidae) and also in Macaca arctoides. A comparison of baculum length relative to body weight was made in 34 species for which detailed information on copulatory behavi or was available. The presence of an elongated baculum was shown to correlate with copulatory patterns involving prolonged intromission and/or the maintenance of intromission during the postejaculatory interval (eg, Galago crassicaudatus, Loris tardigradus, M, arctoides). The evolutionary significance of these observations is discussed and it is suggeted that similar copulatory patterns may occur in species with elongated bacula (eg, Daubentonia, Perodicticus) for which behavioral data are lacking at present. The same hypothesis also applies to an extinct adapid primate which possessed a very large baculum.     

Phylogenetic Relationships of the Macaques (Cercopithecidae: Macaca), Inferred from Mitochondrial DNA Sequences

To study the phylogenetic relationships of the macaques, five gene fragments were sequenced from 40 individuals of eight species: Macaca mulatta, M. cyclopis, M. fascicularis, M. arctoides, M. assamensis, M. thibetana, M. silenus, and M. leonina. In addition, sequences of M. sylvanus were obtained from Genbank. A baboon was used as the outgroup. The phylogenetic trees were constructed using maximum-parsimony and Bayesian methods. Because five gene fragments were from the mitochondrial genome and were inherited as a single entity without recombination, we combined the five genes into a single analysis. The parsimony bootstrap proportions we obtained were higher than those from earlier studies based on the combined mtDNA dataset. Excluding M. arctoides, our results are generally consistent with the classification of Delson (1980). Our phylogenetic analyses agree with earlier studies suggesting that the mitochondrial lineages of M. arctoides share a close evolutionary relationship with the mitochondrial lineages of the fascicularis group of macaques (and M. fascicularis, specifically). M. mulatta (with respect to M. cyclopis), M. assamensis assamensis (with respect to M. thibetana), and M. leonina (with respect to M. silenus) are paraphyletic based on our analysis of mitochondrial genes.     

Influence of the large‐Z effect during contact between butterfly sister species

Recently diverged butterfly populations in North America have been found to exhibit high levels of divergence on the Z chromosome relative to autosomes, as measured by fixation index, Fst. The pattern of divergence appears to result from accumulation of incompatible alleles, obstructing introgression on the Z chromosome in hybrids (i.e., the large‐Z effect); however, it is unknown whether this mechanism is sufficient to explain the data. Here, we simulate the effects of hybrid incompatibility on interbreeding butterfly populations using a model in which populations accumulate cross‐incompatible alleles in allopatry prior to contact. We compute statistics for introgression and population divergence during contact between model populations and compare our results to those for 15 pairs of butterfly species interbreeding along a suture zone in central Texas. Time scales for allopatry and contact in the model are scaled to glacial and interglacial periods during which real populations evolved in isolation and contact. We find that the data for butterflies are explained well by an otherwise neutral model under slow fusion conditions. In particular, levels of divergence on the Z chromosome increase when interacting clusters of genes are closely linked, consistent with clusters of functionally related genes in butterfly genomes.     

Vocal communication of wild bonnet macaques (Macaca radiata)

Field observations and spectrographic analyses of sound recordings of South Indian bonnet macaques revealed a vocal repertoire of at least 25 basic patterns. The repertoire consists of well separated sound classes and acoustic categories connected by structural intergradation. Besides structural variations within and between different elements of the repertoire, the vocal system ofMacaca radiata is characterized by regular combinations of particular basic patterns. These combinations occurred not only between calls of similar structure and function but also between calls usually emitted in entirely different social contexts. According to the qualitative analysis, sex-specific asymmetries of the vocal behaviour were less pronounced than age-dependent characteristics. The comparison of clear call vocalizations ofMacaca radiata andM. fuscata revealed significant species-specific differences on the structural and the behavioural level. Evaluations of the structural features of alarm calls of various macaque species imply marked differences between members of thefascicularis group andsinica group on one hand and thesilenus group andarctoides group on the other.