Saturday, 5 January 2013

Chromosomal rearrangements and human-chimp speciation

Mol Biol Evol (2012) doi: 10.1093/molbev/mss272

Recombination Rates and Genomic Shuffling in Human and Chimpanzee—A New Twist in the Chromosomal Speciation Theory

Marta Farré et al.

A long-standing question in evolutionary biology concerns the effect of recombination in shaping the genomic architecture of organisms and, in particular, how this impacts the speciation process. Despite efforts employed in the last decade, the role of chromosomal reorganizations in the human–chimpanzee speciation process remains unresolved. Through whole-genome comparisons, we have analyzed the genome-wide impact of genomic shuffling in the distribution of human recombination rates during the human–chimpanzee speciation process. We have constructed a highly refined map of the reorganizations and evolutionary breakpoint regions in the human and chimpanzee genomes based on orthologous genes and genome sequence alignments. The analysis of the most recent human and chimpanzee recombination maps inferred from genome-wide single-nucleotide polymorphism data revealed that the standardized recombination rate was significantly lower in rearranged than in collinear chromosomes. In fact, rearranged chromosomes presented significantly lower recombination rates than chromosomes that have been maintained since the ancestor of great apes, and this was related with the lineage in which they become fixed. Importantly, inverted regions had lower recombination rates than collinear and noninverted regions, independently of the effect of centromeres. Our observations have implications for the chromosomal speciation theory, providing new evidences for the contribution of inversions in suppressing recombination in mammals.

Link

Bulging modern human foreheads

AJPA DOI: 10.1002/ajpa.22202

Geometric variation of the frontal squama in the genus homo: Frontal bulging and the origin of modern human morphology

Emiliano Bruner et al.

The majority of studies of frontal bone morphology in paleoanthropology have analyzed the frontal squama and the browridge as a single unit, mixing information from different functional elements. Taking into account that the bulging of the frontal bone is often described as a species-specific trait of Homo sapiens, in this article we analyze variation in the midsagittal profile of the genus Homo, focusing on the frontal squama alone, using landmark-based superimpositions and principal components analysis. Our results demonstrate that anatomically modern humans are definitely separated from extinct human taxa on the basis of frontal bulging. However, there is minor overlap among these groups, indicating that it is necessary to exercise caution when using this trait alone to make taxonomic inferences on individual specimens. Early modern humans do not show differences with recent modern humans, and “transitional” individuals such as Jebel Irhoud 1, Maba, and Florisbad, show modern-like frontal squama morphology. The bulging of the frontal squama in modern humans may represent a structural consequence of more general cranial changes, or it could be a response to changes in the morphology of the underlying prefrontal brain elements. A subtle difference between Neandertals and the Afro-European Middle Pleistocene Homo sample is associated with flattening at bregma in the former group, a result that merits further investigation.

Link

Friday, 4 January 2013

Deep whole-genome sequencing of 100 Malays

The 1000 Genomes Project is the largest collection of full human genomes currently available, but most of its 2.5k samples have been sequenced at low coverage. One downside of this is that infrequent variants are often missed. If an individual is polymorphic at some site, then the chance of detecting this polymorphism increases with the number of reads covering that site. If a number of individuals are sampled, then polymorphisms that are common in the population will probably be detected in a few individuals even if a low number of reads is used for each of them; but, if they are infrequent, then they are more likely to be missed. Hence, low-coverage sequencing of population samples will tend to find common variants and will tend to miss less common variants relative to high-coverage sequencing.

This idea is intuitively correct, but the question of the added power of high-coverage sequencing to detect variants can only be addressed by giving the same individuals both low- and high-coverage sequencing. This is the topic of a new paper in AJHG which creates a useful comparison benchmark for the performance of the two types of sequencing methods. High-coverage sequencing may be needed for things like disease studies (because deleterious alleles tend to be low-frequency), or the study of recent human demography (because recent population growth has resulted in an abundance of low-frequency SNPs that have not had enough time to reach a high population frequency yet).

AJHG dx.doi.org/10.1016/j.ajhg.2012.12.005

Deep Whole-Genome Sequencing of 100 Southeast Asian Malays

Lai-Ping Wong et al.


Whole-genome sequencing across multiple samples in a population provides an unprecedented opportunity for comprehensively characterizing the polymorphic variants in the population. Although the 1000 Genomes Project (1KGP) has offered brief insights into the value of population-level sequencing, the low coverage has compromised the ability to confidently detect rare and low-frequency variants. In addition, the composition of populations in the 1KGP is not complete, despite the fact that the study design has been extended to more than 2,500 samples from more than 20 population groups. The Malays are one of the Austronesian groups predominantly present in Southeast Asia and Oceania, and the Singapore Sequencing Malay Project (SSMP) aims to perform deep whole-genome sequencing of 100 healthy Malays. By sequencing at a minimum of 30? coverage, we have illustrated the higher sensitivity at detecting low-frequency and rare variants and the ability to investigate the presence of hotspots of functional mutations. Compared to the low-pass sequencing in the 1KGP, the deeper coverage allows more functional variants to be identified for each person. A comparison of the fidelity of genotype imputation of Malays indicated that a population-specific reference panel, such as the SSMP, outperforms a cosmopolitan panel with larger number of individuals for common SNPs. For lower-frequency (less than 5%) markers, a larger number of individuals might have to be whole-genome sequenced so that the accuracy currently afforded by the 1KGP can be achieved. The SSMP data are expected to be the benchmark for evaluating the value of deep population-level sequencing versus low-pass sequencing, especially in populations that are poorly represented in population-genetics studies.

Link

Thursday, 3 January 2013

Body form variation of prehistoric Jomon (Fukase et al. 2012)

Am J Phys Anthropol DOI: 10.1002/ajpa.22112

Geographic variation in body form of prehistoric Jomon males in the Japanese archipelago: Its ecogeographic implications 

Hitoshi Fukase et al.

Diversity of human body size and shape is often biogeographically interpreted in association with climatic conditions. According to Bergmann's and Allen's rules, populations in regions with a cold climate are expected to display an overall larger body and smaller/shorter extremities than those in warm/hot environments. In the present study, the skeletal limb size and proportions of prehistoric Jomon hunter-gatherers, who extensively inhabited subarctic to subtropical areas in the ancient Japanese archipelago, were examined to evaluate whether or not the inter-regional differences follow such ecogeographic patterns. Results showed that the Jomon intralimb proportions including relative distal limb lengths did not differ significantly among five regions from northern Hokkaido to the southern Okinawa Islands. This suggests a limited co-variability of the intralimb proportions with climate, particularly within genealogically close populations. In contrast, femoral head breadth (associated with body mass) and skeletal limb lengths were found to be significantly and positively correlated with latitude, suggesting a north-south geographical cline in the body size. This gradient therefore comprehensively conforms to Bergmann's rule, and may stem from multiple potential factors such as phylogenetic constraints, microevolutionary adaptation to climatic/geographic conditions during the Jomon period, and nutritional and physiological response during ontogeny. Specifically, the remarkably small-bodied Jomon in the Okinawa Islands can also be explained as an adjustment to subtropical and insular environments. Thus, the findings obtained in this study indicate that Jomon people, while maintaining fundamental intralimb proportions, displayed body size variation in concert with ambient surroundings.


Link

Tuesday, 1 January 2013

Y-chromosome and mtDNA of Henri IV

A recent paper had determined the Y-chromosome haplotype of Louis XVI of France from a handkerchief preserving his blood after his execution. A new study looks at the mummified head of Henri IV, the first Bourbon King of France. Even though only a limited number of Y-STRs were successfully typed, they match those of Louis XVI, who belonged to the not-so-frequent-anymore haplogroup G2a. So, while we cannot be entirely sure that the two Y-chromosomes were related in a genealogical time frame, the evidence is consistent with their known genealogical relationship and with the attribution of the two samples (mummified head/blood) to the respective kings.

Also of interest, Henri IV's mtDNA haplotype:
The majority of the clones generated show an U5b* mtDNA haplotype defined by three nucleotide changes at positions 16239T 16270T 16311C (see Supplementary material). The three HVR1 diagnostic positions were confirmed in two different amplifications of the L16185-H16378 HVR1 fragment, proving that the results are reproducible. This mtDNA haplotype is present so far in one single individual from France (originally published in [10]) in an in-house database of 22,807 published European sequences, and it is absent in all people involved in the laboratory analysis.
 If I followed the trail of ancestry correctly, this matrilineage leads all the way to a Tochter von Egisheim in the 11th century.

Forensic Science International Available online 30 December 2012

Genetic comparison of the head of Henri IV and the presumptive blood from Louis XVI (both Kings of France)

Philippe Charlier et al.

A mummified head was identified in 2010 as belonging to Henri IV, King of France. A putative blood sample from the King Louis XVI preserved into a pyrographically decorated gourd was analyzed in 2011. Both kings are in a direct male-line descent, separated by seven generations. We have retrieved the hypervariable region 1 of the mitochondrial DNA as well as a partial Y-chromosome profile from Henri IV. Five STR loci match the alleles found in Louis XVI, while another locus shows an allele that is just one mutation step apart. Taking into consideration that the partial Y-chromosome profile is extremely rare in modern human databases, we concluded that both males could be paternally related. The likelihood ratio of the two samples belonging to males separated by seven generations (as opposed to unrelated males) was estimated as 246.3, with a 95% confidence interval between 44.2 and 9729. Historically speaking, this forensic DNA data would confirm the identity of the previous Louis XVI sample, and give another positive argument for the authenticity of the head of Henri IV.

Link

Mating between Modern Humans, Neanderthals and other Archaics (Waddell & Tan 2012)

arXiv:1212.6820 [q-bio.GN]

New g%AIC, g%AICc, g%BIC, and Power Divergence Fit Statistics Expose Mating between Modern Humans, Neanderthals and other Archaics

Peter J. Waddell, Xi Tan

The purpose of this article is to look at how information criteria, such as AIC and BIC, relate to the g%SD fit criterion derived in Waddell et al. (2007, 2010a). The g%SD criterion measures the fit of data to model based on a normalized weighted root mean square percentage deviation between the observed data and model estimates of the data, with g%SD = 0 being a perfectly fitting model. However, this criterion may not be adjusting for the number of parameters in the model comprehensively. Thus, its relationship to more traditional measures for maximizing useful information in a model, including AIC and BIC, are examined. This results in an extended set of fit criteria including g%AIC and g%BIC. Further, a broader range of asymptotically most powerful fit criteria of the power divergence family, which includes maximum likelihood (or minimum G^2) and minimum X^2 modeling as special cases, are used to replace the sum of squares fit criterion within the g%SD criterion. Results are illustrated with a set of genetic distances looking particularly at a range of Jewish populations, plus a genomic data set that looks at how Neanderthals and Denisovans are related to each other and modern humans. Evidence that Homo erectus may have left a significant fraction of its genome within the Denisovan is shown to persist with the new modeling criteria.

Link

Tuesday, 4 December 2012

Disentangling the histories of mtDNA haplogroups M1 and U6

mtDNA haplogroups M1 and U6 are often mentioned in terms of Eurasian back-migration in Africa. The former is the only clade of the Asian haplogroup M which occurs in Africa at all; the latter is the only clade of the West Eurasian haplogroup U that does the same. These haplogroups also tend to co-exist in North and East Africa, although they are largely absent in sub-Saharan Africa. Different ideas have been offered for their occurrence, including a "Paleolithic" spread or a more recent one associated with the spread of Afroasiatic languages.

The new paper offers useful new data on this debate. The most important conclusion is that despite their oft-mentioned association, these two haplogroups appear to have distinct histories. One argument for this is their separate geographic distribution:


M1 (on panel A) is much more common in Northeast Africa and the Near East (including the Caucasus), whereas U6 (panel B) is more confined in Africa, and has its stronger peak in NW Africa, being rare in NE Africa.

An interesting aside, is that all the mysterious M1 from the Caucasus belongs to subclade M1a, while the smaller M1b clade tends to co-occur with M1a in other parts of Africa and the Near East. This indicates a founder effect for the origin of Caucasian M1a, but leaves open the issue of the immediate origins of M1. Hopefully it will become possible to place this haplogroup within the broader M phylogeny in the future.

The Bayesian skyline plots also contrast M1 and U6 in terms of their demographic histories:



The authors argue that these histories are inconsistent with either a very early dispersal history with the Dabban industry, as well as a more recent spread with Afroasiatic. From the paper:
The transition from the Middle Palaeolithic to Upper Palaeolithic in North Africa is characterised by the appearance of the “Dabban”, an industry that is restricted to Cyrenaica in northeast Libya and represented at the caves of Hagfet ed Dabba and Haua Fteah [19]. Whilst a techno-typological shift occurred within the Dabban ~33 KYA [19], starker changes in the archaeological record occurred throughout North Africa and Southwest Asia ~23-20 KYA, represented by the widespread appearance of backed bladelet technologies. The appearance of these backed bladelet industries more or less coincides with the timing of the Last Glacial Maximum (LGM) (~23-18 KYA), including: ~21 KYA in Upper Egypt [20]; ~20 KYA at Haua Fteah with the Oranian [21]; the Iberomaurusian expansion in the Jebel Gharbi ~20 KYA [22]; and the first Iberomaurusian at Tamar Hat in Algeria ~20 KYA [23]. The earliest Iberomaurusian sites in Morocco appear to be only slightly younger ~18 KYA [24].
A disassociation of these haplogroups from the UP in North Africa might be consistent with my idea that the UP was in part a cultural revolution that spread not only with people, but often with ideas across a species that already had the "biological machinery" for behavioral modernity and was already established in both Africa and the Near East.

As for the connection to Afroasiatic, the authors detect a linguistic correlation with M1a, which, however, appears too old to have been involved directly in the spread of this language family:
Concerning haplogroup M1 individually, a significant correlation with languages was observed. Furthermore, within M1, it appears that the correlation is mostly due to M1a. However, given the small sample size of M1b, any potential signal correlating with language might not be detectable. Interestingly, M1a has a likely East African origin, but its coalescent age of ~21 KYA still largely predates that of the proto-AA. Maybe a sub-clade of M1a would still give a similar correlation, but there are not sufficient samples to allow splitting M1a into its various sub-clades, and to test for a correlation. Although we found a correlation, limited sample sizes do not allow drawing unambiguous connection between genes and languages. Furthermore, it is also possible that this putative sub-clade of M1 does not testify for the expansion of AA speaking people, but was already present among the people who inhabited the area before the spread of the AA languages.
Personally, I am in favor of an East African origin of Afroasiatic, as this makes sense of various lines of evidence, one of which is the African shift of the "Southwest_Asian" component that is modal in Semitic populations. I envision that M1 was geographically circumscribed in a NE African population after its much earlier arrival from Asia and piggy-backed onto the expansion of Afroasiatic speakers, thus explaining the observed correlation. A good analogy would be with the expansion of, say, haplogroup H in the Americas which piggybacked on the European colonization, even though the coalescence age of H predates the arrival of Europeans in the New World by many millennia.

BMC Evolutionary Biology 2012, 12:234 doi:10.1186/1471-2148-12-234


Divorcing the Late Upper Palaeolithic demographic histories of mtDNA haplogroups M1 and U6 in Africa

Erwan Pennarun et al.

Abstract (provisional)
Background
A Southwest Asian origin and dispersal to North Africa in the Early Upper Palaeolithic era has been inferred in previous studies for mtDNA haplogroups M1 and U6. Both haplogroups have been proposed to show similar geographic patterns and shared demographic histories.

Results
We report here 24 M1 and 33 U6 new complete mtDNA sequences that allow us to refine the existing phylogeny of these haplogroups. The resulting phylogenetic information was used to genotype a further 131 M1 and 91 U6 samples to determine the geographic spread of their sub-clades. No southwest Asian specific clades for M1 or U6 were discovered. U6 and M1 frequencies in North Africa, the Middle East and Europe do not follow similar patterns, and their sub-clade divisions do not appear to be compatible with their shared history reaching back to the Early Upper Palaeolithic. The Bayesian Skyline Plots testify to non-overlapping phases of expansion, and the haplogroups' phylogenies suggest that there are U6 sub-clades that expanded earlier than those in M1. Some M1 and U6 sub-clades could be linked with certain events. For example, U6a1 and M1b, with their coalescent ages of ~20,000-22,000 years ago and earliest inferred expansion in northwest Africa, could coincide with the flourishing of the Iberomaurusian industry, whilst U6b and M1b1 appeared at the time of the Capsian culture.

Conclusions
Our high-resolution phylogenetic dissection of both haplogroups and coalescent time assessments suggest that the extant main branching pattern of both haplogroups arose and diversified in the mid-later Upper Palaeolithic, with some sub-clades concomitantly with the expansion of the Iberomaurusian industry. Carriers of these maternal lineages have been later absorbed into and diversified further during the spread of Afro-Asiatic languages in North and East Africa.

Link