Inferring the recent ancestry of myostatin alleles affecting muscle mass in cattle

Double muscling is an inherited condition in cattle characterised by large increases in muscle mass. Mutations in the myostatin (MSTN) gene, responsible for double muscling, were targeted in this study to estimate the time since the most recent common ancestor (TMRCA) for Q204X (p.Gln204*), E226X (p.Glu226*), 821del11 (c.821del11), E291X (p.Glu291*), C313Y (p.Cys313Tyr) and the more phenotypically moderate F94L (p.Phe94Leu) mutation. Genetic variability was examined in eight regions upstream and downstream of the MSTN locus. The molecular distance of the homozygous region associated with each MSTN allele was used to estimate the TMRCA. Long homozygous segments were associated with the MSTN alleles (mostly > 2 Mb), compared to short segments (130 kb) for cattle wild type at the double muscling and F94L sites. Estimates of time indicated that each MSTN allele had a recent common ancestor (<400 years ago). The results from this study, and the increasing frequency of these MSTN alleles in some cattle breeds, demonstrate recent positive selection.     

Time to the most recent common ancestor and divergence times of populations of common chaffinches (Fringilla coelebs) in Europe and North Africa: insights into Pleistocene refugia and current levels of migration

We analyzed sequences from a 275-bp hypervariable region in the 5' end of the mitochondrial DNA control region in 190 common chaffinches (Fringilla coelebs) from 19 populations in Europe and North Africa, including new samples from Greece and Morocco. Coalescent techniques were applied to estimate the time to the most recent common ancestor (TMRCA) and divergence times of these populations. The first objective of this study was to infer the locations of refugia where chaffinches survived the last glacial episode, and this was achieved by estimating the TMRCA of populations in regions surrounding the Mediterranean that were unglaciated in the late Pleistocene. Although extant populations in Iberia, Corsica, Greece, and North Africa harbor haplotypes that are basal in a phylogenetic tree, this information alone cannot be used to infer that these localities served as refugia, because it is impossible to infer the ages of populations and their divergence times without also considering the population genetic processes of mutation, migration, and drift. Provided we assume the TMRCAs of populations are a reasonable estimate of a population's age, coalescent-based methods place resident populations in Iberia, Corsica, Greece, and North Africa during the time of the last glacial maximum, suggesting these regions served as refugia for the common chaffinch. The second objective was to determine when populations began diverging from each other and to use this as a baseline to estimate current levels of gene flow. Divergence time estimates suggest that European populations began diverging about 60,000 years before present. The relatively recent divergence of populations in North Africa, Italy, and Iberia may explain why classic migration estimates based on equilibrium assumptions are high for these populations. We compare these estimates with nonequilibrium-based estimates and show that the nonequilibrium estimates are consistently lower than the equilibrium estimates.     

A generalized framework to expand incomplete phylogenies using non-molecular phylogenetic information

Aim

The increasing availability of molecular information has lifted our understanding of species evolutionary relationships to unprecedent levels. However, current estimates of the world's biodiversity suggest that about a fifth of all extant species are yet to be described, and we still lack molecular information for many of the known species. Hence, evolutionary biologists will have to tackle phylogenetic uncertainty for a long time to come. This prospect has urged the development of software to expand phylogenies based on non-molecular phylogenetic information, and while the available tools provide some valuable features, major drawbacks persist and some of the proposed solutions are hardly generalizable to any group of organisms.

Innovation

Here, we present a completely generalized and flexible framework to expand incomplete phylogenies. The framework is implemented in the R package “randtip”, a toolkit of functions that was designed to randomly bind phylogenetically uncertain taxa in backbone phylogenies through a fully customizable and automatic procedure that uses taxonomic ranks as a major source of phylogenetic information. Although randtip can generate fully operative phylogenies for any group of organisms using just a list of species and a backbone tree, we stress that the “blind” expansion of phylogenies using “quick-and-dirty” approaches often leads to suboptimal solutions. Thus, we discuss a variety of circumstances that may require customizing simulation parameters beyond default settings to optimally expand the trees, including a detailed step-by-step tutorial that was designed to provide guidelines to non-specialist users.

Main Conclusions

Phylogenetic uncertainty should be tackled with caution, assessing potential pitfalls and opportunities to optimize parameter space prior to launch any simulation. Used judiciously, our framework will help evolutionary biologists to efficiently expand incomplete phylogenies and thereby account for phylogenetic uncertainty in quantitative analyses.     

Population structure and phylogeography of Aphyocypris kikuchii (Oshima) based on mitochondrial DNA variation

Aphyocypris kikuchii is a cyprinid species endemic to northern and eastern Taiwan and is the only primary freshwater fish native east of the Coastal Mountain Range. In total, 92 individuals of A. kikuchii from seven populations in three regions of the island were surveyed for mitochondrial DNA (mtDNA) variation. High haplotype diversity ( h = 0·989) and low nucleotide diversity ( π = 0·009) of mtDNA were detected. Negative values of Tajima's D and unimodal mismatch distributions probably reflect a history of recent demographic expansions from small populations. Three major haplotype clusters displayed geographically non-overlapping distributions, indicating a long-term isolation between regions. Hierarchical analysis of molecular variance showed significant genetic structuring among populations ( Φ _ST= 0·66). Significant haplotype heterogeneity was also detected among populations within regions ( Φ _SC= 0·41, P < 0·001) and among regions ( Φ _CT= 0·43, P < 0·05). Molecular clock estimates of coalescence in the three major mtDNA lineages indicated coalescence in the most recent common ancestor c. 0·11–0·39 million years ago. Haplotypes of cluster B nested as interior nodes in the haplotype network, indicating that migrations from Shueilian (SL) populations to the northern region (cluster A) and to the eastern region (cluster C) may have occurred independently. Lineages A and B + C should be managed as two distinct evolutionarily significant units, while the northern, SL and southern groups should be managed as separate management units.     

单细胞生物进化研究的进步

20世纪60年代,生物大致分为5界的谱系图,　经历了几次翻新,开始提出了线粒体和叶绿体的内共生学说。 由于分子生物学的发展,首先将各种生物的蛋白质的分子进行比较,构建成蛋白质的分子系统树。再转向核糖核酸,将核糖体的小亚单位,作为区别生命类型之间亲缘关系的指标。发现有些嗜极端条件的细菌,它们不同于原核生物也不同于真核生物,是第三种类型的生命形式。因此,在 80年代为生命建立了细菌、古细菌和真核生物三界的系统树。对许许多多单个基因的系统树的分析,又使人们认识到,古细菌与细菌和真核微生物之间以及各个物种之间,显然皆发生过大量的基因交换。对单细胞进化来说,基因除垂直传递外,横向的或叫侧向的基因转移也十分繁多。在一张完全的系统图中,要同时表现几千个不同的基因家族的超联结的种系型式才符合实际。因而,最新版本的系统树是分枝交缠、无主干的。 Abstract:In 1960s,kimgdoms of organisms were charted generally in a five branching form.Later,the endosymbiont hypothesis for the mitochondria and the chloroplast was proposed.The life-form is divided into two forms,the prokaryotes (bacteria) and the eukaryotes.The study of the molecular biology made the progress faster.In 1980s,Woese,CR.asserted that two-domain view of life was no longer true,a three-domain construct,the Bacteria,the Archaea,and the Eukaryotes had to take its place.At first,phylogeny trees based on differences in the amino acid sequences,then among ribosomal RNAs and also nuclear gene from hundreds of microbial species were depicted and many mini phylogenetic trees grouped the species according to their differences in the sequences.It was found that they shared genes between their contemporaries and across the species barriers.At the root of the phylogeny tree,there was not a single common cell,it was replaced by a common ancestral community of primitive cells.Genes transfered rather freely as the transposons swapping between those cells.There was no last universal common ancestor of single cell that could be found in the revised Tree of Life,It was not easy to represent the genealogical patterns of thousands of different families of genes,in one systematic map,therefor there was no trunk at all.     

The Variance of Coalescent Time Estimates from DNA Sequences

A method for estimating the variance of a simple estimate of the time of the most recent common ancestor of a sample of sequences is described. [Reviewing Editor: Dr. Rasmus Nielson]     

PCA and Clustering Reveal Alternate mtDNA Phylogeny of N and M Clades

Phylogenetic trees based on mtDNA polymorphisms are often used to infer the history of recent human migrations. However, there is no consensus on which method to use. Most methods make strong assumptions which may bias the choice of polymorphisms and result in computational complexity which limits the analysis to a few samples/polymorphisms. For example, parsimony minimizes the number of mutations, which biases the results to minimizing homoplasy events. Such biases may miss the global structure of the polymorphisms altogether, with the risk of identifying a "common" polymorphism as ancient without an internal check on whether it either is homoplasic or is identified as ancient because of sampling bias (from oversampling the population with the polymorphism). A signature of this problem is that different methods applied to the same data or the same method applied to different datasets results in different tree topologies. When the results of such analyses are combined, the consensus trees have a low internal branch consensus. We determine human mtDNA phylogeny from 1737 complete sequences using a new, direct method based on principal component analysis (PCA) and unsupervised consensus ensemble clustering. PCA identifies polymorphisms representing robust variations in the data and consensus ensemble clustering creates stable haplogroup clusters. The tree is obtained from the bifurcating network obtained when the data are split into k = 2,3,4,...,kmax clusters, with equal sampling from each haplogroup. Our method assumes only that the data can be clustered into groups based on mutations, is fast, is stable to sample perturbation, uses all significant polymorphisms in the data, works for arbitrary sample sizes, and avoids sample choice and haplogroup size bias. The internal branches of our tree have a 90% consensus accuracy. In conclusion, our tree recreates the standard phylogeny of the N, M, L0/L1, L2, and L3 clades, confirming the African origin of modern humans and showing that the M and N clades arose in almost coincident migrations. However, the N clade haplogroups split along an East-West geographic divide, with a "European R clade" containing the haplogroups H, V, H/V, J, T, and U and a "Eurasian N subclade" including haplogroups B, R5, F, A, N9, I, W, and X. The haplogroup pairs (N9a, N9b) and (M7a, M7b) within N and M are placed in nonnearest locations in agreement with their expected large TMRCA from studies of their migrations into Japan. For comparison, we also construct consensus maximum likelihood, parsimony, neighbor joining, and UPGMA-based trees using the same polymorphisms and show that these methods give consistent results only for the clade tree. For recent branches, the consensus accuracy for these methods is in the range of 1-20%. From a comparison of our haplogroups to two chimp and one bonobo sequences, and assuming a chimp-human coalescent time of 5 million years before present, we find a human mtDNA TMRCA of 206,000 +/- 14,000 years before present.     

Molecular evolution and multilocus sequence typing of 145 strains of SARS-CoV

In this study, we have identified 876 polymorphism sites in 145 complete or partial genomes of SARS-CoV available in the NCBI GenBank. One hundred and seventy-four of these sites existed in two or more SARS-CoV genome sequences. According to the sequence polymorphism, all SARS-CoVs can be divided into three groups: (I) group 1, animal-origin viruses (such as SARS-CoV SZ1, SZ3, SZ13 and SZ16); (II) group 2, all viruses with clinical origin during first epidemic; and (III) group 3, SARS-CoV GD03T0013. According to 10 special loci, group 2 again can be divided into genotypes C and T, which can be further divided into sub-genotypes C1-C4 and T1-T4. Positive Darwinian selections were identified between any pair of these three groups. Genotype C gives neutral selection. Genotype T, however, shows negative selection. By comparing the death rates of SARS patients in the different regions, it was found that the death rate caused by the viruses of the genotype C was lower than that of the genotype T. SARS-CoVs might originate from an unknown ancestor.     

Recent Evolutionary History of HIV-1 Subtype B—Rebuttal

Abstract The history of the HIV-1 B epidemic is the subject of a continuing debate. Did the epidemic start in the 1970s, as it was established based on the epidemiological data, or decades earlier, as it was suggested based on the analysis of nucleotide distances in the env gene? Our study [Lukashov and Goudsmit,J Mol Evol (2002) 54: 680-691] found that the overestimation of the age of the epidemic in the analysis of env sequences was a bias resulting from the non-clock-like evolution at nonsynonymous sites, while the estimates based on synonymous substitutions agreed with the results of epidemiological studies. Besides the principal difference between the evolution of synonymous and nonsynonymous sites, several issues have to be addressed: (i) the onset of the HIV/AIDS epidemic, and not the circulation of the preepidemic viruses, should be taken as the gold standard for the timeline of HIV diversification; (ii) the circulation of ancient, preepidemic, viruses, whether long- or short-term, does not influence the increase of HIV divergence during the epidemic; and (iii) application of the same random latency probability for all viruses, irrespective to their age and distance from the common ancestor, biases the estimation of the age of the epidemic.     

Studies of the ADP/ATP carrier of mitochondria with fluorescent ADP analogue formycin diphosphate

The ADP/ATP carrier was studied by a fluorescent substrate, formycin diphosphate which is the only fluorescent ADP analogue to bind. Its low quantum yield, short decay time and spectral overlap with tryptophan has as yet prevented its wider use.By incorporating fluorescent acceptors of formycin diphosphate fluorescence, anthracene-maleimide and vinylanthracene, into the membrane, these difficulties were circumvented. Only bound formycin diphosphate transfers energy to the probes so that the secondary emission of these probes is a measure for membrane-bound formycin diphosphate.The fluorescent transfer is inhibited by ADP, bongkrekate and carboxy-atractylate whether added before or after incubation of formycin diphosphate showing that only binding to the adenine nucleotide carrier is measured. It also shows directly that the earlier demonstrated ADP fixation by bongkrekate is indeed a displacement into the matrix.The fluorescence decay time of the bound formycin diphosphate is measured as 1.95 ns compared to 0.95 ns of the free formycin diphosphate, indicating that formycin diphosphate is bound at the carrier in a non-polar environment.The depolarization decay time was found to be larger than 15 ns, indicating that carrier-bound formycin diphosphate is immobile within this time period.     

Meeting the challenges of on-host and off-host water balance in blood-feeding arthropods

In this review, we describe water balance requirements of blood-feeding arthropods, particularly contrasting dehydration tolerance during the unfed, off-host state and the challenges of excess water that accompany receipt of the bloodmeal. Most basic water balance characteristics during the off-host stage are applicable to other terrestrial arthropods, as well. A well-coordinated suite of responses enable arthropods to conserve water resources, enhance their desiccation tolerance, and increase their water supplies by employing a diverse array of molecular, structural and behavioral responses. Water loss rates during the off-host phase are particularly useful for generating a scheme to classify vectors according to their habitat requirements for water, thus providing a convenient tool with potential predictive power for defining suitable current and future vector habitats. Blood-feeding elicits an entirely different set of challenges as the vector responds to overhydration by quickly increasing its rate of cuticular water loss and elevating the rate of diuresis to void excess water and condense the bloodmeal. Immature stages that feed on blood normally have a net increase in water content at the end of a blood-feeding cycle, but in adults the water content reverts to the pre-feeding level when the cycle is completed. Common themes are evident in diverse arthropods that feed on blood, particularly the physiological mechanisms used to respond to the sudden influx of water as well as the mechanisms used to counter water shortfalls that are encountered during the non-feeding, off-host state.     

The adenylate energy charge in marine phytoplantkon: The effect of temperature on the physiological state of Skeletonema costatum (Grev.) Cleve

The adenylate energy charge $\text{([ATP] +}\text{1}\text{2}\text{[ADP])}\text{[0ATP+ADP+AMP]}$ was measured in axenic batch cultures of Skeletonema costatum (Grev.) Cleve at 2°, 10°, 15°, 20°, 24° and 30°C. The results suggest that this eurythermal diatom is physiologically capable of adapting to the 28 °C range of temperature with little apparent difference in the potential energy available to the cell. In N-limited continuous cultures at 15 °C, the energy charge values were lower than those observed in batch culture by 0.2, implying nutrient stress may result in decreased intracellular chemical energy. The utilization of the adenylate energy charge as an indicator of physiological state is suggested.     

Effects of the herbicide atrazine on adenine nucleotide levels in Zostera marina L. (eelgrass)

Response of adenine nucleotides (ATP, ADP, AMP) and adenylate energy charge (EC) to atrazine, a triazine herbicide, was evaluated as an indicator of metabolic state in Zostera marina L. (eelgrass), a submerged marine angiosperm. Short-term (6 h) atrazine stress reduced ATP and total adenylates (AT) at both 10 and 100 ppb, but EC remained constant. Net productivity decreased at 100, but not at 10 ppb atrazine over 6 h. Long-term (21 day) atrazine stress was evidenced by growth inhibition and 50% mortality near 100 ppb. EC was reduced at 0.1, 1.0 and 10 ppb atrazine, but ATP and EC increased with physiological response to severe stress (100 ppb) after 21 days. Apparently, ATP and AT decrease over the short-term but rebound over the long-term with severe atrazine stress, increasing beyond control levels before plant death results. Supplementing adenine nucleotide and EC results with more conventional quantitative analyses should afford greater knowledge of physiological response to environmental variation.     

Invasion history and demographic pattern of Cryphonectria hypovirus 1 across European populations of the chestnut blight fungus

We reconstructed the invasion history of the fungal virus Cryphonectria hypovirus 1 (CHV‐1) in Europe, which infects the chestnut blight fungus Cryphonectria parasitica. The pattern of virus evolution was inferred based on nucleotide sequence variation from isolates sampled across a wide area in Europe at different points in time. Phylogeny and time estimates suggested that CHV‐1 was introduced together with its fungal host to Europe and that it rapidly colonized the central range along the south facing slopes of the Alps and the north‐east facing slopes of the Dinaric Alps. These central populations were the source for two waves of simultaneous invasions toward the southern Balkans and Turkey, as indicated by migration rates. Our results showed that the evolutionary scenarios for CHV‐1 and C. parasitica were spatially congruent. As infection with CHV‐1 reduces the pathogenicity of C. parasitica toward the chestnut tree, CHV‐1 invasions of the newly established C. parasitica populations probably prevented the development of devastating chestnut blight epidemics in Europe. We propose that in this, and supposedly in other pathosystems, geographic, vegetation‐related, demographic, economic, and political factors may help explain the correlated invasion pattern of a parasite and its host.     

The Evolution and Structure Prediction of Coiled Coils across All Genomes

Coiled coils are α-helical interactions found in many natural proteins. Various sequence-based coiled-coil predictors are available, but key issues remain: oligomeric state and protein-protein interface prediction and extension to all genomes. We present SpiriCoil (http://supfam.org/SUPERFAMILY/spiricoil), which is based on a novel approach to the coiled-coil prediction problem for coiled coils that fall into known superfamilies: hundreds of hidden Markov models representing coiled-coil-containing domain families. Using whole domains gives the advantage that sequences flanking the coiled coils help. SpiriCoil performs at least as well as existing methods at detecting coiled coils and significantly advances the state of the art for oligomer state prediction. SpiriCoil has been run on over 16 million sequences, including all completely sequenced genomes (more than 1200), and a resulting Web interface supplies data downloads, alignments, scores, oligomeric state classifications, three-dimensional homology models and visualisation. This has allowed, for the first time, a genomewide analysis of coiled-coil evolution. We found that coiled coils have arisen independently de novo well over a hundred times, and these are observed in 16 different oligomeric states. Coiled coils in almost all oligomeric states were present in the last universal common ancestor of life. The vast majority of occasions that individual coiled coils have arisen de novo were before the last universal common ancestor of life; we do, however, observe scattered instances throughout subsequent evolutionary history, mostly in the formation of the eukaryote superkingdom. Coiled coils do not change their oligomeric state over evolution and did not evolve from the rearrangement of existing helices in proteins; coiled coils were forged in unison with the fold of the whole protein.