Genetic evidence for Near-Eastern origins of European cattle.

The limited ranges of the wild progenitors of many of the primary European domestic species point to their origins further east in Anatolia or the fertile crescent. The wild ox (Bos primigenius), however, ranged widely and it is unknown whether it was domesticated within Europe as one feature of a local contribution to the farming economy. Here we examine mitochondrial DNA control-region sequence variation from 392 extant animals sampled from Europe, Africa and the Near East, and compare this with data from four extinct British wild oxen. The ancient sequences cluster tightly in a phylogenetic analysis and are clearly distinct from modern cattle. Network analysis of modern Bos taurus identifies four star-like clusters of haplotypes, with intra-cluster diversities that approximate to that expected from the time depth of domestic history. Notably, one of these clusters predominates in Europe and is one of three encountered at substantial frequency in the Near East. In contrast, African diversity is almost exclusively composed of a separate haplogroup, which is encountered only rarely elsewhere. These data provide strong support for a derived Near-Eastern origin for European cattle.

European Y-chromosomal lineages in Polynesians: a contrast to the population structure revealed by mtDNA.

We have used Y-chromosomal polymorphisms to trace paternal lineages in Polynesians by use of samples previously typed for mtDNA variants. A genealogical approach utilizing hierarchical analysis of eight rare-event biallelic polymorphisms, seven microsatellite loci, and internal structural analysis of the hypervariable minisatellite, MSY1, has been used to define three major paternal-lineage clusters in Polynesians. Two of these clusters, both defined by novel MSY1 modular structures and representing 55% of the Polynesians studied, are also found in coastal Papua New Guinea. Reduced Polynesian diversity, relative to that in Melanesians, is illustrated by the presence of several examples of identical MSY1 codes and microsatellite haplotypes within these lineage clusters in Polynesians. The complete lack of Y chromosomes having the M4 base substitution in Polynesians, despite their prevalence (64%) in Melanesians, may also be a result of the multiple bottleneck events during the colonization of this region of the world. The origin of the M4 mutation has been dated by use of two independent methods based on microsatellite-haplotype and minisatellite-code diversity. Because of the wide confidence limits on the mutation rates of these loci, the M4 mutation cannot be conclusively dated relative to the colonization of Polynesia, 3,000 years ago. The other major lineage cluster found in Polynesians, defined by a base substitution at the 92R7 locus, represents 27% of the Polynesians studied and, most probably, originates in Europe. This is the first Y-chromosomal evidence of major European admixture with indigenous Polynesian populations and contrasts sharply with the picture given by mtDNA evidence.

Phylogeography of mitochondrial DNA in western Europe.

For most of the past century, prehistorians have had to rely on the fossil and archaeological records in order to reconstruct the past. In the last few decades, this evidence has been substantially supplemented from classical human genetics. More recently, phylogenetic analyses of DNA sequences that incorporate geographical information have provided a high-resolution tool for the investigation of prehistoric demographic events, such as founder effects and population expansions. These events can be dated using a molecular clock when the mutation rate and founder haplotypes are known. We have previously applied such methods to sequence data from the mitochondrial DNA control region, to suggest that most extant mitochondrial sequences in western Europe have a local ancestry in the Early Upper Palaeolithic, with a smaller proportion arriving from the Near East in the Neolithic. Here, we describe a cladistic notation for mitochondrial variation and expand upon our earlier analysis to present a more detailed portrait of the European mitochondrial record.

Variable levels of a heteroplasmic point mutation in individual hair roots.

During direct sequencing of the first hypervariable segment of the human mitochondrial control region, we identified one individual with a heteroplasmic point mutation at nt 16,256. We used primer extension to analyze the proportions of each mitochondrial haplotype in peripheral blood, buccal cells, and single hair roots from this individual and from eight members of his maternal lineage. Significant levels of heteroplasmy were found in only three individuals, and, in these cases, the proportions of each haplotype were similar in both blood and buccal cells. From the changes in mitochondrial haplotypes within mother-offspring pairs, we calculated that the most likely size of a mitochondrial bottleneck during development was 1-27 segregating units. However, highly variable levels of heteroplasmy were found in single hair roots, even among roots from the same individual. We analyzed a large number of hair roots from one individual and found that the proportion of one haplotype was within a range of 9% to > 99% in different roots. Roots originating from within a small patch of skin had haplotype proportions as variable as those from different areas of skin.

Heteroplasmic point mutations in the human mtDNA control region.

As part of an investigation of the fixation mechanisms of mtDNA mutations in humans, we sequenced the first hypervariable segment of the control region in 180 twin pairs and found evidence of site heteroplasmy in 4 pairs. Significant levels of two mitochondrial haplotypes differing by a single point mutation were found in two MZ pairs, and within each pair, both members had similar levels of heteroplasmy. Two DZ pairs were found in which the predominant mitochondrial haplotype differed within the pair. We measured proportions of mitochondrial haplotypes within two twin pairs and their maternal relatives, using primer extension. In both maternal lineages, most family members were heteroplasmic, and the proportions of each genotype varied widely in different individuals. We used the changes in haplotype proportions within mother-offspring pairs to calculate the size range of potential bottlenecks in mitochondrial numbers occurring during development of the offspring. In most individuals, the most likely effective bottleneck sizes ranged from 3 to 20 segregating units, though in two individuals a small bottleneck was very unlikely and there was no upper limit on its possible size. We also used the data from this study, together with unpublished data from other populations, to estimate the frequency of site heteroplasmy in normal human populations. From this, we calculated that the rate of mutation and fixation in the first hypervariable segment of the human mtDNA control region is between 1.2 x 10(-6) and 2.7 x 10(-5) per site per generation. This range is in good agreement with published estimates calculated by other methods.

Ancient DNA suggests a recent expansion of European cattle from a diverse wild progenitor species.

A total of 11 Bos primigenius and Bos taurus bones from archaeological sites between 500 and 12000 years old were examined for the presence of DNA. It was possible to amplify and sequence mitochondrial control region DNA extracted from seven of the 11 samples, including two Pleistocene B. primigenius samples. We compared the results with published data by constructing phylogenetic networks. The two B. primigenius samples clustered with the extant B. taurus samples in the networks. The similarity between B. primigenius and modern taurine cattle confirms that these should be considered members of a single species. The sequences obtained from the B. taurus specimens were either identical to the reference sequence for modern European cattle or closely related to it. They included two sequences not previously documented. The network analysis of the ancient data highlights the intermediary nature of the B. primigenius sequences between modern European and African B. taurus and the proximity of the ancient DNA B. taurus sequences to modern European B. taurus. Further analysis of the extant data in the light of the ancient DNA results suggests that a degree of Pleistocene diversity survives in the extant European Bos population that is mainly derived from a more recent population expansion.

Site 73 in hypervariable region II of the human mitochondrial genome and the origin of European populations.

The majority of published human mitochondrial DNA sequence data are confined to hypervariable region I in the control region. By contrast, this paper focuses on a nucleotide site in hypervariable region II. Unlike most non-European populations whose mtDNA sequences have been studied in the literature, the British 'white Caucasian' population has a high level of variation at site 73 (following the site numbering by Anderson et al. 1981). This variation appears to have its origin largely in a mutation from guanine to adenine at that site with an estimated minimum age between 15,000 and 25,000 years. The data of Piercy et al. (1993) suggest that roughly half of the British 'white Caucasian' mitochondrial gene pool is descended from a common maternal ancestor who carried this mutation at site 73. This site also plays a central role in distinguishing the five major European mtDNA clusters identified in Richards et al. (1996). We suggest that the lineages carrying an A at site 73, together with some other lineages, may have their origins in a small founder population which expanded after the last glacial maximum about 20,000 years ago. We conclude that, in addition to region I sequences, site 73 is worth determining in studies of Caucasian populations.

Genetic diversity in the Iberian Peninsula determined from mitochondrial sequence analysis.

We have analysed 302 bp of the first hypervariable region of the mitochondrial D-loop in 271 individuals from different regions of the Iberian Peninsula and 85 individuals from Algeria. The Basque population is significantly different from neighbouring populations in terms of overall levels of diversity. This is because the majority of sequences in the Basques are restricted to the lineage group defined by the CRS (Cambridge Reference Sequence) and its derivatives although, like other Iberian populations, they showed a unimodal distribution of pairwise sequence differences. The timing of divergence of populations within Iberia points to a shared ancestry of all populations in the Upper Palaeolithic. Further genetic subdivision is apparent in Catalonia and Andalusia, with increased genetic diversity in the latter. Lineage diversity comparisons of Iberian populations with European (Tuscan) and North African (Algerian) populations shows the Iberian Peninsula to be more similar to other European populations, although a small number of Iberian lineages can be traced to North Africa.

Mitochondrial portraits of human populations using median networks.

Analysis of variation in the hypervariable region of mitochondrial DNA (mtDNA) has emerged as an important tool for studying human evolution and migration. However, attempts to reconstruct optimal intraspecific mtDNA phylogenies frequently fail because parallel mutation events partly obscure the true evolutionary pathways. This makes it inadvisable to present a single phylogenetic tree at the expense of neglecting equally acceptable ones. As an alternative, we propose a novel network approach for portraying mtDNA relationships. For small sample sizes (< approximately 50), an unmodified median network contains all most parsimonious trees, displays graphically the full information content of the sequence data, and can easily be generated by hand. For larger sample sizes, we reduce the complexity of the network by identifying parallelisms. This reduction procedure is guided by a compatibility argument and an additional source of phylogenetic information: the frequencies of the mitochondrial haplotypes. As a spin-off, our approach can also assist in identifying sequencing errors, which manifest themselves in implausible network substructures. We illustrate the advantages of our approach with several examples from existing data sets.

Length heteroplasmy in the first hypervariable segment of the human mtDNA control region.

The first hypervariable segment of the human mtDNA control region contains a homopolymeric tract of cytosines between nt 16184 and 16193, interrupted at position 16189 by a thymine, according to the Cambridge reference sequence. A variant commonly found in population screening is a T-to-C transition at nt 16189, resulting in an uninterrupted homopolymeric tract. Direct sequencing of individuals with this variant produces a characteristic blurred sequence in nucleotides beyond the tract. Sequencing clones from these individuals revealed that this is caused by high levels of length heteroplasmy in the homopolymeric tract and low levels of length heteroplasmy in the four adenines following the tract. We have developed a rapid method involving densitometry of sequencing gels to quantify the relative proportions of different length variants present in an individual. We have used this to study the proportions of length variants in individuals from three twin pairs and two maternal lineages. While unrelated individuals usually have different proportions of length variants, all maternally related individuals studied have the same proportions, even if they are only distantly related. It is not obvious how identical heteroplasmic profiles are maintained in maternally related individuals, but some possible mechanisms are suggested.

The calcium binding properties and molecular organization of epidermal growth factor-like domains in human fibrillin-1.

Human fibrillin-1 is a 350-kDa glycoprotein found in 10-nm connective tissue microfibrils. Mutations in the gene encoding this protein cause the Marfan syndrome, a disease characterized by cardiovascular, ocular, and skeletal abnormalities. Fibrillin-1 has a modular structure that includes 47 epidermal growth factor-like (EGF-like) domains, 43 of which contain a consensus sequence associated with calcium binding. A mutation causing an Asn-2144 --> Ser amino acid change in one of the potential calcium binding residues has been described in a patient with the Marfan syndrome. We have chemically synthesized a wild-type EGF-like domain (residues 2126-2165 of human fibrillin-1) and a mutant EGF-like domain containing the Asn-2144 --> Ser amino acid change and measured calcium binding to each using 1H-NMR spectroscopy. The wild-type domain binds calcium with a similar affinity to isolated EGF-like domains from coagulation factors IX and X; however, the mutant domain exhibits > 5-fold reduction in affinity. Rotary shadowing of fibrillin-containing microfibrils, isolated from dermal fibroblast cultures obtained from the Marfan patient, shows that the mutation does not prevent assembly of fibrillin into microfibrils but does alter the appearance of the interbead region. We have modeled a region of fibrillin-1 (residues 2126-2331) encompassing five calcium binding EGF-like domains, using data derived from the recently determined crystal structure of a calcium binding EGF-like domain from human factor IX. Our model suggests that these fibrillin-1 EGF-like domains adopt a helical arrangement stabilized by calcium and that defective calcium binding to a single EGF-like domain results in distortion of the helix. We propose a mechanism for the interaction of contiguous arrays of calcium binding EGF-like domains within the microfibril.

Genetic mapping of a locus for multiple epiphyseal dysplasia (EDM2) to a region of chromosome 1 containing a type IX collagen gene.

Multiple epiphyseal dysplasia (MED) is a dominantly inherited chondrodysplasia characterized by mild short stature and early-onset osteoarthrosis. Some forms of MED clinically resemble another chondrodysplasia phenotype, the mild form of pseudoachondroplasia (PSACH). On the basis of their clinical similarities as well as similar ultrastructural and biochemical features in cartilage from some patients, it has been proposed that MED and PSACH belong to a single bone-dysplasia family. Recently, both mild and severe PSACH as well as a form of MED have been linked to the same interval on chromosome 19, suggesting that they may be allelic disorders. Linkage studies with the chromosome 19 markers were carried out in a large family with MED and excluded the previously identified interval. Using this family, we have identified an MED locus on the short arm of chromosome 1, in a region containing the gene (COL9A2) that encodes the alpha 2 chain of type IX collagen, a structural component of the cartilage extracellular matrix.

A new missense mutation of fibrillin in a patient with Marfan syndrome.

A patient with Marfan syndrome was shown to be heterozygous for a G to A transition at nucleotide 3952 of the FBNI gene. This would result in a cysteine to tyrosine substitution at amino acid 1223 in the fibrillin protein.

Mutation screening by a combination of biotin-SSCP and direct sequencing.

We have developed a mutation detection strategy that combines single strand conformational polymorphism (SSCP) analysis of one strand of a double-stranded amplification product with direct sequencing of the other. Using this strategy, which we find economical of both time and resources, we have identified a G to A transition, which substitutes a serine for glycine residue at position 862 in the major helix of the alpha 1 chain of Type I collagen. We use this mutation, which causes a lethal form of osteogenesis imperfecta, to illustrate the technique.

Tandem competitive polymerase chain reaction (TC-PCR): a method for determining ratios of RNA and DNA templates.

A sensitive and accurate method for determining the ratios of RNA and DNA templates by polymerase chain reaction (PCR) is presented. A common competitor containing tandemly arranged internal standards differing from the target template by the presence of different restriction enzyme sites is coamplified with the target templates under identical conditions. Products from each template and internal standard are identified by the band pattern after digestion with the restriction enzyme. As the amount of the common competitor is kept constant for all target templates, the ratio of PCR products from the templates reflects their ratio in the reaction mix before amplification. The method was used to study the relative abundance of mRNA for the pro-alpha1 and pro-alpha2 chains of type I collagen and for estimating disturbances of normal ratio in the inherited bone disorder, osteogenesis imperfecta.

Segregation analysis of the structural genes of the major fibrillar collagens provides further evidence of molecular heterogeneity in type II Ehlers-Danlos syndrome.

Type II Ehlers-Danlos syndrome (EDS) is one of a group of disorders characterized by striking abnormalities of the soft connective tissues. The major fibrillar collagens (types I and III) found in these tissues have important stress-bearing functions and abnormal collagen could therefore account for the clinical features of this condition. We have used a number of restriction site dimorphisms, tightly linked to the structural genes of type I collagen (COL1A1 COL1A2) and type III collagen (COL3A1), to investigate the segregation of corresponding alleles in three pedigrees in which type II EDS was clearly inherited as a dominant trait. Discordant segregation of all three collagen genes was seen in a large pedigree that included 17 affected individuals with the typical phenotype of type II EDS. Thus mutations in neither type I nor type III collagen genes were responsible for the disease in this family. In a second small pedigree discordant segregation of the disease with both type I collagen loci was observed while the concordant segregation seen at COL3A1 could easily have arisen by chance (P = 0.5). The third pedigree was uninformative at all three collagen loci because of inability to discriminate between the parental alleles. These results suggest that there may be molecular heterogeneity of type II EDS since abnormalities of type I collagen have been described in other individuals phenotypically similar to those in our study.