Molecular evidence to suggest the origin of a colonization: Drosophila subobscura in America

Abstract The recent colonization of America by Drosophila subobscura represents a great opportunity for evolutionary biology studies. Knowledge of the populations from which the colonization started would provide an understanding of how genetic composition changed during adaptation to the new environment. Thus, a 793 nucleotide fragment of the Odh (Octanol dehydrogenase) gene was sequenced in 66 chromosomal lines from Barcelona (western Mediterranean) and in 66 from Mt. Parnes (Greece, eastern Mediterranean). No sequence of Odh fragment in Barcelona or Mt. Parnes was identical to any of those previously detected in America. However, an Odh sequence from Barcelona differed in only one nucleotide from another found in American populations. In both cases, the chromosomal lines presented the same inversion: O7, and the Odh gene was located within this inversion. This evidence suggests a possible western Mediterranean origin for the colonization. Finally, the molecular and inversion data indicate that the colonization was not characterized by multiple reintroductions.

Schmidtea mediterranea phylogeography: an old species surviving on a few mediterranean islands?

Schmidtea mediterranea (Platyhelminthes, Tricladida, Continenticola) is found in scattered localities on a few islands and in coastal areas of the western Mediterranean. Although S. mediterranea is the object of many regeneration studies, little is known about its evolutionary history. Its present distribution has been proposed to stem from the fragmentation and migration of the Corsica-Sardinia microplate during the formation of the western Mediterranean basin, which implies an ancient origin for the species. To test this hypothesis, we obtained a large number of samples from across its distribution area. Using known and new molecular markers and, for the first time in planarians, a molecular clock, we analysed the genetic variability and demographic parameters within the species and between its sexual and asexual populations to estimate when they diverged. Results: A total of 2 kb from three markers (COI, CYB and a nuclear intron N13) was amplified from ~200 specimens. Molecular data clustered the studied populations into three groups that correspond to the west, central and southeastern geographical locations of the current distribution of S. mediterranea. Mitochondrial genes show low haplotype and nucleotide diversity within populations but demonstrate higher values when all individuals are considered. The nuclear marker shows higher values of genetic diversity than the mitochondrial genes at the population level, but asexual populations present lower variability than the sexual ones. Neutrality tests are significant for some populations. Phylogenetic and dating analyses show the three groups to be monophyletic, with the west group being the basal group. The time when the diversification of the species occurred is between ~20 and ~4 mya, although the asexual nature of the western populations could have affected the dating analyses. Conclusions: S. mediterranea is an old species that is sparsely distributed in a harsh habitat, which is probably the consequence of the migration of the Corsica-Sardinia block. This species probably adapted to temperate climates in the middle of a changing Mediterranean climate that eventually became dry and hot. These data also suggest that in the mainland localities of Europe and Africa, sexual individuals of S. mediterranea are being replaced by asexual individuals that are either conspecific or are from other species that are better adapted to the Mediterranean climate.

Ecological and historical drivers of diversification in the fly genus Chiastocheta Pokorny.

Coevolution is among the main forces shaping the biodiversity on Earth. In Eurasia, one of the best-known plant-insect interactions showing highly coevolved features involves the fly genus Chiastocheta and its host-plant Trollius. Although this system has been widely studied from an ecological point of view, the phylogenetic relationships and biogeographic history of the flies have remained little investigated. In this integrative study, we aim to test the monophyly of the five Chiastocheta eco-morphological groups, defined by Pellmyr in 1992, by inferring a mitochondrial phylogeny. We further apply a new approach to assess the effect of (i) different molecular substitution rates and (ii) phylogenetic uncertainty on the inference of the spatio-temporal evolution of the group. From a taxonomic point of view, we demonstrate that only two of Pellmyr's groups (rotundiventris and dentifera) are phylogenetically supported, the other species appearing para- or polyphyletic. We also identify the position of C. lophota, which was not included in previous surveys. From a spatio-temporal perspective, we show that the genus arose during the Pliocene in Europe. Our results also indicate that at least four large-scale dispersal events are required to explain the current distribution of Chiastocheta. Moreover, each dispersal to or from Asia is associated with a host-shift and seems to correspond to an increase in speciation rates. Finally, we highlight the correlation between diversification and climatic fluctuations, which indicate that the cycles of global cooling over the last million years had an influence on the radiation of the group.

Is the Gibraltar strait a barrier to gene flow for the bat Myotis myotis (Chiroptera: Vespertilionidae)?

Because of their role in limiting gene flow, geographical barriers like mountains or seas often coincide with intraspecific genetic discontinuities. Although the Strait of Gibraltar represents such a potential barrier for both plants and animals, few studies have been conducted on its impact on gene flow. Here we test this effect on a bat species (Myotis myotis) which is apparently distributed on both sides of the strait. Six colonies of 20 Myotis myotis each were sampled in southern Spain and northern Morocco along a linear transect of 1350 km. Results based on six nuclear microsatellite loci reveal no significant population structure within regions, but a complete isolation between bats sampled on each side of the strait. Variability at 600 bp of a mitochondrial gene (cytochrome b) confirms the existence of two genetically distinct and perfectly segregating clades, which diverged several million years ago. Despite the narrowness of the Gibraltar Strait (14 km), these molecular data suggest that neither males, nor females from either region have ever reproduced on the opposite side of the strait. Comparisons of molecular divergence with bats from a closely related species (M. blythii) suggest that the North African clade is possibly a distinct taxon warranting full species rank. We provisionally refer to it as Myotis cf punicus Felten 1977, but a definitive systematic understanding of the whole Mouse-eared bat species complex awaits further genetic sampling, especially in the Eastern Mediterranean areas.

Molecular signaling in zebrafish development and the vertebrate phylotypic period

During development vertebrate embryos pass through a stage where their morphology is most conserved between species, the phylotypic period (approximately the pharyngula). To explain the resistance to evolutionary changes of this period, one hypothesis suggests that it is characterized by a high level of interactions. Based on this hypothesis, we examined protein-protein interactions, signal transduction cascades and miRNAs over the course of zebrafish development, and the conservation of expression of these genes in mouse development. We also investigated the characteristics of genes highly expressed before or during the presumed phylotypic period. We show that while there is a high diversity of interactions during the phylotypic period (protein-DNA, RNA-RNA, cell-cell, and between tissues), which is well conserved with mouse, there is no clear difference with later, more morphologically divergent, stages. We propose that the phylotypic period may rather be the expression at the morphological level of strong conservation of molecular processes earlier in development.

Comprehensive molecular profiling of lung adenocarcinoma.

Adenocarcinoma of the lung is the leading cause of cancer death worldwide. Here we report molecular profiling of 230 resected lung adenocarcinomas using messenger RNA, microRNA and DNA sequencing integrated with copy number, methylation and proteomic analyses. High rates of somatic mutation were seen (mean 8.9 mutations per megabase). Eighteen genes were statistically significantly mutated, including RIT1 activating mutations and newly described loss-of-function MGA mutations which are mutually exclusive with focal MYC amplification. EGFR mutations were more frequent in female patients, whereas mutations in RBM10 were more common in males. Aberrations in NF1, MET, ERBB2 and RIT1 occurred in 13% of cases and were enriched in samples otherwise lacking an activated oncogene, suggesting a driver role for these events in certain tumours. DNA and mRNA sequence from the same tumour highlighted splicing alterations driven by somatic genomic changes, including exon 14 skipping in MET mRNA in 4% of cases. MAPK and PI(3)K pathway activity, when measured at the protein level, was explained by known mutations in only a fraction of cases, suggesting additional, unexplained mechanisms of pathway activation. These data establish a foundation for classification and further investigations of lung adenocarcinoma molecular pathogenesis.

A single point mutation in the pore region of the epithelial Na+ channel changes ion selectivity by modifying molecular sieving.

The epithelial Na+ channel (ENaC) belongs to a new class of channel proteins called the ENaC/DEG superfamily involved in epithelial Na+ transport, mechanotransduction, and neurotransmission. The role of ENaC in Na+ homeostasis and in the control of blood pressure has been demonstrated recently by the identification of mutations in ENaC beta and gamma subunits causing hypertension. The function of ENaC in Na+ reabsorption depends critically on its ability to discriminate between Na+ and other ions like K+ or Ca2+. ENaC is virtually impermeant to K+ ions, and the molecular basis for its high ionic selectivity is largely unknown. We have identified a conserved Ser residue in the second transmembrane domain of the ENaC alpha subunit (alphaS589), which when mutated allows larger ions such as K+, Rb+, Cs+, and divalent cations to pass through the channel. The relative ion permeability of each of the alphaS589 mutants is related inversely to the ionic radius of the permeant ion, indicating that alphaS589 mutations increase the molecular cutoff of the channel by modifying the pore geometry at the selectivity filter. Proper geometry of the pore is required to tightly accommodate Na+ and Li+ ions and to exclude larger cations. We provide evidence that ENaC discriminates between cations mainly on the basis of their size and the energy of dehydration.

The hourglass and the early conservation models--co-existing patterns of developmental constraints in vertebrates.

Developmental constraints have been postulated to limit the space of feasible phenotypes and thus shape animal evolution. These constraints have been suggested to be the strongest during either early or mid-embryogenesis, which corresponds to the early conservation model or the hourglass model, respectively. Conflicting results have been reported, but in recent studies of animal transcriptomes the hourglass model has been favored. Studies usually report descriptive statistics calculated for all genes over all developmental time points. This introduces dependencies between the sets of compared genes and may lead to biased results. Here we overcome this problem using an alternative modular analysis. We used the Iterative Signature Algorithm to identify distinct modules of genes co-expressed specifically in consecutive stages of zebrafish development. We then performed a detailed comparison of several gene properties between modules, allowing for a less biased and more powerful analysis. Notably, our analysis corroborated the hourglass pattern at the regulatory level, with sequences of regulatory regions being most conserved for genes expressed in mid-development but not at the level of gene sequence, age, or expression, in contrast to some previous studies. The early conservation model was supported with gene duplication and birth that were the most rare for genes expressed in early development. Finally, for all gene properties, we observed the least conservation for genes expressed in late development or adult, consistent with both models. Overall, with the modular approach, we showed that different levels of molecular evolution follow different patterns of developmental constraints. Thus both models are valid, but with respect to different genomic features.

Additional data for nuclear DNA give new insights into the phylogenetic position of Sorex granarius within the Sorex araneus group.

Many species contain genetic lineages that are phylogenetically intermixed with those of other species. In the Sorex araneus group, previous results based on mtDNA and Y chromosome sequence data showed an incongruent position of Sorex granarius within this group. In this study, we explored the relationship between species within the S. araneus group, aiming to resolve the particular position of S. granarius. In this context, we sequenced a total of 2447 base pairs (bp) of X-linked and nuclear genes from 47 individuals of the S. araneus group. The same taxa were also analyzed within a Bayesian framework with nine autosomal microsatellites. These analyses revealed that all markers apart from mtDNA showed similar patterns, suggesting that the problematic position of S. granarius is best explained by an incongruent behavior by mtDNA. Given their close phylogenetic relationship and their close geographic distribution, the most likely explanation for this pattern is past mtDNA introgression from S. araneus race Carlit to S. granarius.

Genetic isolation of insular populations of the Maghrebian bat, Myotis punicus, in the Mediterranean Basin

We investigate the population genetic structure of the Maghrebian bat, Myotis punicus, between the mainland and islands to assess the island colonization pattern and current gene flow between nearby islands and within the mainland. Location North Africa and the Mediterranean islands of Corsica and Sardinia. Methods We sequenced part of the control region (HVII) of 79 bats across 11 colonies. The phylogeographical pattern was assessed by analysing molecular diversity indices, examining differentiation among populations and estimating divergence time. In addition, we genotyped 182 bats across 10 colonies at seven microsatellite loci. We used analysis of molecular variance and a Bayesian approach to infer nuclear population structure. Finally, we estimated sex-specific dispersal between Corsica and Sardinia. Results Mitochondrial analyses indicated that colonies between Corsica, Sardinia and North Africa are highly differentiated. Within islands there was no difference between colonies, while at the continental level Moroccan and Tunisian populations were highly differentiated. Analyses with seven microsatellite loci showed a similar pattern. The sole difference was the lack of nuclear differentiation between populations in North Africa, suggesting a male-biased dispersal over the continental area. The divergence time of Sardinian and Corsican populations was estimated to date back to the early and mid-Pleistocene. Main conclusions Island colonization by the Maghrebian bats seems to have occurred in a stepping-stone manner and certainly pre-dated human colonization. Currently, open water seems to prevent exchange of bats between the two islands, despite their ability to fly and the narrowness of the strait of Bonifacio. Corsican and Sardinian populations are thus currently isolated from any continental gene pool and must therefore be considered as different evolutionarily significant units (ESU).

Integrative molecular bioinformatics study of human adrenocortical tumors : microRNA, tissue-specific target prediction, and pathway analysis.

MicroRNAs (miRs) are involved in the pathogenesis of several neoplasms; however, there are no data on their expression patterns and possible roles in adrenocortical tumors. Our objective was to study adrenocortical tumors by an integrative bioinformatics analysis involving miR and transcriptomics profiling, pathway analysis, and a novel, tissue-specific miR target prediction approach. Thirty-six tissue samples including normal adrenocortical tissues, benign adenomas, and adrenocortical carcinomas (ACC) were studied by simultaneous miR and mRNA profiling. A novel data-processing software was used to identify all predicted miR-mRNA interactions retrieved from PicTar, TargetScan, and miRBase. Tissue-specific target prediction was achieved by filtering out mRNAs with undetectable expression and searching for mRNA targets with inverse expression alterations as their regulatory miRs. Target sets and significant microarray data were subjected to Ingenuity Pathway Analysis. Six miRs with significantly different expression were found. miR-184 and miR-503 showed significantly higher, whereas miR-511 and miR-214 showed significantly lower expression in ACCs than in other groups. Expression of miR-210 was significantly lower in cortisol-secreting adenomas than in ACCs. By calculating the difference between dCT(miR-511) and dCT(miR-503) (delta cycle threshold), ACCs could be distinguished from benign adenomas with high sensitivity and specificity. Pathway analysis revealed the possible involvement of G2/M checkpoint damage in ACC pathogenesis. To our knowledge, this is the first report describing miR expression patterns and pathway analysis in sporadic adrenocortical tumors. miR biomarkers may be helpful for the diagnosis of adrenocortical malignancy. This tissue-specific target prediction approach may be used in other tumors too.

Molecular epidemiology of clonal diploids: a quick overview and a short DIY (do it yourself) notice.

In this short review we report the basic notions needed for understanding the population genetics of clonal diploids. We focus on the consequences of clonality on the distribution of genetic diversity within individuals, between individuals and between populations. We then summarise how to detect clonality in mainly sexual populations, conversely, how to detect sexuality in mainly clonal populations and also how genetic differentiation between populations is affected by clonality in diploids. This information is then used for building recipes on how to analyse and interpret genetic polymorphism data in molecular epidemiology studies of clonal diploids.

Molecular clock is involved in predictive circadian adjustment of renal function.

Renal excretion of water and major electrolytes exhibits a significant circadian rhythm. This functional periodicity is believed to result, at least in part, from circadian changes in secretion/reabsorption capacities of the distal nephron and collecting ducts. Here, we studied the molecular mechanisms underlying circadian rhythms in the distal nephron segments, i.e., distal convoluted tubule (DCT) and connecting tubule (CNT) and the cortical collecting duct (CCD). Temporal expression analysis performed on microdissected mouse DCT/CNT or CCD revealed a marked circadian rhythmicity in the expression of a large number of genes crucially involved in various homeostatic functions of the kidney. This analysis also revealed that both DCT/CNT and CCD possess an intrinsic circadian timing system characterized by robust oscillations in the expression of circadian core clock genes (clock, bma11, npas2, per, cry, nr1d1) and clock-controlled Par bZip transcriptional factors dbp, hlf, and tef. The clock knockout mice or mice devoid of dbp/hlf/tef (triple knockout) exhibit significant changes in renal expression of several key regulators of water or sodium balance (vasopressin V2 receptor, aquaporin-2, aquaporin-4, alphaENaC). Functionally, the loss of clock leads to a complex phenotype characterized by partial diabetes insipidus, dysregulation of sodium excretion rhythms, and a significant decrease in blood pressure. Collectively, this study uncovers a major role of molecular clock in renal function.

RAE-1, a novel PHR binding protein, is required for axon termination and synapse formation in Caenorhabditis elegans.

Previous studies in Caenorhabditis elegans showed that RPM-1 (Regulator of Presynaptic Morphology-1) regulates axon termination and synapse formation. To understand the mechanism of how rpm-1 functions, we have used mass spectrometry to identify RPM-1 binding proteins, and have identified RAE-1 (RNA Export protein-1) as an evolutionarily conserved binding partner. We define a RAE-1 binding region in RPM-1, and show that this binding interaction is conserved and also occurs between Rae1 and the human ortholog of RPM-1 called Pam (protein associated with Myc). rae-1 loss of function causes similar axon and synapse defects, and synergizes genetically with two other RPM-1 binding proteins, GLO-4 and FSN-1. Further, we show that RAE-1 colocalizes with RPM-1 in neurons, and that rae-1 functions downstream of rpm-1. These studies establish a novel postmitotic function for rae-1 in neuronal development.

Plastid genome characterisation in Brassica and Brassicaceae using a new set of nine SSRs.

We report a new set of nine primer pairs specifically developed for amplification of Brassica plastid SSR markers. The wide utility of these markers is demonstrated for haplotype identification and detection of polymorphism in B. napus, B. nigra, B. oleracea, B. rapa and in related genera Arabidopsis, Camelina, Raphanus and Sinapis. Eleven gene regions (ndhB-rps7 spacer, rbcL-accD spacer, rpl16 intron, rps16 intron, atpB-rbcL spacer, trnE-trnT spacer, trnL intron, trnL-trnF spacer, trnM-atpE spacer, trnR-rpoC2 spacer, ycf3-psaA spacer) were sequenced from a range of Brassica and related genera for SSR detection and primer design. Other sequences were obtained from GenBank/EMBL. Eight out of nine selected SSR loci showed polymorphism when amplified using the new primers and a combined analysis detected variation within and between Brassica species, with the number of alleles detected per locus ranging from 5 (loci MF-6, MF-1) to 11 (locus MF-7). The combined SSR data were used in a neighbour-joining analysis (SMM, D (DM) distances) to group the samples based on the presence and absence of alleles. The analysis was generally able to separate plastid types into taxon-specific groups. Multi-allelic haplotypes were plotted onto the neighbour joining tree. A total number of 28 haplotypes were detected and these differentiated 22 of the 41 accessions screened from all other accessions. None of these haplotypes was shared by more than one species and some were not characteristic of their predicted type. We interpret our results with respect to taxon differentiation, hybridisation and introgression patterns relating to the 'Triangle of U'.