African-Derived Mitochondria in South American Native Cattle Breeds (Bos taurus): Evidence of a New Taurine Mitochondrial Lineage
|
|
- Howard Wade
- 5 years ago
- Views:
Transcription
1 African-Derived Mitochondria in South American Native Cattle Breeds (Bos taurus): Evidence of a New Taurine Mitochondrial Lineage M. M. Miretti, H. A. Pereira, Jr., M. A. Poli, E. P. B. Contel, and J. A. Ferro This article reports the nucleotide diversity within the control region of 42 mitochondrial chromosomes belonging to five South American native cattle breeds (Bos taurus). Analysis of these data in conjunction with B. taurus and B. indicus sequences from Africa, Europe, the Near East, India, and Japan allowed the recognition of eight new mitochondrial haplotypes and their relative positions in a phylogenetic network. The structure of genetic variation among different hypothetical groupings was tested through the molecular variance decomposition, which was best explained by haplotype group components. Haplotypes surveyed were classified as European-related and African-related. Unexpectedly, two haplotypes within the African cluster were more divergent from the African consensus than the latter from the European consensus. A neighbor-joining tree shows the position of two haplotypes compared to European/African mitochondrial lineage splitting. This different and putatively ancestral mitochondrial lineage (AA) is supported by the calibration of sequence divergence based on the Bos Bison separation. The European/African mitochondria divergence might be subsequent (67,100 years before present) to that between AA and Africans (84,700 years before present), also preceding domestication times. These genetic data could reflect the haplotype distribution of Iberian cattle five centuries ago. From the Departmento de Tecnologia, FCAV, UNESP, Via de Acesso Prof. P. D. Castellane km 5, Jaboticabal, SP, Brazil (Miretti, Ferro, and Pereira); Departmento de Genética, FMRP, USP, Av. Bandeirantes 3900, , Ribeirão Preto, SP, Brazil (Miretti and Contel); and Instituto de Genética, CNIA-INTA, Castelar, Argentina (Poli). The authors would like to thank F. Holgado, Andre Carvalho, and Dr. Junqueira for allowing blood collection of an Argentinean Creole, Caracu, and Mocho Nacional population, respectively; Dr. F. Meirelles for providing Nellore DNA; Dr. M. Naves for valuable comments; Dr. D. Bradley for providing the Portuguese cattle mtdna sequences, Dr. J. Sereno, CPAP/ EMBRAPA, and Dr. M. Lara, Pantaneiro and Curraleiro blood samples. This work was funded by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) and Fundação de Apoio ao Ensino, Pesquisa e Assitência do HC da FMRP, USP (FAEPA). Address correspondence to M. M. Miretti at the address above, or jesus@ fcav.unesp.br. Ó 2002 The American Genetic Association 93: The comprehension of phenotypic differences and genetic variations among extant cattle world populations provided clues about their origin, evolutionary history, and in particular their relationship with early human civilizations. Wild aurochs (Bos primigenius), spread over the whole Near East and Europe in the post-pleistocene period and domesticated by Neolithic civilizations in Western Asia, have been considered the common ancestors of all modern cattle breeds (Epstein and Mason 1984; Payne 1991). The analysis of European B. primigenius DNA samples from tissues dated 12,000 years before present (YBP) showed that B. primigenius sequences were phylogenetically closer to extant B. taurus than to B. indicus (Bailey et al. 1996). The extent of divergence between zebuine and taurine cattle, based on variations within the control region of mitochondrial DNA (mtdna) between Bos and Bison, was estimated as 740,000 YBP (Loftus et al. 1994), well before domestication times. This was interpreted as evidence for a predomestic and separate origin of zebuine (humped) and taurine (humpless) cattle population ancestors. Taurine mitochondrial diversity can be grouped within two major mitochondrial haplotype clusters. There is a single highly frequent haplotype in each group representing European and African taurine cattle breeds, namely the African consensus (Afcons) and the European consensus (Eucons), respectively (Bradley et al. 1996). Both mitochondrial types differ in three substitutions within a 240 bp DNA fragment of the control region and are linked to a number of variants through one or a few mutation steps. European extant cattle, including the Iberian peninsula breeds, are presumably descendants of the Near East B. primigenius populations; nevertheless, mitochondrial haplotypes of African origin have also been found in Portuguese cattle breeds (Cymbron et al. 1999). Consequently Iberian cattle appear to have received influence from eastern Sahara B. taurus descents before the spreading of B. indicus in Africa (Mac- Hugh et al. 1997), which is consistent with the demographic history of the Iberian peninsula. Given that all the 323
2 American native cattle have their roots in the Iberian cattle (B. taurus), one would expect to find mitochondria of European as well as of African origin among them. However, current genetic and historical data do not permit us to discern whether those animals first introduced in South America were exclusively European taurus or if African taurus were already present. In any case, mtdna information from South American native breeds could tell us if the Eucons and the Afcons have a similar distribution pattern in America. During the first two decades of the colonization of America, Europeans introduced a relatively small number of animals (Ne;200) (Rabasa 1993; Rouse 1977). Cattle spread throughout South America 50 years later (Salazar and Cardozo 1986), eventually suffering population contraction and expansion periods (Giberti 1974; Hernández 1881). Given this scenario, one would expect to observe a reduced number of mitochondrial haplotypes shared by South American native cattle breeds. With the aim to uncover the source of mitochondrial origin and the nature of mtdna diversity, we examined the nucleotide variation within the mitochondrial control region in Argentinean and Brazilian native cattle breeds as representatives of the South American descendants of the Iberian cattle. We also investigated the contribution of European and African cattle breeds to the formation of the native cattle herds in these two countries. Materials and Methods Animals and DNA Extraction Genomic DNA was extracted from blood or semen samples from 42 animals according to Sambrook et al. (1995) and Zadworny and Kuhnlein (1990) and was used for polymerase chain reaction (PCR) amplification and sequencing. These 42 animals represent the Argentinean Creole (AC, n 5 12) and four Brazilian native cattle breeds: Caracu (CC, n 5 8), Curraleiro (CU, n 5 6), Mocho Nacional (MN, n 5 6), and Pantaneiro (PP, n 5 10). Unrelated AC males were carefully selected based on Breed Association records and pedigree information to represent different AC strains distributed throughout Argentina. Animals from the Brazilian breeds were chosen according to herd records, avoiding relatives except in PP, where no reproductive records are available due to the management nature of this breed (Mazza et al. 1994). mtdna Amplification and Sequencing A 1.39 kb DNA fragment covering the whole mitochondrial control region was amplified by PCR using 200 ng of bovine genomic DNA, 10 pmol each primer (59- TTCCGACCACTCAGCCAATG-39 and 59- CCTAGAGGGCATTCTCACTGGG-39), 2.5 U Taq DNA polymerase, 200 lm of each dntp, 1.5 mm MgCl 2, and 13 reaction buffer. Amplification reactions were performed in a PE9600 cycler (Perkin-Elmer Cetus, Norwalk, CT) starting with a 3 min denaturation step at 948C, followed by 30 cycles of 1 min at 948C, 45 s at 608C, 1 min at 728C, and a 4 min final extension period at 728C. A total of ng of purified PCR product was used as template for sequencing reactions with the following primers: 59-TTCCGACCACT- CAGCCAATG-39, 59-TGCTGGTGCTCAA- GATGC-39, and 59-GCTCGTGATCTAATG- GTAAG-39. A consensus sequence of approximately 800 bp for each animal resulted from contig assembling of sequence reads in both strands. Point mutations were confirmed through chromatogram inspection and by at least two sequence reads covering these positions. PCR and sequencing primers were designed based on the alignment of all Bos mtdna sequences recovered from the GenBank. However, as no B. indicus coding sequence bordering the control region was available, we obtained the nucleotide sequence of CYTB thr trna pro trna and phe trna 12rRNA genes on the upstream and downstream side of the control region, respectively, from a B. indicus animal (Nellore). Nucleotide positions at the mtdna control region were numbered according to the B. taurus mtdna reference sequence (Anderson et al. 1982) (accession no. NC001567). Sequence Analysis mtdna nucleotide sequences obtained in this work (GenBank accession nos. AF AF517828) were aligned to B. taurus, B. indicus, and Bison bison sequences recovered from the GenBank (accession nos. AB AB003801, AF AF016063, AF AF016097, AF AF022924, AF AF034446, AF AF336748, AF AF083354, AF AF162485, BBU12946 BBU12959, L27714 L27715, L27720 L27723, L27728 L27733, L27736 L27737, U51806 U51842, U87633 U87650, U87893, U92230 U92244) using CLUSTALX software (Thompson et al. 1997) totaling 545 Bos mtdna sequences. Positions in the alignment showing gaps were excluded from the analysis. Interhaplotypic distances were estimated using the substitution model of Tamura and Nei (1993), allowing rate heterogeneity among sites. The value for the rate heterogeneity parameter a of the gamma distribution was estimated from the 240 bp ( of the reference sequence) fragment of the control region mtdna sequences from B. taurus and B. indicus using the PUZZLE program (Strimmer and von Haeseler 1996) according to recommendations of Meyer et al. (1999). The corrected average population pairwise sequence difference [p 5 p xy 2 (p x 1 p y )/2] was obtained as implemented by the ARLEQUIN software (Schneider et al. 2000). Analysis of molecular variance (AMOVA), pairwise F ST distance estimates, and the construction of a minimum spanning network (MSN) were also performed using ARLEQUIN. Additional phylogenies were constructed using the NEIGHBOR program incorporated in the PHYLIP package (Felsenstein 1993). Results Variation in mtdna Nucleotide Sequence The nucleotide sequence variation within the first 620 bp of the control region has been examined in 42 animals representing five native cattle breeds from Argentina and Brazil (Figure 1). Sequence alignments revealed the presence of nine haplotypes, eight of them not previously reported. Only haplotype EA1 matched perfectly with the reference sequence (Eucons), which was one of the two most frequently observed in this study (38.09%). The other highly represented haplotype (AA1; 33.33%) shared an identical sequence state with the African consensus (Afcons; Figure 1) at the three differentiating positions (16050, 16113, 16255) that split B. taurus mtdna sequences into two main clusters, the European taurine and the African taurine mitochondria, as described by Bradley et al. (1996). Among the seven remaining haplotypes, four are variations of the European consensus (EA2, EA3, EA4, EA5) and three are closer to the African consensus (AA2, AA3, AA4). The haplotype distribution varied according to the breeds as shown in Figure 1. For instance, Africanrelated haplotypes were absent in CU animals and were found in only one AC sample, the haplotype AA2. AA2 is four 324 The Journal of Heredity 2002:93(5)
3 mutational events away from the haplotype AA1 frequently observed in most of the Brazilian breeds studied in this work, but differs in two substitutions from the African consensus. The remaining AC animals only bore mitochondria from European origin, as did all the CU animals. We did not detect the presence of any B. indicus haplotype in the sample assayed; all the mtdna sequences clustered within the taurine branch in an neighbor-joining tree (Figure 2). We have also sequenced the corresponding mtdna fragment of one B. indicus individual from an imported Indian Nellore strain as PCR and sequencing control (Nell; Figure 1). Genetic Structure The genetic structure of the population was investigated through the AMOVA as implemented by the ARLEQUIN software (Schneider et al. 2000). It performs the analysis of variance (ANOVA) of the gene frequencies taking into account compositional differences among haplotypes. The hierarchical analysis decomposes the total variance into covariance components due to interindividual and interpopulational differences. Accordingly, a particular genetic structure can be tested defining alternative groups of populations. A global AMOVA was performed, considering an mtdna control region fragment of 240 bp of our 42 sequences added to those recovered from the GenBank, namely 209 European, 103 African, 48 Anatolian, 37 Middle East taurine breeds, 33 Japanese Black cattle (B. taurus), and 25 from Indian B. indicus breeds. More importantly, a further 49 mtdna sequences from six Portuguese cattle breeds (Cymbron et al. 1999) were incorporated into the analysis separately. Results indicate that 65.03% of the variance was accounted for by the eight continental subdivisions, 1.75% due to an among-breeds, within-continents component and 33.22% within breeds (Table 1, first column). In order to closely examine the partition of the variance components within the B. taurus mitochondrial lineages, a new round of analyses were carried out excluding the Indian breeds. Consequently the percentage of variation accounted for among continents dropped to 32.15%, increasing the amount of variance due to differences within breeds (64.68%) (Table 1, second column). A portion of the intrabreed variation observed after dropping the B. indicus sequences could result in part from Figure 1. mtdna sequence variations observed in 42 individuals from one Argentinean and four Brazilian native cattle breeds. All mitochondrial control region sequences obtained in this work, added to the one representing the African consensus (Afcons; Bradley et al. 1996), were aligned to the reference sequence of Anderson et al. (1982) (Eucons). Nucleotide positions were also numbered according to Eucons. Identity with the reference sequence is denoted by a dash, variations in the base letter involved in the substitution and deletions by a colon (:). The first two letters of the sequence identification, on the left-hand side, define the breed name (AC, Argentinean Creole; CC, Caracu; CU, Curraleiro; MN, Mocho Nacional; PP, Pantaneiro). Haplotype identification is given on the right side. Nell is the complete control region nucleotide sequence from a B. indicus representative animal of the Nellore breed. Miretti et al African-Derived Mitochondria 325
4 whereas 17 and 11 mtdna sequences pertained to the African-related group. AMOVA of this haplotype clustering significantly reduced the within-breed variation from 64.68% to 39.08% (Table 1) but raised the variance contribution of the among-groups component to 44.06%. Figure 2. Phylogenetic relationships (unrooted neighbor-joining tree) for B. taurus lineages inferred from interhaplotypic distances estimated from mtdna control region sequences ( ). The shaded areas circumscribe haplotypes belonging to the two major taurine lineages (African taurine and European taurine) and the Near Eastern cluster. Note the divergent branch conducting to the AA1 AA4 haplotypes. The AA3 haplotype is assumed to be AA1 when only the 240 bp fragment is considered (Figure 1). B. indicus is represented by the Nell sequence in Figure 1. sequence differences between Africanand European-related haplotypes within South American and Portuguese breeds. To investigate whether this subdivision actually influenced the variance distribution, and to what extent, haplotypes from Argentinean, Brazilian, and Portuguese breeds sharing the characteristic African consensus sequence state at the three nucleotide positions T16050, C16113, and C16255 (TCC) were pooled apart from those closely related to the European consensus sequence (CTT). To perform this new analysis, all sequences were aligned and classified accordingly. The clustering was inspected by the neighbor-joining tree (not shown) and by the haplotype-sharing option of the ARLEQUIN software. Mitotypes from 23 and 38 samples from the South American and Portuguese breeds, respectively, presented European-related haplotypes, Divergence Times Between Cattle Populations and Haplotypes The average corrected pairwise sequence difference between populations was estimated under the Tamura and Nei (1993) model of rate heterogeneity. The a value obtained from these data was 0.30 for the 240 bp fragment. It is assumed, based on paleontological evidence, that cattle and bison diverged at least 1 million years before present, which allows the calibration of mtdna sequence diversity. Differences of approximately 10% in the mutation rate between Bison B. taurus (Bb-Bt) and Bison B. indicus (Bb-Bi) lineages were observed. For instance, it reached 8.92% representing 89,230 years when applying the time depth estimation procedure described in Bradley et al. (1996). Hence the estimation of the age of the most recent common ancestor (MRCA) was carried out, averaging the rate of change between branches according to the relative rate method of Kumar and Hedges (1998) and Ingman et al. (2000). We then applied two other methods in order to check for consistency. First, the method used by Bradley et al. (1996) resulted in estimates similar to those previously obtained by Loftus et al. (1994) for the MRCA age for Bi-Bt (770,000 years, p ), and for Eu-Af divergence as well (61,500 years). Second, following the procedure of Vigilant et al. (1991), 88 substitutions (86 Ts, 2 Tv) were found within a 240 bp region in the alignment of 544 Bos mtdna control region sequences. As four Tv were observed between cattle and bison, the nucleotide differences accumulated within this sequence stretch, since their divergence 1 million years ago was 35.83%. This corresponds to a substitution rate of 1 Ts in 5,814 years, considering both mitochondrial lineages (1 Ts/ 11,628 years, one lineage rate) or to a B. taurus:b. indicus divergence time of 500,000 years. Finally, the age of the Bi-Bt MRCA based on the whole mitochondrial chromosome nucleotide sequence excluding the control region was 580,000 YBP (Miretti M et al., unpublished results). Once the calibration of the D-loop clock proved to be consistent regarding 326 The Journal of Heredity 2002:93(5)
5 the Bi-Bt separation, we proceeded with further time divergence estimations between B. taurus subdivisions. Thus, if Argentinean and Brazilian native cattle breeds (AB) are considered as one population, the estimated divergence time from the extant European, African, and Portuguese cattle is 32,900 YBP, 28,200 YBP, and 19,000 YBP, respectively, given a mutation rate l However, the reliability of these estimates can be questioned as AB, as well as P, is in fact an admixture of highly divergent haplotype derivatives from distant mitochondrial resources. Nevertheless, more confident information can be obtained from gene-tree analysis instead of a population-based analysis. Therefore we implemented haplotype clustering of B. taurus mtdna sequences in Near Eastern, African, and European taurus-derived haplotypes depending on the sequence state at the five positions as defined earlier (C16057 A16189 C16255, T16050 C16113 C16255, C16050 T16113 T16255, respectively). The AB haplotypes were thus separated into African derived (AA) and European derived (EA), a criterion that was also applied to Portuguese mtdna sequences which were thus subdivided into Portuguese-African (PAf)- and Portuguese- European (PEu)-related haplotypes. From the mean number of pairwise differences (p), comparable divergence times were found between EA-PEu and EA-Eu (2,200 and 2,400 years), and EA-Af and EA-PAf (58,700 and 67,300 years), respectively. Meanwhile, AA-EA presented higher divergence dates (157,000 years), even older than the Eu-Af separation and of similar magnitude to AA-Eu and AA-PEu (170,000 and 168,000 years). Anatolian and Middle Eastern cattle were indistinguishable under a population-based analysis (population average corrected p years), but their intrapopulation nucleotide diversity are among the widest (population average corrected p and , corresponding to 74,500 and 93,300 years, respectively). Anatolian and Middle Eastern high population diversity might result from either a considerably larger effective population size or longer genetic history, which is in line with the argument of Troy et al. (2001) in that sequences of the center of origin are expected to retain more ancestral variation and show higher haplotypic and nucleotide diversity. Even though this reasoning certainly can not be extended to explain why AB showed the Table 1. Analysis of molecular variance of cattle mtdna sequences Populations highest within-population diversity (population average corrected p corresponding to 98,700 years). Portuguese cattle also present substantial within-population diversity (population average corrected p corresponding to 68,800 years). Phylogenetic Tree Construction The neighbor-joining tree based on interhaplotypic distances (Figure 2) shows the phylogenetic relationship among B. taurus, B. indicus, and Bison mtdna control region sequences. Within the B. taurus lineage, three clades of sequences can be differentiated corresponding to the African, European, and Near Eastern taurine cattle. Haplotypes observed in Argentinean and Brazilian native cattle are placed within European (EA2, EA3, EA4, EA5) and African (AA1, AA2, AA3, AA4) sequence clusters. Note the extent of divergence in the lineage conducting to one of the most frequently observed haplotypes identified in this survey (AA1). Pairwise F ST values have been adopted to demonstrate the level of genetic distinction between populations (Bradley et al. 1996; Mannen et al. 1998; Roca et al. 2001) and can be used to estimate genetic distance over short time periods (Slatkin 1995). Figure 3 shows phylogenetic reconstruction based on population and haplotype group pairwise F ST values for the 240 bp fragment. All comparisons were significantly different (P, 10 5 ) except differentiation between the An and ME populations, and the PEu and Eu haplotype groups. Both dendrograms clearly show the distinction between the two major clades, B. indicus and B. taurus, as reported by Loftus et al. (1994). However, the branching pattern within B. taurus is markedly different between Haplotypes 1 B. indicus No B. indicus 1 B. indicus No B. indicus Among continents a Among populations within continent Within populations AMOVA testing of two different genetic structures: Populations considers a geographical distribution of haplotypes (American, European, Portuguese, African, Anatolian, Middle Eastern, Japanese, and Indian breeds). Haplotypes data were classified into haplotype groups according to their sequence state at nucleotide positions diagnostic of taurine lineages. The African-related haplotype group includes all those haplotypes sharing T16050 C16113 C16255, whereas the European-related haplotype group includes haplotypes sharing C16050 T16113 T To the Near Eastern-related haplotype group belong haplotypes sharing the following sequence state: 16057C, 16189A, and 16255C. 1 B. indicus and No B. indicus indicate the inclusion and exclusion of B. indicus sequences in the analysis. a Among African-related (AA, PAf, Af), European-related (EA, PEu, Eu), and Near Eastern-related haplotype groups (NE) and Japanese Black cattle (Jp1, Jp2). these two approaches, essentially contrasting distinctive features of population trees and gene trees. In the populationbased analysis, European (Eu) and Japanese Black (Jp) and Portuguese (P), Anatolian (An), and Middle Eastern (ME) cattle populations constitute two subgroups sharing clusters identified as European taurine cattle, corroborating the findings of Cymbron et al. (1999) and Mannen et al. (1998). Instead, African taurine breeds (Af) and Argentinean Brazilian native cattle (AB) are represented as distinct terminal branches of the neighbor-joining tree (Figure 3a). In the haplotype group analysis (gene tree), EA and PEu haplotypes were grouped within the European taurine haplotype cluster, constituted by most of the European sequences and a subdivision of the Japanese Black population (Jp2). On the other side, PAf haplotypes grouped together with the African taurine sequences (Af). Unexpectedly the AA subset of haplotypes observed in Brazilian native cattle constituted a distinct terminal branch excluded from both clusters (Figure 3b). Inspection of sequence alignment led to the observation of further sequencestate differences in four positions when comparing these AA haplotypes to Af and PAf sequences. Notably, the Africanrelated Argentinean Brazilian haplotypes, except AA2, shared the nucleotides CCTA at positions 16053, 16122, 16139, and (Figure 1), but this was not observed in mtdna sequences from Africa and Portugal. Conversely TTCG was found in Af and PAf, and in the European consensus as well. Consequently these haplotypes (AA1, AA3, AA4) differ from other African-derived mtdna sequences in at least these four positions, disregarding the other three Miretti et al African-Derived Mitochondria 327
6 Figure 3. Neighbor-joining tree based on pairwise F ST distances. (a) Population-based analysis. Note that the Argentinean Brazilian (AB) and African (Af) populations are separately grouped from the European (Eu), Japanese Black (Jp), Portuguese (P), Anatolian (An), and Middle Eastern (ME) cattle mtdna haplotypes, both within the taurine cluster. (b) Haplotype group analysis. Haplotype groups were classified into African-related, European-related, and Near Eastern-related according to their sequence state at positions 16050, 16057, 16113, 16185, and AA and PAf are African-derived haplotypes observed in South American and Portuguese cattle breeds, respectively; Jp1 and Jp2, Japanese Black; EA and PEu, European-derived haplotypes observed in South American and Portuguese breeds, respectively; Bi, B. indicus. substitutions that separate them from the European consensus. To focus on the relationship among mtdna haplotypes, a MSN was constructed based on substitutions observed in the 240 bp fragment of 520 B. taurus mtdna sequences (Figure 4). Two predominant and putatively separate ancestral mitochondrial types, defined as the European consensus (Eucons) and the African taurine consensus (Afcons) by Bradley et al. (1996), can be seen. Each sequence forms the center of radiation of a number of variants linked to it by mostly one or only a few substitutions. The most frequently observed African-related Argentinean Brazilian haplotype, AA1, is connected to the African consensus through four substitutions (large dark circle). Also note that the distance, in terms of substitution number, linking AA1 and Afcons is greater than the separation between Afcons and Eucons. Despite the fact that hypothetical relationships of the primary B. taurus haplotypes proposed by Troy et al. (2001) include unsampled intermediate nodes, our MSN displayed other equally parsimonious paths connecting the Eucons and Afcons. In these alternatives, which differ in the way that substitutions separating both clusters are ordered, haplotypes representing all the intermediate nodes were clearly identified, providing consistency to our MSN. Discussion We have assessed the nucleotide diversity in the control region of the mtdna in Argentinean Creole and four Brazilian native cattle breeds. This is the first report on mtdna concerning South American native cattle breeds, which allowed analysis of these data in conjunction with those from European, Anatolian, Middle Eastern, African, Japanese, and Indian cattle. All nine haplotypes described here are of B. taurus origin. The absence of B. indicus mitochondria was unexpected given that there is evidence, based on Y chromosome dimorphism (Britto 1998), of zebuine introgression into some Brazilian native breeds (Lara et al. 1997; Mazza et al. 1994). This result points toward a male-mediated introgression of zebu genes. Conversely Giovambattista et al. (2000) and Sinópoli (1993) observed only submetacentric Y chromosomes in different Argentinean Creole populations. We amplified and sequenced the mtdna control region of a B. indicus animal as a positive control, and primers were also designed to avoid any known polymorphic site. For instance, we found two substitution sites within the pro trna gene at positions (B. taurus C, B. indicus T) and (B. taurus A, B. indicus G), a region commonly used to place amplification primers. Regarding B. taurus mtdna haplotypes, we showed that information is lost when data analysis is only performed considering haplotypes as they occur in breeds (population-based phylogeny and AMO- VA; Figure 3 and Table 1). Under this assumption, estimates of nucleotide diversity and divergence times based only on the corrected MSPD between geographical subdivisions (populations) lead to a partial, and possibly wrong, interpretation of results due to the admixture of mtdna from both African and European taurine mitochondrial lineages within South American and Portuguese native cattle breeds. This reduces the divergence time between populations, inflates the variance component within breeds (Table 1), and distorts the branching pattern of dendrograms constructed with population pairwise F ST genetic distance (Figure 3). Indeed, the shorter separation time between the Argentinean Brazilian native cattle and Portuguese breeds (AB-P, 19,000 YBP) in comparison with its divergence from European (32,900 YBP) and African (28,200 YBP) breeds can be explained by the mitochondrial African component observed within the South American and Portuguese breeds. We then classified mtdna sequences into Near Eastern, African, and European taurine mitochondrial haplotype clusters (Figures 1 and 2), and subsequently subdivided them according to their geographic origin. Divergence estimates based on haplotype group comparisons were more reliable than population-based analysis. Regarding the European mitochondrial lineage of the taurine cluster, we estimate a divergence time of 2,400 YBP between European-related South American haplotypes and European haplotypes (EA-Eu). The distribution of the EA haplotypes within the breeds studied seems to represent the haplotype distribution described for European cattle, with Eucons as predominant and some closely related haplotypes. These haplotypes could have been introduced into America with either Portuguese or Spanish cattle, in spite of the fact that, disregarding EA1, they were absent in the Portuguese cattle breeds already studied (Cymbron et al. 1999) and that no mtdna information is available from Spanish cattle presumably introduced into the Americas. A different picture was found concerning the African-related Argentinean Brazilian haplotypes (AA). In the neighborjoining tree constructed with pairwise F ST distances (Figure 3b), the AA haplo- 328 The Journal of Heredity 2002:93(5)
7 type groups do not share any cluster with either the African (Af) or the African-related haplotypes of Portuguese cattle (PAf). Instead, it occupies an ancestral position regarding the separation of B. taurus mitochondrial lineages. In addition, divergence time estimations support the ancestrality of the AA haplotype group, which might have separated from all the other African-derived haplotypes 84,700 YBP, and 170,000 YBP from the European-derived haplotypes. Considering haplotypes separately, AA1, AA3, and AA4 were classified as Africanrelated due to their nucleotide composition TCC at positions 16050, 16113, and However, variations at four other positions 16053, 16122, 16139, and distinguished them from the two major taurine mitochondrial haplotype lineages: the European and the African taurine (Figures 1, 2, and 4). Notably Eucons and Afcons are more closely related to each other than Afcons to AA1. In this sense European and African mitochondrial lineages diverged 67,100 YBP, whereas African and AA showed older separation times (84,700 YBP). Pairwise haplotype group divergence between AA-PAf (80,100 years) is also consistent with ancestral separation of the mitochondrial lineage that gave rise to AA1, AA3, and AA4 haplotypes. Two other Pleistocene B. primigenius haplotypes are also distantly related to Eucons and Afcons (Bailey et al. 1996), but in this case substitutions involve other positions. Finally, age estimates should be cautiously interpreted, as they explicitly depend on the mutation rate considered ( ). We take time divergence in years as a relative measure of divergence between populations and haplotype groups, but the corresponding p- value for each age estimate can be obtained given the mutation rate, thus allowing comparison with other studies. The separation of African and European mitochondrial lineages should have occurred prior to the African and European cattle population expansion events 9,000 YBP and 5,000 YBP, respectively (Bailey et al. 1996). The profound population growth, which is apparent in the MSN topology, could be a consequence of the domestication process. Based on these observations, Bradley et al. (1996) suggested the existence of two independent domestication events involving two different strains of taurine progenitors. In this sense, if European and African mitochondrial types are indeed relicts Figure 4. Minimum spanning network (MSN) of taurine mtdna haplotypes. The MSN was constructed based on substitutions within the control region (240 bp) of mtdna haplotypes. Haplotypes having CTT at positions 16050, 16113, and are placed within the European consensus circle (Eucons), whereas those presenting TCC at these positions were classified as African consensus (Afcons) sequences. Haplotypes can be viewed as internal or terminal nodes, represented by circles interconnected by lines that correspond to a mutational step. The area of the circles is proportional to the haplotype frequency, except that the Eucons circle size should be twice as big, and short cross-hatched lines indicate further mutational events separating two haplotypes. Not all the alternative paths are shown and not all the haplotypes are represented. of temporally and spatially separate domestic origins (Bradley et al. 1996), it is tempting to interpret our genetic data as evidence of a further taurine subdivision that could have diverged from the taurine branch even before the European-African mitochondria split. This supposedly former taurine mitochondrial lineage, represented by AA1, AA3, and AA4 haplotypes found in South American cattle, was absent in African and Portuguese cattle surveys, even in the comprehensive cattle mtdna diversity study recently published (Troy et al. 2001), which precludes any inference about its geographical origin. Then, and for further analysis, it might be more appropriate to consider these haplotypes as part of another major mitochondrial lineage within B. taurus. It is also feasible that the high proportion of AA1 in South American native cattle, while absent in others, may not reflect European or African taurine history, but the evolutionary process during its prosperous American life. The AA1 haplotype was frequent in most of the Brazilian native breeds, which could be a sign of their introduction by Portuguese conquest through the Brazilian coastal line. But no consistency can be found in the Portuguese cattle study (Cymbron et al. 1999), since, although they reported haplotypes of African origin, the most closely related to AA1 differs in five substitutions. The African-related haplotype found in Argentinean Creole, AA2, belongs to a mitochondrial lineage different from those observed in Brazilian breeds and was also absent in Portuguese cattle. This could represent a hint for further investigations of native cattle breeds of other South American countries. Even though mtdna analysis should be extended to other African and European breeds, specifically Spanish cattle, before any conclusion can be drawn, our results might be mirroring a picture of the mitochondrial lineage distribution among cattle breeds in the Iberian peninsula five centuries ago. Since eight of the nine haplotypes found here have not been described in European and African cattle populations, it can be concluded that these haplotypes might have been lost, or alternatively, that further studies are required. References Anderson S, de Bruijn MHL, Coulson AR, Eperon IC, Sanger F, and Young IG, Complete sequence of bovine mitochondrial DNA. J Mol Biol 156: Bailey JF, Richards MB, Macaulay VA, Colson IB, James IT, Bradley DG, Hedges RE, and Sykes BC, Ancient DNA suggests a recent expansion of European cattle from a diverse wild progenitor species. Proc R Soc Lond B 263: Bradley DG, MacHugh DE, Cunningham P, and Loftus RT, Mitochondrial diversity and the origin of Miretti et al African-Derived Mitochondria 329
8 African and European cattle. Proc Natl Acad Sci USA 93: Britto CMC, Citogenética do Gado Pé Duro. Teresina: EDUFPI. Cymbron T, Loftus RT, Malheiro MI, and Bradley DG, Mitochondrial sequence variation suggests an African influence in Portuguese cattle. Proc R Soc Lond B 266: Epstein H and Mason IL, Cattle. In: Evolution of domesticated animals (Mason IL, ed). London: Longman; Felsenstein J, PHYLIP (phylogeny inference package). Seattle: University of Washington. Giberti H, Historia Económica de la Ganadería Argentina. Buenos Aires: Solar/Hachete. Giovambattista G, Ripoli MV, De Luca JC, Mirol PM, Liron JP, and Dulout FN, Male-mediated introgression of Bos indicus genes into Argentine and Bolivian Creole cattle breeds. Anim Genet 31: Hernández J, Instrucción del Estanciero, 1st ed. [reprinted 1995]. Buenos Aires: Roca. Ingman M, Kaessmann H, Pääbo S, and Gyllesten U, Mitochondrial genome variation and the origin of modern humans. Nature 408: Kumar S and Hedges B, A molecular timescale for vertebrate evolution. Nature 392: Lara MA, Sereno JR, Mazza MC, and Contel EPB, Investigação de variabilidade genética em bovinos Pantaneiros através de polimorfismos de proteínas. Annals of the 34th Annual Meeting of the Brazilian Society for the Zootechnics, Juiz de Fora, MG: Loftus RT, MacHugh DE, Bradley DG, Sharp PM, and Cunningham P, Evidence for two independent domestications of cattle. Proc Natl Acad Sci USA 91: MacHugh DG, Shriver MD, Loftus RT, Cunningham P, and Bradley DG, Microsatellite DNA variation and the evolution, domestication and phylogeography of taurine and zebu cattle (Bos taurus and Bos indicus). Genetics 146: Mannen M, Tsuji S, Loftus RT, and Bradley DG, Mitochondrial DNA variation and evolution of Japanese Black cattle (Bos taurus). Genetics 150: Mazza M, Mazza C, Sereno J, Santos S, and Pellegrin A, Ethnobiology and conservation of Pantaneiro cattle in Brazil. Brasilia: EMBRAPA-CPAP/SPI. Meyer S, Weiss G, and von Haeseler A, Pattern of nucleotide substitution and rate heterogeneity in the hypervariable regions I and II of human mtdna. Genetics 152: Payne WJA, 1991 Cattle genetic resources. In: World Animal Series B7 (Hickman CG, ed). Amsterdam: Elsevier Science; Rabasa SL, El bovino Criollo en los distintos países de América. In: Ganado Bovino Criollo, 3rd ed. Buenos Aires: Orientación Gráfica; Roca AL, Georgiadis N, Pecon-Slattery J, and O Brien SJ, Genetic evidence for two species of elephant in Africa. Science 293: Rouse JE, The Criollo: Spanish cattle in the Americas. Norman: University of Oklahoma Press. Salazar J and Cardozo A, Desarrollo del ganado criollo en América Latina: resumen histórico y distribución actual. In: Ganado Bovino Criollo, 1st ed. Buenos Aires: Orientación Gráfica; Sambrook J, Fritsch EF, and Maniatis T, Molecular cloning: a laboratory manual, 2nd ed. New York: Cold Spring Harbor Laboratory Press. Schneider S, Roessli D, and Excoffier L, ARLEQUIN version 2.0: a software for population genetics data analysis. Geneva, Switzerland: University of Geneva. Sinópoli LF, Estudios citogenéticos en la raza bovina Criolla Argentina, Bos taurus taurus. Posadas, Argentina: FCEQyN, UNaM. Slatkin, M, A measure of population subdivision based on microsatellite allele frequencies. Genetics 139: Strimmer K and von Haeseler A, Quartet puzzling: a quartet maximum likelihood method for reconstructing tree topologies. Mol Biol Evol 13: Tamura K and Nei M, Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol 10: Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, and Higgins DG, The CLUSTALX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 24: Troy CS, MacHugh DE, Bailey JF, Magee DA, Loftus RT, Cunningham P, Chamberlain AT, Sykes B, and Bradley DG, Genetic evidence for Near-Eastern origins of European cattle. Nature 410: Vigilant L, Stoneking M, Harpending H, Hawkes K, and Wilson A, African populations and the evolution of human mitochondrial DNA. Science 253: Zadworny D and Kuhnlein U, The identification of the kappa-casein genotype in Holstein dairy cattle using the polymerase chain reaction. Theor Appl Genet 80: Received April 17, 2001 Accepted August 8, 2002 Corresponding Editor: James Womack 330 The Journal of Heredity 2002:93(5)
Mitochondrial DNA Variation and Genetic Relationships in Japanese and Korean Cattle
1394 Asian-Aust. J. Anim. Sci. Vol. 19, No. 10 : 1394-1398 October 2006 www.ajas.info Mitochondrial DNA Variation and Genetic Relationships in Japanese and Korean Cattle S. Sasazaki*, S. Odahara, C. Hiura
More informationAsian-Aust. J. Anim. Sci. Vol. 22, No. 7 : July 2009
955 Asian-Aust. J. Anim. Sci. Vol. 22, No. 7 : 955-961 July 2009 www.ajas.info Determination of Phylogenetic Relationships of Turkish Native Cattle Breeds with Other Cattle Breeds Using Mitochondrial DNA
More informationCATTLE have had an important but incompletely but whether its domestication led to North East Asian
Copyright 1998 by the Genetics Society of America Mitochondrial DNA Variation and Evolution of Japanese Black Cattle (Bos taurus) H. Mannen,* S. Tsuji,* R. T. Loftus and D. G. Bradley *Faculty of Agriculture,
More informationIntraspecific gene genealogies: trees grafting into networks
Intraspecific gene genealogies: trees grafting into networks by David Posada & Keith A. Crandall Kessy Abarenkov Tartu, 2004 Article describes: Population genetics principles Intraspecific genetic variation
More informationLecture 11 Friday, October 21, 2011
Lecture 11 Friday, October 21, 2011 Phylogenetic tree (phylogeny) Darwin and classification: In the Origin, Darwin said that descent from a common ancestral species could explain why the Linnaean system
More informationQ1) Explain how background selection and genetic hitchhiking could explain the positive correlation between genetic diversity and recombination rate.
OEB 242 Exam Practice Problems Answer Key Q1) Explain how background selection and genetic hitchhiking could explain the positive correlation between genetic diversity and recombination rate. First, recall
More informationC3020 Molecular Evolution. Exercises #3: Phylogenetics
C3020 Molecular Evolution Exercises #3: Phylogenetics Consider the following sequences for five taxa 1-5 and the known outgroup O, which has the ancestral states (note that sequence 3 has changed from
More informationMitochondrial control region genetic diversity and maternal ancestry of a Brangus-Ibage cattle populations
Research Article Genetics and Molecular Biology, 28, 1, 60-66 (2005) Copyright by the Brazilian Society of Genetics. Printed in Brazil www.sbg.org.br Mitochondrial control region genetic diversity and
More informationA (short) introduction to phylogenetics
A (short) introduction to phylogenetics Thibaut Jombart, Marie-Pauline Beugin MRC Centre for Outbreak Analysis and Modelling Imperial College London Genetic data analysis with PR Statistics, Millport Field
More informationthat of Phylotree.org, mtdna tree Build 1756 (Supplementary TableS2). is resulted in 78 individuals allocated to the hg B4a1a1 and three individuals to hg Q. e control region (nps 57372 and nps 1602416526)
More informationBINF6201/8201. Molecular phylogenetic methods
BINF60/80 Molecular phylogenetic methods 0-7-06 Phylogenetics Ø According to the evolutionary theory, all life forms on this planet are related to one another by descent. Ø Traditionally, phylogenetics
More informationHomework Assignment, Evolutionary Systems Biology, Spring Homework Part I: Phylogenetics:
Homework Assignment, Evolutionary Systems Biology, Spring 2009. Homework Part I: Phylogenetics: Introduction. The objective of this assignment is to understand the basics of phylogenetic relationships
More informationMOLECULAR PHYLOGENY AND GENETIC DIVERSITY ANALYSIS. Masatoshi Nei"
MOLECULAR PHYLOGENY AND GENETIC DIVERSITY ANALYSIS Masatoshi Nei" Abstract: Phylogenetic trees: Recent advances in statistical methods for phylogenetic reconstruction and genetic diversity analysis were
More informationPhylogenetic relationship among S. castellii, S. cerevisiae and C. glabrata.
Supplementary Note S2 Phylogenetic relationship among S. castellii, S. cerevisiae and C. glabrata. Phylogenetic trees reconstructed by a variety of methods from either single-copy orthologous loci (Class
More informationApplications of Genetics to Conservation Biology
Applications of Genetics to Conservation Biology Molecular Taxonomy Populations, Gene Flow, Phylogeography Relatedness - Kinship, Paternity, Individual ID Conservation Biology Population biology Physiology
More informationEvolution Problem Drill 10: Human Evolution
Evolution Problem Drill 10: Human Evolution Question No. 1 of 10 Question 1. Which of the following statements is true regarding the human phylogenetic relationship with the African great apes? Question
More informationMolecular phylogeny How to infer phylogenetic trees using molecular sequences
Molecular phylogeny How to infer phylogenetic trees using molecular sequences ore Samuelsson Nov 2009 Applications of phylogenetic methods Reconstruction of evolutionary history / Resolving taxonomy issues
More informationDissecting the Matrilineal Components of Tongjiang Cattle from Southwest China
Biochem Genet (2008) 46:206 215 DOI 10.1007/s10528-008-9144-z Dissecting the Matrilineal Components of Tongjiang Cattle from Southwest China Shi-Yi Chen Æ Yi-Ping Liu Æ Wei Wang Æ Cheng-Zhong Gao Æ Yong-Gang
More informationDr. Amira A. AL-Hosary
Phylogenetic analysis Amira A. AL-Hosary PhD of infectious diseases Department of Animal Medicine (Infectious Diseases) Faculty of Veterinary Medicine Assiut University-Egypt Phylogenetic Basics: Biological
More informationMolecular phylogeny How to infer phylogenetic trees using molecular sequences
Molecular phylogeny How to infer phylogenetic trees using molecular sequences ore Samuelsson Nov 200 Applications of phylogenetic methods Reconstruction of evolutionary history / Resolving taxonomy issues
More informationRapid speciation following recent host shift in the plant pathogenic fungus Rhynchosporium
Rapid speciation following recent host shift in the plant pathogenic fungus Rhynchosporium Tiziana Vonlanthen, Laurin Müller 27.10.15 1 Second paper: Origin and Domestication of the Fungal Wheat Pathogen
More informationPhylogeny and Molecular Evolution. Introduction
Phylogeny and Molecular Evolution Introduction 1 2/62 3/62 Credit Serafim Batzoglou (UPGMA slides) http://www.stanford.edu/class/cs262/slides Notes by Nir Friedman, Dan Geiger, Shlomo Moran, Ron Shamir,
More informationEstimating Evolutionary Trees. Phylogenetic Methods
Estimating Evolutionary Trees v if the data are consistent with infinite sites then all methods should yield the same tree v it gets more complicated when there is homoplasy, i.e., parallel or convergent
More informationAmira A. AL-Hosary PhD of infectious diseases Department of Animal Medicine (Infectious Diseases) Faculty of Veterinary Medicine Assiut
Amira A. AL-Hosary PhD of infectious diseases Department of Animal Medicine (Infectious Diseases) Faculty of Veterinary Medicine Assiut University-Egypt Phylogenetic analysis Phylogenetic Basics: Biological
More informationLetter to the Editor. Department of Biology, Arizona State University
Letter to the Editor Traditional Phylogenetic Reconstruction Methods Reconstruct Shallow and Deep Evolutionary Relationships Equally Well Michael S. Rosenberg and Sudhir Kumar Department of Biology, Arizona
More informationMichael Yaffe Lecture #5 (((A,B)C)D) Database Searching & Molecular Phylogenetics A B C D B C D
7.91 Lecture #5 Database Searching & Molecular Phylogenetics Michael Yaffe B C D B C D (((,B)C)D) Outline Distance Matrix Methods Neighbor-Joining Method and Related Neighbor Methods Maximum Likelihood
More informationEffects of Gap Open and Gap Extension Penalties
Brigham Young University BYU ScholarsArchive All Faculty Publications 200-10-01 Effects of Gap Open and Gap Extension Penalties Hyrum Carroll hyrumcarroll@gmail.com Mark J. Clement clement@cs.byu.edu See
More informationMolecular phylogeny - Using molecular sequences to infer evolutionary relationships. Tore Samuelsson Feb 2016
Molecular phylogeny - Using molecular sequences to infer evolutionary relationships Tore Samuelsson Feb 2016 Molecular phylogeny is being used in the identification and characterization of new pathogens,
More informationPGA: A Program for Genome Annotation by Comparative Analysis of. Maximum Likelihood Phylogenies of Genes and Species
PGA: A Program for Genome Annotation by Comparative Analysis of Maximum Likelihood Phylogenies of Genes and Species Paulo Bandiera-Paiva 1 and Marcelo R.S. Briones 2 1 Departmento de Informática em Saúde
More informationPHYLOGENY AND SYSTEMATICS
AP BIOLOGY EVOLUTION/HEREDITY UNIT Unit 1 Part 11 Chapter 26 Activity #15 NAME DATE PERIOD PHYLOGENY AND SYSTEMATICS PHYLOGENY Evolutionary history of species or group of related species SYSTEMATICS Study
More informationGenetic diversity and phylogenetic relationship using mitochondrial ND5 partial sequences in three cattle breeds
Mal. J. Anim. Sci. 14:7-12 (211) Genetic diversity and phylogenetic relationship using mitochondrial ND5 partial sequences in three cattle breeds Aslinda, K. 1,4 *, Wan Somarny, W.M.Z. 1, Mat Tasol, S.
More informationInDel 3-5. InDel 8-9. InDel 3-5. InDel 8-9. InDel InDel 8-9
Lecture 5 Alignment I. Introduction. For sequence data, the process of generating an alignment establishes positional homologies; that is, alignment provides the identification of homologous phylogenetic
More informationA Phylogenetic Network Construction due to Constrained Recombination
A Phylogenetic Network Construction due to Constrained Recombination Mohd. Abdul Hai Zahid Research Scholar Research Supervisors: Dr. R.C. Joshi Dr. Ankush Mittal Department of Electronics and Computer
More informationAgricultural University
, April 2011 p : 8-16 ISSN : 0853-811 Vol16 No.1 PERFORMANCE COMPARISON BETWEEN KIMURA 2-PARAMETERS AND JUKES-CANTOR MODEL IN CONSTRUCTING PHYLOGENETIC TREE OF NEIGHBOUR JOINING Hendra Prasetya 1, Asep
More informationUnique variations of SRY gene result in distinct patrilineal phylogeny of Capra hircus and other domestic Bovidae*
Animal Science Papers and Reports vol. 30 (2013), no. 3, 219-227 Institute of Genetics and Animal Breeding, Jastrzębiec, Poland Unique variations of SRY gene result in distinct patrilineal phylogeny of
More informationPhylogenetic Tree Reconstruction
I519 Introduction to Bioinformatics, 2011 Phylogenetic Tree Reconstruction Yuzhen Ye (yye@indiana.edu) School of Informatics & Computing, IUB Evolution theory Speciation Evolution of new organisms is driven
More informationPOPULATION GENETICS Winter 2005 Lecture 17 Molecular phylogenetics
POPULATION GENETICS Winter 2005 Lecture 17 Molecular phylogenetics - in deriving a phylogeny our goal is simply to reconstruct the historical relationships between a group of taxa. - before we review the
More informationBig Idea #1: The process of evolution drives the diversity and unity of life
BIG IDEA! Big Idea #1: The process of evolution drives the diversity and unity of life Key Terms for this section: emigration phenotype adaptation evolution phylogenetic tree adaptive radiation fertility
More informationY Chromosome Haplogroup R1b-V88: Biogeographical Evidence for an Iberian Origin
Y Chromosome Haplogroup R1b-V88: Biogeographical Evidence for an Iberian Origin Michael R. Maglio Abstract The status on the origins of haplogroup R1b remains split between the Iberian Peninsula prior
More informationAlgorithms in Bioinformatics
Algorithms in Bioinformatics Sami Khuri Department of Computer Science San José State University San José, California, USA khuri@cs.sjsu.edu www.cs.sjsu.edu/faculty/khuri Distance Methods Character Methods
More informationOrganizing Life s Diversity
17 Organizing Life s Diversity section 2 Modern Classification Classification systems have changed over time as information has increased. What You ll Learn species concepts methods to reveal phylogeny
More informationMETHODS FOR DETERMINING PHYLOGENY. In Chapter 11, we discovered that classifying organisms into groups was, and still is, a difficult task.
Chapter 12 (Strikberger) Molecular Phylogenies and Evolution METHODS FOR DETERMINING PHYLOGENY In Chapter 11, we discovered that classifying organisms into groups was, and still is, a difficult task. Modern
More information"Nothing in biology makes sense except in the light of evolution Theodosius Dobzhansky
MOLECULAR PHYLOGENY "Nothing in biology makes sense except in the light of evolution Theodosius Dobzhansky EVOLUTION - theory that groups of organisms change over time so that descendeants differ structurally
More informationPhylogenetic Analysis. Han Liang, Ph.D. Assistant Professor of Bioinformatics and Computational Biology UT MD Anderson Cancer Center
Phylogenetic Analysis Han Liang, Ph.D. Assistant Professor of Bioinformatics and Computational Biology UT MD Anderson Cancer Center Outline Basic Concepts Tree Construction Methods Distance-based methods
More informationReading for Lecture 13 Release v10
Reading for Lecture 13 Release v10 Christopher Lee November 15, 2011 Contents 1 Evolutionary Trees i 1.1 Evolution as a Markov Process...................................... ii 1.2 Rooted vs. Unrooted Trees........................................
More informationPhylogeny and Molecular Evolution. Introduction
Phylogeny and Molecular Evolution Introduction 1 Credit Serafim Batzoglou (UPGMA slides) http://www.stanford.edu/class/cs262/slides Notes by Nir Friedman, Dan Geiger, Shlomo Moran, Ron Shamir, Sagi Snir,
More informationOMICS Journals are welcoming Submissions
OMICS Journals are welcoming Submissions OMICS International welcomes submissions that are original and technically so as to serve both the developing world and developed countries in the best possible
More informationSupporting Information
Supporting Information Das et al. 10.1073/pnas.1302500110 < SP >< LRRNT > < LRR1 > < LRRV1 > < LRRV2 Pm-VLRC M G F V V A L L V L G A W C G S C S A Q - R Q R A C V E A G K S D V C I C S S A T D S S P E
More informationA comparison of two popular statistical methods for estimating the time to most recent common ancestor (TMRCA) from a sample of DNA sequences
Indian Academy of Sciences A comparison of two popular statistical methods for estimating the time to most recent common ancestor (TMRCA) from a sample of DNA sequences ANALABHA BASU and PARTHA P. MAJUMDER*
More informationThe Phylogenetic Handbook
The Phylogenetic Handbook A Practical Approach to DNA and Protein Phylogeny Edited by Marco Salemi University of California, Irvine and Katholieke Universiteit Leuven, Belgium and Anne-Mieke Vandamme Rega
More informationGENETICS - CLUTCH CH.22 EVOLUTIONARY GENETICS.
!! www.clutchprep.com CONCEPT: OVERVIEW OF EVOLUTION Evolution is a process through which variation in individuals makes it more likely for them to survive and reproduce There are principles to the theory
More informationPiecing It Together. 1) The envelope contains puzzle pieces for 5 vertebrate embryos in 3 different stages of
Piecing It Together 1) The envelope contains puzzle pieces for 5 vertebrate embryos in 3 different stages of development. Lay out the pieces so that you have matched up each animal name card with its 3
More informationTHEORY. Based on sequence Length According to the length of sequence being compared it is of following two types
Exp 11- THEORY Sequence Alignment is a process of aligning two sequences to achieve maximum levels of identity between them. This help to derive functional, structural and evolutionary relationships between
More information31/10/2012. Human Evolution. Cytochrome c DNA tree
Human Evolution Cytochrome c DNA tree 1 Human Evolution! Primate phylogeny! Primates branched off other mammalian lineages ~65 mya (mya = million years ago) Two types of monkeys within lineage 1. New World
More informationSUMMARY KEYWORDS. Mitochondrial DNA, sequencing, genetic characterization, bluefin tuna, Mediterranean Sea, Atlantic Ocean
MITOCHONDRIAL GENETIC CHARACTERIZATION OF BLUEFIN TUNA (Thunnus thynnus) FROM THREE MEDITERRANEAN (LIBYA, MALTA, TUNISIA); AND ONE ATLANTIC LOCATIONS (GULF OF CADIZ). J. Viñas ζ 1; C. Pla 1; M.Y. Tawil
More informationCladistics and Bioinformatics Questions 2013
AP Biology Name Cladistics and Bioinformatics Questions 2013 1. The following table shows the percentage similarity in sequences of nucleotides from a homologous gene derived from five different species
More informationVariances of the Average Numbers of Nucleotide Substitutions Within and Between Populations
Variances of the Average Numbers of Nucleotide Substitutions Within and Between Populations Masatoshi Nei and Li Jin Center for Demographic and Population Genetics, Graduate School of Biomedical Sciences,
More informationWhat is Phylogenetics
What is Phylogenetics Phylogenetics is the area of research concerned with finding the genetic connections and relationships between species. The basic idea is to compare specific characters (features)
More informationUoN, CAS, DBSC BIOL102 lecture notes by: Dr. Mustafa A. Mansi. The Phylogenetic Systematics (Phylogeny and Systematics)
- Phylogeny? - Systematics? The Phylogenetic Systematics (Phylogeny and Systematics) - Phylogenetic systematics? Connection between phylogeny and classification. - Phylogenetic systematics informs the
More informationFields connected to Phylogeography Microevolutionary disciplines Ethology Demography Population genetics
Stephen A. Roussos Fields connected to Phylogeography Microevolutionary disciplines Ethology Demography Population genetics Macrevolutionary disciplines Historical geography Paleontology Phylogenetic biology
More informationPhylogenetic Analysis
Phylogenetic Analysis Aristotle Through classification, one might discover the essence and purpose of species. Nelson & Platnick (1981) Systematics and Biogeography Carl Linnaeus Swedish botanist (1700s)
More informationBioinformatics tools for phylogeny and visualization. Yanbin Yin
Bioinformatics tools for phylogeny and visualization Yanbin Yin 1 Homework assignment 5 1. Take the MAFFT alignment http://cys.bios.niu.edu/yyin/teach/pbb/purdue.cellwall.list.lignin.f a.aln as input and
More informationChapter 26: Phylogeny and the Tree of Life Phylogenies Show Evolutionary Relationships
Chapter 26: Phylogeny and the Tree of Life You Must Know The taxonomic categories and how they indicate relatedness. How systematics is used to develop phylogenetic trees. How to construct a phylogenetic
More informationConstructing Evolutionary/Phylogenetic Trees
Constructing Evolutionary/Phylogenetic Trees 2 broad categories: istance-based methods Ultrametric Additive: UPGMA Transformed istance Neighbor-Joining Character-based Maximum Parsimony Maximum Likelihood
More informationLecture 13: Population Structure. October 8, 2012
Lecture 13: Population Structure October 8, 2012 Last Time Effective population size calculations Historical importance of drift: shifting balance or noise? Population structure Today Course feedback The
More informationGenetic Drift in Human Evolution
Genetic Drift in Human Evolution (Part 2 of 2) 1 Ecology and Evolutionary Biology Center for Computational Molecular Biology Brown University Outline Introduction to genetic drift Modeling genetic drift
More information8/23/2014. Phylogeny and the Tree of Life
Phylogeny and the Tree of Life Chapter 26 Objectives Explain the following characteristics of the Linnaean system of classification: a. binomial nomenclature b. hierarchical classification List the major
More informationPopulation Structure
Ch 4: Population Subdivision Population Structure v most natural populations exist across a landscape (or seascape) that is more or less divided into areas of suitable habitat v to the extent that populations
More informationHow should we organize the diversity of animal life?
How should we organize the diversity of animal life? The difference between Taxonomy Linneaus, and Cladistics Darwin What are phylogenies? How do we read them? How do we estimate them? Classification (Taxonomy)
More informationEVOLUTION change in populations over time
EVOLUTION change in populations over time HISTORY ideas that shaped the current theory James Hutton (1785) proposes that Earth is shaped by geological forces that took place over extremely long periods
More informationPhylogenetic Trees. Phylogenetic Trees Five. Phylogeny: Inference Tool. Phylogeny Terminology. Picture of Last Quagga. Importance of Phylogeny 5.
Five Sami Khuri Department of Computer Science San José State University San José, California, USA sami.khuri@sjsu.edu v Distance Methods v Character Methods v Molecular Clock v UPGMA v Maximum Parsimony
More informationProcesses of Evolution
15 Processes of Evolution Forces of Evolution Concept 15.4 Selection Can Be Stabilizing, Directional, or Disruptive Natural selection can act on quantitative traits in three ways: Stabilizing selection
More informationBustamante et al., Supplementary Nature Manuscript # 1 out of 9 Information #
Bustamante et al., Supplementary Nature Manuscript # 1 out of 9 Details of PRF Methodology In the Poisson Random Field PRF) model, it is assumed that non-synonymous mutations at a given gene are either
More informationNon-independence in Statistical Tests for Discrete Cross-species Data
J. theor. Biol. (1997) 188, 507514 Non-independence in Statistical Tests for Discrete Cross-species Data ALAN GRAFEN* AND MARK RIDLEY * St. John s College, Oxford OX1 3JP, and the Department of Zoology,
More informationExamples of Phylogenetic Reconstruction
Examples of Phylogenetic Reconstruction 1. HIV transmission Recently, an HIV-positive Florida dentist was suspected of having transmitted the HIV virus to his dental patients. Although a number of his
More informationI. Short Answer Questions DO ALL QUESTIONS
EVOLUTION 313 FINAL EXAM Part 1 Saturday, 7 May 2005 page 1 I. Short Answer Questions DO ALL QUESTIONS SAQ #1. Please state and BRIEFLY explain the major objectives of this course in evolution. Recall
More informationPrimate Diversity & Human Evolution (Outline)
Primate Diversity & Human Evolution (Outline) 1. Source of evidence for evolutionary relatedness of organisms 2. Primates features and function 3. Classification of primates and representative species
More informationEVOLUTION change in populations over time
EVOLUTION change in populations over time HISTORY ideas that shaped the current theory James Hutton & Charles Lyell proposes that Earth is shaped by geological forces that took place over extremely long
More informationLeber s Hereditary Optic Neuropathy
Leber s Hereditary Optic Neuropathy Dear Editor: It is well known that the majority of Leber s hereditary optic neuropathy (LHON) cases was caused by 3 mtdna primary mutations (m.3460g A, m.11778g A, and
More informationName: Hour: Teacher: ROZEMA. Inheritance & Mutations Connected to Speciation
Name: Hour: Teacher: ROZEMA Inheritance & Mutations Connected to Speciation Let s Review What We Already Know: What Have We Learned? Lesson 26: PI 1 (Projected Image) - Human Karyotype (image from https://en.wikipedia.org/wiki/karyotype#/media/file:nhgri_human_male_karyotype.png)
More informationMtDNA profiles and associated haplogroups
A systematic approach to an old problem Anita Brandstätter 1, Alexander Röck 2, Arne Dür 2, Walther Parson 1 1 Institute of Legal Medicine Innsbruck Medical University 2 Institute of Mathematics University
More informationBayesian Inference using Markov Chain Monte Carlo in Phylogenetic Studies
Bayesian Inference using Markov Chain Monte Carlo in Phylogenetic Studies 1 What is phylogeny? Essay written for the course in Markov Chains 2004 Torbjörn Karfunkel Phylogeny is the evolutionary development
More informationThe Origin of Species
The Origin of Species Introduction A species can be defined as a group of organisms whose members can breed and produce fertile offspring, but who do not produce fertile offspring with members of other
More informationSequence analysis and comparison
The aim with sequence identification: Sequence analysis and comparison Marjolein Thunnissen Lund September 2012 Is there any known protein sequence that is homologous to mine? Are there any other species
More informationGENETICS OF MODERN HUMAN ORIGINS AND DIVERSITY
Annu. Rev. Anthropol. 1998. 27:1 23 Copyright 1998 by Annual Reviews. All rights reserved GENETICS OF MODERN HUMAN ORIGINS AND DIVERSITY John H. Relethford Department of Anthropology, State University
More informationEVOLUTIONARY DISTANCES
EVOLUTIONARY DISTANCES FROM STRINGS TO TREES Luca Bortolussi 1 1 Dipartimento di Matematica ed Informatica Università degli studi di Trieste luca@dmi.units.it Trieste, 14 th November 2007 OUTLINE 1 STRINGS:
More informationO 3 O 4 O 5. q 3. q 4. Transition
Hidden Markov Models Hidden Markov models (HMM) were developed in the early part of the 1970 s and at that time mostly applied in the area of computerized speech recognition. They are first described in
More informatione.g. population: 500, two alleles: Red (R) and White (r). Total: 1000 genes for flower color in the population
The Evolution of Populations What is Evolution? A change over time in the genetic composition of a population Human evolution The gene pool Is the total aggregate of genes for a particular trait in a population
More informationRobust demographic inference from genomic and SNP data
Robust demographic inference from genomic and SNP data Laurent Excoffier Isabelle Duperret, Emilia Huerta-Sanchez, Matthieu Foll, Vitor Sousa, Isabel Alves Computational and Molecular Population Genetics
More information9/30/11. Evolution theory. Phylogenetic Tree Reconstruction. Phylogenetic trees (binary trees) Phylogeny (phylogenetic tree)
I9 Introduction to Bioinformatics, 0 Phylogenetic ree Reconstruction Yuzhen Ye (yye@indiana.edu) School of Informatics & omputing, IUB Evolution theory Speciation Evolution of new organisms is driven by
More informationAdditive distances. w(e), where P ij is the path in T from i to j. Then the matrix [D ij ] is said to be additive.
Additive distances Let T be a tree on leaf set S and let w : E R + be an edge-weighting of T, and assume T has no nodes of degree two. Let D ij = e P ij w(e), where P ij is the path in T from i to j. Then
More informationPhylogenetic inference
Phylogenetic inference Bas E. Dutilh Systems Biology: Bioinformatic Data Analysis Utrecht University, March 7 th 016 After this lecture, you can discuss (dis-) advantages of different information types
More informationHow to read and make phylogenetic trees Zuzana Starostová
How to read and make phylogenetic trees Zuzana Starostová How to make phylogenetic trees? Workflow: obtain DNA sequence quality check sequence alignment calculating genetic distances phylogeny estimation
More informationSupplementary Information for Hurst et al.: Causes of trends of amino acid gain and loss
Supplementary Information for Hurst et al.: Causes of trends of amino acid gain and loss Methods Identification of orthologues, alignment and evolutionary distances A preliminary set of orthologues was
More informationPopulation Genetics of Modern Human Evolution
Population Genetics of Modern Human Evolution John H Relethford, State University of New York College at Oneonta, Oneonta, New York, USA Rosalind M Harding, University of Oxford, Oxford, UK The aim of
More informationCREATING PHYLOGENETIC TREES FROM DNA SEQUENCES
INTRODUCTION CREATING PHYLOGENETIC TREES FROM DNA SEQUENCES This worksheet complements the Click and Learn developed in conjunction with the 2011 Holiday Lectures on Science, Bones, Stones, and Genes:
More informationConcepts and Methods in Molecular Divergence Time Estimation
Concepts and Methods in Molecular Divergence Time Estimation 26 November 2012 Prashant P. Sharma American Museum of Natural History Overview 1. Why do we date trees? 2. The molecular clock 3. Local clocks
More informationLikelihood Ratio Tests for Detecting Positive Selection and Application to Primate Lysozyme Evolution
Likelihood Ratio Tests for Detecting Positive Selection and Application to Primate Lysozyme Evolution Ziheng Yang Department of Biology, University College, London An excess of nonsynonymous substitutions
More informationGenetic Admixture in the Late Pleistocene
AMERICAN JOURNAL OF PHYSICAL ANTHROPOLOGY 1OO:l-5 (1996) Genetic Admixture in the Late Pleistocene ELISABETH J. MANDERSCHEID AND ALAN R. ROGERS Department of Anthropology, University of Utah, Salt Lake
More informationA PARSIMONY APPROACH TO ANALYSIS OF HUMAN SEGMENTAL DUPLICATIONS
A PARSIMONY APPROACH TO ANALYSIS OF HUMAN SEGMENTAL DUPLICATIONS CRYSTAL L. KAHN and BENJAMIN J. RAPHAEL Box 1910, Brown University Department of Computer Science & Center for Computational Molecular Biology
More information