Partitioning of genetic (Rapd) variability among sexes and populations of the Barn Owl (Tyto alba) in Europe

The white Barn Owl subspecies (Tyto alba alba) is found in southern Europe and the reddish-brown subspecies (T a. guttata) in northern and eastern Europe. In central Europe, the two subspecies interbreed producing a large range of phenotypic variants. Because of the different ratios of the subspecies in different geographic regions, we predict that genetic variation should be greater in Switzerland than in Hungary. We tested this hypothesis by measuring genetic variation with the RAPD method. As predicted, the genetic differentiation within a Swiss population of Barn Owls was significantly greater than the variation within a Hungarian population. This suggests that gene flow is greater in central Europe than at the eastern limit of the Barn Owl distribution in Hungary. In both countries genetic variation was more pronounced in females than in males. As in other birds, this is probably because female Barn Owls are less philopatric than males. The number of migrants between Hungary and Switzerland is ca. 1 individual per generation; if calculated separately for the sexes, then 0.525 for males and ca. I for females (Nm values). The difference in the number of migrants between genders again is likely a consequence of higher male philopatry. The sexual differentiation is greater in the Swiss population than in the Hungarian and the genetic substructuring of the populations of the species is substantial. The reason for the considerable population substructuring could be the nonmigratory behavior and socially monogamous pairing of the species, as well as the geographical barriers (Alps) between the populations examined.

A multilevel analysis on the determinants of regional health care expenditure. A note

We apply a multilevel hierarchical model to explore whether anaggregation fallacy exists in estimating the income elasticity of healthexpenditure by ignoring the regional composition of national healthexpenditure figures. We use data for 110 regions in eight OECD countriesin 1997: Australia, Canada, France, Germany, Italy, Spain, Sweden andUnited Kingdom. In doing this we have tried to identify two sources ofrandom variation: within countries and between-countries. Our resultsshow that: 1- Variability between countries amounts to (SD) 0.5433, andjust 13% of that can be attributed to income elasticity and the remaining87% to autonomous health expenditure; 2- Within countries, variabilityamounts to (SD) 1.0249; and 3- The intra-class correlation is 0.5300. Weconclude that we have to take into account the degree of fiscaldecentralisation within countries in estimating income elasticity ofhealth expenditure. Two reasons lie behind this: a) where there isdecentralisation to the regions, policies aimed at emulating diversitytend to increase national health care expenditure; and b) without fiscaldecentralisation, central monitoring of finance tends to reduce regionaldiversity and therefore decrease national health expenditure.

Influence of learning on range expansion and adaptation to novel habitats.

Learning has been postulated to 'drive' evolution, but its influence on adaptive evolution in heterogeneous environments has not been formally examined. We used a spatially explicit individual-based model to study the effect of learning on the expansion and adaptation of a species to a novel habitat. Fitness was mediated by a behavioural trait (resource preference), which in turn was determined by both the genotype and learning. Our findings indicate that learning substantially increases the range of parameters under which the species expands and adapts to the novel habitat, particularly if the two habitats are separated by a sharp ecotone (rather than a gradient). However, for a broad range of parameters, learning reduces the degree of genetically-based local adaptation following the expansion and facilitates maintenance of genetic variation within local populations. Thus, in heterogeneous environments learning may facilitate evolutionary range expansions and maintenance of the potential of local populations to respond to subsequent environmental changes.

Reduced IFNλ4 activity is associated with improved HCV clearance and reduced expression of interferon-stimulated genes

Hepatitis C virus (HCV) infections are the major cause of chronic liver disease, cirrhosis and hepatocellular carcinoma worldwide. Both spontaneous and treatment-induced clearance of HCV depend on genetic variation within the interferon-lambda locus, but until now no clear causal relationship has been established. Here we demonstrate that an amino-acid substitution in the IFNλ4 protein changing a proline at position 70 to a serine (P70S) substantially alters its antiviral activity. Patients harbouring the impaired IFNλ4-S70 variant display lower interferon-stimulated gene (ISG) expression levels, better treatment response rates and better spontaneous clearance rates, compared with patients coding for the fully active IFNλ4-P70 variant. Altogether, these data provide evidence supporting a role for the active IFNλ4 protein as the driver of high hepatic ISG expression as well as the cause of poor HCV clearance.

Oscillatory magnetic tweezers based on ferromagnetic beads and simple coaxial coils

We report the design and validation of simple magnetic tweezers for oscillating ferromagnetic beads in the piconewton and nanometer scales. The system is based on a single pair of coaxial coils operating in two sequential modes: permanent magnetization of the beads through a large and brief pulse of magnetic field and generation of magnetic gradients to produce uniaxial oscillatory forces. By using this two step method, the magnetic moment of the beads remains constant during measurements. Therefore, the applied force can be computed and varies linearly with the driving signal. No feedback control is required to produce well defined force oscillations over a wide bandwidth. The design of the coils was optimized to obtain high magnetic fields (280 mT) and gradients (2 T/m) with high homogeneity (5% variation) within the sample. The magnetic tweezers were implemented in an inverted optical microscope with a videomicroscopy-based multiparticle tracking system. The apparatus was validated with 4.5 ¿m magnetite beads obtaining forces up to ~2 pN and subnanometer resolution. The applicability of the device includes microrheology of biopolymer and cell cytoplasm, molecular mechanics, and mechanotransduction in living cells.

Cooperation between the transcription factors p63 and IRF6 is essential to prevent cleft palate in mice.

Cleft palate is a common congenital disorder that affects up to 1 in 2,500 live human births and results in considerable morbidity to affected individuals and their families. The etiology of cleft palate is complex, with both genetic and environmental factors implicated. Mutations in the transcription factor-encoding genes p63 and interferon regulatory factor 6 (IRF6) have individually been identified as causes of cleft palate; however, a relationship between the key transcription factors p63 and IRF6 has not been determined. Here, we used both mouse models and human primary keratinocytes from patients with cleft palate to demonstrate that IRF6 and p63 interact epistatically during development of the secondary palate. Mice simultaneously carrying a heterozygous deletion of p63 and the Irf6 knockin mutation R84C, which causes cleft palate in humans, displayed ectodermal abnormalities that led to cleft palate. Furthermore, we showed that p63 transactivated IRF6 by binding to an upstream enhancer element; genetic variation within this enhancer element is associated with increased susceptibility to cleft lip. Our findings therefore identify p63 as a key regulatory molecule during palate development and provide a mechanism for the cooperative role of p63 and IRF6 in orofacial development in mice and humans.

The Lithium, Boron and Beryllium content of serpentinized peridotites from ODP Leg 209 (Sites 1272A and 1274A): Implications for lithium and boron budgets of oceanic lithosphere

Despite the key importance of altered oceanic mantle as a repository and carrier of light elements (B, Li, and Be) to depth, its inventory of these elements has hardly been explored and quantified. In order to constrain the systematics and budget of these elements we have studied samples of highly serpentinized (>50%) spinel harzburgite drilled at the Mid-Atlantic Ridge (Fifteen-Twenty Fracture zone, ODP Leg 209, Sites 1272A and 1274A). In-situ analysis by secondary ion mass spectrometry reveals that the B, Li and Be contents of mantle minerals (olivine, orthopyroxene, and clinopyroxene) remain unchanged during serpentinization. B and Li abundances largely correspond to those of unaltered mantle minerals whereas Be is close to the detection limit. The Li contents of clinopyroxene are slightly higher (0.44-2.8 mu g g(-1)) compared to unaltered mantle clinopyroxene, and olivine and clinopyroxene show an inverse Li partitioning compared to literature data. These findings along with textural observations and major element composition obtained from microprobe analysis suggest reaction of the peridotites with a mafic silicate melt before serpentinization. Serpentine minerals are enriched in B (most values between 10 and 100 mu g g(-1)), depleted in Li (most values below I mu g g(-1)) compared to the primary phases, with considerable variation within and between samples. Be is at the detection limit. Analysis of whole rock samples by prompt gamma activation shows that serpentinization tends to increase B (10.4-65.0 mu g g(-1)), H2O and Cl contents and to lower Li contents (0.07-3.37 mu g g(-1)) of peridotites, implying that-contrary to alteration of oceanic crust-B is fractionated from Li and that the B and Li inventory should depend essentially on rock-water ratios. Based on our results and on literature data, we calculate the inventory of B and Li contained in the oceanic lithosphere, and its partitioning between crust and mantle as a function of plate characteristics. We model four cases, an ODP Leg 209-type lithosphere with almost no igneous crust, and a Semail-type lithosphere with a thick igneous crust, both at I and 75 Ma, respectively. The results show that the Li contents of the oceanic lithosphere are highly variable (17-307 kg in a column of I m x I m x thickness of the lithosphere (kg/col)). They are controlled by the primary mantle phases and by altered crust, whereas the B contents (25-904 kg/col) depend entirely on serpentinization. In all cases, large quantities of B reside in the uppermost part of the plate and could hence be easily liberated during slab dehydration. The most prominent input of Li into subduction zones is to be expected from Semail-type lithosphere because most of the Li is stored at shallow levels in the plate. Subducting an ODP Leg 209-type lithosphere would mean only very little Li contribution from the slab. Serpentinized mantle thus plays an important role in B recycling in subduction zones, but it is of lesser importance for Li. (C) 2008 Elsevier Ltd. All rights reserved.

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'.

The adaptive function of melanin-based coloration in the tawny owl (Strix aluco) ; the proximate role of the melanocortin system

Colour polymorphism is common in wild population. One of the main questioning of evolutionary biologists is to understand how different colour variants could have evolved and be maintained in fluctuating environments, a selective process that forces individuals to constantly adapt their strategies in order to survive. This issue is particularly true for traits that are genetically inherited. Natural selection erodes genotypes with lowest fitness (less adapted), reducing in turn global genetic variation within population. In this context, the study of the evolution and maintenance of melanin- based coloration is relevant since inter-individual variation in the deposition of these pigments is common in animal and plant kingdoms and is under strong genetic control. In this thesis, I focus on the specific case of the tawny owl (Strix aluco), a species displaying continuous variation in reddish pheomelanin-based coloration. Interestingly, empirical studies highlighted covariations between melanin-based coloration and important behavioural, physiological and life history traits. Recently, a genetic model pointed out the melanocortin system and their pleiotropic effects as a potential regulator of these covariations. Accordingly, this PhD thesis further investigates colour-specific behavioural, physiological, or life history strategies, while examining the proximate mechanisms underlying these reaction norms. We found that differently coloured tawny owls differently resolve fundamental trade-off between offspring number and quality (Chapter 1), light melanic individuals producing many low- quality offspring and dark, melanic ones producing few high-quality offspring. These reproductive strategies are likely to induce alternative physiological constraints. Indeed, we demonstrated that light melanic individuals produced higher levels of reactive oxygen species (ROS, Chapter 2), but also expressed higher levels of antioxidant (GSH, Chapters 2 & 3). Interestingly, we showed that light melanic breeding females could modulate their POMC prohormone levels according to the environmental conditions, while dark reddish ones produced constant levels of this prohormone {Chapter 4). Finally, we highlighted colour-specific patterns of prohormone convertase 1 (PCI) gene expression (Chapter 5), an enzyme responsible for POMC prohormone processing to ACTH and a- MSH, for instance. Altogether, these results provide strong evidence of colour-specific strategies, light and melanic tawny owls better coping with stressful and relaxed environments, respectively. Variation in melanin-based coloration is likely to be maintained by the heterogeneity of our study area and strong environmental stochasticity within and between years, these process favouring differently coloured tawny owls at different periods of time. From a proximate point of view, this PhD thesis supports the hypothesis that covariations between phenotypic traits and melanin-based coloration stems from the melanocortin system, especially the fundamental role of POMC gene expression and its processing to melanocortin peptides. - Le polymorphisme de couleur est une variation phénotypique très fréquente dans la nature. En biologie évolutive, une des problématiques clés est donc de comprendre comment différent morphes de couleur peuvent être apparus et maintenus au cours du temps dans des environnements aussi variables que les nôtres, surtout que ces fluctuations forcent ces morphes à s'adapter constamment pour assurer leur survie. Cette thématique est particulièrement réelle lorsque les variations phénotypiques sont héréditaires et donc sous forte influence génétique. La sélection naturelle a en effet le pouvoir d'éroder rapidement la variation génétique en éliminant les génotypes mal adaptés. Dans ce sens, l'étude de l'évolution, et de la maintenance de la coloration mélanique est donc tout à fait pertinente car la variation de coloration entre individus est très répandue à travers les règnes animal et végétal et sous forte influence génétique. Dans cette thèse, je me suis concentré sur le cas spécifique de la chouette hulotte (Strix aluco), une espèce présentant une variation continue dans la déposition de pigments pheomélaniques roux. De précédentes études ont déjà montré que cette variation de coloration était associée avec des variations de traits comportementaux, physiologiques ou d'histoire de vie. Récemment, une étude a souligné l'importance du système des mélanocortines et de leurs effets pléiotropes dans la régulation de ces covariations. En conséquence, cette thèse de doctoral a pour but d'étudier un peu plus les stratégies comportementales, physiologiques ou d'histoire de vie spécifiques à chaque morphe de couleur, tout en examinant un peu plus les mécanismes proximaux potentiellement à la base de ces normes de réactions. Nous constatons tout d'abord que les morphes de couleurs étaient associés à différentes stratégies dans la résolution de compromis telle que la production de beaucoup de jeunes ou des jeunes de qualité (Chapitre 1). Les morphes gris (dit peu mélaniques) ont tendance à produire beaucoup de jeunes mains de moindre qualité, alors que les morphes roux (dit fortement mélaniques) produisent moins de jeunes mais de meilleure qualité. Ces stratégies sont susceptibles alors d'induire certaines contraintes physiologiques. Par exemple, nous montrons que les morphes gris produisent plus de dérivés réactifs de l'oxygène (ROS, Chapitre 2), mais aussi plus d'antioxydants (GSH, Chapitres 2 & 3). Nous montrons ensuite que les femelles grises ont une plus grande capacité à moduler leur niveau de POMC prohormone dans le sang en fonction des conditions environnementales, alors que les femelles rousses gardent un niveau constant (Chapitre 4). Finalement, nous démontrons que les patterns d'expression du gène codant pour la prohormone convertase 1 varient chez des jeunes issus de parents gris ou roux (Chapitre 5). Ceci est particulièrement intéressant car cette enzyme permet de scinder la POMC prohormone en plusieurs peptides importants tels que l'ACTH ou l'a-MSH. En conclusion, ces résultats démontrent qu'il y a bel et bien des stratégies évolutives différentes entre les morphes de couleurs, les chouettes hulottes grises et rousses étant respectivement plus adaptés à des environnements stressants ou favorables. L'hétérogénéité de notre zone d'étude et la stochasticité environnementale qui caractérise ses habitats pourraient donc agir comme une source de sélection temporelle, laquelle favoriserait les différents morphes de couleurs à diverses périodes. D'un point de vue plus proximale maintenant, cette thèse de doctorat soutient l'hypothèse que les covariations observées entre la coloration mélanique et des traits phénotypiques importants sont modulées par les effets pléiotropes du système des mélanocortines, et met en avant le rôle prépondérant que pourrait jouer l'expression du gène POMC et sa post traduction en mélanocortines.

Origin and evolution of homologous repeated sequences in the mitochondrial DNA control region of shrews.

The complete mitochondrial DNA (mtDNA) control region was amplified and directly sequenced in two species of shrew, Crocidura russula and Sorex araneus (Insectivora, Mammalia). The general organization is similar to that found in other mammals: a central conserved region surrounded by two more variable domains. However, we have found in shrews the simultaneous presence of arrays of tandem repeats in potential locations where repeats tend to occur separately in other mammalian species. These locations correspond to regions which are associated with a possible interruption of the replication processes, either at the end of the three-stranded D-loop structure or toward the end of the heavy-strand replication. In the left domain the repeated sequences (R1 repeats) are 78 bp long, whereas in the right domain the repeats are 12 bp long in C. russula and 14 bp long in S. araneus (R2 repeats). Variation in the copy number of these repeated sequences results in mtDNA control region length differences. Southern blot analysis indicates that level of heteroplasmy (more than one mtDNA form within an individual) differs between species. A comparative study of the R2 repeats in 12 additional species representing three shrew subfamilies provides useful indications for the understanding of the origin and the evolution of these homologous tandemly repeated sequences. An asymmetry in the distribution of variants within the arrays, as well as the constant occurrence of shorter repeated sequences flanking only one side of the R2 arrays, could be related to asymmetry in the replication of each strand of the mtDNA molecule. The pattern of sequence and length variation within and between species, together with the capability of the arrays to form stable secondary structures, suggests that the dominant mechanism involved in the evolution of these arrays in unidirectional replication slippage.

Extreme population differentiation in vulnerable slavemaking ant with fragmented distribution

Understanding levels of population differentiation and inbreeding are important issues in conservation biology, especially for social Hymenoptera with fragmented and small population sizes. Isolated populations are more vulnerable to genetic loss and extinction than those with extended continuous distributions. However, small populations are not always a consequence of a recent reduction of their habitat. Thus, determining the history of population isolation and current patterns of genetic variation of a species is crucial for its conservation. Rossomyrmex minuchae is a slave-making ant with patchy distribution in South Eastern Spain and is classified as vulnerable by the IUCN. In contrast, the other three known species of the genus are presumed to show more uniform distributions. Here we investigate the genetic diversity and population structure of R. minuchae and compare it with that found in two other species of the genus: R. anatolicus and R. quandratinodum. We conclude that although genetic diversity of R. minuchae is low, there is no evidence of a recent bottleneck, suggesting a gradual and natural fragmentation process. We also show extreme population differentiation at nuclear and mitochondrial markers, and isolation by distance at a local scale. Despite some evidence for inbreeding and low genetic variation within populations, we found almost no diploid males, a finding which contrasts with that expected in inbred Hymenoptera with single locus complementary sex determination. This could mean that sex is determined by another mechanism. We argue that continued low population size means that detrimental effects of inbreeding and low genetic variation are likely in the future. We suggest that a policy of artificial gene flow aimed at increasing within population variation is considered as a management option.

Measurement of the mechanical power of walking by satellite positioning system (GPS).

PURPOSE: This descriptive article illustrates the application of Global Positioning System (GPS) professional receivers in the field of locomotion studies. The technological challenge was to assess the external mechanical work in outdoor walking. METHODS: Five subjects walked five times during 5 min on an athletic track at different imposed stride frequency (from 70-130 steps x min(-1)). A differential GPS system (carrier phase analysis) measured the variation of the position of the trunk at 5 Hz. A portable indirect calorimeter recorded breath-by-breath energy expenditure. RESULTS: For a walking speed of 1.05 +/- 0.11 m x s(-1), the vertical lift of the trunk (43 +/- 14 mm) induced a power of 46.0 +/- 20.4 W. The average speed variation per step (0.15 +/- 0.03 m x s(-1)) produced a kinetic power of 16.9 +/- 7.2 W. As compared with commonly admitted values, the energy exchange (recovery) between the two energy components was low (39.1 +/- 10.0%), which induced an overestimated mechanical power (38.9 +/- 18.3 W or 0.60 W x kg(-1) body mass) and a high net mechanical efficiency (26.9 +/- 5.8%). CONCLUSION: We assumed that the cause of the overestimation was an unwanted oscillation of the GPS antenna. It is concluded that GPS (in phase mode) is now able to record small body movements during human locomotion, and constitutes a promising tool for gait analysis of outdoor unrestrained walking. However, the design of the receiver and the antenna must be adapted to human experiments and a thorough validation study remains to be conducted.

Dugesia sicula (Platyhelminthes, Tricladida): the colonizing success of an asexual Planarian

Background Dugesia sicula is the only species of its genus not presenting an endemic or restricted distribution within the Mediterranean area. It mostly comprises fissiparous populations (asexual reproduction by body division and regeneration), most likely sexually sterile, and characterized by an extremely low genetic diversity interpreted as the consequence of a recent anthropic expansion. However, its fissiparous reproduction can result in an apparent lack of diversity within the species, since genetic variation within individuals can be as large as between them because most individuals within a population are clones. We have estimated haplotype and nucleotide diversity of cytochrome oxidase I within and among individuals along the species distribution of a broad sample of D. sicula, including asexual and the two only sexual populations known today; and predicted its potential distribution based on climatic variables. Our aim was to determine the centre of colonisation origin, whether the populations are recent, and whether the species is expanding. Results The species presents 3 most frequent haplotypes, differing in a maximum of 11 base pairs. As expected from their fissiparous mode of reproduction, in half of all the analysed localities many individuals have multiple heteroplasmic haplotypes. The distribution of haplotypes is not geographically structured; however, the distribution of haplotypes and heteroplasmic populations shows higher diversity in the central Mediterranean region. The potential distribution predicted by climatic variables based modelling shows a preference for coastal areas and fits well with the observed data. Conclusions The distribution and frequency of the most frequent haplotypes and the presence of heteroplasmic individuals allow us to gain an understanding of the recent history of the species, together with previous knowledge on its phylogenetic relationships and age: The species most probably originated in Africa and dispersed through the central Mediterranean. After one or multiple populations became triploid and fissiparous, the species colonized the Mediterranean basin, likely both by its own means and helped by human activities. Its present distribution practically fulfils its potential distribution as modelled with climatic variables. Its prevalence in coastal regions with higher water temperatures predicts a likely future expansion to northern and more interior areas following the increase in temperatures due to climate change.

Supergenes and complex phenotypes.

Understanding the molecular underpinnings of evolutionary adaptations is a central focus of modern evolutionary biology. Recent studies have uncovered a panoply of complex phenotypes, including locally adapted ecotypes and cryptic morphs, divergent social behaviours in birds and insects, as well as alternative metabolic pathways in plants and fungi, that are regulated by clusters of tightly linked loci. These 'supergenes' segregate as stable polymorphisms within or between natural populations and influence ecologically relevant traits. Some supergenes may span entire chromosomes, because selection for reduced recombination between a supergene and a nearby locus providing additional benefits can lead to locus expansions with dynamics similar to those known for sex chromosomes. In addition to allowing for the co-segregation of adaptive variation within species, supergenes may facilitate the spread of complex phenotypes across species boundaries. Application of new genomic methods is likely to lead to the discovery of many additional supergenes in a broad range of organisms and reveal similar genetic architectures for convergently evolved phenotypes.

Hidden diversity in honey bee gut symbionts detected by single-cell genomics.

Microbial communities in animal guts are composed of diverse, specialized bacterial species, but little is known about how gut bacteria diversify to produce genetically and ecologically distinct entities. The gut microbiota of the honey bee, Apis mellifera, presents a useful model, because it consists of a small number of characteristic bacterial species, each showing signs of diversification. Here, we used single-cell genomics to study the variation within two species of the bee gut microbiota: Gilliamella apicola and Snodgrassella alvi. For both species, our analyses revealed extensive variation in intraspecific divergence of protein-coding genes but uniformly high levels of 16S rRNA similarity. In both species, the divergence of 16S rRNA loci appears to have been curtailed by frequent recombination within populations, while other genomic regions have continuously diverged. Furthermore, gene repertoires differ markedly among strains in both species, implying distinct metabolic capabilities. Our results show that, despite minimal divergence at 16S rRNA genes, in situ diversification occurs within gut communities and generates bacterial lineages with distinct ecological niches. Therefore, important dimensions of microbial diversity are not evident from analyses of 16S rRNA, and single cell genomics has potential to elucidate processes of bacterial diversification.