Site-Specific Evolutionary Rate Inference: Taking Phylogenetic Uncertainty into Account

The evolutionary rate at an amino acid site is indicative of how conserved this site is and, in turn, allows evaluating the importance of this site in maintaining the structure/function of the protein. When evolutionary rates are estimated, one must reconstruct the phylogenetic tree describing the evolutionary relationship among the sequences under study. However, if the inferred phylogenetic tree is incorrect, it can lead to erroneous site-specific rate estimates. Here we describe a novel Bayesian method that uses Markov chain Monte Carlo methodology to integrate over the space of all possible trees and model parameters. By doing so, the method considers alternative evolutionary scenarios weighted by their posterior probabilities. We show that this comprehensive evolutionary approach is superior over methods that are based on only a single tree. We illustrate the potential of our algorithm by analyzing the conservation pattern of the potassium channel protein family.Itay Mayrose, Amir Mitchell contributed equal. Reviewing Editor : Dr. Nicolas Galtier     

Wing-Click Sounds of Heliconius cydno alithea (Nymphalidae: Heliconiinae) Butterflies

Field-collected Heliconius cydno Doubleday females were observed producing audible wing clicks during encounters between conspecifics in greenhouses in a large insectary during the day and at roosting time. Occasionally, these females also were observed producing sounds in aggressive encounters with females of a close relative, H. erato (L). However, the wing-clicks were not observed subsequently from first-generation adults born in the greenhouses. The sounds were produced in short trains of 3–10 wing-clicks at the rate of 10 clicks/s. The individual clicks had a mean duration of 1.48 ms and a broad frequency spectrum, with a peak near 1075 Hz. This peak lies near the 1200-Hz frequency of maximal sensitivity measured previously for auditory neurons of H. erato. The production of these previously unreported sounds suggests that wing clicks may play a role in both intra- and interspecific communication among Heliconius species.     

Quantifying Phylogenetic Beta Diversity: Distinguishing between 'True' Turnover of Lineages and Phylogenetic Diversity Gradients

The evolutionary dissimilarity between communities (phylogenetic beta diversity PBD) has been increasingly explored by ecologists and biogeographers to assess the relative roles of ecological and evolutionary processes in structuring natural communities. Among PBD measures, the PhyloSor and UniFrac indices have been widely used to assess the level of turnover of lineages over geographical and environmental gradients. However, these indices can be considered as 'broad-sense' measures of phylogenetic turnover as they incorporate different aspects of differences in evolutionary history between communities that may be attributable to phylogenetic diversity gradients. In the present study, we extend an additive partitioning framework proposed for compositional beta diversity to PBD. Specifically, we decomposed the PhyloSor and UniFrac indices into two separate components accounting for 'true' phylogenetic turnover and phylogenetic diversity gradients, respectively. We illustrated the relevance of this framework using simple theoretical and archetypal examples, as well as an empirical study based on coral reef fish communities. Overall, our results suggest that using PhyloSor and UniFrac may greatly over-estimate the level of spatial turnover of lineages if the two compared communities show contrasting levels of phylogenetic diversity. We therefore recommend that future studies use the 'true' phylogenetic turnover component of these indices when the studied communities encompass a large phylogenetic diversity gradient.     

Behavioral Adaptations to Pollen-Feeding in Heliconius Butterflies (Nymphalidae, Heliconiinae): An Experiment Using Lantana Flowers

Butterflies in the genera Heliconius and Laparus obtain fitness-related benefits from using amino acids derived from pollen. These butterflies have morphological features of the proboscis that facilitate pollen-feeding. Here we investigate behavioral characteristics potentially involved in pollen-gathering. Analysis of four behavioral characters showed that pollen-feeding species manipulate Lantana flowers faster and more thoroughly compared to non-pollen-feeding relatives. Although this suggests that pollen-feeding species are potentially more efficient in harvesting pollen, every butterfly tested successfully removed pollen from Lantana and non-pollen-feeding butterflies generally extracted larger amounts of pollen than Heliconius and Laparus. Morphological characteristics of the proboscis and the production of abundant fluid exudate help keep pollen attached in the proboscis for long periods of time—possibly key to Heliconius' and Laparus' ability to obtain amino acids from pollen. Our results, in concert with those of previous morphological analysis, indicate that behavioral and structural attributes associated with pollen-feeding in Heliconius and Laparus are subtle modifications of widespread butterfly characteristics and raise the question why other butterflies do not use pollen in their diet.     

陕西延安市区蝶类多样性初探

延安市区有蝶类6科25属36种,其中凤蝶科2属3种;粉蝶科5属8种;眼蝶科3属5种;蛱蝶科7属10种;灰蝶科7属9种;弄蝶科1属1种。陕西省新记录5种。区系分析表明,该地区以广布种为优势种。     

Quantifying the Diversification of Hepatitis C Virus (HCV) during Primary Infection: Estimates of the In Vivo Mutation Rate

Hepatitis C virus (HCV) is present in the host with multiple variants generated by its error prone RNA-dependent RNA polymerase. Little is known about the initial viral diversification and the viral life cycle processes that influence diversity. We studied the diversification of HCV during acute infection in 17 plasma donors, with frequent sampling early in infection. To analyze these data, we developed a new stochastic model of the HCV life cycle. We found that the accumulation of mutations is surprisingly slow: at 30 days, the viral population on average is still 46% identical to its transmitted viral genome. Fitting the model to the sequence data, we estimate the median in vivo viral mutation rate is 2.5×10⁻⁵ mutations per nucleotide per genome replication (range 1.6–6.2×10⁻⁵), about 5-fold lower than previous estimates. To confirm these results we analyzed the frequency of stop codons (N = 10) among all possible non-sense mutation targets (M = 898,335), and found a mutation rate of 2.8–3.2×10⁻⁵, consistent with the estimate from the dynamical model. The slow accumulation of mutations is consistent with slow turnover of infected cells and replication complexes within infected cells. This slow turnover is also inferred from the viral load kinetics. Our estimated mutation rate, which is similar to that of other RNA viruses (e.g., HIV and influenza), is also compatible with the accumulation of substitutions seen in HCV at the population level. Our model identifies the relevant processes (long-lived cells and slow turnover of replication complexes) and parameters involved in determining the rate of HCV diversification.     

Inferring the Dynamics of Diversification: A Coalescent Approach

Recent analyses of the fossil record and molecular phylogenies suggest that there are fundamental limits to biodiversity, possibly arising from constraints in the availability of space, resources, or ecological niches. Under this hypothesis, speciation rates decay over time and biodiversity eventually saturates, with new species emerging only when others are driven to extinction. This view of macro-evolution contradicts an alternative hypothesis that biodiversity is unbounded, with species ever accumulating as they find new niches to occupy. These contrasting theories of biodiversity dynamics yield fundamentally different explanations for the disparity in species richness across taxa and regions. Here, we test whether speciation rates have decayed or remained constant over time, and whether biodiversity is saturated or still expanding. We first derive a general likelihood expression for internode distances in a phylogeny, based on the well-known coalescent process from population genetics. This expression accounts for either time-constant or time-variable rates, time-constant or time-variable diversity, and completely or incompletely sampled phylogenies. We then compare the performance of different diversification scenarios in explaining a set of 289 phylogenies representing amphibians, arthropods, birds, mammals, mollusks, and flowering plants. Our results indicate that speciation rates typically decay over time, but that diversity is still expanding at present. The evidence for expanding-diversity models suggests that an upper limit to biodiversity has not yet been reached, or that no such limit exists.     

Morphological Characters Are Compatible with Mitogenomic Data in Resolving the Phylogeny of Nymphalid Butterflies (Lepidoptera: Papilionoidea: Nymphalidae)

Nymphalidae is the largest family of butterflies with their phylogenetic relationships not adequately approached to date. The mitochondrial genomes (mitogenomes) of 11 new nymphalid species were reported and a comparative mitogenomic analysis was conducted together with other 22 available nymphalid mitogenomes. A phylogenetic analysis of the 33 species from all 13 currently recognized nymphalid subfamilies was done based on the mitogenomic data set with three Lycaenidae species as the outgroups. The mitogenome comparison showed that the eleven new mitogenomes were similar with those of other butterflies in gene content and order. The reconstructed phylogenetic trees reveal that the nymphalids are made up of five major clades (the nymphaline, heliconiine, satyrine, danaine and libytheine clades), with sister relationship between subfamilies Cyrestinae and Biblidinae, and most likely between subfamilies Morphinae and Satyrinae. This whole mitogenome-based phylogeny is generally congruent with those of former studies based on nuclear-gene and mitogenomic analyses, but differs considerably from the result of morphological cladistic analysis, such as the basal position of Libytheinae in morpho-phylogeny is not confirmed in molecular studies. However, we found that the mitogenomic phylogeny established herein is compatible with selected morphological characters (including developmental and adult morpho-characters).     

利用线粒体COII基因序列对中国尾蛱蝶属系统分化的研究(鳞翅目:蛱蝶科)

利用线粒体C0II基因序列分析法,研究了我国尾蛱蝶属5种蝴蝶的系统分化.结果表明,在尾蛱蝶属5个种的12个样品中,405 bp长的COII片段有11.4%的位点为多态性位点,大部分的碱基改变是转换.各物种内不同个体间的差异明显小于不同物种间的差异,种内个体间的差异一般为0.5%～1.5%,各物种间的差异绝大多数在4%以上.利用最大似然性法构建的尾蛱蝶属聚类关系图显示,尾蛱蝶属蝴蝶分为两大分支,一支包括大二尾蛱蝶、二尾蛱蝶和忘忧尾蛱蝶,另外一个分支包括窄斑凤尾蛱蝶和黑凤尾蛱蝶聚在一起.在大二尾蛱蝶、二尾蛱蝶和忘忧尾蛱蝶这一分支中,大二尾蛱蝶和忘忧尾蛱蝶的亲缘关系较近,而二尾蛱蝶较远.这些分子系统学的结果均与形态学的结果相一致,是对形态分类的有力支持.     

湖南省万佛山风景区蝴蝶多样性考察初报

2005年6月对湖南省万佛山风景区的蝴蝶进行了调查,记录到蝴蝶10科、41属、84种。通过应用Shannon-Weiner指数、Pielou指数和KA指数,对农田、水域、林地和居民区4种不同生境的调查及多样性分析,表明万佛山风景区是一个良好的蝴蝶生息地,应予以保护。     

Using Lifetime Risk Estimates in Personal Genomic Profiles: Estimation of Uncertainty

Personal genome tests are now offered direct-to-consumer (DTC) via genetic variants identified by genome-wide association studies (GWAS) for common diseases. Tests report risk estimates (age-specific and lifetime) for various diseases based on genotypes at multiple loci. However, uncertainty surrounding such risk estimates has not been systematically investigated. With breast cancer as an example, we examined the combined effect of uncertainties in population incidence rates, genotype frequency, effect sizes, and models of joint effects among genetic variants on lifetime risk estimates. We performed simulations to estimate lifetime breast cancer risk for carriers and noncarriers of genetic variants. We derived population-based cancer incidence rates from Surveillance, Epidemiology, and End Results (SEER) Program and comparative international data. We used data for non-Hispanic white women from 2003 to 2005. We derived genotype frequencies and effect sizes from published GWAS and meta-analyses. For a single genetic variant in FGFR2 gene (rs2981582), combination of uncertainty in these parameters produced risk estimates where upper and lower 95% simulation intervals differed by more than 3-fold. Difference in population incidence rates was the largest contributor to variation in risk estimates. For a panel of five genetic variants, estimated lifetime risk of developing breast cancer before age 80 for a woman that carried all risk variants ranged from 6.1% to 21%, depending on assumptions of additive or multiplicative joint effects and breast cancer incidence rates. Epidemiologic parameters involved in computation of disease risk have substantial uncertainty, and cumulative uncertainty should be properly recognized. Reliance on point estimates alone could be seriously misleading.     

The Evolution of Cultural Diversity: A Phylogenetic Approach

The Evolution of Cultural Diversity:. Phylogenetic Approach . Ruth Mace, Clare J. Holden, and Stephen Shennan, eds. Oxford, U.K.: Cavendish Publishing, 2005. 291 pp.     

贵州赤水桫椤自然保护区蝶类组成及多样性格局

2008年5～10月期间,对贵州省赤水桫椤自然保护区蝶类资源进行了调查.共采集到蝶类465只,结合文献,计113种,分别隶属于11科66属.从自然保护区3条沟的物种组成和刚捕最来看,幺站沟的种数和网捕量最多,其次是板桥沟,金沙沟种数最少.对保护区的幺站沟、金沙沟和板桥沟的蝶类多样性和均匀度分析表明,幺站沟的物种多样性指数(H')最高,板桥沟的均匀度指数(J)最高.3条沟蝶类之间的相似性系数低.     

Phylogenetic Diversity Analysis of Subterranean Hot Springs in Iceland

Geothermal energy has been harnessed and used for domestic heating in Iceland. In wells that are typically drilled to a depth of 1,500 to 2,000 m, the temperature of the source water is 50 to 130°C. The bottoms of the boreholes can therefore be regarded as subterranean hot springs and provide a unique opportunity to study the subterranean biosphere. Large volumes of geothermal fluid from five wells and a mixture of geothermal water from 50 geothermal wells (hot tap water) were sampled and concentrated through a 0.2-μm-pore-size filter. Cells were observed in wells RG-39 (91.4°C) and MG-18 (71.8°C) and in hot tap water (76°C), but no cells were detected in wells SN-4, SN-5 (95 to 117°C), and RV-5 (130°C). Archaea and Bacteria were detected by whole-cell fluorescent in situ hybridization. DNAs were extracted from the biomass, and small-subunit rRNA genes (16S rDNAs) were amplified by PCR using primers specific for the Archaea and Bacteria domains. The PCR products were cloned and sequenced. The sequence analysis showed 11 new operational taxonomic units (OTUs) out of 14, 3 of which were affiliated with known surface OTUs. Samples from RG-39 and hot tap water were inoculated into enrichment media and incubated at 65 and 85°C. Growth was observed only in media based on geothermal water. 16S rDNA analysis showed enrichments dominated with Desulfurococcales relatives. Two strains belonging to Desulfurococcus mobilis and to the Thermus/Deinococcus group were isolated from borehole RG-39. The results indicate that subsurface volcanic zones are an environment that provides a rich subsurface for novel thermophiles.     

Phylogenetic Patterns of Geographical and Ecological Diversification in the Subgenus Drosophila

Colonisation of new geographic regions and/or of new ecological resources can result in rapid species diversification into the new ecological niches available. Members of the subgenus Drosophila are distributed across the globe and show a large diversity of ecological niches. Furthermore, taxonomic classification of Drosophila includes the rank radiation, which refers to closely related species groups. Nevertheless, it has never been tested if these taxonomic radiations correspond to evolutionary radiations. Here we present a study of the patterns of diversification of Drosophila to test for increased diversification rates in relation to the geographic and ecological diversification processes. For this, we have estimated and dated a phylogeny of 218 species belonging to the major species groups of the subgenus. The obtained phylogenies are largely consistent with previous studies and indicate that the major groups appeared during the Oligocene/Miocene transition or early Miocene, characterized by a trend of climate warming with brief periods of glaciation. Ancestral reconstruction of geographic ranges and ecological resource use suggest at least two dispersals to the Neotropics from the ancestral Asiatic tropical disribution, and several transitions to specialized ecological resource use (mycophagous and cactophilic). Colonisation of new geographic regions and/or of new ecological resources can result in rapid species diversification into the new ecological niches available. However, diversification analyses show no significant support for adaptive radiations as a result of geographic dispersal or ecological resource shift. Also, cactophily has not resulted in an increase in the diversification rate of the repleta and related groups. It is thus concluded that the taxonomic radiations do not correspond to adaptive radiations.     

The Climatic Niche Diversity of Malagasy Primates: A Phylogenetic Perspective

Background

Numerous researchers have posited that there should be a strong negative relationship between the evolutionary distance among species and their ecological similarity. Alternative evidence suggests that members of adaptive radiations should display no relationship between divergence time and ecological similarity because rapid evolution results in near-simultaneous speciation early in the clade''s history. In this paper, we performed the first investigation of ecological diversity in a phylogenetic context using a mammalian adaptive radiation, the Malagasy primates.

Methodology/Principal Findings

We collected data for 43 extant species including: 1) 1064 species by locality samples, 2) GIS climate data for each sampling locality, and 3) the phylogenetic relationships of the species. We calculated the niche space of each species by summarizing the climatic variation at localities of known occurrence. Climate data from all species occurrences at all sites were entered into a principal components analysis. We calculated the mean value of the first two PCA axes, representing rainfall and temperature diversity, for each species. We calculated the K statistic using the Physig program for Matlab to examine how well the climatic niche space of species was correlated with phylogeny.

Conclusions/Significance

We found that there was little relationship between the phylogenetic distance of Malagasy primates and their rainfall and temperature niche space, i.e., closely related species tend to occupy different climatic niches. Furthermore, several species from different genera converged on a similar climatic niche. These results have important implications for the evolution of ecological diversity, and the long-term survival of these endangered species.     

Effects of OTU Clustering and PCR Artifacts on Microbial Diversity Estimates

Next-generation sequencing has increased the coverage of microbial diversity surveys by orders of magnitude, but differentiating artifacts from rare environmental sequences remains a challenge. Clustering 16S rRNA sequences into operational taxonomic units (OTUs) organizes sequence data into groups of 97 % identity, helping to reduce data volumes and avoid analyzing sequencing artifacts by grouping them with real sequences. Here, we analyze sequence abundance distributions across environmental samples and show that 16S rRNA sequences of >99 % identity can represent functionally distinct microorganisms, rendering OTU clustering problematic when the goal is an accurate analysis of organism distribution. Strict postsequencing quality control (QC) filters eliminated the most prevalent artifacts without clustering. Further experiments proved that DNA polymerase errors in polymerase chain reaction (PCR) generate a significant number of substitution errors, most of which pass QC filters. Based on our findings, we recommend minimizing the number of PCR cycles in DNA library preparation and applying strict postsequencing QC filters to reduce the most prevalent artifacts while maintaining a high level of accuracy in diversity estimates. We further recommend correlating rare and abundant sequences across environmental samples, rather than clustering into OTUs, to identify remaining sequence artifacts without losing the resolution afforded by high-throughput sequencing.     

湖南乌云界蝴蝶物种多样性及区系特点研究

2008～2010年,采用路线法对湖南省乌云界国家级自然保护区的蝶类物种多样性进行了为期近3年的调查,共记录蝴蝶147种,隶属于10科94属。蛱蝶科Nymphalidae Swainson,1827种类最多,弄蝶科Hesperiidae Latreille,1809的多样性指数最大。凤蝶属Papilio Linnaeus,1758对属多样性指数贡献值相对较大。记录国家保护种4种(宽尾凤蝶Agehana elwesi、中华虎凤蝶Luehdorfia chinensis、箭环蝶Stichophthalma howqua、大伞弄蝶Bibasis miraculata),其中中华虎凤蝶属于国家二级重点保护动物;湖南省新纪录种3种(金斑剑凤蝶Pazala alebion、报喜斑粉蝶Delias pasithoe和泰环蛱蝶Neptis thestias);宽边黄粉蝶Eurema hecabe采集数量最多(69头),是该保护区优势种群。蝶类昆虫区系主要以东洋成分为主,兼有古北和广布成分。     

Dynamics of host plant use and species diversity in Polygonia butterflies (Nymphalidae)

The ability of insects to utilize different host plants has been suggested to be a dynamic and transient phase. During or after this phase, species can shift to novel host plants or respecialize on ancestral ones. Expanding the range of host plants might also be a factor leading to higher levels of net speciation rates. In this paper, we have studied the possible importance of host plant range for diversification in the genus Polygonia (Nymphalidae, Nymphalini). We have compared species richness between sistergroups in order to find out if there are any differences in number of species between clades including species that utilize only the ancestral host plants ('urticalean rosids') and their sisterclades with a broader (or in some cases potentially broader) host plant repertoire. Four comparisons could be made, and although these are not all phylogenetically or statistically independent, all showed clades including butterfly species using other or additional host plants than the urticalean rosids to be more species-rich than their sisterclade restricted to the ancestral host plants. These results are consistent with the theory that expansions in host plant range are involved in the process of diversification in butterflies and other phytophagous insects, in line with the general theory that plasticity may drive speciation.