Search

The traits involved in sexual selection, such as male secondary sexual characteristics and female mate choice, often co-evolve which can promote population differentiation. However, the genetic architecture of these phenotypes can influence their evolvability and thereby affect the divergence of species. The extraordinary diversity of East African cichlid fishes is often attributed to strong sexual selection and thus this system provides an excellent model to test predictions regarding the genetic architecture of sexually selected traits that contribute to reproductive isolation. In particular, theory predicts that rapid speciation is facilitated when male sexual traits and female mating preferences are controlled by a limited number of linked genes. However, few studies have examined the genetic basis of male secondary sexual traits and female mating preferences in cichlids and none have investigated the genetic architecture of both jointly. In this study, we artificially hybridized a pair of behaviorally isolated cichlid fishes from Lake Malawi and quantified both melanistic color pattern and female mate choice. We investigated the genetic architecture of both phenotypes using quantitative genetic analyses. Our results suggest that 1) many non-additively acting genetic factors influence melanistic color patterns, 2) female mate choice may be controlled by a minimum of 1–2 non-additive genetic factors, and 3) F2 female mate choice is not influenced by male courting effort. Furthermore, a joint analysis of color pattern and female mate choice indicates that the genes underlying these two traits are unlikely to be physically linked. These results suggest that reproductive isolation may evolve rapidly owing to the few genetic factors underlying female mate choice. Hence, female mate choice likely played an important role in the unparalleled speciation of East African cichlid fish. [ABSTRACT FROM AUTHOR]

Genetic drift and selection are ubiquitous evolutionary forces acting to shape genetic variation in populations. While their relative importance has been well studied in plants and animals, less is known about their relative importance in fungal pathogens. Because agro-ecosystems are more homogeneous environments than natural ecosystems, stabilizing selection may play a stronger role than genetic drift or diversifying selection in shaping genetic variation among populations of fungal pathogens in agro-ecosystems. We tested this hypothesis by conducting a QST/FST analysis using agricultural populations of the barley pathogen Rhynchosporium commune. Population divergence for eight quantitative traits (QST) was compared with divergence at eight neutral microsatellite loci (FST) for 126 pathogen strains originating from nine globally distributed field populations to infer the effects of genetic drift and types of selection acting on each trait. Our analyses indicated that five of the eight traits had QST values significantly lower than FST, consistent with stabilizing selection, whereas one trait, growth under heat stress (22°C), showed evidence of diversifying selection and local adaptation (QST>FST). Estimates of heritability were high for all traits (means ranging between 0.55–0.84), and average heritability across traits was negatively correlated with microsatellite gene diversity. Some trait pairs were genetically correlated and there was significant evidence for a trade-off between spore size and spore number, and between melanization and growth under benign temperature. Our findings indicate that many ecologically and agriculturally important traits are under stabilizing selection in R. commune and that high within-population genetic variation is maintained for these traits. [ABSTRACT FROM AUTHOR]

Cystatins are a family of inhibitors of cysteine peptidases that comprises the salivary cystatins (D and S-type cystatins) and cystatin C. These cystatins are encoded by a multigene family (CST3, CST5, CST4, CST1 and CST2) organized in tandem in the human genome. Their presence and functional importance in human saliva has been reported, however the distribution of these proteins in other mammals is still unclear. Here, we performed a proteomic analysis of the saliva of several mammals and studied the evolution of this multigene family. The proteomic analysis detected S-type cystatins (S, SA, and SN) in human saliva and cystatin D in rat saliva. The evolutionary analysis showed that the cystatin C encoding gene is present in species of the most representative mammalian groups, i.e. Artiodactyla, Rodentia, Lagomorpha, Carnivora and Primates. On the other hand, D and S-type cystatins are mainly retrieved from Primates, and especially the evolution of S-type cystatins seems to be a dynamic process as seen in Pongo abelii genome where several copies of CST1-like gene (cystatin SN) were found. In Rodents, a group of cystatins previously identified as D and S has also evolved. Despite the high divergence of the amino acid sequence, their position in the phylogenetic tree and their genome organization suggests a common origin with those of the Primates. These results suggest that the D and S type cystatins have emerged before the mammalian radiation and were retained only in Primates and Rodents. Although the mechanisms driving the evolution of cystatins are unknown, it seems to be a dynamic process with several gene duplications evolving according to the birth-and-death model of evolution. The factors that led to the appearance of a group of saliva-specific cystatins in Primates and its rapid evolution remain undetermined, but may be associated with an adaptive advantage. [ABSTRACT FROM AUTHOR]

A new horseshoe crab species, Limulus darwini, is described from the uppermost Jurassic (ca. 148 Ma) near-shore sediments of the Kcynia Formation, central Poland. The only extant species Limulus polyphemus (Linnaeus) inhabits brackish-marine, shallow water environments of the east coast of the United States. Here it is shown that there are no important morphological differences between the Kcynia Formation specimens and extant juvenile representatives of the genus Limulus. The palaeoecological setting inhabited by the new species and the trophic relationships of extant horseshoe crabs are discussed in an attempt to determine the potential range of food items ingested by these Mesozoic xiphosurans. In this paper we propose the adoption of a new term stabilomorphism, this being: an effect of a specific formula of adaptative strategy among organisms whose taxonomic status does not exceed genus-level. A high effectiveness of adaptation significantly reduces the need for differentiated phenotypic variants in response to environmental changes and provides for long-term evolutionary success. [ABSTRACT FROM AUTHOR]

To understand how landscape characteristics affect gene flow in species with diverging ecological traits, it is important to analyze taxonomically related sympatric species in the same landscape using identical methods. Here, we present such a comparative landscape genetic study involving three closely related Hesperid butterflies of the genus Thymelicus that represent a gradient of diverging ecological traits. We analyzed landscape effects on their gene flow by deriving inter-population connectivity estimates based on different species distribution models (SDMs), which were calculated from multiple landscape parameters. We then used SDM output maps to calculate circuit-theoretic connectivity estimates and statistically compared these estimates to actual genetic differentiation in each species. We based our inferences on two different analytical methods and two metrics of genetic differentiation. Results indicate that land use patterns influence population connectivity in the least mobile specialist T. acteon. In contrast, populations of the highly mobile generalist T. lineola were panmictic, lacking any landscape related effect on genetic differentiation. In the species with ecological traits in between those of the congeners, T. sylvestris, climate has a strong impact on inter-population connectivity. However, the relative importance of different landscape factors for connectivity varies when using different metrics of genetic differentiation in this species. Our results show that closely related species representing a gradient of ecological traits also show genetic structures and landscape genetic relationships that gradually change from a geographical macro- to micro-scale. Thus, the type and magnitude of landscape effects on gene flow can differ strongly even among closely related species inhabiting the same landscape, and depend on their relative degree of specialization. In addition, the use of different genetic differentiation metrics makes it possible to detect recent changes in the relative importance of landscape factors affecting gene flow, which likely change as a result of contemporary habitat alterations. [ABSTRACT FROM AUTHOR]

The evolution of the flower is commonly thought to have spurred angiosperm diversification. Similarly, particular floral traits might have promoted diversification within specific angiosperm clades. We hypothesize that traits promoting the precise positional transfer of pollen between flowers might promote diversification. In particular, precise pollen transfer might produce partial reproductive isolation that facilitates adaptive divergence between parapatric populations differing in their reproductive-organ positions. We investigate this hypothesis with an individual-based model of pollen transfer dynamics associated with heterostyly, a floral syndrome that depends on precise pollen transfer. Our model shows that precise pollen transfer can cause sexual selection leading to divergence in reproductive-organ positions between populations served by different pollinators, pleiotropically causing an increase in reproductive isolation through a “magic trait” mechanism. Furthermore, this increased reproductive isolation facilitates adaptive divergence between the populations in an unlinked, ecologically selected trait. In a different pollination scenario, however, precise pollen transfer causes a decrease in adaptive divergence by promoting asymmetric gene flow. Our results highlight the idea that magic traits are not “magic” in isolation; in particular, the effect size of magic traits in speciation depends on the external environment, and also on other traits that modify the strength of the magic trait's influence on non-random mating. Overall, we show that the evolutionary consequences of pollen transfer dynamics can depend strongly on the available pollinator fauna and on the morphological fit between flowers and pollinators. Furthermore, our results illustrate the potential importance of even weak reproductive isolating barriers in facilitating adaptive divergence. [ABSTRACT FROM AUTHOR]

Stimulated by retinoic acid 6 (STRA6) is the receptor for retinol binding protein and is relevant for the transport of retinol to specific sites such as the eye. The adaptive evolution mechanism that vertebrates have occupied nearly every habitat available on earth and adopted various lifestyles associated with different light conditions and visual challenges, as well as their role in development and adaptation is thus far unknown. In this work, we have investigated different aspects of vertebrate STRA6 evolution and used molecular evolutionary analyses to detect evidence of vertebrate adaptation to the lightless habitat. Free-ratio model revealed significant rate shifts immediately after the species divergence. The amino acid sites detected to be under positive selection are within the extracellular loops of STRA6 protein. Branch-site model A test revealed that STRA6 has undergone positive selection in the different phyla of mammalian except for the branch of rodent. The results suggest that interactions between different light environments and host may be driving adaptive change in STRA6 by competition between species. In support of this, we found that altered functional constraints may take place at some amino acid residues after speciation. We suggest that STRA6 has undergone adaptive evolution in different branch of vertebrate relation to habitat environment. [ABSTRACT FROM AUTHOR]

Availability of genome-wide gene expression datasets provides the opportunity to study gene expression across different organisms under a plethora of experimental conditions. In our previous work, we developed an algorithm called COMODO (COnserved MODules across Organisms) that identifies conserved expression modules between two species. In the present study, we expanded COMODO to detect the co-expression conservation across three organisms by adapting the statistics behind it. We applied COMODO to study expression conservation/divergence between Escherichia coli, Salmonella enterica, and Bacillus subtilis. We observed that some parts of the regulatory interaction networks were conserved between E. coli and S. enterica especially in the regulon of local regulators. However, such conservation was not observed between the regulatory interaction networks of B. subtilis and the two other species. We found co-expression conservation on a number of genes involved in quorum sensing, but almost no conservation for genes involved in pathogenicity across E. coli and S. enterica which could partially explain their different lifestyles. We concluded that despite their different lifestyles, no significant rewiring have occurred at the level of local regulons involved for instance, and notable conservation can be detected in signaling pathways and stress sensing in the phylogenetically close species S. enterica and E. coli. Moreover, conservation of local regulons seems to depend on the evolutionary time of divergence across species disappearing at larger distances as shown by the comparison with B. subtilis. Global regulons follow a different trend and show major rewiring even at the limited evolutionary distance that separates E. coli and S. enterica. [ABSTRACT FROM AUTHOR]

The Mexican Sheartail (Doricha eliza), an endangered hummingbird, is endemic to Mexico where two populations have a disjunct distribution. One population is distributed along the northern tip of the Yucatan Peninsula whereas the other is mostly restricted to central Veracruz. Despite their disjunct distribution, previous work has failed to detect morphological or behavioral differences between these populations. Here we use variation in morphology, mtDNA and nuDNA sequences to determine the degree of morphological and molecular divergence between populations, their divergence time, and historical demography. We use species distribution modeling and niche divergence tests to infer the relative roles of vicariance and dispersal in driving divergence in the genus. Our Bayesian and maximum likelihood phylogenetic analyses revealed that Doricha eliza populations form a monophyletic clade and support their sister relationship with D. enicura. We found marked genetic differentiation, with reciprocal monophyly of haplotypes and highly restricted gene flow, supporting a history of isolation over the last 120,000 years. Genetic divergence between populations is consistent with the lack of overlap in environmental space and slight morphological differences between males. Our findings indicate that the divergence of the Veracruz and Yucatan populations is best explained by a combination of a short period of isolation exacerbated by subsequent divergence in climate conditions, and that rather than vicariance, the two isolated ranges of D. eliza are the product of recent colonization and divergence in isolation. [ABSTRACT FROM AUTHOR]

The Mexican Sheartail (Doricha eliza), an endangered hummingbird, is endemic to Mexico where two populations have a disjunct distribution. One population is distributed along the northern tip of the Yucatan Peninsula whereas the other is mostly restricted to central Veracruz. Despite their disjunct distribution, previous work has failed to detect morphological or behavioral differences between these populations. Here we use variation in morphology, mtDNA and nuDNA sequences to determine the degree of morphological and molecular divergence between populations, their divergence time, and historical demography. We use species distribution modeling and niche divergence tests to infer the relative roles of vicariance and dispersal in driving divergence in the genus. Our Bayesian and maximum likelihood phylogenetic analyses revealed that Doricha eliza populations form a monophyletic clade and support their sister relationship with D. enicura. We found marked genetic differentiation, with reciprocal monophyly of haplotypes and highly restricted gene flow, supporting a history of isolation over the last 120,000 years. Genetic divergence between populations is consistent with the lack of overlap in environmental space and slight morphological differences between males. Our findings indicate that the divergence of the Veracruz and Yucatan populations is best explained by a combination of a short period of isolation exacerbated by subsequent divergence in climate conditions, and that rather than vicariance, the two isolated ranges of D. eliza are the product of recent colonization and divergence in isolation. [ABSTRACT FROM AUTHOR]

Plants have two principal defense mechanisms to decrease fitness losses to herbivory: tolerance, the ability to compensate fitness after damage, and resistance, the ability to avoid damage. Variation in intensity of herbivory among populations should result in variation in plant defense levels if tolerance and resistance are associated with costs. Yet little is known about how levels of tolerance are related to resistance and attack intensity in the field, and about the costs of tolerance. In this study, we used information about tolerance and resistance against larval herbivory by the butterfly Anthocharis cardamines under controlled conditions together with information about damage in the field for a large set of populations of the perennial plant Cardamine pratensis. Plant tolerance was estimated in a common garden experiment where plants were subjected to a combination of larval herbivory and clipping. We found no evidence of that the proportion of damage that was caused by larval feeding vs. clipping influenced plant responses. Damage treatments had a negative effect on the three measured fitness components and also resulted in an earlier flowering in the year after the attack. Tolerance was related to attack intensity in the population of origin, i.e. plants from populations with higher attack intensity were more likely to flower in the year following damage. However, we found no evidence of a relationship between tolerance and resistance. These results indicate that herbivory drives the evolution for increased tolerance, and that changes in tolerance are not linked to changes in resistance. We suggest that the simultaneous study of tolerance, attack intensity in the field and resistance constitutes a powerful tool to understand how plant strategies to avoid negative effects of herbivore damage evolve. [ABSTRACT FROM AUTHOR]

Accurate estimates of species richness are essential to macroecological and macroevolutionary research, as well as to the effective management and conservation of biodiversity. The resolution of taxonomic relationships is therefore of vital importance. While molecular methods have revolutionised taxonomy, contemporary species delimitation requires an integrative, multi-disciplinary approach. Despite boasting a remarkably high level of endemism, the avifauna of the Sulawesi region of Indonesia remains poorly studied. Previous studies of avian diversity in Sulawesi have focussed predominantly on phenotypic characteristics, thus potentially overlooking any genetically distinct lineages. Grey-sided Flowerpecker Dicaeum celebicum populations from the Wakatobi archipelago were originally described as a separate species from those on nearby mainland Sulawesi. However, for reasons that remain unknown, the Wakatobi populations were reclassified as a subspecies of the mainland form. Combining estimates of genetic divergence with phylogenetic and morphological analyses, we reassessed the status of Wakatobi populations. Our results describe the Wakatobi populations as a separate species to those on mainland Sulawesi; reproductively isolated, genetically and morphologically distinct. We therefore recommend the reclassification of these populations to their original status of Dicaeum kuehni and propose the vernacular name ‘Wakatobi Flowerpecker’. In consideration of our findings and the lack of integrative ornithological research within the Sulawesi region, we believe species richness and avian endemism within the region are underestimated. [ABSTRACT FROM AUTHOR]

Distinguishing between hybrid introgression and incomplete lineage sorting causing incongruence among gene trees in that they exhibit topological differences requires application of statistical approaches that are based on biologically relevant models. Such study is especially challenging in hybrid systems, where usual vectors mediating interspecific gene transfers - hybrids with Mendelian heredity - are absent or unknown. Here we study a complex of hybridizing species, which are known to produce clonal hybrids, to discover how one of the species, Cobitis tanaitica, has achieved a pattern of mito-nuclear mosaic genome over the whole geographic range. We appplied three distinct methods, including the method using solely the information on gene tree topologies, and found that the contrasting mito-nuclear signal might not have resulted from the retention of ancestral polymorphism. Instead, we found two signs of hybridization events related to C. tanaitica; one concerning nuclear gene flow and the other suggested mitochondrial capture. Interestingly, clonal inheritance (gynogenesis) of contemporary hybrids prevents genomic introgressions and non-clonal hybrids are either absent or too rare to be detected among European Cobitis. Our analyses therefore suggest that introgressive hybridizations are rather old episodes, mediated by previously existing hybrids whose inheritance was not entirely clonal. Cobitis complex thus supports the view that the type of resulting hybrids depends on a level of genomic divergence between sexual species. [ABSTRACT FROM AUTHOR]

Plants have two principal defense mechanisms to decrease fitness losses to herbivory: tolerance, the ability to compensate fitness after damage, and resistance, the ability to avoid damage. Variation in intensity of herbivory among populations should result in variation in plant defense levels if tolerance and resistance are associated with costs. Yet little is known about how levels of tolerance are related to resistance and attack intensity in the field, and about the costs of tolerance. In this study, we used information about tolerance and resistance against larval herbivory by the butterfly Anthocharis cardamines under controlled conditions together with information about damage in the field for a large set of populations of the perennial plant Cardamine pratensis. Plant tolerance was estimated in a common garden experiment where plants were subjected to a combination of larval herbivory and clipping. We found no evidence of that the proportion of damage that was caused by larval feeding vs. clipping influenced plant responses. Damage treatments had a negative effect on the three measured fitness components and also resulted in an earlier flowering in the year after the attack. Tolerance was related to attack intensity in the population of origin, i.e. plants from populations with higher attack intensity were more likely to flower in the year following damage. However, we found no evidence of a relationship between tolerance and resistance. These results indicate that herbivory drives the evolution for increased tolerance, and that changes in tolerance are not linked to changes in resistance. We suggest that the simultaneous study of tolerance, attack intensity in the field and resistance constitutes a powerful tool to understand how plant strategies to avoid negative effects of herbivore damage evolve. [ABSTRACT FROM AUTHOR]

Congracilaria babae was first reported as a red alga parasitic on the thallus of Gracilaria salicornia based on Japanese materials. It was circumscribed to have deep spermatangial cavities, coloration similar to its host and the absence of rhizoids. We observed a parasitic red alga with morphological and anatomical features suggestive of C. babae on a Hydropuntia species collected from Sabah, East Malaysia. We addressed the taxonomic affinities of the parasite growing on Hydropuntia sp. based on the DNA sequence of molecular markers from the nuclear, mitochondrial and plastid genomes (nuclear ITS region, mitochondrial cox1 gene and plastid rbcL gene). Phylogenetic analyses based on all genetic markers also implied the monophyly of the parasite from Hydropuntia sp. and C. babae, suggesting their conspecificity. The parasite from Hydropuntia sp. has a DNA signature characteristic to C. babae in having plastid rbcL gene sequence identical to G. salicornia. C. babae is likely to have evolved directly from G. salicornia and subsequently radiated onto a secondary host Hydropuntia sp. We also recommend the transfer of C. babae to the genus Gracilaria and propose a new combination, G. babae, based on the anatomical observations and molecular data. [ABSTRACT FROM AUTHOR]

The primate α-/θ-defensin multigene family encodes versatile endogenous cationic and amphipathic peptides that have broad-spectrum antibacterial, antifungal and antiviral activity. Although previous studies have reported that α-/θ-defensin (DEFA/DEFT) genes are under birth-and-death evolution with frequent duplication and rapid evolution, the phylogenetic relationships of the primate DEFA/DEFT genes; the genetic bases for the existence of similar antimicrobial spectra among closely related species; and the evolutionary processes involved in the emergence of cyclic θ-defensins in Old World monkeys and their subsequent loss of function in humans, chimpanzees and gorillas require further investigation. In this study, the DEFA/DEFT gene repertoires from primate and treeshrew were collected, followed by detailed phylogenetic, sequence and structure, selection pressure and comparative genomics analyses. All treeshrew, prosimian and simian DEFA/DEFT genes are grouped into two major clades, which are tissue-specific for enteric and myeloid defensins in simians. The simian enteric and myeloid α-defensins are classified into six functional gene clusters with diverged sequences, variable structures, altered functional constraints and different selection pressures, which likely reflect the antimicrobial spectra among closely related species. Species-specific duplication or pseudogenization within each simian cluster implies that the antimicrobial spectrum is ever-shifting, most likely challenged by the ever-changing pathogen environment. The DEFT evolved from the myeloid DEFA8. The prosegment of θ-defensin is detected with adaptive changes coevolving with the new protein fold of mature peptide, coincident with the importance of the prosegment for the correct folding of the mature peptide. Lastly, a less-is-hitchhiking hypothesis was proposed as a possible explanation for the expansion of pseudogene DEFTP and the loss of functional DEFT, where the gain or loss of the hitchhiker is determined by its adjacent driver gene during the birth-and-death evolutionary process. [ABSTRACT FROM AUTHOR]

The functional divergence of transcriptional factors is critical in the evolution of transcriptional regulation. However, the mechanism of functional divergence among these factors remains unclear. Here, we performed an evolutionary analysis for positive selection in members of the myogenic regulatory factor (MRF) gene family of vertebrates. We selected 153 complete vertebrate MRF nucleotide sequences from our analyses, which revealed substantial evidence of positive selection. Here, we show that sites under positive selection were more frequently detected and identified from the genes encoding the myogenic differentiation factors (MyoG and Myf6) than the genes encoding myogenic determination factors (Myf5 and MyoD). Additionally, the functional divergence within the myogenic determination factors or differentiation factors was also under positive selection pressure. The positive selection sites were more frequently detected from MyoG and MyoD than Myf6 and Myf5, respectively. Amino acid residues under positive selection were identified mainly in their transcription activation domains and on the surface of protein three-dimensional structures. These data suggest that the functional gain and divergence of myogenic regulatory factors were driven by distinct positive selection of their transcription activation domains, whereas the function of the DNA binding domains was conserved in evolution. Our study evaluated the mechanism of functional divergence of the transcriptional regulation factors within a family, whereby the functions of their transcription activation domains diverged under positive selection during evolution. [ABSTRACT FROM AUTHOR]

Through the study of the phylogeographic structure and demographic history of the common goby, Pomatoschistus microps, the influence of Quaternary climatic changes on the evolutionary history of coastal and marine fishes is investigated. Because of its sedentary life cycle in Mediterranean lagoons, it is also a good model to study more specifically if the formation of lagoons during the Holocene had an impact on population structure and demography. Mitochondrial sequences of Northeastern Atlantic and Western Mediterranean specimens were used for phylogenetic reconstructions as well as divergence time estimates, demographic history and population structure analyses. Pomatoschistus microps was a highly supported monophyletic clade including four lineages. It may have appeared 77,000 yr ago, and the divergence of its lineages likely occurred shortly thereafter (between 61,000 and 54,000 yr). Most lineages had polytomic topologies, low nucleotide diversity and demographic analyses providing evidence of population expansion. Each lineage was characterized by a large number of private haplotypes. Most haplotypes found in Mediterranean localities were endemic, and one was dominant. Complex reticulated relationships connecting North European, Atlantic and Mediterranean haplotypes were observed. Moderate to high population structure was underlined. Contrary to previous published studies, no significant differentiation was observed between Atlantic and Mediterranean populations, indicating that the Gibraltar Strait is not a phylogeographic break for P. microps. Indeed, molecular dating combined with the tree topologies, phylogeographic and demographic analyses as well as high haplotype diversity underline a recent and rapid population divergence during the last glacial. However, population structure indicates that differentiation is an ongoing process. From an ancestral population trapped in the Atlantic, this goby colonized first northern Europe and later the Mediterranean Sea. Shared haplotypes could have dispersed in the western Mediterranean basin before the lagoon formation, while most private haplotypes, evidencing a recent isolation, probably diverged in lagoons after their closure. [ABSTRACT FROM AUTHOR]

Human respiratory syncytial virus (HRSV) has three surface glycoproteins: small hydrophobic (SH), attachment (G) and fusion (F), encoded by three consecutive genes (SH-G-F). A 270-nt fragment of the G gene is used to genotype HRSV isolates. This study genotyped and investigated the variability of the gene and amino acid sequences of the three surface proteins of HRSV strains collected from 1987 to 2005 from one center. Sixty original clinical isolates and 5 prototype strains were analyzed. Sequences containing SH, F and G genes were generated, and multiple alignments and phylogenetic trees were analyzed. Genetic variability by protein domains comparing virus genotypes was assessed. Complete sequences of the SH-G-F genes were obtained for all 65 samples: HRSV-A = 35; HRSV-B = 30. In group A strains, genotypes GA5 and GA2 were predominant. For HRSV-B strains, the genotype GB4 was predominant from 1992 to 1994 and only genotype BA viruses were detected in 2004–2005. Different genetic variability at nucleotide level was detected between the genes, with G gene being the most variable and the highest variability detected in the 270-nt G fragment that is frequently used to genotype the virus. High variability (>10%) was also detected in the signal peptide and transmembrane domains of the F gene of HRSV A strains. Variability among the HRSV strains resulting in non-synonymous changes was detected in hypervariable domains of G protein, the signal peptide of the F protein, a not previously defined domain in the F protein, and the antigenic site Ø in the pre-fusion F. Divergent trends were observed between HRSV -A and -B groups for some functional domains. A diverse population of HRSV -A and -B genotypes circulated in Houston during an 18 year period. We hypothesize that diverse sequence variation of the surface protein genes provide HRSV strains a survival advantage in a partially immune-protected community. [ABSTRACT FROM AUTHOR]

Most circulating strains of Human enterovirus 71 (EV-A71) have been classified primarily into three genogroups (A to C) on the basis of genetic divergence between the 1D gene, which encodes the VP1 capsid protein. The aim of the present study was to provide further insights into the diversity of the EV-A71 genogroups following the recent description of highly divergent isolates, in particular those from African countries, including Madagascar. We classified recent EV-A71 isolates by a large comparison of 3,346 VP1 nucleotidic sequences collected from GenBank. Analysis of genetic distances and phylogenetic investigations indicated that some recently-reported isolates did not fall into the genogroups A-C and clustered into three additional genogroups, including one Indian genogroup (genogroup D) and 2 African ones (E and F). Our Bayesian phylogenetic analysis provided consistent data showing that the genogroup D isolates share a recent common ancestor with the members of genogroup E, while the isolates of genogroup F evolved from a recent common ancestor shared with the members of the genogroup B. Our results reveal the wide diversity that exists among EV-A71 isolates and suggest that the number of circulating genogroups is probably underestimated, particularly in developing countries where EV-A71 epidemiology has been poorly studied. [ABSTRACT FROM AUTHOR]

Evolutionary molecular studies of island radiations may lead to insights in the role of vicariance, founder events, population size and drift in the processes of population differentiation. We evaluate the degree of population genetic differentiation and fixation of the Canary Islands blue tit subspecies complex using microsatellite markers and aim to get insights in the population history using coalescence based methods. The Canary Island populations were strongly genetically differentiated and had reduced diversity with pronounced fixation including many private alleles. In population structure models, the relationship between the central island populations (La Gomera, Tenerife and Gran Canaria) and El Hierro was difficult to disentangle whereas the two European populations showed consistent clustering, the two eastern islands (Fuerteventura and Lanzarote) and Morocco weak clustering, and La Palma a consistent unique lineage. Coalescence based models suggested that the European mainland forms an outgroup to the Afrocanarian population, a split between the western island group (La Palma and El Hierro) and the central island group, and recent splits between the three central islands, and between the two eastern islands and Morocco, respectively. It is clear that strong genetic drift and low level of concurrent gene flow among populations have shaped complex allelic patterns of fixation and skewed frequencies over the archipelago. However, understanding the population history remains challenging; in particular, the pattern of extreme divergence with low genetic diversity and yet unique genetic material in the Canary Island system requires an explanation. A potential scenario is population contractions of a historically large and genetically variable Afrocanarian population, with vicariance and drift following in the wake. The suggestion from sequence-based analyses of a Pleistocene extinction of a substantial part of North Africa and a Pleistocene/Holocene eastward re-colonisation of western North Africa from the Canaries remains possible. [ABSTRACT FROM AUTHOR]

Rapid phenotypic adaptation is critical for populations facing environmental changes and can be facilitated by phenotypic plasticity in the selected traits. Whereas recurrent environmental fluctuations can favour the maintenance or de novo evolution of plasticity, strong selection is hypothesized to decrease plasticity or even fix the trait (genetic assimilation). Despite advances in the theoretical understanding of the impact of plasticity on diversification processes, comparatively little empirical data of populations undergoing diversification mediated by plasticity are available. Here we use the planktonic freshwater copepod Acanthodiaptomus denticornis from two lakes as model system to study UV stress responses of two phenotypically different populations under laboratory conditions. Our study reveals heritable lake- and sex-specific differences of behaviour, physiological plasticity, and mortality. We discuss specific selective scenarios causing these differences and argue that phenotypic plasticity will be higher when selection pressure is moderate, but will decrease or even be lost under stronger pressure. [ABSTRACT FROM AUTHOR]

Sauria is the crown-group of Diapsida and is subdivided into Lepidosauromorpha and Archosauromorpha, comprising a high percentage of the diversity of living and fossil tetrapods. The split between lepidosauromorphs and archosauromorphs (the crocodile-lizard, or bird-lizard, divergence) is considered one of the key calibration points for molecular analyses of tetrapod phylogeny. Saurians have a very rich Mesozoic and Cenozoic fossil record, but their late Paleozoic (Permian) record is problematic. Several Permian specimens have been referred to Sauria, but the phylogenetic affinity of some of these records remains questionable. We reexamine and review all of these specimens here, providing new data on early saurian evolution including osteohistology, and present a new morphological phylogenetic dataset. We support previous studies that find that no valid Permian record for Lepidosauromorpha, and we also reject some of the previous referrals of Permian specimens to Archosauromorpha. The most informative Permian archosauromorph is Protorosaurus speneri from the middle Late Permian of Western Europe. A historically problematic specimen from the Late Permian of Tanzania is redescribed and reidentified as a new genus and species of basal archosauromorph: Aenigmastropheus parringtoni. The supposed protorosaur Eorasaurus olsoni from the Late Permian of Russia is recovered among Archosauriformes and may be the oldest known member of the group but the phylogenetic support for this position is low. The assignment of Archosaurus rossicus from the latest Permian of Russia to the archosauromorph clade Proterosuchidae is supported. Our revision suggests a minimum fossil calibration date for the crocodile-lizard split of 254.7 Ma. The occurrences of basal archosauromorphs in the northern (30°N) and southern (55°S) parts of Pangea imply a wider paleobiogeographic distribution for the group during the Late Permian than previously appreciated. Early archosauromorph growth strategies appear to be more diverse than previously suggested based on new data on the osteohistology of Aenigmastropheus. [ABSTRACT FROM AUTHOR]

Contemporary patterns of genetic diversity and population connectivity within species can be influenced by both historical and contemporary barriers to gene flow. In the marine environment, present day oceanographic features such as currents, fronts and upwelling systems can influence dispersal of eggs/larvae and/juveniles/adults, shaping population substructuring. The Benguela Current system in the southeastern Atlantic is one of the oldest upwelling systems in the world, and provides a unique opportunity to investigate the relative influence of contemporary and historical mechanisms shaping the evolutionary history of warm-temperate fish species. Using the genetic variation in the mitochondrial DNA Control Region and eight nuclear microsatellite DNA loci, we identified the presence of two highly divergent populations in a vagile and warm-temperate fish species, Atractoscion aequidens, across the Benguela region. The geographical distributions of the two populations, on either side of the perennial upwelling cell, suggest a strong correlation between the oceanographic features of the system and the breakdown of gene flow within this species. Genetic divergence (mtDNA φST = 0.902, microsatellite FST = 0.055: probability of genetic homogeneity for either marker = p<0.001), absence of migrants (less than 1% per generation) between populations and coalescent estimates of time since most recent common ancestor suggest that the establishment of the main oceanographic features of the system (2 million years ago), particularly the strengthening and position of the perennial upwelling cell, is the most likely mechanism behind the observed isolation. Concordance between mitochondrial and nuclear genetic markers indicates that isolation and divergence of the northern and southern Benguela populations of A. aequidens occurred deep in the past and has continued to the present day. These findings suggest that the Benguela Current system may constitute an ancient and impermeable barrier to gene flow for warm-temperate fish species. [ABSTRACT FROM AUTHOR]

The Rhodnius Pacific group is composed of three species: Rhodnius pallescens, R. colombiensis and R. ecuadoriensis, which are considered important vectors of trypanosomes (Trypanosoma cruzi and T. rangeli) infecting humans. This group is considered as a recent trans-Andean lineage derived from the widespread distributed sister taxa R. pictipes during the later uplift of northern Andes mountain range. The widest spread species R. pallescens may be a complex of two divergent lineages with different chromosomal attributes and a particular biogeographical distribution across Central America and Colombia with several southern populations in Colombia occupying the same sylvatic habitat as its sister species R. colombiensis. Although the taxonomy of Rhodnius Pacific group has been well studied, the unresolved phylogenetic and systematic issues are the target of this paper. Here we explore the molecular phylogeography of this species group analyzing two mitochondrial (ND4 and cyt b) and one nuclear (D2 region of ribosomal 28S gene) gene sequences. The molecular analyses suggest an early divergence of the species R. ecuadoriensis and R. colombiensis, followed by a recent expansion of R. pallescens lineages. The phylogenetic relationship between sympatric R. pallescens Colombian lineage and R. colombiensis was further explored using wing morphometry, DNA genome size measurements, and by analyzing chromosomal behavior of hybrids progeny obtained from experimental crosses. Our results suggest that the diversification of the two R. pallescens lineages was mainly influenced by biogeographical events such as (i) the emergence of the Panama Isthmus, while the origin and divergence of R. colombiensis was associated with (ii) the development of particular genetic and chromosomal features that act as isolation mechanisms from its sister species R. pallescens (Colombian lineage). These findings provide new insights into the evolution of the Rhodnius Pacific group and the underlying biological processes that occurred during its divergence. [ABSTRACT FROM AUTHOR]

Divergence of acoustic signals in a geographic scale results from diverse evolutionary forces acting in parallel and affecting directly inter-male vocal interactions among disjunct populations. Pleurodema thaul is a frog having an extensive latitudinal distribution in Chile along which males' advertisement calls exhibit an important variation. Using the playback paradigm we studied the evoked vocal responses of males of three populations of P. thaul in Chile, from northern, central and southern distribution. In each population, males were stimulated with standard synthetic calls having the acoustic structure of local and foreign populations. Males of both northern and central populations displayed strong vocal responses when were confronted with the synthetic call of their own populations, giving weaker responses to the call of the southern population. The southern population gave stronger responses to calls of the northern population than to the local call. Furthermore, males in all populations were stimulated with synthetic calls for which the dominant frequency, pulse rate and modulation depth were varied parametrically. Individuals from the northern and central populations gave lower responses to a synthetic call devoid of amplitude modulation relative to stimuli containing modulation depths between 30–100%, whereas the southern population responded similarly to all stimuli in this series. Geographic variation in the evoked vocal responses of males of P. thaul underlines the importance of inter-male interactions in driving the divergence of the acoustic traits and contributes evidence for a role of intra-sexual selection in the evolution of the sound communication system of this anuran. [ABSTRACT FROM AUTHOR]

Phylogeographic studies of Holarctic birds are challenging because they involve vast geographic scale, complex glacial history, extensive phenotypic variation, and heterogeneous taxonomic treatment across countries, all of which require large sample sizes. Knowledge about the quality of phylogeographic information provided by different loci is crucial for study design. We use sequences of one mtDNA gene, one sex-linked intron, and one autosomal intron to elucidate large scale phylogeographic patterns in the Holarctic lark genus Eremophila. The mtDNA ND2 gene identified six geographically, ecologically, and phenotypically concordant clades in the Palearctic that diverged in the Early - Middle Pleistocene and suggested paraphyly of the horned lark (E. alpestris) with respect to the Temminck's lark (E. bilopha). In the Nearctic, ND2 identified five subclades which diverged in the Late Pleistocene. They overlapped geographically and were not concordant phenotypically or ecologically. Nuclear alleles provided little information on geographic structuring of genetic variation in horned larks beyond supporting the monophyly of Eremophila and paraphyly of the horned lark. Multilocus species trees based on two nuclear or all three loci provided poor support for haplogroups identified by mtDNA. The node ages calculated using mtDNA were consistent with the available paleontological data, whereas individual nuclear loci and multilocus species trees appeared to underestimate node ages. We argue that mtDNA is capable of discovering independent evolutionary units within avian taxa and can provide a reasonable phylogeographic hypothesis when geographic scale, geologic history, and phenotypic variation in the study system are too complex for proposing reasonable a priori hypotheses required for multilocus methods. Finally, we suggest splitting the currently recognized horned lark into five Palearctic and one Nearctic species. [ABSTRACT FROM AUTHOR]

The skin secretion of many amphibians contains peptides that are able to kill a broad range of microorganisms (antimicrobial peptides: AMPs) and potentially play a role in innate immune defense. Similar to the toxin arsenals of various animals, amphibian AMP repertoires typically show major structural variation, and previous studies have suggested that this may be the result of diversifying selection in adaptation to a diverse spectrum of pathogens. Here we report on transcriptome analyses that indicate a very different pattern in the dwarf clawed frog H. boettgeri. Our analyses reveal a diverse set of transcripts containing two to six tandem repeats, together encoding 14 distinct peptides. Five of these have recently been identified as AMPs, while three more are shown here to potently inhibit the growth of gram-negative bacteria, including multi-drug resistant strains of the medically important Pseudomonas aeruginosa. Although the number of predicted peptides is similar to the numbers of related AMPs in Xenopus and Silurana frog species, they show significantly lower structural variation. Selection analyses confirm that, in contrast to the AMPs of other amphibians, the H. boettgeri peptides did not evolve under diversifying selection. Instead, the low sequence variation among tandem repeats resulted from purifying selection, recent duplication and/or concerted gene evolution. Our study demonstrates that defense peptide repertoires of closely related taxa, after diverging from each other, may evolve under differential selective regimes, leading to contrasting patterns of structural diversity. [ABSTRACT FROM AUTHOR]

The analysis of genes in evolutionarily distant but morphologically similar species is of major importance to unravel the changes in genomes over millions of years, which led to gene silencing and functional diversification. We report the analysis of Wnt8a gene expression in the medakafish and provide a detailed comparison to other vertebrates. In all teleosts analyzed there are two paralogous Wnt8a copies. These show largely overlapping expression in the early developing zebrafish embryo, an evolutionarily distant relative of medaka. In contrast to zebrafish, we find that both maternal and zygotic expression of particularly one Wnt8a paralog has diverged in medaka. While Wnt8a1 expression is mostly conserved at early embryonic stages, the expression of Wnt8a2 differs markedly. In addition, both genes are distinctly expressed during organogenesis unlike the zebrafish homologs, which may hint at the emergence of functional diversification of Wnt8a ligands during evolution. [ABSTRACT FROM AUTHOR]

It is likely that the strength of selection acting upon a mutation varies through time due to changes in the environment. However, most population genetic theory assumes that the strength of selection remains constant. Here we investigate the consequences of fluctuating selection pressures on the quantification of adaptive evolution using McDonald-Kreitman (MK) style approaches. In agreement with previous work, we show that fluctuating selection can generate evidence of adaptive evolution even when the expected strength of selection on a mutation is zero. However, we also find that the mutations, which contribute to both polymorphism and divergence tend, on average, to be positively selected during their lifetime, under fluctuating selection models. This is because mutations that fluctuate, by chance, to positive selected values, tend to reach higher frequencies in the population than those that fluctuate towards negative values. Hence the evidence of positive adaptive evolution detected under a fluctuating selection model by MK type approaches is genuine since fixed mutations tend to be advantageous on average during their lifetime. Never-the-less we show that methods tend to underestimate the rate of adaptive evolution when selection fluctuates. [ABSTRACT FROM AUTHOR]

Large-scale phylogenies provide a valuable source to study background diversification rates and investigate if the rates have changed over time. Unfortunately most large-scale, dated phylogenies are sparsely sampled (fewer than 5% of the described species) and taxon sampling is not uniform. Instead, taxa are frequently sampled to obtain at least one representative per subgroup (e.g. family) and thus to maximize diversity (diversified sampling). So far, such complications have been ignored, potentially biasing the conclusions that have been reached. In this study I derive the likelihood of a birth-death process with non-constant (time-dependent) diversification rates and diversified taxon sampling. Using simulations I test if the true parameters and the sampling method can be recovered when the trees are small or medium sized (fewer than 200 taxa). The results show that the diversification rates can be inferred and the estimates are unbiased for large trees but are biased for small trees (fewer than 50 taxa). Furthermore, model selection by means of Akaike's Information Criterion favors the true model if the true rates differ sufficiently from alternative models (e.g. the birth-death model is recovered if the extinction rate is large and compared to a pure-birth model). Finally, I applied six different diversification rate models – ranging from a constant-rate pure birth process to a decreasing speciation rate birth-death process but excluding any rate shift models – on three large-scale empirical phylogenies (ants, mammals and snakes with respectively 149, 164 and 41 sampled species). All three phylogenies were constructed by diversified taxon sampling, as stated by the authors. However only the snake phylogeny supported diversified taxon sampling. Moreover, a parametric bootstrap test revealed that none of the tested models provided a good fit to the observed data. The model assumptions, such as homogeneous rates across species or no rate shifts, appear to be violated. [ABSTRACT FROM AUTHOR]

Phylogeographic studies of Holarctic birds are challenging because they involve vast geographic scale, complex glacial history, extensive phenotypic variation, and heterogeneous taxonomic treatment across countries, all of which require large sample sizes. Knowledge about the quality of phylogeographic information provided by different loci is crucial for study design. We use sequences of one mtDNA gene, one sex-linked intron, and one autosomal intron to elucidate large scale phylogeographic patterns in the Holarctic lark genus Eremophila. The mtDNA ND2 gene identified six geographically, ecologically, and phenotypically concordant clades in the Palearctic that diverged in the Early - Middle Pleistocene and suggested paraphyly of the horned lark (E. alpestris) with respect to the Temminck's lark (E. bilopha). In the Nearctic, ND2 identified five subclades which diverged in the Late Pleistocene. They overlapped geographically and were not concordant phenotypically or ecologically. Nuclear alleles provided little information on geographic structuring of genetic variation in horned larks beyond supporting the monophyly of Eremophila and paraphyly of the horned lark. Multilocus species trees based on two nuclear or all three loci provided poor support for haplogroups identified by mtDNA. The node ages calculated using mtDNA were consistent with the available paleontological data, whereas individual nuclear loci and multilocus species trees appeared to underestimate node ages. We argue that mtDNA is capable of discovering independent evolutionary units within avian taxa and can provide a reasonable phylogeographic hypothesis when geographic scale, geologic history, and phenotypic variation in the study system are too complex for proposing reasonable a priori hypotheses required for multilocus methods. Finally, we suggest splitting the currently recognized horned lark into five Palearctic and one Nearctic species. [ABSTRACT FROM AUTHOR]

Protein families might evolve paralogous functions on their common tertiary scaffold in two ways. First, the locations of functionally-important sites might be “hard-wired” into the structure, with novel functions evolved by altering the amino acid (e.g. Ala vs Ser) at these positions. Alternatively, the tertiary scaffold might be adaptable, accommodating a unique set of functionally important sites for each paralogous function. To discriminate between these possibilities, we compared the set of functionally important sites in the six largest paralogous subfamilies of the LacI/GalR transcription repressor family. LacI/GalR paralogs share a common tertiary structure, but have low sequence identity (≤30%), and regulate a variety of metabolic processes. Functionally important positions were identified by conservation and co-evolutionary sequence analyses. Results showed that conserved positions use a mixture of the “hard-wired” and “accommodating” scaffold frameworks, but that the co-evolution networks were highly dissimilar between any pair of subfamilies. Therefore, the tertiary structure can accommodate multiple networks of functionally important positions. This possibility should be included when designing and interpreting sequence analyses of other protein families. Software implementing conservation and co-evolution analyses is available at https://sourceforge.net/projects/coevolutils/. [ABSTRACT FROM AUTHOR]

Local adaptation towards divergent ecological conditions often results in genetic differentiation and adaptive phenotypic divergence. To illuminate the ecological distinctiveness of the schizothoracine fish, we studied a Gymnocypris species complex consisting of three morphs distributed across four bodies of water (the Yellow River, Lake Qinghai, the Ganzi River and Lake Keluke) in the Northeast Tibetan Plateau. We used a combination of mitochondrial (16S rRNA and Cyt b) and nuclear (RAG-2) genetic sequences to investigate the phylogeography of these morphs based on a sample of 277 specimens. Analysis of gill rakers allowed for mapping of phenotypic trajectories along the phylogeny. The phylogenetic and morphological analyses showed that the three sparsely rakered morphs were present at two extremes of the phylogenetic tree: the Yellow River morphs were located at the basal phylogenetic split, and the Lake Keluke and Ganzi River morphs at the peak, with the densely rakered Lake Qinghai morphs located between these two extremes. Age estimation further indicated that the sparsely rakered morphs constituted the oldest and youngest lineages, whereas the densely rakered morph was assigned to an intermediate-age lineage. These results are most compatible with the process of evolutionary convergence or reversal. Disruptive natural selection due to divergent habitats and dietary preferences is likely the driving force behind the formation of new morphs, and the similarities between their phenotypes may be attributable to the similarities between their forms of niche tracking associated with food acquisition. This study provides the first genetic evidence for the occurrence of convergence or reversal in the schizothoracine fish of the Tibetan Plateau at small temporal scales. [ABSTRACT FROM AUTHOR]

Patterns of intraspecific geographic variation of signaling systems provide insight into the microevolutionary processes driving phenotypic divergence. The acoustic calls of bats are sensitive to diverse evolutionary forces, but processes that shape call variation are largely unexplored. In China, Rhinolophus ferrumequinum displays a diverse call frequency and inhabits a heterogeneous landscape, presenting an excellent opportunity for this kind of research. We quantified geographic variation in resting frequency (RF) of echolocation calls, estimated genetic structure and phylogeny of R. ferrumequinum populations, and combined this with climatic factors to test three hypotheses to explain acoustic variation: genetic drift, cultural drift, and local adaptation. Our results demonstrated significant regional divergence in frequency and phylogeny among the bat populations in China's northeast (NE), central-east (CE) and southwest (SW) regions. The CE region had higher frequencies than the NE and SW regions. Drivers of RF divergence were estimated in the entire range and just the CE/NE region (since these two regions form a clade). In both cases, RF divergence was not correlated with mtDNA or nDNA genetic distance, but was significantly correlated with geographic distance and mean annual temperature, indicating cultural drift and ecological selection pressures are likely important in shaping RF divergence among different regions in China. [ABSTRACT FROM AUTHOR]

The present study investigated the genetic diversity, population structure, FST outliers, and extent and pattern of linkage disequilibrium in five populations of Keteleeria davidiana var. formosana, which is listed as a critically endangered species by the Council of Agriculture, Taiwan. Twelve amplified fragment length polymorphism primer pairs generated a total of 465 markers, of which 83.74% on average were polymorphic across populations, with a mean Nei’s genetic diversity of 0.233 and a low level of genetic differentiation (approximately 6%) based on the total dataset. Linkage disequilibrium and HICKORY analyses suggested recent population bottlenecks and inbreeding in K. davidiana var. formosana. Both STRUCTURE and BAPS observed extensive admixture of individual genotypes among populations based on the total dataset in various clustering scenarios, which probably resulted from incomplete lineage sorting of ancestral variation rather than a high rate of recent gene flow. Our results based on outlier analysis revealed generally high levels of genetic differentiation and suggest that divergent selection arising from environmental variation has been driven by differences in temperature, precipitation, and humidity. Identification of ecologically associated outliers among environmentally disparate populations further support divergent selection and potential local adaptation. [ABSTRACT FROM AUTHOR]

Phytophthora and other oomycetes secrete a large number of putative host cytoplasmic effectors with conserved FLAK motifs following signal peptides, termed crinkling and necrosis inducing proteins (CRN), or Crinkler. Here, we first investigated the evolutionary patterns and mechanisms of CRN effectors in Phytophthora sojae and compared them to two other Phytophthora species. The genes encoding CRN effectors could be divided into 45 orthologous gene groups (OGG), and most OGGs unequally distributed in the three species, in which each underwent large number of gene gains or losses, indicating that the CRN genes expanded after species evolution in Phytophthora and evolved through pathoadaptation. The 134 expanded genes in P. sojae encoded family proteins including 82 functional genes and expressed at higher levels while the other 68 genes encoding orphan proteins were less expressed and contained 50 pseudogenes. Furthermore, we demonstrated that most expanded genes underwent gene duplication or/and fragment recombination. Three different mechanisms that drove gene duplication or recombination were identified. Finally, the expanded CRN effectors exhibited varying pathogenic functions, including induction of programmed cell death (PCD) and suppression of PCD through PAMP-triggered immunity or/and effector-triggered immunity. Overall, these results suggest that gene duplication and fragment recombination may be two mechanisms that drive the expansion and neofunctionalization of the CRN family in P. sojae, which aids in understanding the roles of CRN effectors within each oomycete pathogen. [ABSTRACT FROM AUTHOR]

Estimated phylogenies of evolutionarily diverse taxa will be well supported and more likely to be historically accurate when the analysis contains large amounts of data–many genes sequenced across many taxa. Inferring such phylogenies for non-model organisms is challenging given limited resources for whole-genome sequencing. We take advantage of genomic data from a single species to test the limits of hybridization-based enrichment of hundreds of exons across frog species that diverged up to 250 million years ago. Enrichment success for a given species depends greatly on the divergence time between it and the reference species, and the resulting alignment contains a significant proportion of missing data. However, our alignment generates a well-supported phylogeny of frogs, suggesting that this technique is a practical solution towards resolving relationships across deep evolutionary time. [ABSTRACT FROM AUTHOR]

Phytophthora and other oomycetes secrete a large number of putative host cytoplasmic effectors with conserved FLAK motifs following signal peptides, termed crinkling and necrosis inducing proteins (CRN), or Crinkler. Here, we first investigated the evolutionary patterns and mechanisms of CRN effectors in Phytophthora sojae and compared them to two other Phytophthora species. The genes encoding CRN effectors could be divided into 45 orthologous gene groups (OGG), and most OGGs unequally distributed in the three species, in which each underwent large number of gene gains or losses, indicating that the CRN genes expanded after species evolution in Phytophthora and evolved through pathoadaptation. The 134 expanded genes in P. sojae encoded family proteins including 82 functional genes and expressed at higher levels while the other 68 genes encoding orphan proteins were less expressed and contained 50 pseudogenes. Furthermore, we demonstrated that most expanded genes underwent gene duplication or/and fragment recombination. Three different mechanisms that drove gene duplication or recombination were identified. Finally, the expanded CRN effectors exhibited varying pathogenic functions, including induction of programmed cell death (PCD) and suppression of PCD through PAMP-triggered immunity or/and effector-triggered immunity. Overall, these results suggest that gene duplication and fragment recombination may be two mechanisms that drive the expansion and neofunctionalization of the CRN family in P. sojae, which aids in understanding the roles of CRN effectors within each oomycete pathogen. [ABSTRACT FROM AUTHOR]

The complete sequences of three mitochondrial genomes from the land snail Cornu aspersum were determined. The mitogenome has a length of 14050 bp, and it encodes 13 protein-coding genes, 22 transfer RNA genes and two ribosomal RNA genes. It also includes nine small intergene spacers, and a large AT-rich intergenic spacer. The intra-specific divergence analysis revealed that COX1 has the lower genetic differentiation, while the most divergent genes were NADH1, NADH3 and NADH4. With the exception of Euhadra herklotsi, the structural comparisons showed the same gene order within the family Helicidae, and nearly identical gene organization to that found in order Pulmonata. Phylogenetic reconstruction recovered Basommatophora as polyphyletic group, whereas Eupulmonata and Pulmonata as paraphyletic groups. Bayesian and Maximum Likelihood analyses showed that C. aspersum is a close relative of Cepaea nemoralis, and with the other Helicidae species form a sister group of Albinaria caerulea, supporting the monophyly of the Stylommatophora clade. [ABSTRACT FROM AUTHOR]

Serine protease inhibitors (SERPINs) are a superfamily of highly conserved proteins that play a key role in controlling the activity of proteases in diverse biological processes. The SERPIN cluster located at the 14q32.1 region includes the gene coding for SERPINA1, and a highly homologous sequence, SERPINA2, which was originally thought to be a pseudogene. We have previously shown that SERPINA2 is expressed in different tissues, namely leukocytes and testes, suggesting that it is a functional SERPIN. To investigate the function of SERPINA2, we used HeLa cells stably transduced with the different variants of SERPINA2 and SERPINA1 (M1, S and Z) and leukocytes as the in vivo model. We identified SERPINA2 as a 52 kDa intracellular glycoprotein, which is localized at the endoplasmic reticulum (ER), independently of the variant analyzed. SERPINA2 is not significantly regulated by proteasome, proposing that ER localization is not due to misfolding. Specific features of SERPINA2 include the absence of insoluble aggregates and the insignificant response to cell stress, suggesting that it is a non-polymerogenic protein with divergent activity of SERPINA1. Using phylogenetic analysis, we propose an origin of SERPINA2 in the crown of primates, and we unveiled the overall conservation of SERPINA2 and A1. Nonetheless, few SERPINA2 residues seem to have evolved faster, contributing to the emergence of a new advantageous function, possibly as a chymotrypsin-like SERPIN. Herein, we present evidences that SERPINA2 is an active gene, coding for an ER-resident protein, which may act as substrate or adjuvant of ER-chaperones. [ABSTRACT FROM AUTHOR]

Cartilaginous fishes, divided into Holocephali (chimaeras) and Elasmoblanchii (sharks, rays and skates), occupy a key phylogenetic position among extant vertebrates in reconstructing their evolutionary processes. Their accurate evolutionary time scale is indispensable for better understanding of the relationship between phenotypic and molecular evolution of cartilaginous fishes. However, our current knowledge on the time scale of cartilaginous fish evolution largely relies on estimates using mitochondrial DNA sequences. In this study, making the best use of the still partial, but large-scale sequencing data of cartilaginous fish species, we estimate the divergence times between the major cartilaginous fish lineages employing nuclear genes. By rigorous orthology assessment based on available genomic and transcriptomic sequence resources for cartilaginous fishes, we selected 20 protein-coding genes in the nuclear genome, spanning 2973 amino acid residues. Our analysis based on the Bayesian inference resulted in the mean divergence time of 421 Ma, the late Silurian, for the Holocephali-Elasmobranchii split, and 306 Ma, the late Carboniferous, for the split between sharks and rays/skates. By applying these results and other documented divergence times, we measured the relative evolutionary rate of the Hox A cluster sequences in the cartilaginous fish lineages, which resulted in a lower substitution rate with a factor of at least 2.4 in comparison to tetrapod lineages. The obtained time scale enables mapping phenotypic and molecular changes in a quantitative framework. It is of great interest to corroborate the less derived nature of cartilaginous fish at the molecular level as a genome-wide phenomenon. [ABSTRACT FROM AUTHOR]

Arctic charr (Salvelinus alpinus) is a highly polymorphic species and in Lake Thingvallavatn, Iceland, four phenotypic morphs have evolved. These differences in morphology, especially in craniofacial structures are already apparent during embryonic development, indicating that genes important in the formation of the craniofacial features are expressed differentially between the morphs. In order to generate tools to examine these expression differences in Arctic charr, the aim of the present study was to identify reference genes for quantitative real-time PCR (qPCR). The specific aim was to select reference genes which are able to detect very small expression differences among different morphs. We selected twelve candidate reference genes from the literature, identified corresponding charr sequences using data derived from transcriptome sequencing (RNA-seq) and examined their expression using qPCR. Many of the candidate reference genes were found to be stably expressed, yet their quality-rank as reference genes varied considerably depending on the type of analysis used. In addition to commonly used software for reference gene validation, we used classical statistics to evaluate expression profiles avoiding a bias for reference genes with similar expression patterns (co-regulation). Based on these analyses we chose three reference genes, ACTB, UB2L3 and IF5A1 for further evaluation. Their consistency was assessed in an expression study of three known craniofacially expressed genes, sparc (or osteonectin), matrix metalloprotease 2 (mmp2) and sox9 (sex-determining region Y box 9 protein) using qPCR in embryo heads derived from four charr groups at three developmental time points. The three reference genes were found to be very suitable for studying expression differences between the morphotypes, enabling robust detection of small relative expression changes during charr development. Further, the results showed that sparc and mmp2 are differentially expressed in embryos of different Arctic charr morphotypes. [ABSTRACT FROM AUTHOR]

EPF1-EPF2 and EPFL9/Stomagen act antagonistically in regulating leaf stomatal density. The aim of this study was to elucidate the evolutionary functional divergence of EPF/EPFL family genes. Phylogenetic analyses showed that AtEPFL9/Stomagen-like genes are conserved only in vascular plants and are closely related to AtEPF1/EPF2-like genes. Modeling showed that EPF/EPFL peptides share a common 3D structure that is constituted of a scaffold and loop. Molecular dynamics simulation suggested that AtEPF1/EPF2-like peptides form an additional disulfide bond in their loop regions and show greater flexibility in these regions than AtEPFL9/Stomagen-like peptides. This study uncovered the evolutionary relationship and the conformational divergence of proteins encoded by the EPF/EPFL family genes. [ABSTRACT FROM AUTHOR]

Serine protease inhibitors (SERPINs) are a superfamily of highly conserved proteins that play a key role in controlling the activity of proteases in diverse biological processes. The SERPIN cluster located at the 14q32.1 region includes the gene coding for SERPINA1, and a highly homologous sequence, SERPINA2, which was originally thought to be a pseudogene. We have previously shown that SERPINA2 is expressed in different tissues, namely leukocytes and testes, suggesting that it is a functional SERPIN. To investigate the function of SERPINA2, we used HeLa cells stably transduced with the different variants of SERPINA2 and SERPINA1 (M1, S and Z) and leukocytes as the in vivo model. We identified SERPINA2 as a 52 kDa intracellular glycoprotein, which is localized at the endoplasmic reticulum (ER), independently of the variant analyzed. SERPINA2 is not significantly regulated by proteasome, proposing that ER localization is not due to misfolding. Specific features of SERPINA2 include the absence of insoluble aggregates and the insignificant response to cell stress, suggesting that it is a non-polymerogenic protein with divergent activity of SERPINA1. Using phylogenetic analysis, we propose an origin of SERPINA2 in the crown of primates, and we unveiled the overall conservation of SERPINA2 and A1. Nonetheless, few SERPINA2 residues seem to have evolved faster, contributing to the emergence of a new advantageous function, possibly as a chymotrypsin-like SERPIN. Herein, we present evidences that SERPINA2 is an active gene, coding for an ER-resident protein, which may act as substrate or adjuvant of ER-chaperones. [ABSTRACT FROM AUTHOR]

Usutu virus (USUV) is a mosquito-borne flavivirus, belonging to the Japanese encephalitis antigenic complex, that circulates among mosquitoes and birds. We describe and analyze the complete genome sequence of the first USUV strain isolated from an immunocompromised patient with neuroinvasive disease. This USUV isolate showed an overall nucleotide identity of 99% and 96%, respectively, with the genomes of isolates from Europe and Africa. Comparison of the human USUV complete polyprotein sequence with bird-derived strains, showed two unique amino acid substitutions. In particular, one substitution (S595G) was situated in the DIII domain of the viral Envelope protein that is recognized by flavivirus neutralizing antibodies. An additional amino acid substitution (D3425E) was identified in the RNA-dependent RNA polymerase (RdRp) domain of the NS5 protein. This substitution is remarkable since E3425 is highly conserved among the other USUV isolates that were not associated with human infection. However, a similar substitution was observed in Japanese encephalitis and in West Nile viruses isolated from humans. Phylogenetic analysis of the human USUV strain revealed a close relationship with an Italian strain isolated in 2009. Analysis of synonymous nucleotide substitutions (SNSs) among the different USUV genomes showed a specific evolutionary divergence among different countries. In addition, 15 SNSs were identified as unique in the human isolate. We also identified four specific nucleotide substitutions in the 5′ and 3′ untranslated regions (UTRs) in the human isolate that were not present in the other USUV sequences. Our analyses provide the basis for further experimental studies aimed at defining the effective role of these mutations in the USUV genome, their potential role in the development of viral variants pathogenic for humans and their evolution and dispersal out of Africa. [ABSTRACT FROM AUTHOR]

Alternative splicing of genes generates novel mRNAs, leading to the evolution of new functional proteins. Cholecystokinin (CCK) induces the release of pancreatic enzymes and the contraction of the gallbladder to promote the digestion of fat and proteins. CCK activates two G-protein-coupled receptors, CCKA and CCKB. Here, we showed that a CCKsv (splicing variant), originated de novo during Catarrhini evolution by including a portion of intronic sequence of the CCK gene, encodes novel C-terminal peptide sequence followed by a new poly-adenylation signal. CCKsv is expressed in many human tissues and likely a secreted peptide retaining the original signal peptide and the N-terminal proteolytic processing signal, together with novel C-terminal sequences. Although CCKsv cannot activate CCK receptors, it partially inhibits the CRE- or SRF-driven reporter activities stimulated by wide type CCK-8 mediated by both CCK receptors. Co-treatment with CCKsv also partially antagonizes Ewing tumor cell growth stimulated by CCK-8. Our study provides an example of new peptide hormone antagonist evolution in primates. [ABSTRACT FROM AUTHOR]

The southern and tropical southeastern regions of the Yunnan Province in southwestern China have similar monsoonal climates and lowland tropical rain forest vegetations. The floras of both regions are dominated by tropical floristic elements (78.3% in southern Yunnan and 68.83% in southeastern Yunnan), and both belong to the Indo-Malaysian flora at the northern margin of tropical Asia. However, some temperate East Asian characteristic families are well represented in the flora of tropical southeastern Yunnan, while families characteristic of tropical Asia are well represented in the flora of southern Yunnan. Additionally, there are 14 mainly east Asian families in tropical southeastern Yunnan that are not found in southern Yunnan. Although the two regions share 80% of their genera, 237 genera are restricted to southern Yunnan, and 349 genera to tropical southeastern Yunnan. Furthermore, 57 genera with an East Asian distribution, 53 genera with a North temperate distribution, 22 genera endemic to China, and 17 genera with an East Asia and North America disjunct distribution are found only in tropical southeastern Yunnan. The flora of tropical southeastern Yunnan is more closely related to Eastern Asian flora, while the flora of southern Yunnan is more closely related to Indo-Malaysian flora. The divergence of the flora is well supported by the geological history of the region; the flora of tropical southeastern Yunnan was mainly derived from the South China Geoblock, while the southern Yunnan flora derived from the Shan-Thai Geoblock. [ABSTRACT FROM AUTHOR]

Divergent selection and local adaptation are responsible for many phenotypic differences between populations, potentially leading to speciation through the evolution of reproductive barriers. Here we evaluated the morphometric divergence among west European populations of Reed Bunting in order to determine the extent of local adaptation relative to two important selection pressures often associated with speciation in birds: migration and diet. We show that, as expected by theory, migratory E. s. schoeniclus had longer and more pointed wings and a slightly smaller body mass than the resident subspecies, with the exception of E. s. lusitanica, which despite having rounder wings was the smallest of all subspecies. Tail length, however, did not vary according to the expectation (shorter tails in migrants) probably because it is strongly correlated with wing length and might take longer to evolve. E. s. witherbyi, which feed on insects hiding inside reed stems during the winter, had a very thick, stubby bill. In contrast, northern populations, which feed on seeds, had thinner bills. Despite being much smaller, the southern E. s. lusitanica had a significantly thicker, longer bill than migratory E. s. schoeniclus, whereas birds from the UK population had significantly shorter, thinner bills. Geometric morphometric analyses revealed that the southern subspecies have a more convex culmen than E. s. schoeniclus, and E. s. lusitanica differs from the nominate subspecies in bill shape to a greater extent than in linear bill measurements, especially in males. Birds with a more convex culmen are thought to exert a greater strength at the bill tip, which is in agreement with their feeding technique. Overall, the three subspecies occurring in Western Europe differ in a variety of traits following the patterns predicted from their migratory and foraging behaviours, strongly suggesting that these birds have became locally adapted through natural selection. [ABSTRACT FROM AUTHOR]