Search

Medicago truncatula Nod Factor Perception (MtNFP) plays a role in both the Rhizobium-Legume (RL) symbiosis and plant immunity, and evidence suggests that the immune-related function of MtNFP is relevant for symbiosis. To better understand these roles of MtNFP, we sought to identify new interacting partners. We screened a yeast-2-hybrid cDNA library from Aphanomyces euteiches infected and noninfected M. truncatula roots. The M. truncatula leucine-rich repeat (LRR) receptor-like kinase SUPPRESSOR OF BIR1 (MtSOBIR1) was identified as an interactor of MtNFP and was characterised for kinase activity, and potential roles in symbiosis and plant immunity. We showed that the kinase domain of MtSOBIR1 is active and can transphosphorylate the pseudo-kinase domain of MtNFP. MtSOBIR1 could functionally complement Atsobir1 and Nbsobir1/sobir1-like mutants for defence activation, and Mtsobir1 mutants were defective in immune responses to A. euteiches. For symbiosis, we showed that Mtsobir1 mutant plants had both a strong, early infection defect and defects in the defence suppression in nodules, and both effects were plant genotype- and rhizobial strain-specific. This work highlights a conserved function for MtSOBIR1 in activating defence responses to pathogen attack, and potentially novel symbiotic functions of downregulating defence in association with the control of symbiotic specificity., (© 2024 John Wiley & Sons Ltd.)

Cannabis sativa has a long history of domestication both for its bioactive compounds and its fibers. This has produced hundreds of varieties, usually characterized in the literature by chemotypes, with Δ 9 -THC and CBD content as the main markers. However, chemotyping could also be done based on minor compounds (phytocannabinoids and others). In this work, a workflow, which we propose to name cannabinomics, combines mass spectrometry of the whole metabolome and statistical analysis to help differentiate C. sativa varieties and deciphering their characteristic markers. By applying this cannabinomics approach to the data obtained from 20 varieties of C. sativa (classically classified as chemotype I, II, or III), we compared the results with those obtained by a targeted quantification of 11 phytocannabinoids. Cannabinomics can be considered as a complementary tool for phenotyping and genotyping, allowing the identification of minor compounds playing a key role as markers of differentiation., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Lipo-chitooligosaccharides (LCOs) were originally found as symbiotic signals called Nod Factors (Nod-LCOs) controlling the nodulation of legumes by rhizobia. More recently, LCOs were also found in symbiotic fungi and, more surprisingly, very widely in the kingdom Fungi, including in saprophytic and pathogenic fungi. The LCO-V(C18:1, fucosylated/methyl fucosylated), hereafter called Fung-LCOs, are the LCO structures most commonly found in fungi. This raises the question of how legume plants such as Medicago truncatula can discriminate between Nod-LCOs and Fung-LCOs. To address this question, we performed a genome-wide association study on 173 natural accessions of M. truncatula, using a root branching phenotype and a newly developed local score approach. Both Nod-LCOs and Fung-LCOs stimulated root branching in most accessions, but the root responses to these two types of LCO molecules were not correlated. In addition, the heritability of the root response was higher for Nod-LCOs than for Fung-LCOs. We identified 123 loci for Nod-LCO and 71 for Fung-LCO responses, of which only one was common. This suggests that Nod-LCOs and Fung-LCOs both control root branching but use different molecular mechanisms. The tighter genetic constraint of the root response to Fung-LCOs possibly reflects the ancestral origin of the biological activity of these molecules., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Experimental Biology. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

The quest for signatures of selection using single nucleotide polymorphism (SNP) data has proven efficient to uncover genes involved in conserved and/or adaptive molecular functions, but none of the statistical methods were designed to identify interacting alleles as targets of selective processes. Here, we propose a statistical test aimed at detecting epistatic selection, based on a linkage disequilibrium (LD) measure accounting for population structure and heterogeneous relatedness between individuals. SNP-based ([Formula: see text]) and window-based ([Formula: see text]) statistics fit a Student distribution, allowing to test the significance of correlation coefficients. As a proof of concept, we use SNP data from the Medicago truncatula symbiotic legume plant and uncover a previously unknown gene coadaptation between the MtSUNN (Super Numeric Nodule) receptor and the MtCLE02 (CLAVATA3-Like) signaling peptide. We also provide experimental evidence supporting a MtSUNN-dependent negative role of MtCLE02 in symbiotic root nodulation. Using human HGDP-CEPH SNP data, our new statistical test uncovers strong LD between SLC24A5 (skin pigmentation) and EDAR (hairs, teeth, sweat glands development) world-wide, which persists after correction for population structure and relatedness in Central South Asian populations. This result suggests that epistatic selection or coselection could have contributed to the phenotypic make-up in some human populations. Applying this approach to genome-wide SNP data will facilitate the identification of coadapted gene networks in model or non-model organisms.

Plant immunity is often negatively impacted by heat stress. However, the underlying molecular mechanisms remain poorly characterized. Based on a genome-wide association mapping approach, this study aims to identify in Arabidopsis thaliana the genetic bases of robust resistance mechanisms to the devastating pathogen Ralstonia solanacearum under heat stress. A local mapping population was phenotyped against the R. solanacearum GMI1000 strain at 27 and 30 °C. To obtain a precise description of the genetic architecture underlying natural variation of quantitative disease resistance (QDR), we applied a genome-wide local score analysis. Alongside an extensive genetic variation found in this local population at both temperatures, we observed a playful dynamics of quantitative trait loci along the infection stages. In addition, a complex genetic network of interacting loci could be detected at 30 °C. As a first step to investigate the underlying molecular mechanisms, the atypical meiotic cyclin SOLO DANCERS gene was validated by a reverse genetic approach as involved in QDR to R. solanacearum at 30 °C. In the context of climate change, the complex genetic architecture underlying QDR under heat stress in a local mapping population revealed candidate genes with diverse molecular functions., (© 2020 The Authors. Molecular Plant Pathology published by British Society for Plant Pathology and John Wiley & Sons Ltd.)

Fabeae legumes such as pea and faba bean form symbiotic nodules with a large diversity of soil Rhizobium leguminosarum symbiovar viciae (Rlv) bacteria. However, bacteria competitive to form root nodules (CFN) are generally not the most efficient to fix dinitrogen, resulting in a decrease in legume crop yields. Here, we investigate differential selection by host plants on the diversity of Rlv. A large collection of Rlv was collected by nodule trapping with pea and faba bean from soils at five European sites. Representative genomes were sequenced. In parallel, diversity and abundance of Rlv were estimated directly in these soils using metabarcoding. The CFN of isolates was measured with both legume hosts. Pea/faba bean CFN were associated to Rlv genomic regions. Variations of bacterial pea and/or faba bean CFN explained the differential abundance of Rlv genotypes in pea and faba bean nodules. No evidence was found for genetic association between CFN and variations in the core genome, but variations in specific regions of the nod locus, as well as in other plasmid loci, were associated with differences in CFN. These findings shed light on the genetic control of CFN in Rlv and emphasise the importance of host plants in controlling Rhizobium diversity., (© 2019 The Authors. New Phytologist © 2019 New Phytologist Trust.)

Quantitative trait loci (QTL) with small effects, which are pervasive in quantitative phenotypic variation, are difficult to detect in genome-wide association studies (GWAS). To improve their detection, we propose to use a local score approach that accounts for the surrounding signal due to linkage disequilibrium, by accumulating association signals from contiguous single markers. Simulations revealed that, in a GWAS context with high marker density, the local score approach outperforms single SNP p-value-based tests for detecting minor QTL (heritability of 5-10%) and is competitive with regard to alternative methods, which also aggregate p-values. Using more than five million SNPs, this approach was applied to identify loci involved in Quantitative Disease Resistance (QDR) to different isolates of the plant root rot pathogen Aphanomyces euteiches, from a GWAS performed on a collection of 174 accessions of the model legume Medicago truncatula. We refined the position of a previously reported major locus, underlying MYB/NB-ARC/tyrosine kinase candidate genes conferring resistance to two closely related A. euteiches isolates belonging to pea pathotype I. We also discovered a diversity of minor resistance QTL, not detected using p-value-based tests, some of which being putatively shared in response to pea (pathotype I and III) and/or alfalfa (race 1 and 2) isolates. Candidate genes underlying these QTL suggest pathogen effector recognition and plant proteasome as key functions associated with M. truncatula resistance to A. euteiches. GWAS on any organism can benefit from the local score approach to uncover many weak-effect QTL.

Plant -specific lysin-motif receptor-like kinases (LysM-RLKs) are implicated in the perception of N-acetyl glucosamine-containing compounds, some of which are important signal molecules in plant-microbe interactions. Among these, both lipo-chitooligosaccharides (LCOs) and chitooligosaccharides (COs) are proposed as arbuscular mycorrhizal (AM) fungal symbiotic signals. COs can also activate plant defence, although there are scarce data about CO production by pathogens, especially nonfungal pathogens. We tested Medicago truncatula mutants in the LysM-RLK MtLYK9 for their abilities to interact with the AM fungus Rhizophagus irregularis and the oomycete pathogen Aphanomyces euteiches. This prompted us to analyse whether A. euteiches can produce COs. Compared with wild-type plants, Mtlyk9 mutants had fewer infection events and were less colonised by the AM fungus. By contrast, Mtlyk9 mutants were more heavily infected by A. euteiches and showed more disease symptoms. Aphanomyces euteiches was also shown to produce short COs, mainly CO II, but also CO III and CO IV, and traces of CO V, both ex planta and in planta. MtLYK9 thus has a dual role in plant immunity and the AM symbiosis, which raises questions about the functioning and the ancestral origins of such a receptor protein., (© 2019 The Authors. New Phytologist © 2019 New Phytologist Trust.)

Natural variations help in identifying genetic mechanisms of morphologically and developmentally complex traits. Mountainous habitats provide an altitudinal gradient where one species encounters different abiotic conditions. We report the study of 341 individuals of Arabidopsis thaliana derived from 30 natural populations not belonging to the 1001 genomes, collected at increasing altitudes, between 200 and 1800 m in the Pyrenees. Class III peroxidases and ribosomal RNA sequences were used as markers to determine the putative genetic relationships among these populations along their altitudinal gradient. Using Bayesian-based statistics and phylogenetic analyses, these Pyrenean populations appear with significant divergence from the other regional accessions from 1001 genome (i.e., from north Spain or south France). Individuals of these populations exhibited varying phenotypic changes, when grown at sub-optimal temperature (22 vs. 15°C). These phenotypic variations under controlled conditions reflected intraspecific morphological variations. This study could bring new information regarding the west European population structure of A. thaliana and its phenotypic variations at different temperatures. The integrative analysis combining genetic, phenotypic variation and environmental datasets is used to analyze the acclimation of population in response to temperature changes. Regarding their geographical proximity and environmental diversity, these populations represent a tool of choice for studying plant response to temperature variation., Highlights: -Studying the natural diversity of A. thaliana in the Pyrenees mountains helps to understand European population structure and to evaluate the phenotypic trait variation in response to climate change.

The roots of most land plants are colonized by symbiotic arbuscular mycorrhiza (AM) fungi. To facilitate this symbiosis, plant genomes encode a set of genes required for microbial perception and accommodation. However, the extent to which infection by filamentous root pathogens also relies on some of these genes remains an open question. Here, we used genome-wide association mapping to identify genes contributing to colonization of Medicago truncatula roots by the pathogenic oomycete Phytophthora palmivora. Single-nucleotide polymorphism (SNP) markers most significantly associated with plant colonization response were identified upstream of RAD1, which encodes a GRAS transcription regulator first negatively implicated in root nodule symbiosis and recently identified as a positive regulator of AM symbiosis. RAD1 transcript levels are up-regulated both in response to AM fungus and, to a lower extent, in infected tissues by P. palmivora where its expression is restricted to root cortex cells proximal to pathogen hyphae. Reverse genetics showed that reduction of RAD1 transcript levels as well as a rad1 mutant are impaired in their full colonization by AM fungi as well as by P. palmivora. Thus, the importance of RAD1 extends beyond symbiotic interactions, suggesting a general involvement in M. truncatula microbe-induced root development and interactions with unrelated beneficial and detrimental filamentous microbes., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Experimental Biology.)

Plant NF-Y transcription factors control a wide array of biological functions enabling appropriate reproductive and developmental processes as well as adaptation to various abiotic and biotic environments. In Medicago truncatula , MtNF-YA1 was previously identified as a key determinant for nodule development and establishment of rhizobial symbiosis. Here, we highlight a new role for this protein in compatibility to Aphanomyces euteiches , a root pathogenic oomycete. The Mtnf-ya1-1 mutant plants showed better survival rate, reduced symptoms, and increased development of their root apparatus as compared to their wild-type (WT) background A17. MtNF-YA-1 was specifically up-regulated by A. euteiches in F83005.5, a highly susceptible natural accession of M. truncatula while transcript level remained stable in A17, which is partially resistant. The role of MtNF-YA1 in F83005.5 susceptibility was further documented by reducing MtNF-YA1 expression either by overexpression of the miR169q, a microRNA targeting MtNF-YA1 , or by RNAi approaches leading to a strong enhancement in the resistance of this susceptible line. Comparative analysis of the transcriptome of WT and Mtnf-ya1-1 led to the identification of 1509 differentially expressed genes. Among those, almost 36 defense-related genes were constitutively expressed in Mtnf-ya1-1 , while 20 genes linked to hormonal pathways were repressed. In summary, we revealed an unexpected dual role for this symbiotic transcription factor as a key player in the compatibility mechanisms to a pathogen.

Local climatic conditions likely constitute an important selective pressure on genes underlying important fitness-related traits such as flowering time, and in many species, flowering phenology and climatic gradients strongly covary. To test whether climate shapes the genetic variation on flowering time genes and to identify candidate flowering genes involved in the adaptation to environmental heterogeneity, we used a large Medicago truncatula core collection to examine the association between nucleotide polymorphisms at 224 candidate genes and both climate variables and flowering phenotypes. Unlike genome-wide studies, candidate gene approaches are expected to enrich for the number of meaningful trait associations because they specifically target genes that are known to affect the trait of interest. We found that flowering time mediates adaptation to climatic conditions mainly by variation at genes located upstream in the flowering pathways, close to the environmental stimuli. Variables related to the annual precipitation regime reflected selective constraints on flowering time genes better than the other variables tested (temperature, altitude, latitude or longitude). By comparing phenotype and climate associations, we identified 12 flowering genes as the most promising candidates responsible for phenological adaptation to climate. Four of these genes were located in the known flowering time QTL region on chromosome 7. However, climate and flowering associations also highlighted largely distinct gene sets, suggesting different genetic architectures for adaptation to climate and flowering onset., (© 2016 John Wiley & Sons Ltd.)

To gain an insight into the molecular mechanisms of quantitative disease resistance in Medicago truncatula to the root-infecting oomycete Aphanomyces euteiches, we selected two near-isogenic lines (NILs), NR and NS, partially resistant and susceptible, respectively, differing in the allelic state of the quantitative resistance locus (QRL) prAe1 (partially resistant to A. euteiches 1). Complementary molecular and cytological phenotyping methods showed that prAe1 alone confers quantitative resistance to A. euteiches. Root and stem tissues were colonized in NS plants and 80% of NS plants died by 21 days post-inoculation (dpi). In contrast, A. euteiches mycelium was restricted to the root cortex and the spread of symptoms was arrested in aerial parts of NR plants. A transcriptome analysis performed at 0, 1 and 6 dpi identified 1198 differentially expressed genes (DEGs) between NR and NS lines. More than 87% of the DEGs were significantly more expressed in NR. The highest number of DEGs was found in control conditions, with 723 genes over-expressed in NR versus 85 in NS. Genes belonging to secondary metabolism, pathogenesis-related (PR) proteins and kinases were significantly enriched. The significant role of the flavonoid pathway in resistance was corroborated by the detection of larger amounts of flavonoids in NR roots and the inhibition of A. euteiches zoospore germination by 2'-O-methyl-isoliquiritigenin, a compound synthesized by enzymes specifically induced in NR. Our study revealed that prAe1-dependent resistance relies mainly on the constitutive expression of defence-related pathways and signalling elements, which can be re-amplified in later time points of the infection., (© 2015 BSPP AND JOHN WILEY & SONS LTD.)

Medicago truncatula is a model legume species used to investigate plant-microorganism interactions, notably root symbioses. Massive population genomic and transcriptomic data now available for this species open the way for a comprehensive investigation of genomic variations associated with adaptation of M. truncatula to its environment. Here we performed a fine-scale genome scan of selective sweep signatures in M. truncatula using more than 15 million single nucleotide polymorphisms identified on 283 accessions from two populations (Circum and Far West), and exploited annotation and published transcriptomic data to identify biological processes associated with molecular adaptation. We identified 58 swept genomic regions with a 15 kb average length and comprising 3.3 gene models on average. The unimodal sweep state probability distribution in these regions enabled us to focus on the best single candidate gene per region. We detected two unambiguous species-wide selective sweeps, one of which appears to underlie morphological adaptation. Population genomic analyses of the remaining 56 sweep signatures indicate that sweeps identified in the Far West population are less population-specific and probably more ancient than those identified in the Circum population. Functional annotation revealed a predominance of immunity-related adaptations in the Circum population. Transcriptomic data from accessions of the Far West population allowed inference of four clusters of coregulated genes putatively involved in the adaptive control of symbiotic carbon flow and nodule senescence, as well as in other root adaptations upon infection with soil microorganisms. We demonstrate that molecular adaptations in M. truncatula were primarily triggered by selective pressures from root-associated microorganisms., (© The Author 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Understanding how immune genetic variation is shaped by selective and neutral processes in wild populations is of prime importance in both evolutionary biology and epidemiology. The European roe deer (Capreolus capreolus) has considerably expanded its distribution range these last decades, notably by colonizing agricultural landscapes. This range shift is likely to have led to bottlenecks and increased roe deer exposure to a new range of pathogens that until recently predominantly infected humans and domestic fauna. We therefore investigated the historical and contemporary forces that have shaped variability in a panel of genes involved in innate and acquired immunity in roe deer, including Mhc-Drb and genes encoding cytokines or toll-like receptors (TLRs). Together, our results suggest that genetic drift is the main contemporary evolutionary force shaping immunogenetic variation within populations. However, in contrast to the classical view, we found that some innate immune genes involved in micropathogen recognition (e.g. Tlrs) continue to evolve dynamically in roe deer in response to pathogen-mediated positive selection. Most studied Tlrs (Tlr2, Tlr4 and Tlr5) had similarly high levels of amino acid diversity in the three studied populations including one recently established in southwestern France that showed a clear signature of genetic bottleneck. Tlr2 implicated in the recognition of Gram-positive bacteria in domestic ungulates, showed strong evidence of balancing selection. The high immunogenetic variation revealed here implies that roe deer are able to cope with a wide spectrum of pathogens and to respond rapidly to emerging infectious diseases., (© 2015 John Wiley & Sons Ltd.)

• The use of quantitative disease resistance (QDR) is a promising strategy for promoting durable resistance to plant pathogens, but genes involved in QDR are largely unknown. To identify genetic components and accelerate improvement of QDR in legumes to the root pathogen Aphanomyces euteiches, we took advantage of both the recently generated massive genomic data for Medicago truncatula and natural variation of this model legume. • A high-density (≈5.1 million single nucleotide polymorphisms (SNPs)) genome-wide association study (GWAS) was performed with both in vitro and glasshouse phenotyping data collected for 179 lines. • GWAS identified several candidate genes and pinpointed two independent major loci on the top of chromosome 3 that were detected in both phenotyping methods. Candidate SNPs in the most significant locus (σ(A)²= 23%) were in the promoter and coding regions of an F-box protein coding gene. Subsequent qRT-PCR and bioinformatic analyses performed on 20 lines demonstrated that resistance is associated with mutations directly affecting the interaction domain of the F-box protein rather than gene expression. • These results refine the position of previously identified QTL to specific candidate genes, suggest potential molecular mechanisms, and identify new loci explaining QDR against A. euteiches., (© 2013 The Authors. New Phytologist © 2013 New Phytologist Trust.)

Plant LysM proteins control the perception of microbial-derived N-acetylglucosamine compounds for the establishment of symbiosis or activation of plant immunity. This raises questions about how plants, and notably legumes, can differentiate friends and foes using similar molecular actors and whether any receptors can intervene in both symbiosis and resistance. To study this question, nfp and lyk3 LysM-receptor like kinase mutants of Medicago truncatula that are affected in the early steps of nodulation, were analysed following inoculation with Aphanomyces euteiches, a root oomycete. The role of NFP in this interaction was further analysed by overexpression of NFP and by transcriptome analyses. nfp, but not lyk3, mutants were significantly more susceptible than wildtype plants to A. euteiches, whereas NFP overexpression increased resistance. Transcriptome analyses on A. euteiches inoculation showed that mutation in the NFP gene led to significant changes in the expression of c. 500 genes, notably involved in cell dynamic processes previously associated with resistance to pathogen penetration. nfp mutants also showed an increased susceptibility to the fungus Colletotrichum trifolii. These results demonstrate that NFP intervenes in M. truncatula immunity, suggesting an unsuspected role for NFP in the perception of pathogenic signals., (© 2013 The Authors. New Phytologist © 2013 New Phytologist Trust.)

Detecting genetic signatures of selection is of great interest for many research issues. Common approaches to separate selective from neutral processes focus on the variance of F(ST) across loci, as does the original Lewontin and Krakauer (LK) test. Modern developments aim to minimize the false positive rate and to increase the power, by accounting for complex demographic structures. Another stimulating goal is to develop straightforward parametric and computationally tractable tests to deal with massive SNP data sets. Here, we propose an extension of the original LK statistic (T(LK)), named T(F-LK), that uses a phylogenetic estimation of the population's kinship (F) matrix, thus accounting for historical branching and heterogeneity of genetic drift. Using forward simulations of single-nucleotide polymorphisms (SNPs) data under neutrality and selection, we confirm the relative robustness of the LK statistic (T(LK)) to complex demographic history but we show that T(F-LK) is more powerful in most cases. This new statistic outperforms also a multinomial-Dirichlet-based model [estimation with Markov chain Monte Carlo (MCMC)], when historical branching occurs. Overall, T(F-LK) detects 15-35% more selected SNPs than T(LK) for low type I errors (P < 0.001). Also, simulations show that T(LK) and T(F-LK) follow a chi-square distribution provided the ancestral allele frequencies are not too extreme, suggesting the possible use of the chi-square distribution for evaluating significance. The empirical distribution of T(F-LK) can be derived using simulations conditioned on the estimated F matrix. We apply this new test to pig breeds SNP data and pinpoint outliers using T(F-LK), otherwise undetected using the less powerful T(LK) statistic. This new test represents one solution for compromise between advanced SNP genetic data acquisition and outlier analyses.

Rhesus macaque (Macaca mulatta) and long-tailed macaque (Macaca fascicularis) are the 2 most commonly used primate model species in biomedical sciences. Although morphological studies have revealed a weak hybridization at the interspecific contact zone, in the north of Indochina, a molecular study has suggested an ancient introgression from rhesus to long-tailed macaque into the Indo-Chinese peninsula. However, the gene flow between these 2 taxa has never been quantified using genetic data and theoretical models. In this study, we have examined genetic variation within and between the parapatric Chinese rhesus macaque and Indo-Chinese long-tailed macaque populations, using 13 autosomal, 5 sex-linked microsatellite loci and mitochondrial DNA sequence data. From these data, we assessed genetic structure and estimated gene flow using a Bayesian clustering approach and the "Isolation with Migration" model. Our results reveal a weak interspecific genetic differentiation at both autosomal and sex-linked loci, suggesting large population sizes and/or gene flow between populations. According to the Bayesian clustering, Chinese rhesus macaque is a highly homogeneous gene pool that contributes strongly to the current Indo-Chinese long-tailed macaque genetic makeup, whether or not current admixture is assumed. Coalescent simulations, which integrated the characteristics of the loci, pointed out 1) a higher effective population size in rhesus macaque, 2) no mitochondrial gene flow, and 3) unilateral and male-mediated nuclear gene flow of approximately 10 migrants per generation from rhesus to long-tailed macaque. These patterns of genetic structure and gene flow suggest extensive ancient introgression from Chinese rhesus macaque into the Indo-Chinese long-tailed macaque population.

Many genetic studies on catarrhines use microsatellite markers that were isolated from human DNA. A large number of these markers have been characterized in the great apes, macaques and baboons. However, there are few or no markers available for other members of this group. In this study, an extensive literature search was performed to find microsatellite markers that had been successfully amplified across a range of catarrhine species. These conserved loci can provide a valuable starting point for characterizing loci in other catarrhines. Finally, microsatellite markers were tested in a range of species that are not well represented in the literature., (Copyright 2009 S. Karger AG, Basel.)

In humans, the single polymorphic B locus of the major histocompatibility complex is linked to the microsatellite MIB. In rhesus macaques, however, haplotypes are characterized by the presence of unique combinations of multiple B genes, which may display different levels of polymorphism. The aim of the study was to shed light on the evolutionary history of this highly complex region. First, the robustness of the microsatellite MIB-linked to almost half of the B genes in rhesus macaques (Mamu-B)-for accurate B haplotyping was studied. Based on the physical map of an established haplotype comprising 7 MIB loci, each located next to a certain Mamu-B gene, two MIB loci, MIB1 and MIB6, were investigated in a panel of MHC homozygous monkeys. MIB1 revealed a complex genotyping pattern, whereas MIB6 analysis resulted in the detection of one or no amplicon. Both patterns are specific for a given B haplotype, show Mendelian segregation, and even allow a more precise haplotype definition than do traditional typing methods. Second, a search was performed for retroelements that may have played a role in duplication processes as observed in the macaque B region. This resulted in the description of two types of duplicons. One basic unit comprises an expressed Mamu-B gene, adjacent to an HERV16 copy closely linked to MIB. The second type of duplicon comprises a Mamu-B (pseudo)gene, linked to a truncated HERV16 structure lacking its MIB segment. Such truncation seems to coincide with the loss of B gene transcription. Subsequent to the duplication processes, recombination between MIB and Mamu-B loci appears to have occurred, resulting in a hyperplastic B region. Thus, analysis of MIB in addition to B loci allows deciphering of the compound evolutionary history of the class I B region in Old World monkeys.

Background: In sharp contrast to humans and great apes, the expanded Mhc-B region of rhesus and cynomolgus macaques is characterized by the presence of differential numbers and unique combinations of polymorphic class I B genes per haplotype. The MIB microsatellite is closely linked to the single class I B gene in human and in some great apes studied. The physical map of the Mhc of a heterozygous rhesus monkey provides unique material to analyze MIB and Mamu-B copy number variation and then allows one to decipher the compound evolutionary history of this region in primate species., Results: In silico research pinpointed 12 MIB copies (duplicons), most of which are associated with expressed B-genes that cluster in a separate clade in the phylogenetic tree. Generic primers tested on homozygous rhesus and pedigreed cynomolgus macaques allowed the identification of eight to eleven MIB copies per individual. The number of MIB copies present per haplotype varies from a minimum of three to six in cynomolgus macaques and from five to eight copies in rhesus macaques. Phylogenetic analyses highlight a strong transpecific sharing of MIB duplicons. Using the physical map, we observed that, similar to MIB duplicons, highly divergent Mamu-B genes can be present on the same haplotype. Haplotype variation as reflected by the copy number variation of class I B loci is best explained by recombination events, which are found to occur between MIBs and Mamu-B., Conclusion: The data suggest the existence of highly divergent MIB and Mamu-B lineages on a given haplotype, as well as variable MIB and B copy numbers and configurations, at least in rhesus macaque. Recombination seems to occur between MIB and Mamu-B loci, and the resulting haplotypic plasticity at the individual level may be a strategy to better cope with pathogens. Therefore, evolutionary inferences based on the multiplicated MIB loci but also other markers close to B-genes appear to be promising for the study of B-region organization and evolution in primates.

We studied the mitochondrial DNA (mtDNA) polymorphism of 304 Macaca fascicularis fascicularis (M. f. fascicularis) individuals, representative of 4 cynomolgus macaque populations (Indochina, Indonesia, Philippines, and Mauritius). By sequencing a 590-bp fragment in the hypervariable II region of the D-loop region, we defined 70 haplotypes. The homologous region was also characterized in 22 Chinese Macaca mulatta and 2 Macaca sylvanus. The phylogenetic analysis confirms the monophyly of M. f. fascicularis and defines 2 haplotype groups inside the M. f. fascicularis clade: one "insular," encompassing 6 Philippines, 2 Mauritius, and 31 Indonesian haplotypes, the other "continental" that contains all Indochinese and 6 Indonesian haplotypes. Continental and insular group divergence time was estimated to be approximately 10(6) years before present (BP). Among Indonesian haplotypes, some have a continental origin. This suggests either direct migration from mainland to Indonesia or that remnant lineages from an ancient population genetically close to the mainland (i.e., in the Sunda Shelf, <550 000 years BP) were subsequently brought southward to Indonesia. The low nucleotide diversity in the Philippines population suggests a bottleneck following colonization by Indonesian individuals, around 110 000 years BP. mtDNA and further observations of nuclear genetic data corroborate the mixed origin (Indonesian/continental) hypothesis of Mauritius individuals and a population bottleneck.

Chimpanzee populations are diminishing as a consequence of human activities, and as a result this species is now endangered. In the context of conservation programmes, genetic data can add vital information, for instance on the genetic diversity and structure of threatened populations. Single nucleotide polymorphisms (SNP) are biallelic markers that are widely used in human molecular studies and can be implemented in efficient microarray systems. This technology offers the potential of robust, multiplexed SNP genotyping at low reagent cost in other organisms than humans, but it is not commonly used yet in wild population studies. Here, we describe the characterization of new SNPs in Y-chromosomal intronic regions in chimpanzees and also identify SNPs from mitochondrial genes, with the aim of developing a microarray system that permits the simultaneous study of both paternal and maternal lineages. Our system consists of 42 SNPs for the Y chromosome and 45 SNPs for the mitochondrial genome. We demonstrate the applicability of this microarray in a captive population where genotypes accurately reflected its large pedigree. Two wild-living populations were also analysed and the results show that the microarray will be a useful tool alongside microsatellite markers, since it supplies complementary information about population structure and ecology. SNP genotyping using microarray technology, therefore, is a promising approach and may become an essential tool in conservation genetics to help in the management and study of captive and wild-living populations. Moreover, microarrays that combine SNPs from different genomic regions could replace microsatellite typing in the future., (© 2007 The Authors.)

Chimpanzees experienced a reduction of the allelic repertoire at the major histocompatibility complex (MHC) class I A and B loci, which may have been caused by a retrovirus belonging to the simian immunodeficiency virus (SIV) family. Extended MHC haplotypes were defined in a pedigreed chimpanzee colony. Comparison of genetic variation at microsatellite markers mapping inside and outside the Mhc region was carried out in humans and chimpanzees to investigate the genomic extent of the repertoire reduction. Multilocus demographic analyses underscored that chimpanzees indeed experienced a selective sweep that mainly targeted the chromosomal segment carrying the Mhc class I region. Probably due to genetic linkage, the sweep also affected other polymorphic loci, mapping in the close vicinity of the Mhc class I region genes. Nevertheless, although the allelic repertoire at particular Mhc class I and II loci appears to be limited, naturally occurring recombination events allowed the establishment of haplotype diversity after the sweep. However, recombination did not have sufficient time to erase the signal of the selective sweep.