Incorporation of transition to transversion ratio and nonsense mutations,improves the estimation of the number of synonymous and non-synonymous sites in codons

A common approach to estimate the strength and direction of selection acting on protein coding sequences is to calculate the dN/dS ratio. The method to calculate dN/dS has been widely used by many researchers and many critical reviews have been made on its application after the proposition by Nei and Gojobori in 1986. However, the method is still evolving considering the non-uniform substitution rates and pretermination codons. In our study of SNPs in 586 genes across 156 Escherichia coli strains, synonymous polymorphism in 2-fold degenerate codons were higher in comparison to that in 4-fold degenerate codons, which could be attributed to the difference between transition (Ti) and transversion (Tv) substitution rates where the average rate of a transition is four times more than that of a transversion in general. We considered both the Ti/Tv ratio, and nonsense mutation in pretermination codons, to improve estimates of synonymous (S) and non-synonymous (NS) sites. The accuracy of estimating dN/dS has been improved by considering the Ti/Tv ratio and nonsense substitutions in pretermination codons. We showed that applying the modified approach based on Ti/Tv ratio and pretermination codons results in higher values of dN/dS in 29 common genes of equal reading-frames between E. coli and Salmonella enterica. This study emphasizes the robustness of amino acid composition with varying codon degeneracy, as well as the pretermination codons when calculating dN/dS values.     

Modeling Future Conservation of Hawaiian Honeycreepers by Mosquito Management and Translocation of Disease-Tolerant Amakihi

Avian malaria is an important cause of the decline of endemic Hawaiian honeycreepers. Because of the complexity of this disease system we used a computer model of avian malaria in forest birds to evaluate how two proposed conservation strategies: 1) reduction of habitat for mosquito larvae and 2) establishment of a low-elevation, malaria-tolerant honeycreeper (Hawaii Amakihi) to mid-elevation forests would affect native Hawaiian honeycreeper populations. We evaluated these approaches in mid-elevation forests, where malaria transmission is seasonal and control strategies are more likely to work. Our model suggests the potential benefit of larval habitat reduction depends on the level of malaria transmission, abundance of larval cavities, and the ability to substantially reduce these cavities. Permanent reduction in larval habitat of >80% may be needed to control abundance of infectious mosquitoes and benefit bird populations. Establishment of malaria-tolerant Amakihi in mid-elevation forests increases Amakihi abundance, creates a larger disease reservoir, and increases the abundance of infectious mosquitoes which may negatively impact other honeycreepers. For mid-elevation sites where bird populations are severely affected by avian malaria, malaria-tolerant Amakihi had little impact on other honeycreepers. Both management strategies may benefit native Hawaiian honeycreepers, but benefits depend on specific forest characteristics, the amount of reduction in larval habitat that can be achieved, and how malaria transmission is affected by temperature.     

Clade-specific morphological diversification and adaptive radiation in Hawaiian songbirds.

The Hawaiian honeycreepers are a dramatic example of adaptive radiation but contrast with the four other songbird lineages that successfully colonized the Hawaiian archipelago and failed to undergo similar diversification. To explore the processes that produced the diversity dichotomy in this insular fauna, we compared clade age and morphological diversity between the speciose honeycreepers and the comparatively depauperate Hawaiian thrushes. Mitochondrial-DNA-based genetic distances between these Hawaiian clades and their continental sister taxa indicate that the ancestral thrush colonized the Hawaiian Islands as early as the common ancestor of the honeycreepers. This similar timing of colonization indicates that the marked difference in diversity between the Hawaiian honeycreeper and thrush clades is unlikely to result from differences in these clades' tenures within the archipelago. If time cannot explain the contrasting diversities of these taxa, then an intrinsic, clade-specific trait may have fostered the honeycreeper radiation. As the honeycreepers have diversified most dramatically in morphological characters related to resource utilization, we used principal components analyses of bill characters to compare the magnitudes of morphological variation in the ancestral clades from which the Hawaiian honeycreeper and thrush lineages are derived, the Carduelini and Turdinae respectively. Although the Carduelini share a more recent common ancestor and have a lower species diversity than the Turdinae, these finch-like relatives of the honeycreepers exhibit significantly greater variation in bill morphology than do the continental relatives of the Hawaiian thrushes. The higher magnitude of morphological variation in the non-Hawaiian Carduelini suggests that the honeycreepers fall within a clade exhibiting a generally high evolutionary flexibility in bill morphology. Accordingly, although the magnitude of bill variation among the honeycreepers is similar to that of the entire passerine radiation, this dramatic morphological radiation represents an extreme manifestation of a general clade-specific ability to evolve novel morphologies.     

A Morphological Test of the Monophyly of the Cardueline Finches (Aves: Fringillidae, Carduelinae)

Monophyly of the cardueline finches (Aves: Fringillidae, Carduelinae), an assemblage of 119 species of seed-eating songbirds, has not yet been rigorously tested. To test the hypothesis of cardueline monophyly I examined 37 taxa, among them 1 or 2 representatives from 16 of the 19 cardueline-finch genera as well as 4 Hawaiian honeycreepers. Each taxon was scored for 225 morphological characters, 148 from the skeleton and 77 from the integument. Cladistic analysis of these data found the carduelines to be monophyletic, so long as the honeycreeper Telespiza cantans is considered to be cardueline; as for the other honeycreepers, they form a clade but do not group with either Telespiza or the carduelines. Both cardueline monophyly and the nonmonophyly of honeycreepers have comparatively strong support, with support for cardueline monophyly coming primarily from evolutionary modifications to the head skeleton.     

Two mitochondrial genes under episodic positive selection in subterranean octodontoid rodents

Tuco-tucos (Ctenomys) and related coruros (Spalacopus) are South American subterranean rodents. An energetically demanding lifestyle within the hypoxic, underground atmosphere may change the selective regime on oxidative phosphorylation. We examined whether weak and/or episodic positive directional selection affected the evolution of two mitochondrial genes (COX2, CytB), in a background of purifying selection in these lineages. We estimated rates of synonymous (dS) and non-synonymous (dN) substitutions and found: 1) significantly higher dN/dS ratio in subterranean groups relative to non-subterranean related species, and 2) two codons in each gene under episodic selection: 94 and 277 of COX2 and 269 and 307 of CytB.     

Patterns of variation of the major histocompatibility complex class IIB loci in Chinese goose (<Emphasis Type="Italic">Anser Cygnoides</Emphasis>)

In order to understand the variations of genomic organization of the major histocompatibility complex (MHC) and provide data for the studies on disease resistance of avian species, the MHC class II polymorphism in Chinese Z-goose was investigated for the first time in the present study. Eight alleles, which probably came from different loci, were found in six different geese with only one obvious band in the restriction fragment length polymorphism data. The numbers of nonsynonymous substitutions (dN) in peptide binding region of exon 2 were higher than that of synonymous substitutions (dS), and no stop codons or frameshift mutations were found in this region, indicating that balance selection was in operation, and the sequences are not likely to be pseudogenes. In addition, we successfully obtained five different long MHC class II fragments (about 1,162 bp) in six geese and found that the length of intron 1 was longer than that in chicken and some other birds, but intron 2 seemed to be intermediate in length. The phylogenetic tree appeared to branch in an order consistent with accepted evolutionary pathway. X. Zhou and C. Li contributed equally to this work and shared the co-first authorship.     

Patterns of Spatial use and movement of the Poóuli – a critically endangered Hawaiian honeycreeper

The Poóuli (Melamprosops phaeosoma) is a highly endangered Hawaiian honeycreeper endemic to Maui, and is currently one of the World's rarest birds. Only two wild individuals of this species are now known to exist, and they are restricted to the windward slopes of Haleakala volcano on east Maui. Studies of the more common honeycreepers on the Hawaiian islands describe a diverse array of spatial use and movement patterns, which vary according to specific ecological needs. In contrast, spatial use and movement of Poóuli are very poorly known, despite continuous field monitoring of all three individuals since 1997. We analyzed annual data by breeding season between 1995 and 2001 for three individuals, and eight days of telemetry data derived from the radio-tracking of one individual in 2002, using GIS and conventional methods, to estimate home range size, and to interpret these data alongside those for other honeycreepers. We estimated mean home range sizes of 8.43 hectares (ha) from annual re-sights using kernel analysis and 2.14 ha using minimum convex polygons, and 8.44 and 3.51 ha respectively from telemetry data. Our estimates conform to those derived for other insectivorous honeycreepers, but indicate that Poóuli may forage widely to support their diet of forest snails. Our home range size estimates are compatible with estimates of population density for Poóuli that were derived from field transects between 1975 and 1985, suggesting that such field methods may be a reliable density estimator for rare forest birds. Jim J. Groombridge and Bill Sparklin: Joint first author     

A PCR test for avian malaria in Hawaiian birds

The decline of native Hawaiian forest birds since European contact is attributed to factors ranging from habitat destruction to interactions with introduced species. Remaining populations of Hawaiian honeycreepers (Fringillidae: Drepanidinae) are most abundant and diverse in high elevation refuges above the normal range of disease-carrying mosquitoes. Challenge experiments suggest that honeycreepers are highly susceptible to avian malaria (Plasmodium sp.) but resistance exists in some species. In order to detect low levels of malarial infection and quantify prevalence of Plasmodium in high elevation natural populations of Hawaiian birds, a polymerase chain reaction (PCR) based diagnostic test was developed that identifies rRNA genes of Plasmodium in avian blood samples. Quantitative competitive PCR (QC-PCR) experiments indicate that the detection limit of our test is an order of magnitude greater than that reported for human malaria DNA blot tests. Compared with standard histological methods, the PCR test detected a higher prevalence of diseased birds at mid-elevations. Malaria was detected in three species of native birds living in a high elevation wildlife refuge on the island of Hawaii and in four species from Maui. Our results show that avian malaria is more widespread in Hawaiian forests than previously thought, a finding that has important conservation implications for these threatened species.     

Phylogenetic relationships of Neotropical honeycreepers and the evolution of feeding morphology

Phylogenetic relationships among Neotropical honeycreepers were studied by using cytochrome b sequence data. Although honeycreepers were once placed in their own family (Coerebidae), these species did not form a monophyletic group in any of the phylogenies in this study. Thus, our results indicate that nectar-feeding evolved independently multiple times among these birds. The Neotropical honeycreepers are best considered disparate members of a larger radiation of tanagers and finches, in which bill size and shape, and associated feeding behaviors have changed frequently to fill a variety of niches. Our phylogenies also provide specific taxonomic recommendations on the placement of each honeycreeper genus. The use of the taxon Coerebidae is no longer warranted given the well-supported relationships between different honeycreeper species and a variety of tanagers and finches.     

Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution

An excess of nonsynonymous substitutions over synonymous ones is an important indicator of positive selection at the molecular level. A lineage that underwent Darwinian selection may have a nonsynonymous/synonymous rate ratio (dN/dS) that is different from those of other lineages or greater than one. In this paper, several codon-based likelihood models that allow for variable dN/dS ratios among lineages were developed. They were then used to construct likelihood ratio tests to examine whether the dN/dS ratio is variable among evolutionary lineages, whether the ratio for a few lineages of interest is different from the background ratio for other lineages in the phylogeny, and whether the dN/dS ratio for the lineages of interest is greater than one. The tests were applied to the lysozyme genes of 24 primate species. The dN/dS ratios were found to differ significantly among lineages, indicating that the evolution of primate lysozymes is episodic, which is incompatible with the neutral theory. Maximum- likelihood estimates of parameters suggested that about nine nonsynonymous and zero synonymous nucleotide substitutions occurred in the lineage leading to hominoids, and the dN/dS ratio for that lineage is significantly greater than one. The corresponding estimates for the lineage ancestral to colobine monkeys were nine and one, and the dN/dS ratio for the lineage is not significantly greater than one, although it is significantly higher than the background ratio. The likelihood analysis thus confirmed most, but not all, conclusions Messier and Stewart reached using reconstructed ancestral sequences to estimate synonymous and nonsynonymous rates for different lineages.      

Isolation of polymorphic microsatellite loci in the Hawaii amakihi (Hemignathus virens) and their use in other honeycreeper species

We report the isolation and characterization of five polymorphic microsatellite loci in the Hawaii amakihi, Hemignathus virens, one of the most common native Hawaiian forest birds. These loci exhibit high levels of allelic diversity and heterozygosity in the three honeycreeper species we screened, and promise to be useful in our investigation of differences between honeycreeper populations at different elevations.     

Sequence variability at three MHC loci of finless porpoises (<Emphasis Type="Italic">Neophocaena phocaenoides</Emphasis>)

Major histocompatibility complex (MHC) class II DQB and DRA genes and class I gene of finless porpoises (Neophocaena phocaenoides) were investigated by single-strand conformation polymorphism and sequence analysis. The DRA, DQB, and MHC-I loci each contained 5, 14, and 34 unique sequences, respectively, and considerable sequence variation was found at the MHC-I and DQB loci. Gene duplication was manifested as three to five distinct sequences at each of the DQB and MHC-I loci from some individuals, and these sequences at each of the two loci separately clustered into four groups (cluster A, B, C, and D) based on the phylogenetic trees. Phylogenetic reconstruction revealed a trans-species pattern of evolution. Relatively high rates of non-synonymous (dN) vs synonymous (dS) substitution in the peptide-binding region (PBR) suggested balancing selection for maintaining polymorphisms at the MHC-I and DQB loci. In contrast, one single locus with little sequence variation was detected in the DRA gene, and no non-synonymous substitutions in the PBR indicated no balancing selection on this gene.     

Energetic cost of feeding territories in an Hawaiian honeycreeper

Summary By analysis of time budgets the daily energy expenditure in territorial individuals of a Hawaiian honeycreeper (Vestiaria coccinea, Fam. Drepanididae) were estimated during the nonbreeding season and compared to that of nonterritorial individuals. The mean rise in living costs was 2.3 kcal/24 h or 17% of the nonterritorial energy budget. The most costly territorial behavior was advertisement rather than chasing, and total territorial cost was seen to be little affected by the number of intruders or the size of the territory. These results are compared with data on feeding (nonbreeding) territories of other nectar-feeding birds. The suggestion is made that hummingbirds may be more likely to develop nonbreeding territorial behavior in any set of environmental circumstances than are honeycreepers because of relatively lower total cost of advertisement plus chasing.     

Molecular Evolution of the Genomic RNA of Apple Stem Grooving Capillovirus

The complete genome of the German isolate AC of Apple stem grooving virus (ASGV) was sequenced. It encodes two overlapping open reading frames (ORFs), similarly to previously described ASGV isolates. Two regions of high variability were detected between the ASGV isolates, variable region 1 (V1, from amino acids (aa) 532 to 570), and variable region 2 (V2, from aa 1,583 to 1,868). The phylogenetic analysis of the V1 and V2 regions suggested that the ASGV diversity was structured by host plant species rather than geographical origin. The dN/dS ratio between nonsynonymous and synonymous nucleotide substitution rates varied greatly along the ASGV genome. Most of ORF1 showed predominant negative selection except for the two regions V1 and V2. V1 showed an elevated dN and an average dS when compared to the ORF1 background but no significant positive selection was detected. The V2 region of ORF1 showed an elevated dN and a low dS when compared to the ORF1 background with an average dN/dS????3.0 indicative of positive selection. However, the V2 area includes overlapping ORFs, making the dN/dS estimate biased. Joint estimates of the selection intensity in the different ORFs by a recent method indicated that this region of ORF1 was in fact evolving close to neutrality. This was convergent with previous results showing that introduction of stop codons in this region of ORF1 did not impair plant infection. These data suggest that the elimination of a stop codon caused the overprinting of a novel coding region over the ancestral ORF.     

A Sliding Window-Based Method to Detect Selective Constraints in Protein-Coding Genes and Its Application to RNA Viruses

Here we present a new sliding window-based method specially designed to detect selective constraints in specific regions of a multiple protein-coding sequence alignment. In contrast to previous window-based procedures, our method is based on a nonarbitrary statistical approach to find the appropriate codon-window size to test deviations of synonymous (dS) and nonsynonymous (dN) nucleotide substitutions from the expectation. The probabilities of dN and dS are obtained from simulated data and used to detect significant deviations of dN and dS in a specific window region of the real sequence alignment. The nonsynonymous-to-synonymous rate ratio (w = dN/dS) was used to highlight selective constraints in any window wherein dS or dN was significantly different from the expectation. In these significant windows, w and its variance [V(w)] were calculated and used to test the neutral hypothesis. Computer simulations showed that the method is accurate even for highly divergent sequences. The main advantages of the new method are that it (i) uses a statistically appropriate window size to detect different selective patterns, (ii) is computationally less intensive than maximum likelihood methods, and (iii) detects saturation of synonymous sites, which can give deviations from neutrality. Hence, it allows the analysis of highly divergent sequences and the test of different alternative hypothesis as well. The application of the method to different human immunodeficiency virus type 1 and to foot-and-mouth disease virus genes confirms the action of positive selection on previously described regions as well as on new regions.     

Unmelanized plumage patterns in Old World leaf warblers do not correspond to sequence variation at the melanocortin-1 receptor locus (MC1R)

Evolutionary changes in patterns and coloration of plumage are likely to represent a major mechanism for speciation among birds, yet the molecular basis for such changes remains poorly understood. Recently much attention has focused on the melanocortin-1 receptor (MC1R) as a candidate locus for determining the level and extent of epidermal melanin deposition. We tested the hypothesis that MC1R sequence variation is associated with interspecific variation in unmelanized plumage pattern elements in Old World leaf warblers (genus Phylloscopus). This genus is characterized by a variety of plumage patterns that nonetheless vary along similar lines. Species vary in the presence or absence of pale (unmelanized) pattern elements against a dark background, and these patterns are used in species recognition and courtship. We sequenced most of the MC1R coding region for eight Phylloscopus species, representing the full range of plumage patterns found in this genus. Although MC1R sequence varied among species, this variation was not related to melanin-based plumage variation. Rather, evolution of this locus in these birds appears to be conservative. Ratios of nonsynonymous to synonymous substitutions (dN/dS) were consistently low, suggesting that strong purifying selection has operated at this locus, and likelihood ratio testing revealed no evidence of variable selective pressures among lineages or across codons. Adaptive evolution at MC1R may be constrained by the adaptive importance of plumage pattern elements in this genus.     

Multilocus resolution of phylogeny and timescale in the extant adaptive radiation of Hawaiian honeycreepers

Evolutionary theory has gained tremendous insight from studies of adaptive radiations. High rates of speciation, morphological divergence, and hybridization, combined with low sequence variability, however, have prevented phylogenetic reconstruction for many radiations. The Hawaiian honeycreepers are an exceptional adaptive radiation, with high phenotypic diversity and speciation that occurred within the geologically constrained setting of the Hawaiian Islands. Here we analyze a new data set of 13 nuclear loci and pyrosequencing of mitochondrial genomes that resolves the Hawaiian honeycreeper phylogeny. We show that they are a sister taxon to Eurasian rosefinches (Carpodacus) and probably came to Hawaii from Asia. We use island ages to calibrate DNA substitution rates, which vary substantially among gene regions, and calculate divergence times, showing that the radiation began roughly when the oldest of the current large Hawaiian Islands (Kauai and Niihau) formed, ~5.7 million years ago (mya). We show that most of the lineages that gave rise to distinctive morphologies diverged after Oahu emerged (4.0-3.7 mya) but before the formation of Maui and adjacent islands (2.4-1.9 mya). Thus, the formation of Oahu, and subsequent cycles of colonization and speciation between Kauai and Oahu, played key roles in generating the morphological diversity of the extant honeycreepers.     

Differentiation of Xylella fastidiosa strains via multilocus sequence analysis of environmentally mediated genes (MLSA-E)

Isolates of the plant pathogen Xylella fastidiosa are genetically very similar, but studies on their biological traits have indicated differences in virulence and infection symptomatology. Taxonomic analyses have identified several subspecies, and phylogenetic analyses of housekeeping genes have shown broad host-based genetic differences; however, results are still inconclusive for genetic differentiation of isolates within subspecies. This study employs multilocus sequence analysis of environmentally mediated genes (MLSA-E; genes influenced by environmental factors) to investigate X. fastidiosa relationships and differentiate isolates with low genetic variability. Potential environmentally mediated genes, including host colonization and survival genes related to infection establishment, were identified a priori. The ratio of the rate of nonsynonymous substitutions to the rate of synonymous substitutions (dN/dS) was calculated to select genes that may be under increased positive selection compared to previously studied housekeeping genes. Nine genes were sequenced from 54 X. fastidiosa isolates infecting different host plants across the United States. Results of maximum likelihood (ML) and Bayesian phylogenetic (BP) analyses are in agreement with known X. fastidiosa subspecies clades but show novel within-subspecies differentiation, including geographic differentiation, and provide additional information regarding host-based isolate variation and specificity. dN/dS ratios of environmentally mediated genes, though <1 due to high sequence similarity, are significantly greater than housekeeping gene dN/dS ratios and correlate with increased sequence variability. MLSA-E can more precisely resolve relationships between closely related bacterial strains with low genetic variability, such as X. fastidiosa isolates. Discovering the genetic relationships between X. fastidiosa isolates will provide new insights into the epidemiology of populations of X. fastidiosa, allowing improved disease management in economically important crops.     

Comparative genomic analysis of two avian (quail and chicken) MHC regions

We mapped two different quail Mhc haplotypes and sequenced one of them (haplotype A) for comparative genomic analysis with a previously sequenced haplotype of the chicken Mhc. The quail haplotype A spans 180 kb of genomic sequence, encoding a total of 41 genes compared with only 19 genes within the 92-kb chicken Mhc. Except for two gene families (B30 and tRNA), both species have the same basic set of gene family members that were previously described in the chicken "minimal essential" Mhc. The two Mhc regions have a similar overall organization but differ markedly in that the quail has an expanded number of duplicated genes with 7 class I, 10 class IIB, 4 NK, 6 lectin, and 8 B-G genes. Comparisons between the quail and chicken Mhc class I and class II gene sequences by phylogenetic analysis showed that they were more closely related within species than between species, suggesting that the quail Mhc genes were duplicated after the separation of these two species from their common ancestor. The proteins encoded by the NK and class I genes are known to interact as ligands and receptors, but unlike in the quail and the chicken, the genes encoding these proteins in mammals are found on different chromosomes. The finding of NK-like genes in the quail Mhc strongly suggests an evolutionary connection between the NK C-type lectin-like superfamily and the Mhc, providing support for future studies on the NK, lectin, class I, and class II interaction in birds.