Polymorphic Sequence Research Articles

Genetic Polymorphism and Differentiation of Populations of Sterlet Acipenser ruthenus (Acipenseridae) in the Lower Irtysh and Middle Ob Basins

The article presents data on polymorphism of intermicrosatellite sequences in the sterlet Acipenser ruthenus of the lower reaches of the Irtysh River and the middle reaches of the Ob River. We assessed intra- and interpopulation variability and genetic differentiation of A. ruthenus and revealed a high ISSR polymorphism in the species from the central part of the Ob-Irtysh basin. The proportion of polymorphic amplicons was 0.966, genetic diversity was 0.355, and the average number of alleles per locus was 1.97. The highest polymorphism was typical for the sterlet from the Tobol River at the confluence with the Irtysh River. Genetic differentiation between the sterlet groups of the Irtysh and Ob rivers is well pronounced, the interpopulation component accounts for 42% of variability (Gst = 0.42), gene flow is limited (Nm = 0.67). The sterlet groups inhabiting the Lower Irtysh from the mouth of the Tobol River to the mouth of the Konda River do not differ genetically and form one population (Gst = 0.08–0.12, Nm = 3.76–5.55). The sterlet from the Irtysh River within the Vagay region is genetically different from the other Irtysh samples (Gst = 0.22, Nm = 1.68) and belongs to a different population group. The differentiation between samples of sterlet from the Ob basin is higher than between samples from the Irtysh basin. Groups of sterlet from the Ob River and the Yuganskaya Ob canal are genetically different (Gst = 0.30, Nm = 1.19) and form various subpopulations. Spawning migrations, as well as confinement to wintering pits, play a decisive role in the formation of the sterlet population structure in the studied part of the distribution area. The identified sterlet population groups should be considered as separate units of environmental and economic management.

Kazakhstan tulips: comparative analysis of complete chloroplast genomes of four local and endangered species of the genus Tulipa L.

Species of Tulipa are important ornamental plants used for horticultural purposes in various countries, across Asia, Europe, and North Africa. The present study is the first report on typical features of the complete chloroplast genome sequence of four local and endangered species including T. alberti, T. kaufmanniana, T. greigii, and T. dubia from Kazakhstan using Illumina sequencing technology. The comparative analyses revealed that the complete genomes of four species were highly conserved in terms of total genome size (152. 006 bp - 152. 382 bp), including a pair of inverted repeat regions (26. 330 bp - 26. 371 bp), separated by a large single copy region (82.169 bp - 82,378 bp) and a small copy region (17.172 bp -17.260 bp). Total GC content (36.58-36.62 %), gene number (131), and intron length (540 bp - 2620 bp) of 28 genes. The complete genomes of four species showed nucleotide diversity (π =0,003257). The total number of SSR loci was 159 in T. alberti, 158 in T. kaufmanniana, 174 in T. greigii, and 163 in T. dubia. The result indicated that ten CDS genes, namely rpoC2, cemA, rbcL, rpl36, psbH, rps3, rpl22, ndhF, ycf1, and matK, with effective polymorphic simple sequence repeats (SSRs), high sequence variability (SV) ranging from 2.581 to 6.102, and high nucleotide diversity (Pi) of these loci ranging from 0,004 to 0,010. For all intergenic regions longer than 150 bp, twenty one most variable regions were found with high sequence variability (SV) ranging from 4,848 to 11,862 and high nucleotide diversity (Pi) ranging from 0,01599 to 0,01839. Relative synonymous codon usage (RSCU) analysis was used to identify overrepresented and underrepresented codons for each amino acid. Based on the phylogenic analysis, the sequences clustered into four major groups, reflecting distinct evolutionary lineages corresponding to the subgenera Eriostemons, Tulipa, and Orithyia. Notably, T. greigii was distinctively grouped with species from Orithyia and Eriostemons rather than with other Tulipa species, suggesting a unique evolutionary history potentially shaped by geographical isolation or specific ecological pressures. The complete chloroplast genome of the four Tulipa species provides fundamental information for future research studies, even for designing the high number of available molecular markers.

Genotyping of human isolates from human toxoplasmosis outbreak: Restriction Fragment Length Polymorphism and Sanger Sequencing.

Toxoplasmosis is a parasitic disease caused by infection with Toxoplasma gondii, capable of infecting a wide range of hosts. The parasite exhibits a broad genetic diversity, necessitating genotypic characterization for genotype identification and associations with epidemiological information. Therefore, the Restriction Fragment Length Polymorphism (RFLP) technique is used for characterization. This study aimed to perform genotypic characterization of isolates from pregnant women infected during a human toxoplasmosis outbreak, comparing the RFLP and Sanger Sequencing methodologies. For this purpose, six human isolates were subjected to conventional PCR, Multiplex PCR, Nested PCR, Enzymatic Digestion, and Sanger Sequencing. Additionally, the standard strains GTI (Type I), PTG (Type II), and CTG (Type III) were also subjected to the same techniques described above. Subsequently, the amplified DNA products were compared with the standard strains. As a result, it was observed that Sanger Sequencing provides the same information as RFLP PCR, as well as the possibility of cost reduction for genotypic characterization, and providing greater agility in issuing results. Additionally, Sanger Sequencing of T. gondii isolates allows for detailed evaluation of nucleotide sequences, including the assessment of SNPs and enzymatic restriction sites, which the RFLP technique does not.

Fine mapping of QGPC.caas-7AL for grain protein content in bread wheat.

A major stable QTL, QGPC.caas-7AL, for grain protein content of wheat, was narrowed down to a 1.82-Mb inter on chromosome 7AL, and four candidate genes were predicated. Wheat grain protein content (GPC) is important for end-use quality. Identification of genetic loci for GPC is helpful to create new varieties with good processing quality and nutrients. Zhongmai 578 (ZM578) and Jimai 22 (JM22) are two elite wheat varieties with different contents of GPC. In the present study, 262 recombinant inbred lines (RILs) derived from a cross between ZM578 and JM22 were used to map the GPC with high-density wheat Illumina iSelect 50K single-nucleotide polymorphism (SNP) array. Seven quantitative trait loci (QTLs) were identified for GPC on chromosomes 3AS, 3AL, 3BS, 4AL, 5BS, 5DL and 7AL by inclusive composite interval mapping, designated as QGPC.caas-3AS, QGPC.caas-3AL, QGPC.caas-3BS, QGPC.caas-4AL, QGPC.caas-5BS, QGPC.caas-5DL and QGPC.caas-7AL, respectively. Among these, alleles for increasing GPC at QGPC.caas-3AS, QGPC.caas-3BS, QGPC.caas-4AL and QGPC.caas-7AL loci were contributed by ZM578, whereas those at the other three loci were from JM22. The stable QTL QGPC.caas-7AL was fine mapped to a 1.82-Mb physical interval using secondary populations from six heterozygous recombinant plants obtained by selfing a residual RIL. Four genes were predicted as candidates of QGPC.caas-7AL based on sequence polymorphism and expression patterns. The near-isogenic lines (NILs) with the favorable allele at the QGPC.caas-7AL locus increased Farinograph stability time, Extensograph extension area, extensibility and maximum resistance by 19.6%, 6.3%, 6.0% and 20.3%, respectively. Kompetitive allele-specific PCR (KASP) marker for QGPC.caas-7AL was developed and validated in a diverse panel of 166 Chinese wheat cultivars. These results provide further insight into the genetic basis of GPC, and the fine-mapped QGPC.caas-7AL will be an attractive target for map-based cloning and marker-assisted selection in wheat breeding programs.

Dynamic atlas of histone modifications and gene regulatory networks in endosperm of bread wheat.

Dissecting the genetic basis of seed traits in wheat is impeded by limited genetic polymorphisms and significant variations caused by environmental conditions and seed position in a spikelet. Seed performance is largely determined by endosperm development controlled by spatiotemporal variation in gene activities, which is greatly affected by chromatin status. Here, we map genome-wide dynamic distributions of H3K27me3, H3K4me3 and H3K9ac modifications and profile gene transcription across wheat endosperm development. The combinatorial effects of active and repressive marks ensure spatiotemporal dynamic gene expression, especially for starch biosynthesis. By scanning the transcription factor binding motifs in the ATAC-seq peaks, hub regulators are identified from the regulatory network. In addition, we observe significant correlations between sequence polymorphisms of hub regulators and variations in seed traits in a germplasm population. Thus, the analysis of genomic regulatory activities together with genetic variation provides a robust approach to dissect seed traits in bread wheat.

Cleaved Amplified Polymorphic Sequence Markers in Horticultural Crops: Current Status and Future Perspectives

DNA markers have broad applications, including marker-assisted selection (MAS) for breeding new cultivars. Currently, single nucleotide polymorphisms (SNPs) have become a preferred choice of markers for molecular geneticists and breeders. They offer many advantages, such as high abundance and coverage in the genome, codominant inheritance, locus specificity, and flexibility for high-throughput genotyping/detection formats, and they are relatively inexpensive. The availability of reference genome sequences enables precise identification of candidate genes and SNPs associated with a trait of interest through quantitative trait loci (QTL) mapping and genome-wide association studies (GWAS). Such SNPs can be converted into markers for their application in MAS in crop breeding programs. Cleaved amplified polymorphic sequence (CAPS) markers amplify short genomic sequences around the polymorphic endonuclease restriction site. This review provides insight into the recent advancements made in the development and application of CAPS markers in several horticultural plants. We discussed many new tools that aid faster and more accurate design of CAPS markers from the whole genome resequencing (WGRS) data. The developed CAPS markers offer immense application in germplasm screening and field trials, genomic loci mapping, identifying candidate genes, and MAS of important horticultural traits such as disease resistance, fruit quality and morphology, and genetic purity.

Impacts of Matrix Metalloproteinase-11 Genotypes on Pterygium Risk.

The dysregulation of matrix metalloproteinase (MMP) proteins has been reported to be involved in the etiology of pterygium. However, studies about the role of matrix metalloproteinase-11 (MMP-11) are lacking. This study is the first to examine the genomic role of MMP-11 in pterygium. The genotypes of MMP-11 rs738791, rs2267029, rs738792, and rs28382575 were determined in 140 pterygium cases and 280 non-pterygium controls by utilizing polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) and direct sequencing. The genotypic frequencies of MMP-11 rs738792 TT, CT and CC were 40.0%, 50.7%, and 9.3% in the pterygium group, significantly different from those in the non-pterygium group (55.0%, 37.1%, and 7.9%, respectively; p for trend=0.0139). Specifically, individuals carrying the variants CT and CC had a 1.88- and 1.63-fold odds ratio of pterygium risk [95% confidence interval (CI)=1.22-2.89 and 0.77-3.44, p=0.0054 and 0.2834, respectively]. In the dominant model, individuals carrying CT+CC had significantly higher pterygium risk (odds ratio=1.83, 95%CI=1.21-2.77, p=0.0052). No association was found for other MMP-11 polymorphisms. Allelic analysis showed that MMP-11 rs738792 C allele was significantly associated with pterygium risk (odds ratio=1.48, 95%CI=1.08-2.01, p=0.0169). for the variant alleles of other MMP-11 polymorphisms were not associated with pterygium risk. The MMP-11 rs738792 genotypes can serve as a predictive marker for pterygium risk in Taiwanese. Additionally, elucidating the role of MMP-11 in the pathogenesis of pterygium could inform targeted therapies based on MMP-11 modulation.

Detection of Two Common ACCase Mutations Associated with High Levels of Fenoxaprop-P-Ethyl Resistance in Shortawn Foxtail (Alopecurus aequalis) Using Loop-Mediated Isothermal Amplification

Abstract The resistance to fenoxaprop-P-ethyl, a herbicide that inhibits acetyl-CoA carboxylase (ACCase), has emerged in shortawn foxtail (Alopecurus aequalis Sobol.) since the 1990s, presenting a considerable challenge to wheat (Triticum aestivum L.) production in China. One of the primary mechanisms responsible for this high-level resistance is the presence of mutations at codons 1781, 2041, and 2078 in the ACCase gene. However, the conventional methods used to detect these mutations, such as polymerase chain reaction (PCR) and gene sequencing, are time-consuming and labor-intensive. To address this issue and enable the prompt and effective detection of these common ACCase mutations in A. aequalis, a loop-mediated isothermal amplification (LAMP) strategy was developed. The LAMP assay specifically targets the Ile-1781-Leu and Asp-2078-Gly mutations within the ACCase gene of A. aequalis. Through the optimization of primers, systems, and conditions, the LAMP assay enables rapid differentiation between wild-type individuals and mutants of A. aequalis carrying either of these two mutations. By including SYBR Green I dye in the final reaction mixtures, the target mutation can be visually detected through a noticeable color change that can be observed with the naked eye. It is noteworthy that the sensitivity of the LAMP assay was approximately 104-fold greater than that of conventional PCR methods. Additionally, a derived cleaved amplified polymorphic sequence (dCAPS) assay was established for each mutation to distinguish between homozygous and heterozygous mutants. Overall, the developed LAMP assay could efficiently detect the Ile-1781-Leu and Asp-2078-Gly mutations in the ACCase gene of A. aequalis, offering significant advantages for the monitoring and management of fenoxaprop-P-ethyl resistance.

Active DNA Demethylation Mediated by OsGADD45a2 Regulates Growth, Development, and Blast (Magnaporthe oryzea) Resistance in Rice.

OsGADD45a1, a member of the growth arrest and DNA damage-inducible 45 (GADD45) family in rice, has a newly identified homologue, OsGADD45a2, which differs from OsGADD45a1 in only three amino acids. The role and function of the OsGADD45a2 in DNA demethylation are not well-understood and were investigated in this study. Osgadd45a2 mutants exhibited reduced height, shorter panicle length, fewer grains per panicle, and a lower seed setting rate compared with wild-type plants. Moreover, the results showed that OsGADD45a2 negatively regulates rice blast fungus resistance and exhibited high expression in various tissues. Using the 3000 Rice Genomes Project database, we identified four major haplotypes (each with over 100 cultivars) based on single-nucleotide polymorphisms in the coding sequence of OsGADD45a2. Among these, Hap4 was associated with a significantly greater plant height than Hap1-3, possibly due to a functional alteration of OsGADD45a2 linked to the SNP at position 2614993. In OsGADD45a2 overexpression lines, significant decreases in CG and CHG methylation levels were observed in protein-coding genes, leading to their upregulation. Overall, our findings indicate that OsGADD45a2 acts as a methylation regulator, mediating the expression of genes essential for plant growth and development and blast resistance.

Haplotype analysis of Tongil-type rice varieties through sequence polymorphism analysis of cold tolerance-enhancing genes

Haplotype analysis of Tongil-type rice varieties through sequence polymorphism analysis of cold tolerance-enhancing genes

Chloroplast genome sequencing and divergence analysis of 18 Pyrus species: insights into intron length polymorphisms and evolutionary processes.

Pears constitute an essential temperate crop and are primarily produced through interspecific hybridization owing to self-incompatibility that complicates their breeding history. To address this, we sequenced the complete chloroplast (cp) genomes of 18 Pyrus and one Malus species using the Illumina HiSeq4000 platform. The cp genomes ranged from 159,885 bp to 160,153bp and exhibited a conserved circular DNA structure with an average GC content of 36.5%. Each cp genome contained 127 genes, including 83 protein-coding, 36 tRNA, and 8 rRNA genes. Divergence analysis with mVISTA showed high conservation in the coding regions and notable variations in the non-coding regions. All species shared 17 intron-containing genes, with ycf3 and clpP each having two introns. Five intron-containing genes (ndhB, rpl2, rps12, trnA-UGC, and trnE-UUC) were located in the inverted repeat regions, while trnL-UAA was located in the large single-copy region, with conserved intron lengths across Pomoideae. We identified polymorphic intron sequences in the rpl22, petB, clpP, ndhA, and rps16 genes and designed primers for these regions. Notably, the two Pyrus ussuriensis accessions Doonggeullebae and Cheongdangrori showed intron-length polymorphisms despite being classified as the same species. Phylogenetic analysis of the cp genome sequences revealed two major clusters, indicating distinct maternal lineages and evolutionary origins. This study underscores the importance of cp gene polymorphisms in P. fauriei, P. calleryana, P. ussuriensis, and P. pyrifolia, providing valuable insights into Pyrus evolution as well as aiding in the conservation and breeding of pear germplasm.

Correction for Welsh et al., "Sequence polymorphisms in the reovirus σ1 attachment protein modulate encapsidation efficiency and replication in mice".

Correction for Welsh et al., "Sequence polymorphisms in the reovirus σ1 attachment protein modulate encapsidation efficiency and replication in mice".

Genome-wide discovery of polymorphic simple sequence repeat (SSR) markers through comparative analysis of closely related hami melon cultivars

Genome-wide discovery of polymorphic simple sequence repeat (SSR) markers through comparative analysis of closely related hami melon cultivars

Resolving Ambiguities in Deceased Donor HLA Typing: Next Generation Sequencing’s Role in DPB1 Allele Discrimination and Immunological Implications

Abstract HLA typing is crucial for assessing immunological compatibility between donors and recipients before transplantation. Deceased donor (DD) typing, typically via real-time PCR (rtPCR), provides rapid results vital for organ allocation. rtPCR primarily targets exons encoding the highly polymorphic peptide binding region (PBR) of HLA molecules. rtPCR cannot target all polymorphic sequences. Consequently, several possibilities of allele combinations are provided, leading to ambiguous typing results. As recent evidence suggests antibody responses can target regions in exons outside the PBR, notably in Exon 3 of the DPB1 locus, typing ambiguities can lead to the inability to assess histocompatibility between a donor and recipient. Herein, the aim was to investigate the utility of Next Generation Sequencing (NGS) in resolving DPB1 ambiguities in DD typing. A retrospective analysis of DD typing from February 2022 to January 2024 examined instances where NGS was used to resolve ambiguities reported via rtPCR. Each case was reviewed to understand the reason for resorting to NGS typing. The study focused on DPB1 ambiguities and assessed whether NGS effectively reduced ambiguity and differentiated alleles with polymorphisms outside the PBR. In total, 29 DD cases were found to have NGS testing results. In 14 cases, NGS was specifically performed to resolve DPB1 ambiguities. The ambiguities were significantly reduced from hundreds detected in rtPCR to 4 or fewer possible allele combinations with NGS. Eleven cases demonstrated polymorphism in exons outside the PBRs. Of these, only 2 cases showed no typing ambiguity with NGS. In the remaining 9 cases, though polymorphisms could be resolved in the PBR, ambiguities remained in exon 3, specifically at codon 96. In these cases, NGS could not distinguish lysine (K) from arginine (R). For example, DPB1*835:01 with the 96R motif could not be distinguished from DPB1*34:01with the 96K motif. The current investigation found that NGS cannot resolve all ambiguities at position 96, likely due, at least in part, to persistent cis-trans ambiguities. Our results suggest that although NGS drastically limited the number of ambiguities in 9 of 11 cases, it still could not adequately distinguish all immunologically relevant alleles. Inability to resolve the 96R/K epitope, which has been shown to be immunogenic, could impact solid organ transplantation, potentially leading to rejection in mismatched recipient-donor pairs at this position. Sequencing methods that generate longer reads, such as nanopore, may help resolve some of these cases if they are indeed due to cis-trans ambiguities.

Identification of sequence polymorphism in the D-loop region of mitochondrial DNA as a risk factor for breast cancer.

Mitochondrial DNA (mtDNA) variations affect the efficiency of the electron transport chain and production of reactive oxygen species, contributing to carcinogenesis. The D-loop region of mtDNA has emerged as a variation hotspot region in human neoplasia; however, the potential contribution of these variations in breast cancer risk prediction remains unknown. We investigated the relationship between germline single nucleotide polymorphisms (SNPs) in the entire D-loop region and breast cancer risk in Chinese women. Peripheral blood-isolated mtDNA from 2329 patients with breast cancer and 2328 cancer-free controls was examined for SNPs. In the combined cohort, we used traditional risk factors, susceptibility germline polymorphisms, and logistic regression analysis to evaluate the predictive value of susceptibility variants for breast cancer risk. We calculated the area under the receiver operating characteristic curve (AUC) as a measure. We also measured the content of 8-hydroxy-2'-deoxyguanosine (8-OHdG). Individual polymorphisms SNP573 were significantly associated with breast cancer risk in both the discovery and validation cohorts. In the combined cohort, the AUC of the traditional risk factors was 64.3%; after adding susceptibility variants, the AUC was 64.9% (DeLong test, p = 0.007). 8-OHdG levels were significantly higher in patients with breast cancer than in controls and higher in individuals with SNP573 than in those negative for this variation. Overall, oxidative stress might be associated with the risk of breast cancer, and SNP573 might be associated with oxidative stress. Our results indicate the risk potential of polymorphisms in the D-loop region in breast cancer in Southern China.

Evolutionary Modes of wtf Meiotic Driver Genes in Schizosaccharomyces pombe.

Killer meiotic drivers are a class of selfish genetic elements that bias inheritance in their favor by destroying meiotic progeny that do not carry them. How killer meiotic drivers evolve is not well understood. In the fission yeast, Schizosaccharomyces pombe, the largest gene family, known as the wtf genes, is a killer meiotic driver family that causes intraspecific hybrid sterility. Here, we investigate how wtf genes evolve using long-read-based genome assemblies of 31 distinct S. pombe natural isolates, which encompass the known genetic diversity of S. pombe. Our analysis, involving nearly 1,000 wtf genes in these isolates, yields a comprehensive portrayal of the intraspecific diversity of wtf genes. Leveraging single-nucleotide polymorphisms in adjacent unique sequences, we pinpoint wtf gene-containing loci that have recently undergone gene conversion events and infer their ancestral state. These events include the revival of wtf pseudogenes, lending support to the notion that gene conversion plays a role in preserving this gene family from extinction. Moreover, our investigation reveals that solo long terminal repeats of retrotransposons, frequently found near wtf genes, can act as recombination arms, influencing the upstream regulatory sequences of wtf genes. Additionally, our exploration of the outer boundaries of wtf genes uncovers a previously unrecognized type of directly oriented repeats flanking wtf genes. These repeats may have facilitated the early expansion of the wtf gene family in S. pombe. Our findings enhance the understanding of the mechanisms influencing the evolution of this killer meiotic driver gene family.

Whole-Genome Resequencing Reveals Significant Genetic Differentiation Between Exserohilum turcicum Populations from Maize and Sorghum and Candidate Effector Genes Related to Host Specificity.

Exserohilum turcicum is a devastating fungal pathogen that infects both maize and sorghum, leading to severe leaf diseases of the two crops. According to host specificity, pathogenic isolates of E. turcicum are divided into two formae speciales, namely E. turcicum f. sp. zeae and E. turcicum f. sp. sorghi. To date, the molecular mechanism underlying the host specificity of E. turcicum is marginally known. In this study, the whole genomes of 60 E. turcicum isolates collected from both maize and sorghum were resequenced, which enabled identification of 233,022 single-nucleotide polymorphisms (SNPs) in total. Phylogenetic analysis indicated that all isolates are clustered into four genetic groups that have a close relationship with host source. This observation is validated by the result of principal component analysis. Analysis of population structure revealed that there is obvious genetic differentiation between two populations from maize and sorghum. Further analysis showed that 5,431 SNPs, including 612 nonsynonymous SNPs, are completely co-segregated with the host source. These nonsynonymous SNPs are located in 539 genes, among which 18 genes are predicted to encode secretory proteins, including six putative effector genes named SIX13-like, Ecp6, GH12, GH28-1, GH28-2, and CHP1. Sequence polymorphism analysis revealed various numbers of SNPs in the coding regions of these genes. These findings provide new insights into the molecular basis of host specificity in E. turcicum.

Genetic Variability of the Laxmann’s Shrew (<i>Sorex caecutiens</i> Laxmann, 1788) of the Kolyma River Basin and Chukotka

Nucleotide sequence polymorphism and haplotypic diversity of the cytochrome b gene were studied in samples of the Laxmann’s Shrew Sorex caecutiens koreni from the populations of the Kolyma River basin and Chukotka. 56 cytb-haplotypes differing among themselves by 87 mutations in 84 sites were found. The cytb-haplotypes of individuals of the subspecies S. c. koreni and S. c. macropygmaeus were shown to belong to the same maternal mtDNA lineage, their monophyletic origin from the same ancestral haplotype Scb1, and the presence of identical cytb-haplotypes in both subspecies. The mutual isolation of the shrew populations in the upper basin of the Kolyma River and Chukotka is shown. The proportion of intrapopulation genetic variability is 95.74%, intergroup variability is 5.74%, and intragroup variability is –1.48%, which reflects a very high level of intrapopulation variability of cytb haplotypes in the studied groups of S. c. koreni, their monophyletic origin, and the genetic unity of the Chukchi samples. The values of molecular diversity indices testify to the stability of populations that have a high value of effective abundance for a long time.

Phylogeography of the Sinica Group of Macaques in the Himalayas: Taxonomic and Evolutionary Implications.

Owing to the taxonomic incongruence between the morphological features and genetic relationships of the sinica group of macaques (genus Macaca), the taxonomy of this macaque group has remained inconclusive. We aimed to resolve the taxonomic quandary and improve our understanding of the historical biogeography of the group by including macaque DNA samples from previously unsampled areas in the Himalayas. We sequenced and analyzed three mitochondrial DNA loci [cytochrome b (CYTB), cytochrome oxidase subunit 1 (COI) and D-loop; 2898 bp] for sequence polymorphism, phylogenetics, species delimitation, and ancestral area reconstruction. We confirmed the occurrence of Arunachal macaque (Macaca munzala) on the southern slopes of the Eastern Himalayas in the Xizang Zizhiqu (Tibet Autonomous Region) of China. The results revealed that the sinica group of macaques is a parapatric species group composed of seven distinct species. Phylogenetic and species delimitation analyses revealed that the two previously considered subspecies of Assamese macaques (the eastern subspecies M. assamensis assamensis and the western subspecies M. a. pelops) are two distinct species. The eastern Assamese macaque is a sister species to the Tibetan macaque, whereas the western Assamese macaque and Arunachal macaque are the closest genetic sister species. The sinica group of macaques underwent five vicariance and seven dispersal radiations in the past, which mainly coincided with the Quaternary climatic oscillations between the late Pliocene and the late Pleistocene. By integrating our phylogenetic and ancestral area reconstruction results with findings from previous paleontological and molecular studies, we propose a robust hypothesis about the phylogeography of the sinica group of macaques.

Comparative single-cell multiome identifies evolutionary changes in neural progenitor cells during primate brain development.

Understanding the cellular and genetic mechanisms driving human-specific features of cortical development remains a challenge. We generated a cell-type resolved atlas of transcriptome and chromatin accessibility in the developing macaque and mouse prefrontal cortex (PFC). Comparing with published human data, our findings demonstrate that although the cortex cellular composition is overall conserved across species, progenitor cells show significant evolutionary divergence in cellular properties. Specifically, human neural progenitors exhibit extensive transcriptional rewiring in growth factor and extracellular matrix (ECM) pathways. Expression of the human-specific progenitor marker ITGA2 in the fetal mouse cortex increases the progenitor proliferation and the proportion of upper-layer neurons. These transcriptional divergences are primarily driven by altered activity in the distal regulatory elements. The chromatin regions with human-gained accessibility are enriched with human-specific sequence changes and polymorphisms linked to intelligence and neuropsychiatric disorders. Our results identify evolutionary changes in neural progenitors and putative gene regulatory mechanisms shaping primate brain evolution.