Genetic diversity and population structure in Vicia faba L. landraces and wild related species assessed by nuclear SSRs

Research Article

Nuclear factor I-C links platelet-derived growth factor and transforming growth factor beta1 signaling to skin wound healing progression.

Transforming growth factor beta (TGF-beta) and platelet-derived growth factor A (PDGFAlpha) play a central role in tissue morphogenesis and repair, but their interplay remain poorly understood. The nuclear factor I C (NFI-C) transcription factor has been implicated in TGF-beta signaling, extracellular matrix deposition, and skin appendage pathologies, but a potential role in skin morphogenesis or healing had not been assessed. To evaluate this possibility, we performed a global gene expression analysis in NFI-C(-/-) and wild-type embryonic primary murine fibroblasts. This indicated that NFI-C acts mostly to repress gene expression in response to TGF-beta1. Misregulated genes were prominently overrepresented by regulators of connective tissue inflammation and repair. In vivo skin healing revealed a faster inflammatory stage and wound closure in NFI-C(-/-) mice. Expression of PDGFA and PDGF-receptor alpha were increased in wounds of NFI-C(-/-) mice, explaining the early recruitment of macrophages and fibroblasts. Differentiation of fibroblasts to contractile myofibroblasts was also elevated, providing a rationale for faster wound closure. Taken together with the role of TGF-beta in myofibroblast differentiation, our results imply a central role of NFI-C in the interplay of the two signaling pathways and in regulation of the progression of tissue regeneration.

NR2E3 mutations in enhanced S-cone sensitivity syndrome (ESCS), Goldmann-Favre syndrome (GFS), clumped pigmentary retinal degeneration (CPRD), and retinitis pigmentosa (RP).

NR2E3, also called photoreceptor-specific nuclear receptor (PNR), is a transcription factor of the nuclear hormone receptor superfamily whose expression is uniquely restricted to photoreceptors. There, its physiological activity is essential for proper rod and cone photoreceptor development and maintenance. Thirty-two different mutations in NR2E3 have been identified in either homozygous or compound heterozygous state in the recessively inherited enhanced S-cone sensitivity syndrome (ESCS), Goldmann-Favre syndrome (GFS), and clumped pigmentary retinal degeneration (CPRD). The clinical phenotype common to all these patients is night blindness, rudimental or absent rod function, and hyperfunction of the "blue" S-cones. A single p.G56R mutation is inherited in a dominant manner and causes retinitis pigmentosa (RP). We have established a new locus-specific database for NR2E3 (www.LOVD.nl/eye), containing all reported mutations, polymorphisms, and unclassified sequence variants, including novel ones. A high proportion of mutations are located in the evolutionarily-conserved DNA-binding domains (DBDs) and ligand-binding domains (LBDs) of NR2E3. Based on homology modeling of these NR2E3 domains, we propose a structural localization of mutated residues. The high variability of clinical phenotypes observed in patients affected by NR2E3-linked retinal degenerations may be caused by different disease mechanisms, including absence of DNA-binding, altered interactions with transcriptional coregulators, and differential activity of modifier genes.

Superfamily of steroid nuclear receptors: positive and negative regulators of gene expression.

The nuclear hormone receptor superfamily is characterized by an impressive functional diversity of its members despite a remarkable overall structural unity. A variety of ligands bind specifically to them and these receptors control gene networks that have profound effects on growth, development, and homeostasis. The ligand-receptor complexes recognize transcriptional enhancer DNA sequences, the hormone response elements, resulting in induction or repression of gene activity. The similarity between all these hormone response enhancer elements, as well as between the receptors themselves, indicates a conserved general strategy for the hormonal control of transcription by steroids. The activated receptors bind to responsive promoters and most likely mediate the assembly of stage- and tissue-specific transcription factor complexes that stimulate or inhibit gene expression.

Biogeographic origin and radiation of the Old World crocidurine shrews (Mammalia: Soricidae) inferred from mitochondrial and nuclear genes.

The crocidurine shrews include the most speciose genus of mammals, Crocidura. The origin and evolution of their radiation is, however, poorly understood because of very scant fossil records and a rather conservative external morphology between species. Here, we use an alignment of 3560 base pairs of mitochondrial and nuclear DNA to generate a phylogenetic hypothesis for the evolution of Old World shrews of the subfamily Crocidurinae. These molecular data confirm the monophyly of the speciose African and Eurasian Crocidura, which also includes the fossorial, monotypic genus Diplomesodon. The phylogenetic reconstructions give further credit to a paraphyletic position of Suncus shrews, which are placed into at least two independent clades (one in Africa and sister to Sylvisorex and one in Eurasia), at the base of the Crocidura radiation. Therefore, we recommend restricting the genus Suncus to the Palaearctic and Oriental taxa, and to consider all the African Suncus as Sylvisorex. Using molecular dating and biogeographic reconstruction analyses, we suggest a Palaearctic-Oriental origin for Crocidura dating back to the Upper Miocene (6.8 million years ago) and several subsequent colonisations of the Afrotropical region by independent lineages of Crocidura.

Design and analysis of ChIP-Seq experiments for nuclear factor I DNA-binding proteins

SUMMARY : Eukaryotic DNA interacts with the nuclear proteins using non-covalent ionic interactions. Proteins can recognize specific nucleotide sequences based on the sterical interactions with the DNA and these specific protein-DNA interactions are the basis for many nuclear processes, e.g. gene transcription, chromosomal replication, and recombination. New technology termed ChIP-Seq has been recently developed for the analysis of protein-DNA interactions on a whole genome scale and it is based on immunoprecipitation of chromatin and high-throughput DNA sequencing procedure. ChIP-Seq is a novel technique with a great potential to replace older techniques for mapping of protein-DNA interactions. In this thesis, we bring some new insights into the ChIP-Seq data analysis. First, we point out to some common and so far unknown artifacts of the method. Sequence tag distribution in the genome does not follow uniform distribution and we have found extreme hot-spots of tag accumulation over specific loci in the human and mouse genomes. These artifactual sequence tags accumulations will create false peaks in every ChIP-Seq dataset and we propose different filtering methods to reduce the number of false positives. Next, we propose random sampling as a powerful analytical tool in the ChIP-Seq data analysis that could be used to infer biological knowledge from the massive ChIP-Seq datasets. We created unbiased random sampling algorithm and we used this methodology to reveal some of the important biological properties of Nuclear Factor I DNA binding proteins. Finally, by analyzing the ChIP-Seq data in detail, we revealed that Nuclear Factor I transcription factors mainly act as activators of transcription, and that they are associated with specific chromatin modifications that are markers of open chromatin. We speculate that NFI factors only interact with the DNA wrapped around the nucleosome. We also found multiple loci that indicate possible chromatin barrier activity of NFI proteins, which could suggest the use of NFI binding sequences as chromatin insulators in biotechnology applications. RESUME : L'ADN des eucaryotes interagit avec les protéines nucléaires par des interactions noncovalentes ioniques. Les protéines peuvent reconnaître les séquences nucléotidiques spécifiques basées sur l'interaction stérique avec l'ADN, et des interactions spécifiques contrôlent de nombreux processus nucléaire, p.ex. transcription du gène, la réplication chromosomique, et la recombinaison. Une nouvelle technologie appelée ChIP-Seq a été récemment développée pour l'analyse des interactions protéine-ADN à l'échelle du génome entier et cette approche est basée sur l'immuno-précipitation de la chromatine et sur la procédure de séquençage de l'ADN à haut débit. La nouvelle approche ChIP-Seq a donc un fort potentiel pour remplacer les anciennes techniques de cartographie des interactions protéine-ADN. Dans cette thèse, nous apportons de nouvelles perspectives dans l'analyse des données ChIP-Seq. Tout d'abord, nous avons identifié des artefacts très communs associés à cette méthode qui étaient jusqu'à présent insoupçonnés. La distribution des séquences dans le génome ne suit pas une distribution uniforme et nous avons constaté des positions extrêmes d'accumulation de séquence à des régions spécifiques, des génomes humains et de la souris. Ces accumulations des séquences artéfactuelles créera de faux pics dans toutes les données ChIP-Seq, et nous proposons différentes méthodes de filtrage pour réduire le nombre de faux positifs. Ensuite, nous proposons un nouvel échantillonnage aléatoire comme un outil puissant d'analyse des données ChIP-Seq, ce qui pourraient augmenter l'acquisition de connaissances biologiques à partir des données ChIP-Seq. Nous avons créé un algorithme d'échantillonnage aléatoire et nous avons utilisé cette méthode pour révéler certaines des propriétés biologiques importantes de protéines liant à l'ADN nommés Facteur Nucléaire I (NFI). Enfin, en analysant en détail les données de ChIP-Seq pour la famille de facteurs de transcription nommés Facteur Nucléaire I, nous avons révélé que ces protéines agissent principalement comme des activateurs de transcription, et qu'elles sont associées à des modifications de la chromatine spécifiques qui sont des marqueurs de la chromatine ouverte. Nous pensons que lés facteurs NFI interagir uniquement avec l'ADN enroulé autour du nucléosome. Nous avons également constaté plusieurs régions génomiques qui indiquent une éventuelle activité de barrière chromatinienne des protéines NFI, ce qui pourrait suggérer l'utilisation de séquences de liaison NFI comme séquences isolatrices dans des applications de la biotechnologie.

Towards the Physical Map of the Trypanosoma cruzi Nuclear Genome: Construction of YAC and BAC Libraries of the Reference Clone T. cruzi CL-Brener

Strategies to construct the physical map of the Trypanosoma cruzi nuclear genome have to capitalize on three main advantages of the parasite genome, namely (a) its small size, (b) the fact that all chromosomes can be defined, and many of them can be isolated by pulse field gel electrophoresis, and (c) the fact that simple Southern blots of electrophoretic karyotypes can be used to map sequence tagged sites and expressed sequence tags to chromosomal bands. A major drawback to cope with is the complexity of T. cruzi genetics, that hinders the construction of a comprehensive genetic map. As a first step towards physical mapping, we report the construction and partial characterization of a T. cruzi CL-Brener genomic library in yeast artificial chromosomes (YACs) that consists of 2,770 individual YACs with a mean insert size of 365 kb encompassing around 10 genomic equivalents. Two libraries in bacterial artificial chromosomes (BACs) have been constructed, BACI and BACII. Both libraries represent about three genome equivalents. A third BAC library (BAC III) is being constructed. YACs and BACs are invaluable tools for physical mapping. More generally, they have to be considered as a common resource for research in Chagas disease

Analysis and chromosomal mapping of Leishmania (Leishmania) amazonensis amastigote expressed sequence tags

The characterization of expressed sequence tags (ESTs) generated from a cDNA library of Leishmania (Leishmania) amazonensis amastigotes is described. The sequencing of 93 clones generated new L. (L.) amazonensis ESTs from which 32% are not related to any other sequences in database and 68% presented significant similarities to known genes. The chromosome localization of some L. (L.) amazonensis ESTs was also determined in L. (L.) amazonensis and L. (L.) major. The characterization of these ESTs is suitable for the genome physical mapping, as well as for the identification of genes encoding cysteine proteinases implicated with protective immune responses in leishmaniasis.

Phylogeny and circumscription of Sapindaceae revisited: molecular sequence data, morphology and biogeography support recognition of a new family, Xanthoceraceae

Background and aims Recent studies have adopted a broad definition of Sapindaceae that includes taxa traditionally placed in Aceraceae and Hippocastanaceae, achieving monophyly but yielding a family difficult to characterize and for which no obvious morphological synapomorphy exists. This expanded circumscription was necessitated by the finding that the monotypic, temperate Asian genus Xanthoceras, historically placed in Sapindaceae tribe Harpullieae, is basal within the group. Here we seek to clarify the relationships of Xanthoceras based on phylogenetic analyses using a dataset encompassing nearly 3/4 of sapindaceous genera, comparing the results with information from morphology and biogeography, in particular with respect to the other taxa placed in Harpullieae. We then re-examine the appropriateness of maintaining the current broad, morphologically heterogeneous definition of Sapindaceae and explore the advantages of an alternative family circumscription. Methods Using 243 samples representing 104 of the 142 currently recognized genera of Sapindaceae s. lat. (including all in Harpullieae), sequence data were analyzed for nuclear (ITS) and plastid (matK, rpoB, trnD-trnT, trnK-matK, trnL-trnF and trnS-trnG) markers, adopting the methodology of a recent family-wide study, performing single-gene and total evidence analyses based on maximum likelihood (ML) and maximum parsimony (MP) criteria, and applying heuristic searches developed for large datasets, viz, a new strategy implemented in RAxML (for ML) and the parsimony ratchet (for MP). Bootstrap analyses were performed for each method to test for congruence between markers. Key results Our findings support earlier suggestions that Harpullieae are polyphyletic: Xanthoceras is confirmed as sister to all other sampled taxa of Sapindaceae s. lat.; the remaining members belong to three other clades within Sapindaceae s. lat., two of which correspond respectively to the groups traditionally treated as Aceraceae and Hippocastanaceae, together forming a clade sister to the largely tropical Sapindaceae s. str., which is monophyletic and morphologically coherent provided Xanthoceras is excluded. Conclusion To overcome the difficulties of a broadly circumscribed Sapindaceae, we resurrect the historically recognized temperate families Aceraceae and Hippocastanaceae, and describe a new family, Xanthoceraceae, thus adopting a monophyletic and easily characterized circumscription of Sapindaceae nearly identical to that used for over a century.

Expression and localization of PPARs in the rat ovary during follicular development and the periovulatory period.

PPARs are a family of nuclear hormone receptors involved in various processes that could influence ovarian function. We investigated the cellular localization and expression of PPARs during follicular development in ovarian tissue collected from rats 0, 6, 12, 24, and 48 h post-PMSG. A second group of animals received human CG (hCG) 48 h post-PMSG. Their ovaries were removed 0, 4, 8, 12, and 24 h post-hCG to study the periovulatory period. mRNAs corresponding to the PPAR isotypes (alpha, delta, and gamma) were localized by in situ hybridization. Changes in the levels of mRNA for the PPARs were determined by ribonuclease protection assays. PPAR gamma mRNA was localized primarily to granulosa cells, and levels of expression did not change during follicular development. Four hours post-hCG, levels of mRNA for PPAR gamma decreased (P < 0.05) but not uniformly in all follicles. At 24 h post-hCG, levels of PPAR gamma mRNA were reduced 64%, but some follicles maintained high expression. In contrast, mRNAs for PPAR alpha and delta were located primarily in theca and stroma, and their levels did not change during the intervals studied. To investigate the physiologic significance of PPAR gamma in the ovary, granulosa cells from PMSG-primed rats were cultured for 48 h with prostaglandin J(2) (PGJ(2)) and ciglitazone, PPAR gamma activators. Both compounds increased progesterone and E2 secretion (P < 0.05). These data suggest that PPAR gamma is involved in follicular development, has a negative influence on the luteinization of granulosa cells, and/or regulates the periovulatory shift in steroid production. The more general and steady expression of PPARs alpha and delta indicate that they may play a role in basal ovarian function.

Integrating receptor interactions and dynamics and Interference with endocrine disruptors in the mechanisms of action of PPAR nuclear receptors

Abstract Peroxisome Proliferator-Activated Receptors (PPARs) form a family of three nuclear receptors regulating important cellular and metabolic functions. PPARs control gene expression by directly binding to target promoters as heterodimers with the Retinoid X Receptor (RXR), and their transcriptional activity is enhanced upon activation by natural or pharmacological ligands. The binding of PPAR/RXR heterodimers on target promoters allows the anchoring of a series of coactivators and corepressors involved in promoter remodeling and the recruitment of the transcription machinery. The transcriptional output finally depends on a complex interplay between (i) the respective expression levels of PPARs, RXRs and of other nuclear receptors competing for DNA binding and RXR recruitment, (ii) the availability and the nature of PPAR and RXR ligands, (iii) the expression levels and the nature of the different coactivators and corepressors and (iv) the sequence and the epigenetic status of the promoter. Understanding how all these factors and signals integrate and fine-tune transcription remains a challenge but is necessary to understand the specificity of the physiological functions regulated by PPARs. The work presented herein focuses on the molecular mechanisms of PPAR action and aims at understanding how the interactions and mobility of the receptor modulate transcription in the physiological context of a living cell: Such observations in vivo rely on the use of engineered fluorescent protein chimeras and require the development and the application of complementary imaging techniques such as Fluorescence Recovery After Photobleaching (FRAP), Fluorescence Resonance Energy Transfer (FRET) and Fluorescence Correlation Spectroscopy (FCS). Using such techniques, PPARs are shown to reside solely in the nucleus where they are constitutively associated with RXR but transcriptional activation by ligand binding -does not promote the formation of sub-nuclear structures as observed with other nuclear receptors. In addition, the engagement of unliganded PPARs in large complexes of cofactors in living cells provides a molecular basis for their ligand-independent activity. Ligand binding reduces receptor diffusion by promoting the recruitment of coactivators which further enlarge the size of PPAR complexes to acquire full transcriptional competence. Using these molecular approaches, we deciphered the molecular mechanisms through which phthalates, a class of pollutants from the plastic industry, interfere with PPARγ signaling. Mono-ethyl-hexyl-phthalate (MEHP) binding induces the recruitment of a specific subset of cofactors and translates into the expression of a specific subset of target genes, the transcriptional output being strongly conditioned by the differentiation status of the cell. This selective PPARγ modulation induces limited adipogenic effects in cellular models while exposure to phthalates in animal models leads to protective effects on glucose tolerance and diet-induced obesity. These results demonstrate that phthalates influence lipid and carbohydrate metabolism through complex mechanisms which most likely involve PPARγ but also probably PPARα and PPARß, Altogether, the molecular and physiological demonstration of the interference of pollutants with PPAR action outlines an important role of chemical exposure in metabolic regulations. Résumé Les PPARs (Peroxisome Proliferator-Activated Receptors) forment une famille de récepteurs nucléaires qui régulent des fonctions cellulaires et métaboliques importantes. Les PPARs contrôlent l'expression des gènes en se liant directement à leurs promoteurs sous forme d'hétérodimères avec les récepteurs RXR (Retinoid X Receptor), et leur activité transcriptionnelle est stimulée par la liaison de ligands naturels ou pharmacologiques. L'association des hétérodimères PPAR/RXR avec les promoteurs des gènes cibles permet le recrutement de coactivateurs et de corépresseurs qui vont permettre le remodelage de la chromatine et le recrutement de la machinerie transcriptionnelle. Les actions transcriptionnelles du récepteur dépendent toutefois d'interactions complexes qui sont régulées par (i) le niveau d'expression des PPARs, des RXRs et d'autres récepteurs nucléaires entrant en compétition pour la liaison à l'ADN et l'association avec RXR, (ii) la disponibilité et la nature de ligands de PPAR et de RXR, (iii) les niveaux d'expression et la nature des différents coactivateurs et corépresseurs et (iv) la séquence et le marquage épigénétique des promoteurs. La compréhension des mécanismes qui permettent d'intégrer ces aspects pour assurer une régulation fine de l'activité transcriptionnelle est un défi qu'il est nécessaire de relever pour comprendre la spécificité des fonctions physiologiques régulées par les PPARs. Ce travail concerne l'étude des mécanismes d'action moléculaire des PPARs et vise à mieux comprendre comment les interactions du récepteur avec d'autres protéines ainsi que la mobilité de ce dernier régulent son activité transcriptionnelle dans le contexte physiologique des cellules vivantes. De telles observations reposent sur l'emploi de protéines fusionnées à des protéines fluorescentes ainsi que sur le développement et l'utilisation de techniques d'imagerie complémentaires telles que le FRAP (Fluorescence Recovery After Photobleaching), le FRET (Fluorescence Resonance Energy Transfer) ou la FCS (Fluorescence Corrélation Spectroscopy). En appliquant ces méthodes, nous avons pu montrer que les PPARs résident toujours dans le noyau où ils sont associés de manière constitutive à RXR, mais que l'ajout de ligand n'induit pas la formation de structures sub-nucléaires comme cela a pu être décrit pour d'autres récepteurs nucléaires. De plus, les PPARs sont engagés dans de larges complexes protéiques de cofacteurs en absence de ligand, ce qui procure une explication moléculaire à leur activité ligand-indépendante. La liaison du ligand réduit la vitesse de diffusion du récepteur en induisant le recrutement de coactivateurs qui augmente encore plus la taille des complexes afin d'acquérir un potentiel d'activation maximal. En utilisant ces approches moléculaires, nous avons pu caractériser les mécanismes permettant aux phtalates, une classe de polluants provenant de l'industrie plastique, d'interférer avec PPARγ. La liaison du mono-ethyl-hexyl-phtalate (NERF) à PPARγ induit un recrutement sélectif de cofacteurs, se traduisant par l'induction spécifique d'un sous-ensemble de gènes qui varie en fonction du niveau de différentiation cellulaire. La modulation sélective de PPARγ par le MEHP provoque une adipogenèse modérée dans des modèles cellulaires alors que l'exposition de modèles animaux aux phtalates induit des effets bénéfiques sur la tolérance au glucose et sur le développement de l'obésité. Toutefois, les phtalates ont une action complexe sur le métabolisme glucido-lipidique en faisant intervenir PPARγ mais aussi probablement PPARα et PPARß. Cette démonstration moléculaire et physiologique de l'interférence des polluants avec les récepteurs nucléaires PPAR souligne un rôle important de l'exposition à de tels composés dans les régulations métaboliques.

Localization of peroxisome proliferator-activated receptor alpha (PPARα) and N-acyl phosphatidylethanolamine phospholipase D (NAPE-PLD) in cells expressing the Ca(2+)-binding proteins calbindin, calretinin, and parvalbumin in the adult rat hippocampus.

The N-acylethanolamines (NAEs), oleoylethanolamide (OEA) and palmithylethanolamide (PEA) are known to be endogenous ligands of PPARα receptors, and their presence requires the activation of a specific phospholipase D (NAPE-PLD) associated with intracellular Ca(2+) fluxes. Thus, the identification of a specific population of NAPE-PLD/PPARα-containing neurons that express selective Ca(2+)-binding proteins (CaBPs) may provide a neuroanatomical basis to better understand the PPARα system in the brain. For this purpose, we used double-label immunofluorescence and confocal laser scanning microscopy for the characterization of the co-existence of NAPE-PLD/PPARα and the CaBPs calbindin D28k, calretinin and parvalbumin in the rat hippocampus. PPARα expression was specifically localized in the cell nucleus and, occasionally, in the cytoplasm of the principal cells (dentate granular and CA pyramidal cells) and some non-principal cells of the hippocampus. PPARα was expressed in the calbindin-containing cells of the granular cell layer of the dentate gyrus (DG) and the SP of CA1. These principal PPARα(+)/calbindin(+) cells were closely surrounded by NAPE-PLD(+) fiber varicosities. No pyramidal PPARα(+)/calbindin(+) cells were detected in CA3. Most cells containing parvalbumin expressed both NAPE-PLD and PPARα in the principal layers of the DG and CA1/3. A small number of cells containing PPARα and calretinin was found along the hippocampus. Scattered NAPE-PLD(+)/calretinin(+) cells were specifically detected in CA3. NAPE-PLD(+) puncta surrounded the calretinin(+) cells localized in the principal cells of the DG and CA1. The identification of the hippocampal subpopulations of NAPE-PLD/PPARα-containing neurons that express selective CaBPs should be considered when analyzing the role of NAEs/PPARα-signaling system in the regulation of hippocampal functions.

Genetic consequences of the ice ages on nurseries of the bat Myotis myotis: a mitochondrial and nuclear survey.

Analyses of mitochondrial DNA (mtDNA) control region polymorphism and of variation at 10 nuclear microsatellite loci were used to investigate the mechanisms and genetic consequences of postglacial expansion of Myotis myotis in Europe. Initial sampling consisted of 480 bats genotyped in 24 nursery colonies arranged along a transect of approximately 3000 km. The phylogeographical survey based on mtDNA sequences revealed the existence of major genetic subdivisions across this area, with several suture zones between haplogroups. Such zones of secondary contact were found in the Alps and Rhodopes, whereas other potential barriers to gene flow, like the Pyrenees, did not coincide with genetic discontinuities. Areas of population admixture increased locally the genetic diversity of colonies, which confounded the northward decrease in nucleotide diversity predicted using classical models of postglacial range expansion. However, when analyses were restricted to a subset of 15 nurseries originating from a single presumed glacial refugium, mtDNA polymorphism did indeed support a northwards decrease in diversity. Populations were also highly structured (PhiST = 0.384). Conversely, the same subset of colonies showed no significant latitudinal decrease in microsatellite diversity and much less population structure (FST = 0.010), but pairwise genetic differentiation at these nuclear markers was strongly correlated with increasing geographical distance. Together, this evidence suggests that alleles carried via male bats have maintained enough nuclear gene flow to counteract the effects of recurrent bottlenecks generally associated with recolonization processes. As females are highly philopatric, we argue that the maternally transmitted mtDNA marker better reflects the situation of past, historical gene flow, whereas current levels of gene flow are better reflected by microsatellite markers.