Inbreeding, pedigree size, and the most recent common ancestor of humanity

How many generations ago did the common ancestor of all present-day individuals live, and how does inbreeding affect this estimate? The number of ancestors within family trees determines the timing of the most recent common ancestor of humanity. However, mating is often non-random and inbreeding is ubiquitous in natural populations. Rates of pedigree growth are found for multiple types of inbreeding. This data is then combined with models of global population structure to estimate biparental coalescence times. When pedigrees for regular systems of mating are constructed, the growth rates of inbred populations contain Fibonacci n-step constants. The timing of the most recent common ancestor depends on global population structure, the mean rate of pedigree growth, mean fitness, and current population size. Inbreeding reduces the number of ancestors in a pedigree, pushing back global common ancestry times. These results are consistent with the remarkable findings of previous studies: all humanity shares common ancestry in the recent past.     

Genetic variation and evolutionary demography of Fenneropenaeus chinensis populations, as revealed by the analysis of mitochondrial control region sequences

Genetic variation and evolutionary demography of the shrimp Fenneropenaeus chinensis were investigated using sequence data of the complete mitochondrial control region (CR). Fragments of 993 bp of the CR were sequenced for 93 individuals from five localities over most of the species' range in the Yellow Sea and the Bohai Sea. There were 84 variable sites defining 68 haplotypes. Haplotype diversity levels were very high (0.95 ± 0.03-0.99 ± 0.02) in F. chinensis populations, whereas those of nucleotide diversity were moderate to low (0.66 ± 0.36%-0.84 ± 0.46%). Analysis of molecular variance and conventional population statistics (F(ST) ) revealed no significant genetic structure throughout the range of F. chinensis. Mismatch distribution, estimates of population parameters and neutrality tests revealed that the significant fluctuations and shallow coalescence of mtDNA genealogies observed were coincident with estimated demographic parameters and neutrality tests, in implying important past-population size fluctuations or range expansion. Isolation with Migration (IM) coalescence results suggest that F. chinensis, distributed along the coasts of northern China and the Korean Peninsula (about 1000 km apart), diverged recently, the estimated time-split being 12,800 (7,400-18,600) years ago.     

Estimation of parameters in large offspring number models and ratios of coalescence times

The ratio of singletons to the total number of segregating sites is used to estimate a reproduction parameter in a population model of large offspring numbers without having to jointly estimate the mutation rate. For neutral genetic variation, the ratio of singletons to the total number of segregating sites is equivalent to the ratio of total length of external branches to the total length of the gene genealogy. A multinomial maximum likelihood method that takes into account more frequency classes than just the singletons is developed to estimate the parameter of another large offspring number model. The performance of these methods with regard to sample size, mutation rate, and bias, is investigated by simulation. The expected value of the ratio of the total length of external branches to the total length of the whole tree is, using simulation, shown to decrease for the Kingman coalescent as sample size increases, but can increase or decrease, depending on parameter values, for Λ coalescents. Considering ratios of tree statistics, as opposed to considering lengths of various subtrees separately, can yield better insight into the dynamics of gene genealogies.     

Effect of oceanographic barriers and overfishing on the population genetic structure of the European spiny lobster (Palinurus elephas)

Defining population structure and genetic diversity levels is of the utmost importance for developing efficient conservation strategies. Overfishing has caused mean annual catches of the European spiny lobster (Palinurus elephas) to decrease alarmingly along its distribution area. In this context, there is a need for comprehensive studies aiming to evaluate the genetic health of the exploited populations. The present study is based on a set of ten nuclear markers amplified in 331 individuals from ten different localities covering most of P. elephas distribution area. Samples from Atlantic and Mediterranean basins showed small but significant differences, indicating that P. elephas populations do not behave as a single panmictic unit but form two partially‐overlapping groups. Despite intense overfishing, our dataset did not recover a recent bottleneck signal, and instead showed a large and stable historical effective size. This result could be accounted for by specific life‐history traits (reproduction and longevity) and the limitations of molecular markers in covering recent timescales for nontemporal samples. The findings of the present study emphasize the need to integrate information on effective population sizes and life‐history parameters when evaluating population connectivity levels from genetic data. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 104 , 407–418.     

Using all Gene Families Vastly Expands Data Available for Phylogenomic Inference

Traditionally, single-copy orthologs have been the gold standard in phylogenomics. Most phylogenomic studies identify putative single-copy orthologs using clustering approaches and retain families with a single sequence per species. This limits the amount of data available by excluding larger families. Recent advances have suggested several ways to include data from larger families. For instance, tree-based decomposition methods facilitate the extraction of orthologs from large families. Additionally, several methods for species tree inference are robust to the inclusion of paralogs and could use all of the data from larger families. Here, we explore the effects of using all families for phylogenetic inference by examining relationships among 26 primate species in detail and by analyzing five additional data sets. We compare single-copy families, orthologs extracted using tree-based decomposition approaches, and all families with all data. We explore several species tree inference methods, finding that identical trees are returned across nearly all subsets of the data and methods for primates. The relationships among Platyrrhini remain contentious; however, the species tree inference method matters more than the subset of data used. Using data from larger gene families drastically increases the number of genes available and leads to consistent estimates of branch lengths, nodal certainty and concordance, and inferences of introgression in primates. For the other data sets, topological inferences are consistent whether single-copy families or orthologs extracted using decomposition approaches are analyzed. Using larger gene families is a promising approach to include more data in phylogenomics without sacrificing accuracy, at least when high-quality genomes are available.     

King penguin demography since the last glaciation inferred from genome-wide data

How natural climate cycles, such as past glacial/interglacial patterns, have shaped species distributions at the high-latitude regions of the Southern Hemisphere is still largely unclear. Here, we show how the post-glacial warming following the Last Glacial Maximum (ca 18 000 years ago), allowed the (re)colonization of the fragmented sub-Antarctic habitat by an upper-level marine predator, the king penguin Aptenodytes patagonicus. Using restriction site-associated DNA sequencing and standard mitochondrial data, we tested the behaviour of subsets of anonymous nuclear loci in inferring past demography through coalescent-based and allele frequency spectrum analyses. Our results show that the king penguin population breeding on Crozet archipelago steeply increased in size, closely following the Holocene warming recorded in the Epica Dome C ice core. The following population growth can be explained by a threshold model in which the ecological requirements of this species (year-round ice-free habitat for breeding and access to a major source of food such as the Antarctic Polar Front) were met on Crozet soon after the Pleistocene/Holocene climatic transition.     

Coalescence with Background and Balancing Selection in Systems with Bi- and Uniparental Reproduction: Contrasting Partial Asexuality and Selfing

Uniparental reproduction in diploids, via asexual reproduction or selfing, reduces the independence with which separate loci are transmitted across generations. This is expected to increase the extent to which a neutral marker is affected by selection elsewhere in the genome. Such effects have previously been quantified in coalescent models involving selfing. Here we examine the effects of background selection and balancing selection in diploids capable of both sexual and asexual reproduction (i.e., partial asexuality). We find that the effect of background selection on reducing coalescent time (and effective population size) can be orders of magnitude greater when rates of sex are low than when sex is common. This is because asexuality enhances the effects of background selection through both a recombination effect and a segregation effect. We show that there are several reasons that the strength of background selection differs between systems with partial asexuality and those with comparable levels of uniparental reproduction via selfing. Expectations for reductions in N_e via background selection have been verified using stochastic simulations. In contrast to background selection, balancing selection increases the coalescence time for a linked neutral site. With partial asexuality, the effect of balancing selection is somewhat dependent upon the mode of selection (e.g., heterozygote advantage vs. negative frequency-dependent selection) in a manner that does not apply to selfing. This is because the frequency of heterozygotes, which are required for recombination onto alternative genetic backgrounds, is more dependent on the pattern of selection with partial asexuality than with selfing.     

INFERRING PHYLOGENIES FROM mtDNA VARIATION: MITOCHONDRIAL-GENE TREES VERSUS NUCLEAR-GENE TREES

An accurately resolved gene tree may not be congruent with the species tree because of lineage sorting of ancestral polymorphisms. DNA sequences from the mitochondrially encoded genes (mtDNA) are attractive sources of characters for estimating the phylogenies of recently evolved taxa because mtDNA evolves rapidly, but its utility is limited because the mitochondrial genes are inherited as a single linkage group (haplotype) and provide only one independent estimate of the species tree. In contrast, a set of nuclear genes can be selected from distinct chromosomes, such that each gene tree provides an independent estimate of the species tree. Another aspect of the gene-tree versus species-tree problem, however, favors the use of mtDNA for inferring species trees. For a three-species segment of a phylogeny, the branching order of a gene tree will correspond to that of the species tree if coalescence of the alleles or haplotypes occurred in the internode between the first and second bifurcation. From neutral theory, it is apparent that the probability of coalescence increases as effective population size decreases. Because the mitochondrial genome is maternally inherited and effectively haploid, its effective population size is one-fourth that of a nuclear-autosomal gene. Thus, the mitochondrial-haplotype tree has a substantially higher probability of accurately tracking a short internode than does a nuclear-autosomal-gene tree. When an internode is sufficiently long that the probability that the mitochondrial-haplotype tree will be congruent with the species tree is 0.95, the probability that a nuclear-autosomalgene tree will be congruent is only 0.62. If each of k independently sampled nuclear-gene trees has a probability of congruence with the species tree of 0.62, then a sample of 16 such trees would be required to be as confident of the inference based on the mitochondrial-haplotype tree. A survey of mtDNA-haplotype diversity in 34 species of birds indicates that coalescence is generally very recent, which suggests that coalescence times are typically much shorter than internodal branch lengths of the species tree, and that sorting of mtDNA lineages is not likely to confound the species tree. Hybridization resulting in transfer of mtDNA haplotypes among branches could also result in a haplotype tree that is incongruent with the species tree; if undetected, this could confound the species tree. However, hybridization is usually easy to detect and should be incorporated in the historical narrative of the group, because reticulation, as well as cladistic events, contributed to the evolution of the group.     

Marine biogeography of southern Australia: phylogeographical structure in a temperate sea-star

Aim To test whether marine biogeographical patterns observed at the community level are also important within species. It is postulated that historical hydrogeographic barriers have driven in situ diversification. Location The intertidal and shallow subtidal zones of southern Australia, New Zealand and nearby islands. Australia's temperate marine communities are characterized by a high degree of endemism and show strong biogeographical structure along an east–west axis. Methods Phylogeographical analysis of the widespread asteriid sea‐star Coscinasterias muricata Verrill across southern Australia and New Zealand. Forty‐two samples from 27 locations were included in phylogenetic analyses of mitochondrial (CO1; control region) and nuclear (ITS2) DNA sequences. Results Analysis of mtDNA revealed a deep phylogenetic split within Australian C. muricata, strongly correlated with latitude. ‘Northern’ haplotypes (latitude ≤ 37.6° S, nine sites, 15 samples) were 7.3–9.4% divergent from ‘southern’ haplotypes (latitude ≥ 37.6° S, 19 sites, 27 samples), consistent with late Pliocene separation. Eastern and western representatives of the ‘northern’ clade were 0.5–1.0% divergent, probably reflecting Pleistocene isolation. The ‘southern’ clade of Australia is also represented in New Zealand, indicating Pleistocene oceanic dispersal. Nuclear DNA (ITS2) sequences yielded relatively little phylogenetic resolution, but were generally congruent with mtDNA‐based groupings. Main conclusions The phylogeographical pattern detected within Australian C. muricata closely resembles marine biogeographical groupings proposed on the basis of community and species distributions. Recurring evolutionary patterns may have been driven by the hydrographic history of southern Australia. Specifically, we suggest that Plio‐Pleistocene temperature change and the repeated opening and closure of Bass Strait promoted allopatric divergence and perhaps cryptic speciation in C. muricata.