The variation and evolution of complete human centromeres.

Human centromeres have been traditionally very difficult to sequence and assemble owing to their repetitive nature and large size1. As a result, patterns of human centromeric variation and models for their evolution and function remain incomplete, despite centromeres being among the most rapidly mutating regions2,3. Here, using long-read sequencing, we completely sequenced and assembled all centromeres from a second human genome and compared it to the finished reference genome4,5. We find that the two sets of centromeres show at least a 4.1-fold increase in single-nucleotide variation when compared with their unique flanks and vary up to 3-fold in size. Moreover, we find that 45.8% of centromeric sequence cannot be reliably aligned using standard methods owing to the emergence of new α-satellite higher-order repeats (HORs). DNA methylation and CENP-A chromatin immunoprecipitation experiments show that 26% of the centromeres differ in their kinetochore position by >500 kb. To understand evolutionary change, we selected six chromosomes and sequenced and assembled 31 orthologous centromeres from the common chimpanzee, orangutan and macaque genomes. Comparative analyses reveal a nearly complete turnover of α-satellite HORs, with characteristic idiosyncratic changes in α-satellite HORs for each species. Phylogenetic reconstruction of human haplotypes supports limited to no recombination between the short (p) and long (q) arms across centromeres and reveals that novel α-satellite HORs share a monophyletic origin, providing a strategy to estimate the rate of saltatory amplification and mutation of human centromeric DNA.

Comparative genomics of macaques and integrated insights into genetic variation and population history.

The crab-eating macaques ( Macaca fascicularis ) and rhesus macaques ( M. mulatta ) are widely studied nonhuman primates in biomedical and evolutionary research. Despite their significance, the current understanding of the complex genomic structure in macaques and the differences between species requires substantial improvement. Here, we present a complete genome assembly of a crab-eating macaque and 20 haplotype-resolved macaque assemblies to investigate the complex regions and major genomic differences between species. Segmental duplication in macaques is ∼42% lower, while centromeres are ∼3.7 times longer than those in humans. The characterization of ∼2 Mbp fixed genetic variants and ∼240 Mbp complex loci highlights potential associations with metabolic differences between the two macaque species (e.g., CYP2C76 and EHBP1L1 ). Additionally, hundreds of alternative splicing differences show post-transcriptional regulation divergence between these two species (e.g., PNPO ). We also characterize 91 large-scale genomic differences between macaques and humans at a single-base-pair resolution and highlight their impact on gene regulation in primate evolution (e.g., FOLH1 and PIEZO2 ). Finally, population genetics recapitulates macaque speciation and selective sweeps, highlighting potential genetic basis of reproduction and tail phenotype differences (e.g., STAB1 , SEMA3F , and HOXD13 ). In summary, the integrated analysis of genetic variation and population genetics in macaques greatly enhances our comprehension of lineage-specific phenotypes, adaptation, and primate evolution, thereby improving their biomedical applications in human diseases.

On the Way to the Technological Development of Newly Selected Non-Saccharomyces Yeasts Selected as Innovative Biocontrol Agents in Table Grapes.

Post-harvest decay of fresh table grapes causes considerable annual production losses. The main fungal agents of decay both in pre- and post-harvest are B. cinerea, Penicillium spp., Aspergillus spp., Alternaria spp., and Cladosporium spp. To date, the use of agrochemicals and SO2 are the main methods to control grape molds in pre- and postharvest, respectively. Significant improvements, however, have already been made in to apply innovative and more environmentally sustainable control strategies, such as Biological Control Agents (BCAs), which can reduce disease severity in both pre- and post-harvest. In this study, 31 new non-Saccharomyces yeast strains, isolated from berries of native Apulian table grape genotypes, were tested for their in vivo effectiveness against grey mold of table grapes, resulting in two St. bacillaris ('N22_I1' and 'S13_I3'), one S. diversa ('N22_I3'), one A. pullulans ('OLB_9.1_VL') and one H. uvarum ('OLB_9.1_BR') yeast strains that were marked as efficient and good BCAs. Their mechanisms of action were characterized through in vitro assays, and additional characteristics were evaluated to assess the economic feasibility and viability for future technological employment. Their effectiveness was tested by reducing the working concentration, their antagonistic effect on a wide range of fungal pathogens, their ability to survive in formulations with long shelf life, and their safety to human health.

A Rare Case of Concurrent 2q34q36 Duplication and 2q37 Deletion in a Neonate with Syndromic Features.

Large-scale genomic structural variations can have significant clinical implications, depending on the specific altered genomic region. Briefly, 2q37 microdeletion syndrome is a prevalent subtelomeric deletion disorder characterized by variable-sized deletions. Affected patients exhibit a wide range of clinical manifestations, including short stature, facial dysmorphism, and features of autism spectrum disorder, among others. Conversely, isolated duplications of proximal chromosome 2q are rare and lack a distinct phenotype. In this report, we provide an extensive molecular analysis of a 15-day-old newborn referred for syndromic features. Our analysis reveals an 8.5 Mb microdeletion at 2q37.1, which extends to the telomere, in conjunction with an 8.6 Mb interstitial microduplication at 2q34q36.1. Our findings underscore the prominence of 2q37 terminal deletions as commonly reported genomic anomalies. We compare our patient's phenotype with previously reported cases in the literature to contribute to a more refined classification of 2q37 microdeletion syndrome and assess the potential impact of 2q34q36.1 microduplication. We also investigate multiple hypotheses to clarify the genetic mechanisms responsible for the observed genomic rearrangement.

Structural Variation Evolution at the 15q11-q13 Disease-Associated Locus.

The impact of segmental duplications on human evolution and disease is only just starting to unfold, thanks to advancements in sequencing technologies that allow for their discovery and precise genotyping. The 15q11-q13 locus is a hotspot of recurrent copy number variation associated with Prader-Willi/Angelman syndromes, developmental delay, autism, and epilepsy and is mediated by complex segmental duplications, many of which arose recently during evolution. To gain insight into the instability of this region, we characterized its architecture in human and nonhuman primates, reconstructing the evolutionary history of five different inversions that rearranged the region in different species primarily by accumulation of segmental duplications. Comparative analysis of human and nonhuman primate duplication structures suggests a human-specific gain of directly oriented duplications in the regions flanking the GOLGA cores and HERC segmental duplications, representing potential genomic drivers for the human-specific expansions. The increasing complexity of segmental duplication organization over the course of evolution underlies its association with human susceptibility to recurrent disease-associated rearrangements.

Occurrence of L1M Elements in Chromosomal Rearrangements Associated to Chronic Myeloid Leukemia (CML): Insights from Patient-Specific Breakpoints Characterization.

Chronic myeloid leukemia (CML) is a rare myeloproliferative disorder caused by the reciprocal translocation t(9;22)(q34;q11) in hematopoietic stem cells (HSCs). This chromosomal translocation results in the formation of an extra-short chromosome 22, called a Philadelphia chromosome (Ph), containing the BCR-ABL1 fusion gene responsible for the expression of a constitutively active tyrosine kinase that causes uncontrolled growth and replication of leukemic cells. Mechanisms behind the formation of this chromosomal rearrangement are not well known, even if, as observed in tumors, repetitive DNA may be involved as core elements in chromosomal rearrangements. We have participated in the explorative investigations of the PhilosoPhi34 study to evaluate residual Ph+ cells in patients with negative FISH analysis on CD34+/lin- cells with gDNA qPCR. Using targeted next-generation deep sequencing strategies, we analyzed the genomic region around the t(9;22) translocations of 82 CML patients and one CML cell line and assessed the relevance of interspersed repeat elements at breakpoints (BP). We found a statistically higher presence of LINE elements, in particular belonging to the subfamily L1M, in BP cluster regions of both chromosome 22 and 9 compared to the whole human genome. These data suggest that L1M elements could be potential drivers of t(9;22) translocation leading to the generation of the BCR-ABL1 chimeric gene and the expression of the active BCR-ABL1-controlled tyrosine kinase chimeric protein responsible for CML.

A systematic mapping study on machine learning techniques for the prediction of CRISPR/Cas9 sgRNA target cleavage.

CRISPR/Cas9 technology has greatly accelerated genome engineering research. The CRISPR/Cas9 complex, a bacterial immune response system, is widely adopted for RNA-driven targeted genome editing. The systematic mapping study presented in this paper examines the literature on machine learning (ML) techniques employed in the prediction of CRISPR/Cas9 sgRNA on/off-target cleavage, focusing on improving support in sgRNA design activities and identifying areas currently being researched. This area of research has greatly expanded recently, and we found it appropriate to work on a Systematic Mapping Study (SMS), an investigation that has proven to be an effective secondary study method. Unlike a classic review, in an SMS, no comparison of methods or results is made, while this task can instead be the subject of a systematic literature review that chooses one theme among those highlighted in this SMS. The study is illustrated in this paper. To the best of the authors' knowledge, no other SMS studies have been published on this topic. Fifty-seven papers published in the period 2017-2022 (April, 30) were analyzed. This study reveals that the most widely used ML model is the convolutional neural network (CNN), followed by the feedforward neural network (FNN), while the use of other models is marginal. Other interesting information has emerged, such as the wide availability of both open code and platforms dedicated to supporting the activity of researchers or the fact that there is a clear prevalence of public funds that finance research on this topic.

A high-resolution map of small-scale inversions in the gibbon genome.

Gibbons are the most speciose family of living apes, characterized by a diverse chromosome number and rapid rate of large-scale rearrangements. Here we performed single-cell template strand sequencing (Strand-seq), molecular cytogenetics, and deep in silico analysis of a southern white-cheeked gibbon genome, providing the first comprehensive map of 238 previously hidden small-scale inversions. We determined that more than half are gibbon specific, at least fivefold higher than shown for other primate lineage-specific inversions, with a significantly high number of small heterozygous inversions, suggesting that accelerated evolution of inversions may have played a role in the high sympatric diversity of gibbons. Although the precise mechanisms underlying these inversions are not yet understood, it is clear that segmental duplication-mediated NAHR only accounts for a small fraction of events. Several genomic features, including gene density and repeat (e.g., LINE-1) content, might render these regions more break-prone and susceptible to inversion formation. In the attempt to characterize interspecific variation between southern and northern white-cheeked gibbons, we identify several large assembly errors in the current GGSC Nleu3.0/nomLeu3 reference genome comprising more than 49 megabases of DNA. Finally, we provide a list of 182 candidate genes potentially involved in gibbon diversification and speciation.

Assessing temporal and geographic contacts across the Adriatic Sea through the analysis of genome-wide data from Southern Italy.

Southern Italy was characterised by a complex prehistory that started with different Palaeolithic cultures, later followed by the Neolithization and the demic dispersal from the Pontic-Caspian Steppe during the Bronze Age. Archaeological and historical evidences point to a link between Southern Italians and the Balkans still present in modern times. To shed light on these dynamics, we analysed around 700 South Mediterranean genomes combined with informative ancient DNAs. Our findings revealed high affinities of South-Eastern Italians with modern Eastern Peloponnesians, and a closer affinity of ancient Greek genomes with those from specific regions of South Italy than modern Greek genomes. The higher similarity could be associated with a Bronze Age component ultimately originating from the Caucasus with high Iranian and Anatolian Neolithic ancestries. Furthermore, extremely differentiated allele frequencies among Northern and Southern Italy revealed putatively adapted SNPs in genes involved in alcohol metabolism, nevi features and immunological traits.

Segmental duplications and their variation in a complete human genome.

Despite their importance in disease and evolution, highly identical segmental duplications (SDs) are among the last regions of the human reference genome (GRCh38) to be fully sequenced. Using a complete telomere-to-telomere human genome (T2T-CHM13), we present a comprehensive view of human SD organization. SDs account for nearly one-third of the additional sequence, increasing the genome-wide estimate from 5.4 to 7.0% [218 million base pairs (Mbp)]. An analysis of 268 human genomes shows that 91% of the previously unresolved T2T-CHM13 SD sequence (68.3 Mbp) better represents human copy number variation. Comparing long-read assemblies from human (n = 12) and nonhuman primate (n = 5) genomes, we systematically reconstruct the evolution and structural haplotype diversity of biomedically relevant and duplicated genes. This analysis reveals patterns of structural heterozygosity and evolutionary differences in SD organization between humans and other primates.

Genomic and Personalized Medicine Approaches for Substance Use Disorders (SUDs) Looking at Genome-Wide Association Studies.

Drug addiction, or substance use disorder (SUD), is a chronic, relapsing disorder in which compulsive drug-seeking and drug-taking behaviour persist despite serious negative consequences. Drug abuse represents a problem that deserves great attention from a social point of view, and focuses on the importance of genetic studies to help in understanding the genetic basis of addiction and its medical treatment. Despite the complexity of drug addiction disorders, and the high number of environmental variables playing a role in the onset, recurrence, and duration of the symptoms, several studies have highlighted the non-negligible role of genetics, as demonstrated by heritability and genome-wide association studies. A correlation between the relative risk of addiction to specific substances and heritability has been recently observed, suggesting that neurobiological mechanisms may be, at least in part, inherited. All these observations point towards a scenario where the core neurobiological factors of addiction, involving the reward system, impulsivity, compulsivity, stress, and anxiety response, are transmitted, and therefore, genes and mutations underlying their variation might be detected. In the last few years, the development of new and more efficient sequencing technologies has paved the way for large-scale studies in searching for genetic and epigenetic factors affecting drug addiction disorders and their treatments. These studies have been crucial to pinpoint single nucleotide polymorphisms (SNPs) in genes that affect the reaction to medical treatments. This is critically important to identify pharmacogenomic approaches for substance use disorder, such as OPRM1 SNPs and methadone required doses for maintenance treatment (MMT). Nevertheless, despite the promising results obtained by genome-wide association and pharmacogenomic studies, specific studies related to population genetics diversity are lacking, undermining the overall applicability of the preliminary findings, and thus potentially affecting the portability and the accuracy of the genetic studies. In this review, focusing on cannabis, cocaine and heroin use, we report the state-of-the-art genomics and pharmacogenomics of SUDs, and the possible future perspectives related to medical treatment response in people that ask for assistance in solving drug-related problems.

Evaluating the Impact of Sex-Biased Genetic Admixture in the Americas through the Analysis of Haplotype Data.

A general imbalance in the proportion of disembarked males and females in the Americas has been documented during the Trans-Atlantic Slave Trade and the Colonial Era and, although less prominent, more recently. This imbalance may have left a signature on the genomes of modern-day populations characterised by high levels of admixture. The analysis of the uniparental systems and the evaluation of continental proportion ratio of autosomal and X chromosomes revealed a general sex imbalance towards males for European and females for African and Indigenous American ancestries. However, the consistency and degree of this imbalance are variable, suggesting that other factors, such as cultural and social practices, may have played a role in shaping it. Moreover, very few investigations have evaluated the sex imbalance using haplotype data, containing more critical information than genotypes. Here, we analysed genome-wide data for more than 5000 admixed American individuals to assess the presence, direction and magnitude of sex-biased admixture in the Americas. For this purpose, we applied two haplotype-based approaches, ELAI and NNLS, and we compared them with a genotype-based method, ADMIXTURE. In doing so, besides a general agreement between methods, we unravelled that the post-colonial admixture dynamics show higher complexity than previously described.

Continental-scale genomic analysis suggests shared post-admixture adaptation in the Americas.

American populations are one of the most interesting examples of recently admixed groups, where ancestral components from three major continental human groups (Africans, Eurasians and Native Americans) have admixed within the last 15 generations. Recently, several genetic surveys focusing on thousands of individuals shed light on the geography, chronology and relevance of these events. However, even though gene flow could drive adaptive evolution, it is unclear whether and how natural selection acted on the resulting genetic variation in the Americas. In this study, we analysed the patterns of local ancestry of genomic fragments in genome-wide data for ~ 6000 admixed individuals from 10 American countries. In doing so, we identified regions characterized by a divergent ancestry profile (DAP), in which a significant over or under ancestral representation is evident. Our results highlighted a series of genomic regions with DAPs associated with immune system response and relevant medical traits, with the longest DAP region encompassing the human leukocyte antigen locus. Furthermore, we found that DAP regions are enriched in genes linked to cancer-related traits and autoimmune diseases. Then, analysing the biological impact of these regions, we showed that natural selection could have acted preferentially towards variants located in coding and non-coding transcripts and characterized by a high deleteriousness score. Taken together, our analyses suggest that shared patterns of post admixture adaptation occurred at a continental scale in the Americas, affecting more often functional and impactful genomic variants.

Special Issue: A Tale of Genes and Genomes.

Variability is the source on which selective pressure acts, allowing genome evolution and adaptation [...].

A high-quality bonobo genome refines the analysis of hominid evolution.

The divergence of chimpanzee and bonobo provides one of the few examples of recent hominid speciation1,2. Here we describe a fully annotated, high-quality bonobo genome assembly, which was constructed without guidance from reference genomes by applying a multiplatform genomics approach. We generate a bonobo genome assembly in which more than 98% of genes are completely annotated and 99% of the gaps are closed, including the resolution of about half of the segmental duplications and almost all of the full-length mobile elements. We compare the bonobo genome to those of other great apes1,3-5 and identify more than 5,569 fixed structural variants that specifically distinguish the bonobo and chimpanzee lineages. We focus on genes that have been lost, changed in structure or expanded in the last few million years of bonobo evolution. We produce a high-resolution map of incomplete lineage sorting and estimate that around 5.1% of the human genome is genetically closer to chimpanzee or bonobo and that more than 36.5% of the genome shows incomplete lineage sorting if we consider a deeper phylogeny including gorilla and orangutan. We also show that 26% of the segments of incomplete lineage sorting between human and chimpanzee or human and bonobo are non-randomly distributed and that genes within these clustered segments show significant excess of amino acid replacement compared to the rest of the genome.

Mayer-Rokitansky-Küster-Hauser syndrome with 22q11.21 microduplication: a case report.

BACKGROUND: Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome (Online Mendelian Inheritance in Man [OMIM] #277000) is a congenital condition characterized by the total or partial agenesis of vagina and uterus. Agenesis can be isolated (MRKH 1) or associated with other renal, vertebral or cardiac defects (MRKH 2). CASE PRESENTATION: In this paper, we report a case of a Caucasian patient showing the clinical signs associated with MRKH. Array-based comparative genomic hybridization (a-CGH) analysis revealed a microduplication of approximately 3.01 megabases (Mb) located on the long arm of chromosome 22 (22q11.21). Microduplications affecting the 22q11.21 region have been shown to be associated with MRKH syndrome and Müllerian aplasia. The phenotype of patients with 22q11.2 duplication (OMIM #608363) appears extremely variable, ranging from apparently normal to mild learning difficulties or with multiple defects, sharing features with DiGeorge/velocardiofacial (DGS/VCFS) syndrome. CONCLUSIONS: The altered gene expression together with other genetic, nongenetic, epigenetic or environmental factors can cause the extremely variable phenotype in patients carrying such duplication. Therefore, we can consider MRKH syndrome to be one of the clinical features of DGS/VCFS syndrome.

The structure, function and evolution of a complete human chromosome 8.

The complete assembly of each human chromosome is essential for understanding human biology and evolution1,2. Here we use complementary long-read sequencing technologies to complete the linear assembly of human chromosome 8. Our assembly resolves the sequence of five previously long-standing gaps, including a 2.08-Mb centromeric α-satellite array, a 644-kb copy number polymorphism in the ß-defensin gene cluster that is important for disease risk, and an 863-kb variable number tandem repeat at chromosome 8q21.2 that can function as a neocentromere. We show that the centromeric α-satellite array is generally methylated except for a 73-kb hypomethylated region of diverse higher-order α-satellites enriched with CENP-A nucleosomes, consistent with the location of the kinetochore. In addition, we confirm the overall organization and methylation pattern of the centromere in a diploid human genome. Using a dual long-read sequencing approach, we complete high-quality draft assemblies of the orthologous centromere from chromosome 8 in chimpanzee, orangutan and macaque to reconstruct its evolutionary history. Comparative and phylogenetic analyses show that the higher-order α-satellite structure evolved in the great ape ancestor with a layered symmetry, in which more ancient higher-order repeats locate peripherally to monomeric α-satellites. We estimate that the mutation rate of centromeric satellite DNA is accelerated by more than 2.2-fold compared to the unique portions of the genome, and this acceleration extends into the flanking sequence.

Native Vineyard Non-Saccharomyces Yeasts Used for Biological Control of Botrytis cinerea in Stored Table Grape.

Postharvest spoilage fungi, such as Botrytis cinerea, are considered the main cause of losses of fresh fruit quality and vegetables during storage, distribution, and consumption. The current control strategy is the use of SO2 generator pads whose application is now largely under observation. A high quantity of SO2 can be deleterious for fresh fruits and vegetables and it is not allowed in organic agriculture. For this reason, great attention has been recently focused on identifying Biological Control Agents (BCA) to implement biological approaches devoid of chemicals. In this direction, we carried out our study in isolating five different non-Saccharomyces yeast strains from local vineyards in the South of Italy as possible BCA. We performed both in vitro and in vivo assays in semi-commercial conditions on detached grape berries stored at 0 °C, simulating the temperature normally used during cold storage, and obtained relevant results. We isolated three M. pulcherrima strains and one L. thermotolerans strain able to largely antagonize the development of the B. cinerea, at both in vitro and in vivo conditions. In particular, we detected the ability of the three isolates of M. pulcherrima strains Ale4, N20/006, and Pr7 and the L. thermotolerans strain N10 to completely inhibit (100% in reduction) the mycelial growth of B. cinerea by producing fungistatic compounds. We found, using an extracellular lytic enzymes activity assay, that such activity could be related to lipid hydrolyzation, ß-1,3-glucanase and pectinase activity, and pectinase and protease activity, depending on the yeasts used. Results from our in vitro assays allowed us to hypothesize for M. pulcherrima strains Ale4 and N20/006 a possible combination of both the production of soluble metabolites and volatile organic compounds to antagonize against B. cinerea growth. Moreover, in semi-commercial conditions, the M. pulcherrima strain N20/006 and L. thermotolerans strain N10 showed relevant antagonistic effect also at low concentrations (with a significantly reduction of 'slip skin' incidence of 86.4% and 72.7%, respectively), thus highlighting a peculiar property to use in commercial development for organic agriculture and the handling process.

Color Stabilization of Apulian Red Wines through the Sequential Inoculation of Starmerella bacillaris and Saccharomyces cerevisiae.

Mixed fermentation using Starmerella bacillaris and Saccharomyces cerevisiae has gained attention in recent years due to their ability to modulate the qualitative parameters of enological interest, such as the color intensity and stability of wine. In this study, three of the most important red Apulian varieties were fermented through two pure inoculations of Saccharomyces cerevisiae strains or the sequential inoculation of Saccharomyces cerevisiae after 48 h from Starmerella bacillaris. The evolution of anthocyanin profiles and chromatic characteristics were determined in the produced wines at draining off and after 18 months of bottle aging in order to assess the impact of the different fermentation protocols on the potential color stabilization and shelf-life. The chemical composition analysis showed titratable acidity and ethanol content exhibiting marked differences among wines after fermentation and aging. The 48 h inoculation delay produced wines with higher values of color intensity and color stability. This was ascribed to the increased presence of compounds, such as stable A-type vitisins and reddish/violet ethylidene-bridge flavonol-anthocyanin adducts, in the mixed fermentation. Our results proved that the sequential fermentation of Starmerella bacillaris and Saccharomyces cerevisiae could enhance the chromatic profile as well as the stability of the red wines, thus improving their organoleptic quality.

Single-cell strand sequencing of a macaque genome reveals multiple nested inversions and breakpoint reuse during primate evolution.

Rhesus macaque is an Old World monkey that shared a common ancestor with human ∼25 Myr ago and is an important animal model for human disease studies. A deep understanding of its genetics is therefore required for both biomedical and evolutionary studies. Among structural variants, inversions represent a driving force in speciation and play an important role in disease predisposition. Here we generated a genome-wide map of inversions between human and macaque, combining single-cell strand sequencing with cytogenetics. We identified 375 total inversions between 859 bp and 92 Mbp, increasing by eightfold the number of previously reported inversions. Among these, 19 inversions flanked by segmental duplications overlap with recurrent copy number variants associated with neurocognitive disorders. Evolutionary analyses show that in 17 out of 19 cases, the Hominidae orientation of these disease-associated regions is always derived. This suggests that duplicated sequences likely played a fundamental role in generating inversions in humans and great apes, creating architectures that nowadays predispose these regions to disease-associated genetic instability. Finally, we identified 861 genes mapping at 156 inversions breakpoints, with some showing evidence of differential expression in human and macaque cell lines, thus highlighting candidates that might have contributed to the evolution of species-specific features. This study depicts the most accurate fine-scale map of inversions between human and macaque using a two-pronged integrative approach, such as single-cell strand sequencing and cytogenetics, and represents a valuable resource toward understanding of the biology and evolution of primate species.