Sequence variation, differential expression, and divergent evolution in starch-related genes among accessions of Arabidopsis thaliana

Transitory starch metabolism is a nonlinear and highly regulated process. It originated very early in the evolution of chloroplast-containing cells and is largely based on a mosaic of genes derived from either the eukaryotic host cell or the prokaryotic endosymbiont. Initially located in the cytoplasm, starch metabolism was rewired into plastids in Chloroplastida. Relocation was accompanied by gene duplications that occurred in most starch-related gene families and resulted in subfunctionalization of the respective gene products. Starch-related isozymes were then evolutionary conserved by constraints such as internal starch structure, posttranslational protein import into plastids and interactions with other starch-related proteins. 25 starch-related genes in 26 accessions of Arabidopsis thaliana were sequenced to assess intraspecific diversity, phylogenetic relationships, and modes of selection. Furthermore, sequences derived from additional 80 accessions that are publicly available were analyzed. Diversity varies significantly among the starch-related genes. Starch synthases and phosphorylases exhibit highest nucleotide diversities, while pyrophosphatases and debranching enzymes are most conserved. The gene trees are most compatible with a scenario of extensive recombination, perhaps in a Pleistocene refugium. Most genes are under purifying selection, but disruptive selection was inferred for a few genes/substitutiones. To study transcript levels, leaves were harvested throughout the light period. By quantifying the transcript levels and by analyzing the sequence of the respective accessions, we were able to estimate whether transcript levels are mainly determined by genetic (i.e., accession dependent) or physiological (i.e., time dependent) parameters. We also identified polymorphic sites that putatively affect pattern or the level of transcripts.

CyBorD induction therapy in clinical practice

Cyclophosphamide, bortezomib and dexamethasone (CyBorD) is a highly active three-drug induction regimen for untreated transplant-eligible multiple myeloma patients. Although CyBorD has been evaluated only in the phase 2 setting in a limited number of patients, its high efficacy and ease of administration have led to its widespread use. Given that clinical trial efficacy can overestimate real-life effectiveness, we reviewed our institutional experience with 109 newly diagnosed patients who were treated with CyBorD in a non-clinical trial setting. After a median of four cycles, overall response rate (ORR) and very good partial response rate or better (⩾VGPR) were 95 and 66%, respectively, comparable to phase 2 studies of CyBorD and other three/four-drug induction regimens. All patients subsequently underwent successful stem cell collection and upgraded responses to ORR 98% and ⩾VGPR 79% post transplant. At a median follow-up of 19.8 months after diagnosis, the 2-year OS probability was 95.3% (95%CI: 89-98). The presence of concurrent plasmacytoma at diagnosis was the only prognostic factor predicting poorer survival (HR=5.56; 95%CI: 0.92-33.74; P=0.03). CyBorD was well-tolerated, with no severe peripheral neuropathy and minimal hematologic toxicity. Therefore, CyBorD is a convenient, well-tolerated, highly effective induction regimen in preparation for autologous SCT in real-life clinical practice.

Effects of habitat structure and land-use intensity on the genetic structure of the grasshopper species Chorthippus parallelus

Land-use intensity (LUI) is assumed to impact the genetic structure of organisms. While effects of landscape structure on the genetics of local populations have frequently been analysed, potential effects of variation in LUI on the genetic diversity of local populations have mostly been neglected. In this study, we used six polymorphic microsatellites to analyse the genetic effects of variation in land use in the highly abundant grasshopper Chorthippus parallelus. We sampled a total of 610 individuals at 22 heterogeneous grassland sites in the Hainich-Dün region of Central Germany. For each of these grassland sites we assessed habitat size, LUI (combined index of mowing, grazing and fertilization), and the proportion of grassland adjoining the sampling site and the landscape heterogeneity (the latter two factors within a 500 m buffer zone surrounding each focal site). We found only marginal genetic differentiation among all local populations and no correlation between geographical and genetic distance. Habitat size, LUI and landscape characteristics had only weak effects on most of the parameters of genetic diversity of C. parallelus; only expected heterozygosity and the grasshopper abundances were affected by interacting effects of LUI, habitat size and landscape characteristics. The lack of any strong relationships between LUI, abundance and the genetic structure might be due to large local populations of the species in the landscape, counteracting local differentiation and potential genetic drift effects.

Isolation and characterization of eight microsatellite loci in the brook lamprey Lampetra planeri (Petromyzontiformes) using 454 sequence data

Eight polymorphic microsatellite loci were developed for the brook lamprey Lampetra planeri through 454 sequencing and their usefulness was tested in 45 individuals of both L. planeri and the river lamprey Lampetra fluviatilis. The number of alleles per loci ranged between two and five; the Italian and Irish populations had a mean expected heterozygosity of 0·388 and 0·424 and a mean observed heterozygosity of 0·418 and 0·411, respectively.

De novo assembly and characterization of the skeletal muscle and electric organ transcriptomes of the African weakly electric fish Campylomormyrus compressirostris (Mormyridae, Teleostei)

African weakly electric fishes (Mormyridae) underwent an outstanding adaptive radiation (about 200 species), putatively owing to their ability to communicate through species-specific weak electric signals. The electric organ discharge (EOD) is produced by muscle-derived electrocytes organized in piles to form an electric organ. Despite the importance of this trait as a prezygotic isolation mechanism, genomic resources remained limited. We present here a first draft of the skeletal muscle and electric organ transcriptomes from the weakly electric fish species Campylomormyrus compressirostris, obtained using the Illumina HiSeq2000 sequencing technology. Approximately 6.8 Gbp of cDNA sequence data were produced from both tissues, resulting in 57,268,109 raw reads for the skeletal muscle and 46,934,923 for the electric organ, and assembled de novo into 46,143 and 89,270 contigs, respectively. About 50% of both transcriptomes were annotated after protein databases search. The two transcriptomes show similar profiles in terms of Gene Ontology categories composition. We identified several candidate genes which are likely to play a central role in the production and evolution of the electric signal. For most of these genes, and for many other housekeeping genes, we were able to obtain the complete or partial coding DNA sequences (CDS), which can be used for the development of primers to be utilized in qRT-PCR experiments. We present also the complete mitochondrial genome and compare it to those available from other weakly electric fish species. Additionally, we located 1671 SSR-containing regions with their flanking sites and designed the relative primers. This study establishes a first step in the development of genomic tools aimed at understanding the role of electric communication during speciation.

Comparative analysis of the gonadal transcriptomes of the all-female species Poecilia formosa and its maternal ancestor Poecilia mexicana

Background

The Amazon molly, Poecilia formosa (Teleostei: Poeciliinae) is an unisexual, all-female species. It evolved through the hybridisation of two closely related sexual species and exhibits clonal reproduction by sperm dependent parthenogenesis (or gynogenesis) where the sperm of a parental species is only used to activate embryogenesis of the apomictic, diploid eggs but does not contribute genetic material to the offspring.

Here we provide and describe the first de novo assembled transcriptome of the Amazon molly in comparison with its maternal ancestor, the Atlantic molly Poecilia mexicana. The transcriptome data were produced through sequencing of single end libraries (100 bp) with the Illumina sequencing technique.

Results

83,504,382 reads for the Amazon molly and 81,625,840 for the Atlantic molly were assembled into 127,283 and 78,961 contigs for the Amazon molly and the Atlantic molly, respectively. 63% resp. 57% of the contigs could be annotated with gene ontology terms after sequence similarity comparisons. Furthermore, we were able to identify genes normally involved in reproduction and especially in meiosis also in the transcriptome dataset of the apomictic reproducing Amazon molly.

Conclusions

We assembled and annotated the transcriptome of a non-model organism, the Amazon molly, without a reference genome (de novo). The obtained dataset is a fundamental resource for future research in functional and expression analysis. Also, the presence of 30 meiosis-specific genes within a species where no meiosis is known to take place is remarkable and raises new questions for future research.

De novo assembly and characterization of the skeletal muscle and electric organ transcriptomes of the African weakly electric fish Campylomormyrus compressirostris (Mormyridae, Teleostei)

African weakly electric fishes (Mormyridae) underwent an outstanding adaptive radiation (about 200 species), putatively owing to their ability to communicate through species-specific weak electric signals. The electric organ discharge (EOD) is produced by muscle-derived electrocytes organized in piles to form an electric organ. Despite the importance of this trait as a prezygotic isolation mechanism, genomic resources remained limited. We present here a first draft of the skeletal muscle and electric organ transcriptomes from the weakly electric fish species Campylomormyrus compressirostris, obtained using the Illumina HiSeq2000 sequencing technology. Approximately 6.8 Gbp of cDNA sequence data were produced from both tissues, resulting in 57,268,109 raw reads for the skeletal muscle and 46,934,923 for the electric organ, and assembled de novo into 46,143 and 89,270 contigs, respectively. About 50% of both transcriptomes were annotated after protein databases search. The two transcriptomes show similar profiles in terms of Gene Ontology categories composition. We identified several candidate genes which are likely to play a central role in the production and evolution of the electric signal. For most of these genes, and for many other housekeeping genes, we were able to obtain the complete or partial coding DNA sequences (CDS), which can be used for the development of primers to be utilized in qRT-PCR experiments. We present also the complete mitochondrial genome and compare it to those available from other weakly electric fish species. Additionally, we located 1671 SSR-containing regions with their flanking sites and designed the relative primers. This study establishes a first step in the development of genomic tools aimed at understanding the role of electric communication during speciation.

Intraspecific rearrangement of duplicated mitochondrial control regions in the Luzon Tarictic Hornbill Penelopides manillae (Aves: Bucerotidae)

Philippine hornbills of the genera Aceros and Penelopides (Bucerotidae) are known to possess a large tandemly duplicated fragment in their mitochondrial genome, whose paralogous parts largely evolve in concert. In the present study, we surveyed the two distinguishable duplicated control regions in several individuals of the Luzon Tarictic Hornbill Penelopides manillae, compare their characteristics within and across individuals, and report on an intraspecific mitochondrial gene rearrangement found in one single specimen, i.e., an interchange between the two control regions. To our knowledge, this is the first observation of two distinct mitochondrial genome rearrangements within a bird species. We briefly discuss a possible evolutionary mechanism responsible for this pattern, and highlight potential implications for the application of control region sequences as a marker in population genetics and phylogeography.

Microsatellite usefulness is independent of phylogenetic distance in Tyrant flycatchers (Aves: Tyrannidae): a test using two globally threatened species

Tyrant flycatchers (Aves: Tyrannidae) are endemic to the New World, and many species of this group are threatened or near-threatened at the global level. The aim of this study was to test the 18 microsatellite markers that have been published for other Tyrant flycatchers in the Strange-tailed Tyrant (Alectrurus risora) and the Sharp-tailed Tyrant (Culicivora caudacuta), two endemic species of southern South American grasslands that are classified as vulnerable. We also analyzed the usefulness of loci in relation to phylogenetic distance to the source species. Amplification success was high in both species (77 to 83%) and did not differ between the more closely and more distantly related species to the source species. Polymorphism success was also similar for both species, with 9 and 8 loci being polymorphic, respectively. An increased phylogenetic distance thus does not gradually lead to allelic or locus dropouts, implying that in Tyrant flycatchers, the published loci are useful independent of species relatedness.

Diversification of Ponto-Caspian aquatic fauna: morphology and molecules retrieve congruent evolutionary relationships in Pontogammarus maeoticus (Amphipoda: Pontogammaridae)

The geological history of the Ponto-Caspian region, with alternating cycles of isolation and reconnection among the three main basins (Black and Azov Seas, and the more distant Caspian Sea) as well as between them and the Mediterranean Sea, profoundly affected the diversification of its aquatic fauna, leading to a high degree of endemism. Two alternative hypotheses on the origin of this amazing biodiversity have been proposed, corresponding to phases of allopatric separation of aquatic fauna among sea basins: a Late Miocene origin (10-6 MYA) vs. a more recent Pleistocene ancestry (<2 MYA). Both hypotheses support a vicariant origin of (1) Black+Azov Sea lineages on the one hand, and (2) Caspian Sea lineages on the other. Here, we present a study on the Ponto-Caspian endemic amphipod Pontogammarus maeoticus. We assessed patterns of divergence based on (a) two mitochondrial and one nuclear gene, and (b) a morphometric analysis of 23 morphological traits in 16 populations from South and West Caspian Sea, South Azov Sea and North-West Black Sea. Genetic data indicate a long and independent evolutionary history, dating back from the late Miocene to early Pleistocene (6.6-1.6 MYA), for an unexpected, major split between (i) a Black Sea clade and (ii) a well-supported clade grouping individuals from the Caspian and Azov Seas. Absence of shared haplotypes argues against either recent or human-mediated exchanges between Caspian and Azov Seas. A mismatch distribution analysis supports more stable population demography in the Caspian than in the Black Sea populations. Morphological divergence largely followed patterns of genetic divergence: our analyses grouped samples according to the basin of origin and corroborated the close phylogenetic affinity between Caspian and Azov Sea lineages. Altogether, our results highlight the necessity of careful (group-specific) evaluation of evolutionary trajectories in marine taxa that should certainly not be inferred from the current geographical proximity of sea basins alone.

Genetic connectivity between land and sea: the case of the beachflea Orchestia montagui (Crustacea, Amphipoda, Talitridae) in the Mediterranean Sea

INTRODUCTION

We examined patterns of genetic divergence in 26 Mediterranean populations of the semi-terrestrial beachflea Orchestia montagui using mitochondrial (cytochrome oxidase subunit I), microsatellite (eight loci) and allozymic data. The species typically forms large populations within heaps of dead seagrass leaves stranded on beaches at the waterfront. We adopted a hierarchical geographic sampling to unravel population structure in a species living at the sea-land transition and, hence, likely subjected to dramatically contrasting forces.

RESULTS

Mitochondrial DNA showed historical phylogeographic breaks among Adriatic, Ionian and the remaining basins (Tyrrhenian, Western and Eastern Mediterranean Sea) likely caused by the geological and climatic changes of the Pleistocene. Microsatellites (and to a lesser extent allozymes) detected a further subdivision between and within the Western Mediterranean and the Tyrrhenian Sea due to present-day processes. A pattern of isolation by distance was not detected in any of the analyzed data set.

CONCLUSIONS

We conclude that the population structure of O. montagui is the result of the interplay of two contrasting forces that act on the species population genetic structure. On one hand, the species semi-terrestrial life style would tend to determine the onset of local differences. On the other hand, these differences are partially counter-balanced by passive movements of migrants via rafting on heaps of dead seagrass leaves across sites by sea surface currents. Approximate Bayesian Computations support dispersal at sea as prevalent over terrestrial regionalism.

Deglacial pulses of deep-ocean silicate into the subtropical North Atlantic Ocean

Growing evidence suggests that the low atmospheric CO2 concentration of the ice ages resulted from enhanced storage of CO2 in the ocean interior, largely as a result of changes in the Southern Ocean. Early in the most recent deglaciation, a reduction in North Atlantic overturning circulation seems to have driven CO2 release from the Southern Ocean, but the mechanism connecting the North Atlantic and the Southern Ocean remains unclear. Biogenic opal export in the low-latitude ocean relies on silicate from the underlying thermocline, the concentration of which is affected by the circulation of the ocean interior. Here we report a record of biogenic opal export from a coastal upwelling system off the coast of northwest Africa that shows pronounced opal maxima during each glacial termination over the past 550,000 years. These opal peaks are consistent with a strong deglacial reduction in the formation of silicate-poor glacial North Atlantic intermediate water (GNAIW). The loss of GNAIW allowed mixing with underlying silicate-rich deep water to increase the silicate supply to the surface ocean. An increase in westerly-wind-driven upwelling in the Southern Ocean in response to the North Atlantic change has been proposed to drive the deglacial rise in atmospheric CO2 (refs 3, 4). However, such a circulation change would have accelerated the formation of Antarctic intermediate water and sub-Antarctic mode water, which today have as little silicate as North Atlantic Deep Water and would have thus maintained low silicate concentrations in the Atlantic thermocline. The deglacial opal maxima reported here suggest an alternative mechanism for the deglacial CO2 release. Just as the reduction in GNAIW led to upward silicate transport, it should also have allowed the downward mixing of warm, low-density surface water to reach into the deep ocean. The resulting decrease in the density of the deep Atlantic relative to the Southern Ocean surface promoted Antarctic overturning, which released CO2 to the atmosphere.

Magic trait electric organ discharge (EOD): Dual function of electric signals promotes speciation in African weakly electric fish

A unique evolutionary specialization of African weakly electric fish (Mormyridae) is their ability to produce and perceive electric signals. Mormyrids use their electric organs discharge (EOD) for electrolocation and electrocommunication. Here we discuss the adaptive significance of the EOD in foraging (electric prey detection) in light of recent results demonstrating that mormyrid fish mate assortatively according to EOD waveform characteristics (electric mate choice). Therefore the EOD as a single trait pleiotropically combines natural divergent selection and reproductive isolation. Consequently we postulate the EOD as a "magic trait" promoting the diversification of African weakly electric fish.

Mitochondrial control region I and microsatellite analyses of endangered Philippine hornbill species (Aves; Bucerotidae) detect gene flow between island populations and genetic diversity loss

Background

The Visayan Tarictic Hornbill (Penelopides panini) and the Walden’s Hornbill (Aceros waldeni) are two threatened hornbill species endemic to the western islands of the Visayas that constitute - between Luzon and Mindanao - the central island group of the Philippine archipelago. In order to evaluate their genetic diversity and to support efforts towards their conservation, we analyzed genetic variation in ~ 600 base pairs (bp) of the mitochondrial control region I and at 12–19 nuclear microsatellite loci. The sampling covered extant populations, still occurring only on two islands (P. panini: Panay and Negros, A. waldeni: only Panay), and it was augmented with museum specimens of extinct populations from neighboring islands. For comparison, their less endangered (= more abundant) sister taxa, the Luzon Tarictic Hornbill (P. manillae) from the Luzon and Polillo Islands and the Writhed Hornbill (A. leucocephalus) from Mindanao Island, were also included in the study. We reconstructed the population history of the two Penelopides species and assessed the genetic population structure of the remaining wild populations in all four species.

Results

Mitochondrial and nuclear data concordantly show a clear genetic separation according to the island of origin in both Penelopides species, but also unravel sporadic over-water movements between islands. We found evidence that deforestation in the last century influenced these migratory events. Both classes of markers and the comparison to museum specimens reveal a genetic diversity loss in both Visayan hornbill species, P. panini and A. waldeni, as compared to their more abundant relatives. This might have been caused by local extinction of genetically differentiated populations together with the dramatic decline in the abundance of the extant populations.

Conclusions

We demonstrated a loss in genetic diversity of P. panini and A. waldeni as compared to their sister taxa P. manillae and A. leucocephalus. Because of the low potential for gene flow and population exchange across islands, saving of the remaining birds of almost extinct local populations - be it in the wild or in captivity - is particularly important to preserve the species’ genetic potential.

Recovery of a nearly extinct Galápagos tortoise despite minimal genetic variation

A species of Galápagos tortoise endemic to Española Island was reduced to just 12 females and three males that have been bred in captivity since 1971 and have produced over 1700 offspring now repatriated to the island. Our molecular genetic analyses of juveniles repatriated to and surviving on the island indicate that none of the tortoises sampled in 1994 had hatched on the island versus 3% in 2004 and 24% in 2007, which demonstrates substantial and increasing reproduction in situ once again. This recovery occurred despite the parental population having an estimated effective population size <8 due to a combination of unequal reproductive success of the breeders and nonrandom mating in captivity. These results provide guidelines for adapting breeding regimes in the parental captive population and decreasing inbreeding in the repatriated population. Using simple morphological data scored on the sampled animals, we also show that a strongly heterogeneous distribution of tortoise sizes on Española Island observed today is due to a large variance in the number of animals included in yearly repatriation events performed in the last 40 years. Our study reveals that, at least in the short run, some endangered species can recover dramatically despite a lack of genetic variation and irregular repatriation efforts.

Prevalence of avian haemosporidian parasites and their host fidelity in the central Philippine islands

We examined the prevalence and host fidelity of avian haemosporidian parasites belonging to the genera Haemoproteus, Leucocytozoon and Plasmodium in the central Philippine islands by sampling 23 bird families (42 species). Using species-specific PCR assays of the mitochondrial cytochrome b gene (471 base pairs, bp), we detected infections in 91 of the 215 screened individuals (42%). We also discriminated between single and multiple infections. Thirty-one infected individuals harbored a single Haemoproteus lineage (14%), 18 a single Leucocytozoon lineage (8%) and 12 a single Plasmodium lineage (6%). Of the 215 screened birds, 30 (14%) presented different types of multiple infections. Intrageneric mixed infections were generally more common (18 Haemoproteus/Haemoproteus, 3 Leucocytozoon/Leucocytozoon, and 1 Plasmodium/Plasmodium) than intergeneric mixed infections (7 Haemoproteus/Leucocytozoon and 1 Haemoproteus/Leucocytozoon/Plasmodium). We recovered 81 unique haemosporidian mitochondrial haplotypes. These clustered in three strongly supported monophyletic clades that correspond to the three haemosporidian genera. Related lineages of Haemoproteus and Leucocytozoon were more likely to derive from the same host family than predicted by chance; however, this was not the case for Plasmodium. These results indicate that switches between host families are more likely to occur in Plasmodium. We conclude that Haemoproteus has undergone a recent diversification across well-supported host-family specific clades, while Leucocytozoon shows a longer association with its host(s). This study supports previous evidence of a higher prevalence and stronger host-family specificity of Haemoproteus and Leucocytozoon compared to Plasmodium.

Mitochondrial DNA regionalism and historical demography in the extant populations of Chirocephalus kerkyrensis (Branchiopoda: Anostraca)

BACKGROUND

Mediterranean temporary water bodies are important reservoirs of biodiversity and host a unique assemblage of diapausing aquatic invertebrates. These environments are currently vanishing because of increasing human pressure. Chirocephalus kerkyrensis is a fairy shrimp typical of temporary water bodies in Mediterranean plain forests and has undergone a substantial decline in number of populations in recent years due to habitat loss. We assessed patterns of genetic connectivity and phylogeographic history in the seven extant populations of the species from Albania, Corfu Is. (Greece), Southern and Central Italy.

METHODOLOGY/PRINCIPAL FINDINGS

We analyzed sequence variation at two mitochondrial DNA genes (Cytochrome Oxidase I and 16s rRNA) in all the known populations of C. kerkyrensis. We used multiple phylogenetic, phylogeographic and coalescence-based approaches to assess connectivity and historical demography across the whole distribution range of the species. C. kerkyrensis is genetically subdivided into three main mitochondrial lineages; two of them are geographically localized (Corfu Is. and Central Italy) and one encompasses a wide geographic area (Albania and Southern Italy). Most of the detected genetic variation (≈81%) is apportioned among the aforementioned lineages.

CONCLUSIONS/SIGNIFICANCE

Multiple analyses of mismatch distributions consistently supported both past demographic and spatial expansions with the former predating the latter; demographic expansions were consistently placed during interglacial warm phases of the Pleistocene while spatial expansions were restricted to cold periods. Coalescence methods revealed a scenario of past isolation with low levels of gene flow in line with what is already known for other co-distributed fairy shrimps and suggest drift as the prevailing force in promoting local divergence. We recommend that these evolutionary trajectories should be taken in proper consideration in any effort aimed at protecting Mediterranean temporary water bodies.

A comparison of two methods to assess audience-induced changes in male mate choice

Multidirectional communicative interactions in social networks can have a profound effect on mate choice behavior. Male Atlantic molly Poecilia mexicana exhibit weaker mating preferences when an audience male is presented. This could be a male strategy to reduce sperm competition risk: interacting more equally with different females may be advantageous because rivals might copy mate choice decisions. In line with this hypothesis, a previous study found males to show a strong audience effect when being observed while exercising mate choice, but not when the rival was presented only before the choice tests. Audience effects on mate choice decisions have been quantified in poeciliid fishes using association preference designs, but it remains unknown if patterns found from measuring association times translate into actual mating behavior. Thus, we created five audience treatments simulating different forms of perceived sperm competition risk and determined focal males’ mating preferences by scoring pre-mating (nipping) and mating behavior (gonopodial thrusting). Nipping did not reflect the pattern that was found when association preferences were measured, while a very similar pattern was uncovered in thrusting behavior. The strongest response was observed when the audience could eavesdrop on the focal male’s behavior. A reduction in the strength of focal males’ preferences was also seen after the rival male had an opportunity to mate with the focal male’s preferred mate. In comparison, the reduction of mating preferences in response to an audience was greater when measuring association times than actual mating behavior. While measuring direct sexual interactions between the focal male and both stimulus females not only the male’s motivational state is reflected but also females’ behavior such as avoidance of male sexual harassment.

Exploiting gene families for phylogenomic analysis of myzostomid transcriptome data

BACKGROUND

In trying to understand the evolutionary relationships of organisms, the current flood of sequence data offers great opportunities, but also reveals new challenges with regard to data quality, the selection of data for subsequent analysis, and the automation of steps that were once done manually for single-gene analyses. Even though genome or transcriptome data is available for representatives of most bilaterian phyla, some enigmatic taxa still have an uncertain position in the animal tree of life. This is especially true for myzostomids, a group of symbiotic (or parasitic) protostomes that are either placed with annelids or flatworms.

METHODOLOGY

Based on similarity criteria, Illumina-based transcriptome sequences of one myzostomid were compared to protein sequences of one additional myzostomid and 29 reference metazoa and clustered into gene families. These families were then used to investigate the phylogenetic position of Myzostomida using different approaches: Alignments of 989 sequence families were concatenated, and the resulting superalignment was analyzed under a Maximum Likelihood criterion. We also used all 1,878 gene trees with at least one myzostomid sequence for a supertree approach: the individual gene trees were computed and then reconciled into a species tree using gene tree parsimony.

CONCLUSIONS

Superalignments require strictly orthologous genes, and both the gene selection and the widely varying amount of data available for different taxa in our dataset may cause anomalous placements and low bootstrap support. In contrast, gene tree parsimony is designed to accommodate multilocus gene families and therefore allows a much more comprehensive data set to be analyzed. Results of this supertree approach showed a well-resolved phylogeny, in which myzostomids were part of the annelid radiation, and major bilaterian taxa were found to be monophyletic.

Hidden diversity in diatoms of Kenyan Lake Naivasha: a genetic approach detects temporal variation

This study provides insights into the morphological and genetic diversity in diatoms occurring in core sediments from tropical lakes in Kenya. We developed a genetic survey technique specific for diatoms utilizing a short region (76-67 bp) of the ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit (rbcL) gene as genetic barcode. Our analyses (i) validated the use of rbcL as a barcoding marker for diatoms, applied to sediment samples, (ii) showed a significant correlation between the results obtained by morphological and molecular data and (iii) indicated temporal variation in diatom assemblages on the inter- and intra-specific level. Diatom assemblages from a short core from Lake Naivasha show a drastic shift over the last 200 years, as littoral species (e.g. Navicula) are replaced by more planktonic ones (e.g. Aulacoseira). Within that same period, we detected periodic changes in the respective frequencies of distinct haplotype groups of Navicula, which coincide with wet and dry periods of Lake Naivasha between 1820 and 1938 AD. Our genetic analyses on historical lake sediments revealed inter- and intra-specific variation in diatoms, which is partially hidden behind single morphotypes. The occurrence of particular genetic lineages is probably correlated with environmental factors.

Multiple paternity in different populations of the sailfin molly, Poecilia latipinna

Rates of multiple paternities were investigated in the sailfin molly (Poecilialatipinna), using eight microsatellite loci. Genotyping was performed for offspring and mothers in 40 broods from four allopatric populations from the south-eastern U.S.A. along a geographic stretch of 1200 km in west-east direction and approximately 200 km from north to south. No significant differences regarding rates of multiple paternities were found between populations despite sample populations stemming from ecologically divergent habitats. Even the most conservative statistical approach revealed a minimum of 70% of the broods being sired by at least two males, with an average of 1.80-2.95 putative fathers per brood. Within broods, one male typically sired far more offspring than would be expected under an assumed equal probability of all detected males siring offspring.

Molecular profiling of diatom assemblages in tropical lake sediments using taxon-specific PCR and Denaturing High-Performance Liquid Chromatography (PCR-DHPLC)

Here we present a protocol to genetically detect diatoms in sediments of the Kenyan tropical Lake Naivasha, based on taxon-specific PCR amplification of short fragments (approximately 100 bp) of the small subunit ribosomal (SSU) gene and subsequent separation of species-specific PCR products by PCR-based denaturing high-performance liquid chromatography (DHPLC). An evaluation of amplicons differing in primer specificity to diatoms and length of the fragments amplified demonstrated that the number of different diatom sequence types detected after cloning of the PCR products critically depended on the specificity of the primers to diatoms and the length of the amplified fragments whereby shorter fragments yielded more species of diatoms. The DHPLC was able to discriminate between very short amplicons based on the sequence difference, even if the fragments were of identical length and if the amplicons differed only in a small number of nucleotides. Generally, the method identified the dominant sequence types from mixed amplifications. A comparison with microscopic analysis of the sediment samples revealed that the sequence types identified in the molecular assessment corresponded well with the most dominant species. In summary, the PCR-based DHPLC protocol offers a fast, reliable and cost-efficient possibility to study DNA from sediments and other environmental samples with unknown organismic content, even for very short DNA fragments.

Full mitochondrial genome sequences of two endemic Philippine hornbill species (Aves: Bucerotidae) provide evidence for pervasive mitochondrial DNA recombination

Background

Although nowaday it is broadly accepted that mitochondrial DNA (mtDNA) may undergo recombination, the frequency of such recombination remains controversial. Its estimation is not straightforward, as recombination under homoplasmy (i.e., among identical mt genomes) is likely to be overlooked. In species with tandem duplications of large mtDNA fragments the detection of recombination can be facilitated, as it can lead to gene conversion among duplicates. Although the mechanisms for concerted evolution in mtDNA are not fully understood yet, recombination rates have been estimated from "one per speciation event" down to 850 years or even "during every replication cycle".

Results

Here we present the first complete mt genome of the avian family Bucerotidae, i.e., that of two Philippine hornbills, Aceros waldeni and Penelopides panini. The mt genomes are characterized by a tandemly duplicated region encompassing part of cytochrome b, 3 tRNAs, NADH6, and the control region. The duplicated fragments are identical to each other except for a short section in domain I and for the length of repeat motifs in domain III of the control region. Due to the heteroplasmy with regard to the number of these repeat motifs, there is some size variation in both genomes; with around 21,657 bp (A. waldeni) and 22,737 bp (P. panini), they significantly exceed the hitherto longest known avian mt genomes, that of the albatrosses. We discovered concerted evolution between the duplicated fragments within individuals. The existence of differences between individuals in coding genes as well as in the control region, which are maintained between duplicates, indicates that recombination apparently occurs frequently, i.e., in every generation.

Conclusions

The homogenised duplicates are interspersed by a short fragment which shows no sign of recombination. We hypothesize that this region corresponds to the so-called Replication Fork Barrier (RFB), which has been described from the chicken mitochondrial genome. As this RFB is supposed to halt replication, it offers a potential mechanistic explanation for frequent recombination in mitochondrial genomes.

A gene duplication/loss event in the ribulose-1,5-bisphosphate-carboxylase/oxygenase (rubisco) small subunit gene family among accessions of Arabidopsis thaliana

Rubisco (ribulose-1,5-bisphosphate carboxylase/oxygenase; EC 4.1.1.39), the most abundant protein in nature, catalyzes the assimilation of CO(2) (worldwide about 10(11) t each year) by carboxylation of ribulose-1,5-bisphosphate. It is a hexadecamer consisting of eight large and eight small subunits. Although the Rubisco large subunit (rbcL) is encoded by a single gene on the multicopy chloroplast genome, the Rubisco small subunits (rbcS) are encoded by a family of nuclear genes. In Arabidopsis thaliana, the rbcS gene family comprises four members, that is, rbcS-1a, rbcS-1b, rbcS-2b, and rbcS-3b. We sequenced all Rubisco genes in 26 worldwide distributed A. thaliana accessions. In three of these accessions, we detected a gene duplication/loss event, where rbcS-1b was lost and substituted by a duplicate of rbcS-2b (called rbcS-2b*). By screening 74 additional accessions using a specific polymerase chain reaction assay, we detected five additional accessions with this duplication/loss event. In summary, we found the gene duplication/loss in 8 of 100 A. thaliana accessions, namely, Bch, Bu, Bur, Cvi, Fei, Lm, Sha, and Sorbo. We sequenced an about 1-kb promoter region for all Rubisco genes as well. This analysis revealed that the gene duplication/loss event was associated with promoter alterations (two insertions of 450 and 850 bp, one deletion of 730 bp) in rbcS-2b and a promoter deletion (2.3 kb) in rbcS-2b* in all eight affected accessions. The substitution of rbcS-1b by a duplicate of rbcS-2b (i.e., rbcS-2b*) might be caused by gene conversion. All four Rubisco genes evolve under purifying selection, as expected for central genes of the highly conserved photosystem of green plants. We inferred a single positive selected site, a tyrosine to aspartic acid substitution at position 72 in rbcS-1b. Exactly the same substitution compromises carboxylase activity in the cyanobacterium Anacystis nidulans. In A. thaliana, this substitution is associated with an inferred recombination. Functional implications of the substitution remain to be evaluated.

Locally adapted fish populations maintain small-scale genetic differentiation despite perturbation by a catastrophic flood event

Background

Local adaptation to divergent environmental conditions can promote population genetic differentiation even in the absence of geographic barriers and hence, lead to speciation. Perturbations by catastrophic events, however, can distort such parapatric ecological speciation processes. Here, we asked whether an exceptionally strong flood led to homogenization of gene pools among locally adapted populations of the Atlantic molly (Poecilia mexicana, Poeciliidae) in the Cueva del Azufre system in southern Mexico, where two strong environmental selection factors (darkness within caves and/or presence of toxic H₂S in sulfidic springs) drive the diversification of P. mexicana. Nine nuclear microsatellites as well as heritable female life history traits (both as a proxy for quantitative genetics and for trait divergence) were used as markers to compare genetic differentiation, genetic diversity, and especially population mixing (immigration and emigration) before and after the flood.

Results

Habitat type (i.e., non-sulfidic surface, sulfidic surface, or sulfidic cave), but not geographic distance was the major predictor of genetic differentiation. Before and after the flood, each habitat type harbored a genetically distinct population. Only a weak signal of individual dislocation among ecologically divergent habitat types was uncovered (with the exception of slightly increased dislocation from the Cueva del Azufre into the sulfidic creek, El Azufre). By contrast, several lines of evidence are indicative of increased flood-induced dislocation within the same habitat type, e.g., between different cave chambers of the Cueva del Azufre.

Conclusions

The virtual absence of individual dislocation among ecologically different habitat types indicates strong natural selection against migrants. Thus, our current study exemplifies that ecological speciation in this and other systems, in which extreme environmental factors drive speciation, may be little affected by temporary perturbations, as adaptations to physico-chemical stressors may directly affect the survival probability in divergent habitat types.

Audience effects in the Atlantic molly (Poecilia mexicana)-prudent male mate choice in response to perceived sperm competition risk?

Background

Multidirectional interactions in social networks can have a profound effect on mate choice behavior; e.g., Poecilia mexicana males show weaker expression of mating preferences when being observed by a rival. This may be an adaptation to reduce sperm competition risk, which arises because commonly preferred female phenotypes will receive attention also from surrounding males, and/or because other males can copy the focal male's mate choice. Do P. mexicana males indeed respond to perceived sperm competition risk? We gave males a choice between two females and repeated the tests under one of the following conditions: (1) an empty transparent cylinder was presented (control); (2) another ("audience") male inside the cylinder observed the focal male throughout the 2^nd part, or (3) the audience male was presented only before the tests, but could not eavesdrop during the actual choice tests (non-specific sperm competition risk treatments); (4) the focal male could see a rival male interact sexually with the previously preferred, or (5) with the non-preferred female before the 2^nd part of the tests (specific sperm competition risk treatments).

Results

The strength of individual male preferences declined slightly also during the control treatment (1). However, this decrease was more than two-fold stronger in audience treatment (2), i.e., with non-specific sperm competition risk including the possibility for visual eavesdropping by the audience male. No audience effect was found in treatments (3) and (5), but a weak effect was also observed when the focal male had seen the previously preferred female sexually interact with a rival male (treatment 4; specific sperm competition risk).

Conclusion

When comparing the two 'non-specific sperm competition risk' treatments, a very strong effect was found only when the audience male could actually observe the focal male during mate choice [treatment (2)]. This suggests that focal males indeed attempt to conceal their mating preferences so as to prevent surrounding males from copying their mate choice. When there is no potential for eavesdropping [treatment (3)], non-specific specific sperm competition risk seems to play a minor or no role. Our results also show that P. mexicana males tend to share their mating effort more equally among females when the resource value of their previously preferred mate decreases after mating with a rival male (perceived specific sperm competition risk), but this effect is comparatively weak.

On the phylogenetic position of Myzostomida: can 77 genes get it wrong?

BACKGROUND

Phylogenomic analyses recently became popular to address questions about deep metazoan phylogeny. Ribosomal proteins (RP) dominate many of these analyses or are, in some cases, the only genes included. Despite initial hopes, phylogenomic analyses including tens to hundreds of genes still fail to robustly place many bilaterian taxa.

RESULTS

Using the phylogenetic position of myzostomids as an example, we show that phylogenies derived from RP genes and mitochondrial genes produce incongruent results. Whereas the former support a position within a clade of platyzoan taxa, mitochondrial data recovers an annelid affinity, which is strongly supported by the gene order data and is congruent with morphology. Using hypothesis testing, our RP data significantly rejects the annelids affinity, whereas a platyzoan relationship is significantly rejected by the mitochondrial data.

CONCLUSION

We conclude (i) that reliance of a set of markers belonging to a single class of macromolecular complexes might bias the analysis, and (ii) that concatenation of all available data might introduce conflicting signal into phylogenetic analyses. We therefore strongly recommend testing for data incongruence in phylogenomic analyses. Furthermore, judging all available data, we consider the annelid affinity hypothesis more plausible than a possible platyzoan affinity for myzostomids, and suspect long branch attraction is influencing the RP data. However, this hypothesis needs further confirmation by future analyses.

Occurrence of Can-SINEs and intron sequence evolution supports robust phylogeny of pinniped carnivores and their terrestrial relatives

Investigating the dog genome we found 178965 introns with a moderate length of 200-1000 bp. A screening of these sequences against 23 different repeat libraries to find insertions of short interspersed elements (SINEs) detected 45276 SINEs. Virtually all of these SINEs (98%) belong to the tRNA-derived Can-SINE family. Can-SINEs arose about 55 million years ago before Carnivora split into two basal groups, the Caniformia (dog-like carnivores) and the Feliformia (cat-like carnivores). Genome comparisons of dog and cat recovered 506 putatively informative SINE loci for caniformian phylogeny. In this study we show how to use such genome information of model organisms to research the phylogeny of related non-model species of interest. Investigating a dataset including representatives of all major caniformian lineages, we analysed 24 randomly chosen loci for 22 taxa. All loci were amplifiable and revealed 17 parsimony-informative SINE insertions. The screening for informative SINE insertions yields a large amount of sequence information, in particular of introns, which contain reliable phylogenetic information as well. A phylogenetic analysis of intron- and SINE sequence data provided a statistically robust phylogeny which is congruent with the absence/presence pattern of our SINE markers. This phylogeny strongly supports a sistergroup relationship of Musteloidea and Pinnipedia. Within Pinnipedia, we see strong support from bootstrapping and the presence of a SINE insertion for a sistergroup relationship of the walrus with the Otariidae.

Electrifying love: electric fish use species-specific discharge for mate recognition

Mate choice is mediated by a range of sensory cues, and assortative mating based on these cues can drive reproductive isolation among diverging populations. A specific feature of mormyrid fish, the electric organ discharge (EOD), is used for electrolocation and intraspecific communication. We hypothesized that the EOD also facilitates assortative mating and ultimately promotes prezygotic reproductive isolation in African weakly electric fishes. Our behavioural experiments using live males as well as EOD playback demonstrated that female mate recognition is influenced by EOD signals and that females are attracted to EOD characteristics of conspecific males. The dual function of the EOD for both foraging and social communication (including mate recognition leading to assortative mating) underlines the importance of electric signal differentiation for the divergence of African weakly electric fishes. Thus, the EOD provides an intriguing mechanism promoting trophic divergence and reproductive isolation between two closely related Campylomormyrus species occurring in sympatry in the lower Congo rapids.