On the relationship between an Asian haplotype on chromosome 6 that reduces androstenone levels in boars and the differential expression of SULT2A1 in the testis

Background

Androstenone is one of the major compounds responsible for boar taint, a pronounced urine-like odor produced when cooking boar meat. Several studies have identified quantitative trait loci (QTL) for androstenone level on Sus scrofa chromosome (SSC) 6. For one of the candidate genes in the region SULT2A1, a difference in expression levels in the testis has been shown at the protein and RNA level.

Results

Haplotypes were predicted for the QTL region and their effects were estimated showing that haplotype 1 was consistently related with a lower level, and haplotype 2 with a higher level of androstenone. A recombinant haplotype allowed us to narrow down the QTL region from 3.75 Mbp to 1.94 Mbp. An RNA-seq analysis of the liver and testis revealed six genes that were differentially expressed between homozygotes of haplotypes 1 and 2. Genomic sequences of these differentially expressed genes were checked for variations within potential regulatory regions. We identified one variant located within a CpG island that could affect expression of SULT2A1 gene. An allele-specific expression analysis in the testis did not show differential expression between the alleles of SULT2A1 located on the different haplotypes in heterozygous animals. However a synonymous mutation C166T (SSC6: 49,117,861 bp in Sscrofa 10.2; C/T) was identified within the exon 2 of SULT2A1 for which the haplotype 2 only had the C allele which was higher expressed than the T allele, indicating haplotype-independent allelic-imbalanced expression between the two alleles. A phylogenetic analysis for the 1.94 Mbp region revealed that haplotype 1, associated with low androstenone level, originated from Asia.

Conclusions

Differential expression could be observed for six genes by RNA-seq analysis. No difference in the ratio of C:T expression of SULT2A1 for the haplotypes was found by the allele-specific expression analysis, however, a difference in expression between the C over T allele was found for a variation within SULT2A1, showing that the difference in androstenone levels between the haplotypes is not caused by the SNP in exon 2.     

Mitochondrial DNA and Y‐chromosome diversity in East Adriatic sheep

Variation in mitochondrial DNA (mtDNA) and Y‐chromosome haplotypes was analysed in nine domestic sheep breeds (159 rams) and 21 mouflon ( Ovis musimon) sampled in the East Adriatic. Mitochondrial DNA analyses revealed a high frequency of type B haplotypes, predominantly in European breeds, and a very low frequency of type A haplotypes, which are more frequent in some Asian breeds. Mitochondrial haplotype Hmt‐3 was the most frequent (26.4%), and 37.1%, 20.8% and 7.6% of rams had haplotypes one, two and three mutations remote from Hmt‐3 respectively. In contrast, Y‐chromosome analyses revealed extraordinary paternal allelic richness: HY‐6, 89.3%; HY‐8, 5.0%; HY‐18, 3.1%; HY‐7, 1.3%; and HY‐5, 1.3%. In fact, the number of haplotypes observed is comparable to the number found in Turkish breeds and greater than the number found in European breeds so far. Haplotype HY‐18 (A‐oY1/135‐SRYM18), identified here for the first time, provides a link between the haplotype HY‐12 (A‐oY1/139‐SRYM18) found in a few rams in Turkey and haplotype HY‐9 (A‐oY1/131‐SRYM18) found in one ram in Ethiopia. All mouflons had type B mtDNA haplotypes, including the private haplotype (Hmt‐55), and all were paternally monomorphic for haplotype HY‐6. Our data support a quite homogeneous maternal origin of East Adriatic sheep, which is a characteristic of European breeds. At the same time, the high number of haplotypes found was surprising and intriguing, and it begs for further analysis. Simultaneous analysis of mtDNA and Y‐chromosome information allowed us to detect a large discrepancy between maternal and paternal lineages in some populations. This is most likely the result of breeder efforts to ‘upgrade’ local populations using rams with different paternal origins.     

Haplotype hitchhiking promotes trait coselection in Brassica napus

Local haplotype patterns surrounding densely spaced DNA markers with significant trait associations can reveal information on selective sweeps and genome diversity associated with important crop traits. Relationships between haplotype and phenotype diversity, coupled with analysis of gene content in conserved haplotype blocks, can provide insight into coselection for nonrelated traits. We performed genome‐wide analysis of haplotypes associated with the important physiological and agronomic traits leaf chlorophyll and seed glucosinolate content, respectively, in the major oilseed crop species Brassica napus. A locus on chromosome A01 showed opposite effects on leaf chlorophyll content and seed glucosinolate content, attributed to strong linkage disequilibrium (LD) between orthologues of the chlorophyll biosynthesis genes EARLY LIGHT‐INDUCED PROTEIN and CHLOROPHYLL SYNTHASE, and the glucosinolate synthesis gene ATP SULFURYLASE 1. Another conserved haplotype block, on chromosome A02, contained a number of chlorophyll‐related genes in LD with orthologues of the key glucosinolate biosynthesis genes METHYLTHIOALKYMALATE SYNTHASE‐LIKE 1 and 3. Multigene haplogroups were found to have a significantly greater contribution to variation for chlorophyll content than haplotypes for any single gene, suggesting positive effects of additive locus accumulation. Detailed reanalysis of population substructure revealed a clade of ten related accessions exhibiting high leaf chlorophyll and low seed glucosinolate content. These accessions each carried one of the above‐mentioned haplotypes from A01 or A02, generally in combination with further chlorophyll‐associated haplotypes from chromosomes A05 and/or C05. The phenotypic rather than pleiotropic correlations between leaf chlorophyll content index and seed GSL suggest that LD may have led to inadvertent coselection for these two traits.     

Genome‐wide prediction of age at puberty and reproductive longevity in sows

Traditional selection for sow reproductive longevity is ineffective due to low heritability and late expression of the trait. Incorporation of DNA markers into selection programs is potentially a more practical approach for improving sow lifetime productivity. Using a resource population of crossbred gilts, we explored pleiotropic sources of variation that influence age at puberty and reproductive longevity. Of the traits recorded before breeding, only age at puberty significantly affected the probability that females would produce a first parity litter. The genetic variance explained by 1‐Mb windows of the sow genome, compared across traits, uncovered regions that influence both age at puberty and lifetime number of parities. Allelic variants of SNPs located on SSC5 (27–28 Mb), SSC8 (36–37 Mb) and SSC12 (1.2–2 Mb) exhibited additive effects and were associated with both early expression of puberty and a greater than average number of lifetime parities. Combined analysis of these SNPs showed that an increase in the number of favorable alleles had positive impact on reproductive longevity, increasing number of parities by up to 1.36. The region located on SSC5 harbors non‐synonymous alleles in the arginine vasopressin receptor 1A (AVPR1A) gene, a G‐protein‐coupled receptor associated with social and reproductive behaviors in voles and humans and a candidate for the observed effects. This region is characterized by high levels of linkage disequilibrium in different lines and could be exploited in marker‐assisted selection programs across populations to increase sow reproductive longevity.     

Maternal genealogical patterns of chicken breeds sampled in Europe

The aim of this study was to investigate the maternal genealogical pattern of chicken breeds sampled in Europe. Sequence polymorphisms of 1256 chickens of the hypervariable region (D‐loop) of mitochondrial DNA (mtDNA) were used. Median‐joining networks were constructed to establish evolutionary relationships among mtDNA haplotypes of chickens, which included a wide range of breeds with different origin and history. Chicken breeds which have had their roots in Europe for more than 3000 years were categorized by their founding regions, encompassing Mediterranean type, East European type and Northwest European type. Breeds which were introduced to Europe from Asia since the mid‐19th century were classified as Asian type, and breeds based on crossbreeding between Asian breeds and European breeds were classified as Intermediate type. The last group, Game birds, included fighting birds from Asia. The classification of mtDNA haplotypes was based on Liu et al.'s (2006) nomenclature. Haplogroup E was the predominant clade among the European chicken breeds. The results showed, on average, the highest number of haplotypes, highest haplotype diversity, and highest nucleotide diversity for Asian type breeds, followed by Intermediate type chickens. East European and Northwest European breeds had lower haplotype and nucleotide diversity compared to Mediterranean, Intermediate, Game and Asian type breeds. Results of our study support earlier findings that chicken breeds sampled in Europe have their roots in the Indian subcontinent and East Asia. This is consistent with historical and archaeological evidence of chicken migration routes to Europe.     

Long‐term selection for litter size in swine results in shifts in allelic frequency in regions involved in reproductive processes

High‐density genotype data were analyzed in three lines of swine that express substantial variation in sow fertility to uncover regions of the genome potentially influenced during selection for litter size traits. The experimental lines examined include the Nebraska Index Line (NIL), which has been subjected to long‐term selection for litter size; a control line derived from the same population that founded NIL; and a commercial Duroc × Hampshire (D × H) population, in which no selection for litter size was practiced. Regions of the genome potentially affected by selection for litter size traits in NIL were determined by multiple lines of evidence, including altered allelic frequency compared to the other lines, loss of heterozygosity and relative extended haplotype homozygosity. Additionally, a genome‐wide association study for litter size traits was conducted in a population based on NIL and commercial maternal line genetics. Several genomic regions identified as putative signatures of selection overlapped with QTL for litter size traits. One of these regions, located on SSC2 (13–14 Mb), includes the candidate gene P2X3R, which plays a role in implantation and sustained release of hormones associated with reproductive processes. Sequencing identified synonymous SNPs in P2X3R that are fixed in NIL but polymorphic with nearly equal frequencies in the D × H line, indicating a potential role of P2X3R in sow fertility. These results suggest that data derived from these lines can help to uncover and understand a portion of the genetic variance associated with fertility traits in swine.     

Geographically congruent large‐scale patterns of plastid haplotype variation in the European herbs Silene dioica and S. latifolia (Caryophyllaceae)

The closely related dioecious herbs Silene latifolia and Silene dioica are widespread and predominantly sympatric in Europe. The species are interfertile, but morphologically and ecologically distinct. A study of large‐scale patterns of plastid DNA (polymerase chain reaction–restriction fragment length polymorphism) haplotypes in a sample of 198 populations from most of the European ranges of both species revealed extensive interspecific haplotype sharing. Four of the 28 detected haplotypes were frequent (found in > 40 populations) and widespread. Three of these frequent haplotypes occurred in both species and the geographic distribution of each haplotype was broadly congruent in both species. Each of these three, shared and widespread haplotypes is likely to have colonized central and/or northern Europe after the last glaciation from one or more of refugial areas in southern Europe. Interspecific hybridization and plastid introgression within refugial regions and/or during the early stages of postglacial expansion is the most plausible explanation for the broadly similar distribution patterns of the shared, frequent chloroplast haplotypes in the two species. The fourth frequent, widespread haplotype was absent from S. latifolia and almost entirely restricted to Nordic S. dioica. It is most likely that this haplotype spread into the Nordic countries from a central or northern European source or from a refugial area in Russia. © 2009 The Linnean Society of London, Botanical Journal of the Linnean Society, 2009, 161 , 153–170.     

Genotyping of five single nucleotide polymorphisms in the OCA2 and HERC2 genes associated with blue‐brown eye color in the Japanese population

Human eye color is a polymorphic phenotype influenced by multiple genes. It has recently been reported that three single nucleotide polymorphisms (SNPs) within intron 1 of the OCA2 gene (rs7495174, rs4778241, rs4778138) and two SNPs in intron 86 (rs12913832) and the 3′ UTR region (rs1129038) of the HERC2 gene—located in the upstream of the OCA2 locus —have a high statistical association with human eye color. The present study is the first to examine in detail the genotype and haplotype frequencies for these five SNPs in an Asian (Japanese) population (n = 523) comprising solely brown‐eyed individuals. Comparison of the genotype and haplotype distributions in Japanese with those in African and European subjects revealed significant differences between Japanese and other populations. Analysis of haplotypes consisting of four SNPs at the HERC2‐OCA2 locus (rs12913832/rs7495174/rs4778241/rs4778138) showed that the most frequent haplotype in the Japanese population is A‐GAG (0.568), while the frequency of this haplotype is rather low in the European population, even in the brown‐eyed group (0.167). The haplotype distribution in the Japanese population was significantly different from that in the brown‐eyed European group (F_ST = 0.18915). Copyright © 2009 John Wiley & Sons, Ltd.     

TYR基因外显子1的序列变异

为了分析家猪与野猪的遗传多样性及起源,测定了来自１２个中国地方家猪品种、３个欧洲引进猪品种以及8个中国野猪和2个越南野猪共36个个体的酪氨酸酶基因（TYR）外显子1的序列,共检出６个单核苷酸多态性位点（SNPs）,且这6个位点的变异均为同义突变,根据这些变异可将酪氨酸酶基因DNA序列归结为４种单倍型。结合已发表的数据,构建了简约中介网络图。 在网络图中,单倍型TYR*2主要为欧洲家猪与欧洲野猪和三条亚洲家猪染色体。大部分亚洲家猪和野猪共享单倍型TYR*1,表明这是一个亚洲类型的单倍型;同时也有部分欧洲家猪与野猪携带这一单倍型。 而单倍型TYR*3和TYR*4为本研究检测到的稀有单倍型,这两种单倍型主要由中国家猪与亚洲野猪组成。这种网络图结构支持家猪的欧洲和亚洲独立起源学说,同时也表明相当部分的欧洲家猪品种受到亚洲猪的基因渗透,而少量中国家猪和日本野猪也受到了欧洲猪的基因渗透。     

Paternal phylogeography and genetic diversity of East Asian goats

This study was a first analysis of paternal genetic diversity for extensive Asian domestic goats using SRY gene sequences. Sequencing comparison of the SRY 3′‐untranslated region among 210 Asian goats revealed four haplotypes (Y1A, Y1B, Y2A and Y2B) derived from four variable sites including a novel substitution detected in this study. In Asian goats, the predominant haplotype was Y1A (62%) and second most common was Y2B (30%). Interestingly, the Y2B was a unique East Asian Y chromosomal variant, which differentiates eastern and western Eurasian goats. The SRY geographic distribution in Myanmar and Cambodia indicated predominant the haplotype Y1A in plains areas and a high frequency of Y2B in mountain areas. The results suggest recent genetic infiltration of modern breeds into South‐East Asian goats and an ancestral SRY Y2B haplotype in Asian native goats.     

Untangling the hybrid nature of modern pig genomes: a mosaic derived from biogeographically distinct and highly divergent Sus scrofa populations

The merging of populations after an extended period of isolation and divergence is a common phenomenon, in natural settings as well as due to human interference. Individuals with such hybrid origins contain genomes that essentially form a mosaic of different histories and demographies. Pigs are an excellent model species to study hybridization because European and Asian wild boars diverged ~1.2 Mya, and pigs were domesticated independently in Europe and Asia. During the Industrial Revolution in England, pigs were imported from China to improve the local pigs. This study utilizes the latest genomics tools to identify the origin of haplotypes in European domesticated pigs that are descendant from Asian and European populations. Our results reveal fine‐scale haplotype structure representing different ancient demographic events, as well as a mosaic composition of those distinct histories due to recently introgressed haplotypes in the pig genome. As a consequence, nucleotide diversity in the genome of European domesticated pigs is higher when at least one haplotype of Asian origin is present, and haplotype length correlates negatively with recombination frequency and nucleotide diversity. Another consequence is that the inference of past effective population size is influenced by the background of the haplotypes in an individual, but we demonstrate that by careful sorting based on the origin of haplotypes, both distinct demographic histories can be reconstructed. Future detailed mapping of the genomic distribution of variation will enable a targeted approach to increase genetic diversity of captive and wild populations, thus facilitating conservation efforts in the near future.     

Genomewide analysis of admixture and adaptation in the Africanized honeybee

Genetic exchange by hybridization or admixture can make an important contribution to evolution, and introgression of favourable alleles can facilitate adaptation to new environments. A small number of honeybees (Apis mellifera) with African ancestry were introduced to Brazil ~60 years ago, which dispersed and hybridized with existing managed populations of European origin, quickly spreading across much of the Americas in an example of a massive biological invasion. Here, we analyse whole‐genome sequences of 32 Africanized honeybees sampled from throughout Brazil to study the effect of this process on genome diversity. By comparison with ancestral populations from Europe and Africa, we infer that these samples have 84% African ancestry, with the remainder from western European populations. However, this proportion varies across the genome and we identify signals of positive selection in regions with high European ancestry proportions. These observations are largely driven by one large gene‐rich 1.4‐Mbp segment on chromosome 11 where European haplotypes are present at a significantly elevated frequency and likely confer an adaptive advantage in the Africanized honeybee population. This region has previously been implicated in reproductive traits and foraging behaviour in worker bees. Finally, by analysing the distribution of ancestry tract lengths in the context of the known time of the admixture event, we are able to infer an average generation time of 2.0 years. Our analysis highlights the processes by which populations of mixed genetic ancestry form and adapt to new environments.     

A morphological and genetic analysis of the European bitterling species complex

Bitterling fishes lay their eggs on the gills of living freshwater mussels and are valuable models in behavioural and evolutionary ecology. We used morphological and genetic data to resolve the taxonomic relationships of bitterling in Europe. Previous studies have suggested the European bitterling is either a single species with a wide but discontinuous geographic distribution, or a complex of species. Morphometric and meristic data identified differences between three putative species; with a clear distinction between the eastern Asian Rhodeus sericeus, western European bitterling Rhodeus amarus, and colchian bitterling, Rhodeus colchicus. Polymorphism in the mitochondrial DNA control region was predominantly due to insertion/deletion events, making phylogenetic inference difficult, but the single haplotype found in R. sericeus populations was detected at low frequency (one of 24 individuals) in R. amarus and R. colchicus populations. Eight control region haplotypes were found in R. amarus populations, which were distinct from the two haplotypes in a R. colchicus population. Cytochrome b data produced a phylogeny with strongly‐supported differentiation between a clade of two R. sericeus haplotypes and a clade of six R. amarus/colchicus haplotypes. The star‐like topology of the R. amarus/colchicus haplotypes in a minimum spanning network suggested a rapid radiation in this clade. Our results are consistent with an hypothesis of relatively ancient divergence of R. sericeus from R. amarus/colchicus and more recent and rapid differentiation between R. amarus and R. colchicus. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 337–347.     

Detecting selected haplotype blocks in evolve and resequence experiments

Shifting from the analysis of single nucleotide polymorphisms to the reconstruction of selected haplotypes greatly facilitates the interpretation of evolve and resequence (E&R) experiments. Merging highly correlated hitchhiker SNPs into haplotype blocks reduces thousands of candidates to few selected regions. Current methods of haplotype reconstruction from Pool‐seq data need a variety of data‐specific parameters that are typically defined ad hoc and require haplotype sequences for validation. Here, we introduce haplovalidate, a tool which detects selected haplotypes in Pool‐seq time series data without the need for sequenced haplotypes. Haplovalidate makes data‐driven choices of two key parameters for the clustering procedure, the minimum correlation between SNPs constituting a cluster and the window size. Applying haplovalidate to simulated E&R data reliably detects selected haplotype blocks with low false discovery rates. Importantly, our analyses identified a restriction of the haplotype block‐based approach to describe the genomic architecture of adaptation. We detected a substantial fraction of haplotypes containing multiple selection targets. These blocks were considered as one region of selection and therefore led to underestimation of the number of selection targets. We demonstrate that the separate analysis of earlier time points can significantly increase the separation of selection targets into individual haplotype blocks. We conclude that the analysis of selected haplotype blocks has great potential for the characterization of the adaptive architecture with E&R experiments.     

Geographic pattern of genetic variation in the fox tapeworm Echinococcus multilocularis

Intraspecific genetic variation of Echinococcus multilocularis, the etiologic agent of human alveolar echinococcosis, has been evaluated among 76 geographic isolates from Europe, Asia and North America by using sequence data of mitochondrial and nuclear DNA. Relatively low genetic variation was found only in the mitochondrial DNA sequence consisting of 3 protein-coding genes. Pairwise divergence among the resultant 18 haplotypes ranged from 0.03 to 1.91%. Phylogenetic trees and parsimony network of these haplotypes depicted a geographic division into European, Asian and North American clades, but 1 haplotype from Inner Mongolia was unrelated to other haplotypes. The coexistence of the Asian and North American haplotypes could be seen, particularly on the St. Lawrence Island in the Bering Sea. These data suggest an evolutionary scenario in which distinct parasite populations derived from glacial refugia have been maintained by indigenous host mammals. The nuclear DNA sequence for the immunodominant B cell epitope region of ezrin/radixin/moesin-like protein (elp) was extremely conservative, indicating that the elp antigen is available for immunodiagnosis in any endemic areas.     

Y chromosome haplotype diversity of domestic sheep (Ovis aries) in northern Eurasia

Variation in two SNPs and one microsatellite on the Y chromosome was analyzed in a total of 663 rams representing 59 breeds from a large geographic range in northern Eurasia. SNPA‐oY1 showed the highest allele frequency (91.55%) across the breeds, whereas SNPG‐oY1 was present in only 56 samples. Combined genotypes established seven haplotypes (H4, H5, H6, H7, H8, H12 and H19). H6 dominated in northern Eurasia, and H8 showed the second‐highest frequency. H4, which had been earlier reported to be absent in European breeds, was detected in one European breed (Swiniarka), whereas H7, which had been previously identified to be unique to European breeds, was present in two Chinese breeds (Ninglang Black and Large‐tailed Han), one Buryatian (Transbaikal Finewool) and two Russian breeds (North Caucasus Mutton‐Wool and Kuibyshev). H12, which had been detected only in Turkish breeds, was also found in Chinese breeds in this work. An overall low level of haplotype diversity (median h = 0.1288) was observed across the breeds with relatively higher median values in breeds from the regions neighboring the Near Eastern domestication center of sheep. H6 is the dominant haplotype in northwestern and eastern China, in which the haplotype distribution could be explained by the historical translocations of the H4 and H8 Y chromosomes to China via the Mongol invasions followed by expansions to northwestern and eastern China. Our findings extend previous results of sheep Y chromosomal genetic variability and indicate probably recent paternal gene flows between sheep breeds from distinct major geographic regions.     

Genome‐wide association study reveals regions associated with gestation length in two pig populations

Reproduction traits, such as gestation length (GLE), play an important role in dam line breeding in pigs. The objective of our study was to identify single nucleotide polymorphisms (SNPs) that are associated with GLE in two pig populations. Genotypes and deregressed breeding values were available for 2081 Dutch Landrace‐based (DL) and 2301 Large White‐based (LW) pigs. We identified two QTL regions for GLE, one in each population. For DL, three associated SNPs were detected in one QTL region spanning 0.52 Mbp on Sus scrofa chromosome (SSC) 2. For LW, four associated SNPs were detected in one region of 0.14 Mbp on SSC5. The region on SSC2 contains the heparin‐binding EGF‐like growth factor (HBEGF) gene, which promotes embryo implantation and has been described to be involved in embryo survival throughout gestation. The associated SNP can be used for marker‐assisted selection in the studied populations, and further studies of the HBEGF gene are warranted to investigate its role in GLE.     

Phylogenetic relationships of Eurasian Nuthatches (Sitta europaea Linnaeus, 1758) from the Alborz and Zagros Mountains,Iran

The Eurasian Nuthatch (Sitta europaea Linnaeus, 1758) is a resident bird in the Alborz and Zagros deciduous forests. We investigated the phylogenetic relationships and the taxonomic position of the Eurasian Nuthatch among other separated lineages of Eurasia with the help of blood samples collected from 19 individuals belonging to four populations in the Eastern and Western Alborz, as well as in the Northern and Southern Zagros forests. Genetic variation was then analysed using complete ND2 gene sequence (1041?bp) and phylogenetic analysis was done using Bayesian and maximum likelihood inference. Additionally, a median-joining algorithm was used to reveal the relationships among haplotypes. The results of the phylogenetic and haplo-type network analyses indicated that Eurasian Nuthatch haplotypes from the Alborz and Zagros Mountains form lineages distinct from the Asian, Caucasian and European haplotypes. Furthermore, an analysis of molecular variance (AMOVA) detected significant (P<0.001) genetic structure among the lineages. The Asian, European, Caucasian and Alborz lineages diverged from one another by an uncorrected genetic distance ranging from 0.029 to 0.039, while the Zagros lineage showed a slightly lower genetic divergence from the Caucasian lineage (0.006), but it did not share any haplotype with the Caucasian lineage. Thus, we suggest considering five Conservation Significant Units (CSU) for the Eurasian Nuthatches as the result of used dataset.     

Haplotypes of <Emphasis Type="Italic">Fallopia</Emphasis> introduced into the US

In the US, clonal growth of Fallopia japonica, Fallopia sachalinensis and their hybrid Fallopia x bohemica (Polygonaceae) is prominent, yet sexual reproduction and hybridization contribute to the genetic complexity of swarms. The contribution to this diversity from multiple introductions is unknown. Using 800 bp of the non-coding chloroplast marker accD–rbcL, we compared 21 Japanese haplotypes with 46 US samples from 11 states, 2 Canadian samples, and 6 European samples from 4 countries, in order to investigate if there were repeated introductions from Asia. While most North American and all European haplotypes accessions in our collection matched a single widespread haplotype, we identified 8 other haplotypes. Three haplotypes of F. japonica (including the widespread haplotype) and one F. sachalinensis matched previously identified Japanese haplotypes, supporting the hypothesis of multiple introductions in the US. Five additional US haplotypes were detected once. Four of these differed from Japanese haplotypes by one single nucleotide polymorphism (SNP), possibly indicating a recent in situ change. The fifth haplotype represents a garden cultivar, which differed from all F. japonica haplotypes. It therefore appears that the US genetic diversity of these taxa has three sources: intra-specific reproduction, inter-specific reproduction, and multiple sources of introduction.