Adhesion and protease activity in cell lines from human salivary gland tumors are regulated by the laminin-derived peptide AG73, syndecan-1 and beta 1 integrin

We studied the induction of protease activity by the laminin alpha 1-derived peptide AG73 in cells from adenoid cystic carcinoma (CAC2) and myoepithelioma (M1), respectively a malignant and a benign salivary gland tumors. Laminin alpha 1 chain and MMP9 were immunolocalized in adenoid cystic carcinoma and myoepithelioma in vivo and in vitro. Cells grown inside AG73-enriched laminin-111 exhibited large spaces in the extracellular matrix, suggestive of remodeling. The broad spectrum MMP inhibitor GM6001 decreased spaces induced by AG73 in CAC2 and M I cells. This result strongly suggests that AG73-mediated matrix remodeling involves matrix metalloproteinases. CAC2 and M1 cells cultured on AG73 showed a dose-dependent increase of MMP9 secretion, as detected by zymography. Furthermore, siRNA silencing of MMP9 decreased remodeling in 3D cultures. We searched for AG73 receptors regulating MMP9 activity in our cell lines. CAC2 and M1 cells grown on AG73 exhibited colocalization of syndecan-1 and beta 1 integrin. siRNA knockdown of syndecan-1 expression in these cells resulted in decreased adhesion to AG73 and reduced protease and remodeling activity. We investigated syndecan-1 co-receptors in both cell lines. Silencing beta 1 integrin inhibited adhesion to AG73, matrix remodeling and protease activity. Double-knockdown experiments were carried out to further explore syndecan-1 and beta 1 integrin cooperation. CAC2 cells transfected with both syndecan-1 and beta 1 integrin siRNA oligos showed significant decrease in adhesion to AG73. Simultaneous silencing of receptors also induced a decrease in protease activity. Our results suggest that syndecan-1 and beta 1 integrin signaling downstream of AG73 regulate adhesion and MMP production by CAC2 and M1 cells. (c) 2008 Elsevier B.V./International Society of Matrix Biology. All rights reserved.

Identification of regions within the third FnIII-like domain of the IL-5R alpha involved in IL-5 interaction

Previously, two binding sites for interleukin 5 (IL-5) were identified on the IL-5 receptor alpha chain (IL-5R alpha). They are located within the CD loop of the first fibronectin type III (FnIII)-like domain and the EF loop of the second FnIII-like domain. The first binding site was identified by exploiting the different abilities of human IL-5R alpha (hIL-5R alpha) and mouse IL-5R alpha (mIL-5R alpha) to bind hIL-5. Here we show that ovine IL-5 (oIL-5) has the ability to activate the hIL-5R alpha but not the mIL-5R alpha. By using chimeras of the mIL-5R alpha and hIL-5R alpha we demonstrate that residues within the first and third FnIII-like domains of mIL-5R alpha are responsible for this lack of activity. Furthermore, mutation of residues on hIL-5R alpha to mIL-5R alpha within the predicted DE and FG loop regions of the third FnIII domain reduces oIL-5 activity, These results show that regions of the third FnIII domain of IL-5R alpha are involved in binding, in addition to the regions in domains one and two of the IL-5R alpha that were identified in an earlier study. (C) 2000 Academic Press.

Genome Size, Karyotype Polymorphism and Chromosomal Evolution in Trypanosoma cruzi

Background: The Trypanosoma cruzi genome was sequenced from a hybrid strain (CL Brener). However, high allelic variation and the repetitive nature of the genome have prevented the complete linear sequence of chromosomes being determined. Determining the full complement of chromosomes and establishing syntenic groups will be important in defining the structure of T. cruzi chromosomes. A large amount of information is now available for T. cruzi and Trypanosoma brucei, providing the opportunity to compare and describe the overall patterns of chromosomal evolution in these parasites. Methodology/Principal Findings: The genome sizes, repetitive DNA contents, and the numbers and sizes of chromosomes of nine strains of T. cruzi from four lineages (TcI, TcII, TcV and TcVI) were determined. The genome of the TcI group was statistically smaller than other lineages, with the exception of the TcI isolate Tc1161 (Jose-IMT). Satellite DNA content was correlated with genome size for all isolates, but this was not accompanied by simultaneous amplification of retrotransposons. Regardless of chromosomal polymorphism, large syntenic groups are conserved among T. cruzi lineages. Duplicated chromosome-sized regions were identified and could be retained as paralogous loci, increasing the dosage of several genes. By comparing T. cruzi and T. brucei chromosomes, homologous chromosomal regions in T. brucei were identified. Chromosomes Tb9 and Tb11 of T. brucei share regions of syntenic homology with three and six T. cruzi chromosomal bands, respectively. Conclusions: Despite genome size variation and karyotype polymorphism, T. cruzi lineages exhibit conservation of chromosome structure. Several syntenic groups are conserved among all isolates analyzed in this study. The syntenic regions are larger than expected if rearrangements occur randomly, suggesting that they are conserved owing to positive selection. Mapping of the syntenic regions on T. cruzi chromosomal bands provides evidence for the occurrence of fusion and split events involving T. brucei and T. cruzi chromosomes.

Human Airway Epithelial Cell Responses to Neisseria lactamica and Purified Porin via Toll-Like Receptor 2-Dependent Signaling

The human airway epithelium is constantly exposed to microbial products from colonizing organisms. Regulation of Toll-like receptor (TLR) expression and specific interactions with bacterial ligands is thought to mitigate exacerbation of inflammatory processes induced by the commensal flora in these cells. The genus Neisseria comprises pathogenic and commensal organisms that colonize the human nasopharynx. Neisseria lactamica is not associated with disease, but N. meningitidis occasionally invades the host, causing meningococcal disease and septicemia. Upon colonization of the airway epithelium, specific host cell receptors interact with numerous Neisseria components, including the PorB porin, at the immediate bacterial-host cell interface. This major outer membrane protein is expressed by all Neisseria strains, regardless of pathogenicity, but its amino acid sequence varies among strains, particularly in the surface-exposed regions. The interaction of Neisseria PorB with TLR2 is essential for driving TLR2/TLR1-dependent cellular responses and is thought to occur via the porin`s surface-exposed loop regions. Our studies show that N. lactamica PorB is a TLR2 ligand but its binding specificity for TLR2 is different from that of meningococcal PorB. Furthermore, N. lactamica PorB is a poor inducer of proinflammatory mediators and of TLR2 expression in human airway epithelial cells. These effects are reproduced by whole N. lactamica organisms. Since the responsiveness of human airway epithelial cells to colonizing bacteria is in part regulated via TLR2 expression and signaling, commensal organisms such as N. lactamica would benefit from expressing a product that induces low TLR2-dependent local inflammation, likely delaying or avoiding clearance by the host.

Structural characterization of the urea-unfolded state of Colicin Immunity Protein Im7 and Im9

Dissertação para obtenção do Grau de Mestre em Bioquímica Estrutural e Funcional

A simple RT-PCR-based strategy for screening connexin identity

Vertebrate gap junctions are aggregates of transmembrane channels which are composed of connexin (Cx) proteins encoded by at least fourteen distinct genes in mammals. Since the same Cx type can be expressed in different tissues and more than one Cx type can be expressed by the same cell, the thorough identification of which connexin is in which cell type and how connexin expression changes after experimental manipulation has become quite laborious. Here we describe an efficient, rapid and simple method by which connexin type(s) can be identified in mammalian tissue and cultured cells using endonuclease cleavage of RT-PCR products generated from "multi primers" (sense primer, degenerate oligonucleotide corresponding to a region of the first extracellular domain; antisense primer, degenerate oligonucleotide complementary to the second extracellular domain) that amplify the cytoplasmic loop regions of all known connexins except Cx36. In addition, we provide sequence information on RT-PCR primers used in our laboratory to screen individual connexins and predictions of extension of the "multi primer" method to several human connexins.

Polymorphisme de la capside des papillomavirus appartenant à l’espèce 9 des alphapapillomavirus

Le gène L1 encode pour la protéine majeure de la capside des papillomavirus humains (VPH). L’information relative au polymorphisme de L1 pour les types autres que VPH- 16 est jusqu’ici limitée. Cet ouvrage explore le polymorphisme de L1 en comparant les séquences des types phylogénétiquement apparentés VPH-31, -33, -35 à VPH-16. Des spécimens génitaux recueillis de 732 femmes VIH-séropositives et 323 VIHséronégatives ont été criblés à le recherche d’ADN de VPH par PCR consensus au niveau du gène L1. Les échantillons positifs pour VPH-16 (n=74), -31 (n=74), -33 (n=37) et -35 (n=58) étaient analysés par PCR-séquençage pour la totalité du gène L1. Le nombre de nucléotides substitués pour L1 variait de 19 pour VPH-33 à 52 pour VPH-31. Le rapport du nombre de variantes sur le nombre d’isolats testés était plus élevé pour VPH-31 (56.4%, p=0.05) et VPH-35 (60.3%, p=0.04) comparativement à VPH-16 (40.5%), alors que ce ratio était inférieur pour VPH-33 mais sans différence statistiquement significative (24.3%, p=0.14). La distance entre les variantes était plus grande à l’intérieur des cinq boucles présumément exposées à la surface de la protéine L1 que dans la séquence à l’extérieur (p<0.01) Des variations synonymes étaient observées chez 1.7% (95% CI 1.1- 2.3) des nucléotides intra-boucles et 2.4% (95% CI 1.2-3.7) de ceux extra-boucles. Les variations non-synonymes étaient rencontrées pour 1.8% (95% CI 1.1-2.5) des nucléotides intra-boucles et 0.2% (95% CI 0-0.4) pour les nucléotides extra-boucles. Les ratios dN/dS étaient inférieurs à 1.0 pour les régions extra-boucles et encore davantage pour les régions intra-boucles. Ces résultats suggèrent que les séquences des régions hypervariables de L1 ont été sélectionnées positivement.

Nuclear magnetic resonance structure shows that the severe acute respiratory syndrome coronavirus-unique domain contains a macrodomain fold

The nuclear magnetic resonance (NMR) structure of a central segment of the previously annotated severe acute respiratory syndrome (SARS)-unique domain (SUD-M, for "middle of the SARS-unique domain") in SARS coronavirus (SARS-CoV) nonstructural protein 3 (nsp3) has been determined. SUD-M(513-651) exhibits a macrodomain fold containing the nsp3 residues 528 to 648, and there is a flexibly extended N-terminal tail with the residues 513 to 527 and a C-terminal flexible tail of residues 649 to 651. As a follow-up to this initial result, we also solved the structure of a construct representing only the globular domain of residues 527 to 651 [SUD-M(527-651)]. NMR chemical shift perturbation experiments showed that SUD-M(527-651) binds single-stranded poly(A) and identified the contact area with this RNA on the protein surface, and electrophoretic mobility shift assays then confirmed that SUD-M has higher affinity for purine bases than for pyrimidine bases. In a further search for clues to the function, we found that SUD-M(527-651) has the closest three-dimensional structure homology with another domain of nsp3, the ADP-ribose-1 ''-phosphatase nsp3b, although the two proteins share only 5% sequence identity in the homologous sequence regions. SUD-M(527-651) also shows three-dimensional structure homology with several helicases and nucleoside triphosphate-binding proteins, but it does not contain the motifs of catalytic residues found in these structural homologues. The combined results from NMR screening of potential substrates and the structure-based homology studies now form a basis for more focused investigations on the role of the SARS-unique domain in viral infection.

NMR structure shows that the SARS-unique domain contains a macrodomain fold

The NMR structure of a central segment of the previously annotated "SARS-unique domain" (SUD-M; "middle of the SARS-unique domain") in the SARS coronavirus (SARS-CoV) non-structural protein 3 (nsp3) has been determined. SUD-M(513-651) exhibits a macrodomain fold containing the nsp3-residues 528-648, and there is a flexibly extended N-terminal tail with the residues 513-527 and a C-terminal flexible tail of residues 649-651. As a follow-up to this initial result, we also solved the structure of a construct representing only the globular domain of residues 527-651 [SUD-M(527-651)]. NMR chemical shift perturbation experiments showed that SUD-M(527-651) binds single-stranded poly-A and identified the contact area with this RNA on the protein surface, and electrophoretic mobility shift assays then confirmed that SUD-M has higher affinity for purine bases than for pyrimidine bases. In further search for clues to the function, we found that SUD-M(527-651) has the closest three-dimensional structure homology with another domain of nsp3, the ADP-ribose-1''-phosphatase nsp3b, although the two proteins share only 5% sequence identity in the homologous sequence regions. SUD-M(527-651) also shows 3D structure homology with several helicases and NTP-binding proteins, but it does not contain the motifs of catalytic residues found in these structural homologues. The combined results from NMR screening of potential substrates and the structure-based homology studies now form a basis for more focused investigations on the role of the SARS-unique domain in viral infection.

Chromosome homologies of the highly rearranged karyotypes of four Akodon species (Rodentia, Cricetidae) resolved by reciprocal chromosome painting: the evolution of the lowest diploid number in rodents

Traditionally comparative cytogenetic studies are based mainly on banding patterns. Nevertheless, when dealing with species with highly rearranged genomes, as in Akodon species, or with other highly divergent species, cytogenetic comparisons of banding patterns prove inadequate. Hence, comparative chromosome painting has become the method of choice for genome comparisons at the cytogenetic level since it allows complete chromosome probes of a species to be hybridized in situ onto chromosomes of other species, detecting homologous genomic regions between them. In the present study, we have explored the highly rearranged complements of the Akodon species using reciprocal chromosome painting through species-specific chromosome probes obtained by chromosome sorting. The results revealed complete homology among the complements of Akodon sp. n. (ASP), 2n = 10; Akodon cursor (ACU), 2n = 15; Akodon montensis (AMO), 2n = 24; and Akodon paranaensis (APA), 2n = 44, and extensive chromosome rearrangements have been detected within the species with high precision. Robertsonian and tandem rearrangements, pericentric inversions and/or centromere repositioning, paracentric inversion, translocations, insertions, and breakpoints, where chromosomal rearrangements, seen to be favorable, were observed. Chromosome painting using the APA set of 21 autosomes plus X and Y revealed eight syntenic segments that are shared with A. montensis, A. cursor, and ASP, and one syntenic segment shared by A. montensis and A. cursor plus five exclusive chromosome associations for A. cursor and six for ASP chromosome X, except for the heterochromatin region of ASP X, and even chromosome Y shared complete homology among the species. These data indicate that all those closely related species have experienced a recent extensive process of autosomal rearrangement in which, except for ASP, there is still complete conservation of sex chromosomes homologies.

The endangered species Brycon orbignyanus: genetic analysis and definition of priority areas for conservation

Currently, biodiversity is threatened by several factors often associated with human population growth and the extension of areas occupied by human activity. In particular, freshwater fish fauna is affected by overfishing, deforestation, water pollution, introduction of exotic species and habitat fragmentation promoted by hydroelectric dams, among other environmental impact factors. Several action plans to preserve ichthyofauna biodiversity have been adopted; however, these plans frequently cover only a small number of species, and decisions are often made without strong scientific support. This study aimed to evaluate the genetic aspects of wild groups of Brycon orbignyanus, an endangered fish species, using microsatellites and D-loop regions to identify the genetic structure of the samples and to establish priority areas for conservation based on the genetic patterns of this species. The results indicated that the samples showed levels of genetic variability compatible with others studies with Neotropical fishes. However, the results obtained in the analysis of molecular variance (AMOVA) for microsatellites (F (ST) = 0.258) and D-loop (F (ST) = 0.234) and the interpopulation fixation index revealed that B. orbignyanus was structured in different subpopulations in the La Plata River basin; the areas with better environmental conditions also showed subgroups with higher rates of genetic variability. Future conservation actions addressing these sites should consider two different management units: the complex formed by the Ivinhema River, Upper Parana, Camargo Port and Ilha Grande groups; and the complex formed by the Verde River and Sucuriu River groups.

Ursachen der unterschiedlichen Oligomerisierung von Lichtsammelproteinen

Die verschiedenen Lichtsammelproteine (Lhc-Proteine) höherer Pflanzen unterscheiden sich im Oligomerisierungsverhalten. Im Photosystem II existieren 6 Lhc-Proteine, die entweder die monomeren Lichtsammelkomplexe (LHC) CP24 (Lhcb6), CP26 (Lhcb5) und CP29 (Lhcb4) oder den trimeren LHCII (Lhcb1, Lhcb2 und Lhcb3) bilden. Im Photosystem I sind laut Kristallstruktur vier Lhc-Proteine lokalisiert, die als Heterodimere organisiert vorliegen. Der schwerpunktmäßig untersuchte LHCI-730 setzt sich aus Lhca1 und Lhca4 zusammen, während der LHCI-680 aus Lhca2 und Lhca3 besteht. Das Ziel der Arbeit bestand in der Identifizierung der für das unterschiedliche Oligomerisierungsverhalten verantwortlichen Proteinbereiche und Aminosäuren. Die für diese Arbeit generierten Consensussequenzalignments verschiedener Lhca- und Lhcb-Proteine vieler Arten unterstützen die Folgerungen aus Strukturdaten und anderen Sequenzalignments, dass den LHCs eine gemeinsame Monomerstruktur zu Grunde liegt. Die Helices 1 und 3 weisen weitgehend sehr hohe Sequenzidentitäten auf, während die N- und C-Termini, die zwei Schleifenregionen und die Helix 2 nur schwach konserviert sind. Falls die Bereiche mit hoher Sequenzübereinstimmung für das Zustandekommen ähnlicher monomerer LHC-Strukturen verantwortlich sind, könnten in den schwach konservierten Domänen die Ursachen für das unterschiedliche Oligomerisierungsverhalten lokalisiert sein. Aufgrund dessen wurden die schwach konservierten Domänen des monomerisierenden Lhcb4, des mit dem Lhca1 dimerisierenden Lhca4 und des Trimere bildenden Lhcb1 gegen die entsprechenden Domänen der anderen Proteine ausgetauscht und bezüglich ihres Oligomerisierungsverhaltens untersucht. Im Lhca4 konnten mit der Helix 2 und der stromalen Schleife zwei für eine Heterodimerisierung essentielle Domänen gefunden werden. Im Lhcb1 waren neben dem N-Terminus auch die 2. Helix und die stromale Schleifendomäne unentbehrlich für eine Trimerisierung. Zusätzlich waren Dimerisierung und Trimerisierung bei Austausch der luminalen Schleife beeinträchtigt. Ein geringer Beitrag zur Lhcb1-Trimerisierung konnte auch für den C-Terminus belegt werden. Ein zusätzliches Ziel der Arbeit sollte der Transfer der Oligomerisierungseigenschaften durch umfangreichen Domänentausch von einem auf ein anderes Protein sein. Der Transfer der Fähigkeit zur Dimerbildung durch Substitution gegen essentielle Lhca4-Domänen (50% luminale Schleife, 100% Helix 2 und 100% stromale Schleife) gelang beim Lhcb4, nicht aber beim Lhcb1. Der Transfer der Trimerisierungsfähigkeit auf Lhca4 und Lhcb4 scheiterte. Eine Lhca1-Mutante mit allen für eine Dimerisierung essentiellen Lhca4-Domänen, die durch Interaktion einzelner Moleküle untereinander multimere LHCs bilden sollte, war bereits in ihrer Monomerbildung beeinträchtigt. Eine Übertragung der Oligomerisierungsfähigkeit auf andere Proteine durch massiven Domänentransfer gestaltete sich somit schwierig, da vermutlich im mutierten Protein immer noch ursprüngliche Tertiärstrukturanteile enthalten waren, die nicht mit den transferierten Proteinbestandteilen kompatibel sind. Bei zukünftigen Experimenten zur Klärung der Transferierbarkeit der Oligomerisierungseigenschaft sollten deswegen neben dem unberücksichtigten 1. Teil der luminalen Schleife auch wenig konservierte Aminosäuren in der 1. und 3. Helix Beachtung finden. Ein weiteres Ziel dieser Arbeit war es, die LHCI-730-Dimerisierung im Detail zu untersuchen. Mutationsanalysen bestätigten den von früheren Untersuchungen bekannten Einfluss des Isoleucins 103 und Histidins 99. Letzteres geht möglicherweise durch sein gebundenes Chlorophyll eine Interaktion mit dem Lhca1 ein. Das Phenylalanin 95 stellte sich ebenfalls als ein wichtiger Interaktionspartner heraus und könnte in Wechselwirkung mit einem zwischen Lhca1 und Lhca4 lokalisierten Phosphatidylglycerin treten. Das ebenfalls an der Dimerbildung beteiligte Serin 88 des Lhca4 könnte auf Grund der räumlichen Nähe bei Modellierungen direkt mit dem am C-Terminus des Lhca1 lokalisierten Glycin 190 interagieren. Darüber hinaus wurde ein in der luminalen Lhca4-Schleife lokalisiertes Phenylalanin 84 als Interaktionspartner des Tryptophans 185 im C-Terminus von Lhca1 identifiziert. Der simultane Austausch des Isoleucins 109 und Lysins 110 in der stromalen Schleife des Lhca4, konnte deren Einfluss auf die Dimerisierung belegen. Nachdem bislang an der Dimerbildung beteiligte Aminosäuren am N- und C-Terminus des Lhca1 und Lhca4 identifiziert werden konnten, wurden in dieser Arbeit viele an einer Dimerbildung beteiligten Proteinbereiche und Aminosäuren in der Helix 2 und den Schleifenregionen des Lhca4 identifiziert. Um alle an der Lhca1-Lhca4-Interaktion beteiligten Aminosäuren aufzuklären, müssten durch Mutationsanalysen die in der stromalen Lhca4-Schleife vermuteten Interaktionspartner des für die Dimerisierung wichtigen Tryptophans 4 am N-Terminus von Lhca1 identifiziert, und die in der Helix 3 des Lhca1 vermuteten Interaktionspartner der Helix 2 des Lhca4 ermittelt werden.

Fic domain-catalyzed adenylylation: insight provided by the structural analysis of the type IV secretion system effector BepA

Numerous bacterial pathogens subvert cellular functions of eukaryotic host cells by the injection of effector proteins via dedicated secretion systems. The type IV secretion system (T4SS) effector protein BepA from Bartonella henselae is composed of an N-terminal Fic domain and a C-terminal Bartonella intracellular delivery domain, the latter being responsible for T4SS-mediated translocation into host cells. A proteolysis resistant fragment (residues 10-302) that includes the Fic domain shows autoadenylylation activity and adenylyl transfer onto Hela cell extract proteins as demonstrated by autoradiography on incubation with α-[(32)P]-ATP. Its crystal structure, determined to 2.9-Å resolution by the SeMet-SAD method, exhibits the canonical Fic fold including the HPFxxGNGRxxR signature motif with several elaborations in loop regions and an additional β-rich domain at the C-terminus. On crystal soaking with ATP/Mg(2+), additional electron density indicated the presence of a PP(i) /Mg(2+) moiety, the side product of the adenylylation reaction, in the anion binding nest of the signature motif. On the basis of this information and that of the recent structure of IbpA(Fic2) in complex with the eukaryotic target protein Cdc42, we present a detailed model for the ternary complex of Fic with the two substrates, ATP/Mg(2+) and target tyrosine. The model is consistent with an in-line nucleophilic attack of the deprotonated side-chain hydroxyl group onto the α-phosphorus of the nucleotide to accomplish AMP transfer. Furthermore, a general, sequence-independent mechanism of target positioning through antiparallel β-strand interactions between enzyme and target is suggested.

Structure of neurolysin reveals a deep channel that limits substrate access

The zinc metallopeptidase neurolysin is shown by x-ray crystallography to have large structural elements erected over the active site region that allow substrate access only through a deep narrow channel. This architecture accounts for specialization of this neuropeptidase to small bioactive peptide substrates without bulky secondary and tertiary structures. In addition, modeling studies indicate that the length of a substrate N-terminal to the site of hydrolysis is restricted to approximately 10 residues by the limited size of the active site cavity. Some structural elements of neurolysin, including a five-stranded β-sheet and the two active site helices, are conserved with other metallopeptidases. The connecting loop regions of these elements, however, are much extended in neurolysin, and they, together with other open coil elements, line the active site cavity. These potentially flexible elements may account for the ability of the enzyme to cleave a variety of sequences.

Monophyletic origin and unique dispersal patterns of domestic fowls.

With the aim of elucidating in greater detail the genealogical origin of the present domestic fowls of the world, we have determined mtDNA sequences of the D-loop regions for a total of 21 birds, of which 12 samples belong to red junglefowl (Gallus gallus) comprising three subspecies (six Gallus gallus gallus, three Gallus gallus spadiceus, and three Gallus gallus bankiva) and nine represent diverse domestic breeds (Gallus gallus domesticus). We also sequenced four green junglefowl (Gallus varius), two Lafayette's junglefowl (Gallus lafayettei), and one grey junglefowl (Gallus sonneratii). We then constructed a phylogenetic tree for these birds by the use of nucleotide sequences, choosing the Japanese quail (Coturnix coturnix japonica) as an outgroup. We found that a continental population of G. g. gallus was the real matriarchic origin of all the domestic poultries examined in this study. It is also of particular interest that there were no discernible differences among G. gallus subspecies; G. g. bankiva was a notable exception. This was because G. g. spadiceus and a continental population of G. g. gallus formed a single cluster in the phylogenetic tree. G. g. bankiva, on the other hand, was a distinct entity, thus deserving its subspecies status. It implies that a continental population of G. g. gallus sufficed as the monophyletic ancestor of all domestic breeds. We also discussed a possible significance of the initial dispersal pattern of the present domestic fowls, using the phylogenetic tree.