The genetic basis of floral mechanical isolation between two hummingbird-pollinated Neotropical understorey herbs.

Floral divergence can contribute to reproductive isolation among plant lineages, and thus provides an opportunity to study the genetics of speciation, including the number, effect size, mode of action and interactions of quantitative trait loci (QTL). Moreover, flowers represent suites of functionally interrelated traits, but it is unclear to what extent the phenotypic integration of the flower is underlain by a shared genetic architecture, which could facilitate or constrain correlated evolution of floral traits. Here, we examine the genetic architecture of floral morphological traits involved in an evolutionary switch from bill to forehead pollen placement between two species of hummingbird-pollinated Neotropical understorey herbs that are reproductively isolated by these floral differences. For the majority of traits, we find multiple QTL of relatively small effect spread throughout the genome. We also find substantial colocalization and alignment of effects of QTL underlying different floral traits that function together to promote outcrossing and reduce heterospecific pollen transfer. Our results are consistent with adaptive pleiotropy or linkage of many co-adapted genes, either of which could have facilitated a response to correlated selection and helped to stabilize divergent phenotypes in the face of low levels of hybridization. Moreover, our results indicate that floral mechanical isolation can be consistent with an infinitesimal model of adaptation.

Multilocus phylogeny, species age and biogeography of the Lesser Antillean anoles.

Lesser Antillean anoles provide classic examples of island radiations. A detailed knowledge of their phylogeny and biogeography, in particular how the age of species relate to the ages of their respective islands and the age of their radiation, is essential to elucidate the tempo and mechanisms of these radiations. We conduct a large-scale phylogenetic and phylogeographic investigation of the Lesser Antillean anoles using multiple genetic markers and comprehensive geographic sampling of most species. The multilocus phylogeny gives the first well-supported reconstruction of the interspecific relationships, and the densely sampled phylogeography reveals a highly dynamic system, driven by overseas dispersal, with several alternative post-dispersal colonisation trajectories. These radiations currently occupy both the outer-older (Eocene to Miocene), and the inner-younger (<8mybp), Lesser Antillean arcs. The origin of these radiations corresponds with the age of the ancient outer arc. However, the ages of extant species (compatible with the age of other small terrestrial amniotes) are much younger, about the age of the emergence of the younger arc, or less. The difference between the age of the radiation and the age of the extant species suggests substantial species turnover on older arc islands, most likely through competitive replacement. Although extant anoles are extremely speciose, this may represent only a fraction of their biodiversity over time. While paraphyly enables us to infer several recent colonization events, the absence of the younger arc islands and extant species at the earlier and middle stages of the radiation, does not allow the earlier inter-island colonization to be reliably inferred. Reproductive isolation in allopatry takes a very considerable time (in excess of 8my) and sympatry appears to occur only late in the radiation. The resolved multilocus phylogeny, and relative species age, raise difficulties for some earlier hypotheses regarding size evolution, and provide no evidence for within-island speciation.

Divergence maintained by climatic selection despite recurrent gene flow: a case study of Castanopsis carlesii (Fagaceae).

Local adaptation to different environments has the potential to maintain divergence between populations despite recurrent gene flow and is an important driver for generating biological diversity. In this study, we investigate the role of adaptation in the maintenance of two parapatric varieties of a forest tree. We used sequence variation of chloroplastic DNA and restriction site-associated DNA to investigate the genetic structure of two varieties of Castanopsis carlesii in subtropical China and relate it to climatic variation. We used niche reconstruction methods to investigate niche differentiation between the two varieties and to estimate the past distribution of this species. A deep divergence was observed between the two varieties, but evidence of introgression and genetic admixture was detected in two phenotypically and geographically intermediate populations. Niche reconstruction suggests that the distribution of the two varieties was disjunct during periods of global cooling and that the two varieties occupy significantly different niches. The genetic structure was mainly driven by environmental factors, and 13 outlier loci under divergent selection were correlated with climatic variation. These results suggest that the two varieties evolved in allopatry and came back into secondary contact after the last glacial maximum and that they are an evolutionary example of divergence maintained by climatic selection despite recurrent gene flow.

An assembly and alignment-free method of phylogeny reconstruction from next-generation sequencing data.

BACKGROUND: Next-generation sequencing technologies are rapidly generating whole-genome datasets for an increasing number of organisms. However, phylogenetic reconstruction of genomic data remains difficult because de novo assembly for non-model genomes and multi-genome alignment are challenging. RESULTS: To greatly simplify the analysis, we present an Assembly and Alignment-Free (AAF) method ( https://sourceforge.net/projects/aaf-phylogeny ) that constructs phylogenies directly from unassembled genome sequence data, bypassing both genome assembly and alignment. Using mathematical calculations, models of sequence evolution, and simulated sequencing of published genomes, we address both evolutionary and sampling issues caused by direct reconstruction, including homoplasy, sequencing errors, and incomplete sequencing coverage. From these results, we calculate the statistical properties of the pairwise distances between genomes, allowing us to optimize parameter selection and perform bootstrapping. As a test case with real data, we successfully reconstructed the phylogeny of 12 mammals using raw sequencing reads. We also applied AAF to 21 tropical tree genome datasets with low coverage to demonstrate its effectiveness on non-model organisms. CONCLUSION: Our AAF method opens up phylogenomics for species without an appropriate reference genome or high sequence coverage, and rapidly creates a phylogenetic framework for further analysis of genome structure and diversity among non-model organisms.

Widespread parallel population adaptation to climate variation across a radiation: implications for adaptation to climate change.

Global warming will impact species in a number of ways, and it is important to know the extent to which natural populations can adapt to anthropogenic climate change by natural selection. Parallel microevolution within separate species can demonstrate natural selection, but several studies of homoplasy have not yet revealed examples of widespread parallel evolution in a generic radiation. Taking into account primary phylogeographic divisions, we investigate numerous quantitative traits (size, shape, scalation, colour pattern and hue) in anole radiations from the mountainous Lesser Antillean islands. Adaptation to climatic differences can lead to very pronounced differences between spatially close populations with all studied traits showing some evidence of parallel evolution. Traits from shape, scalation, pattern and hue (particularly the latter) show widespread evolutionary parallels within these species in response to altitudinal climate variation greater than extreme anthropogenic climate change predicted for 2080. This gives strong evidence of the ability to adapt to climate variation by natural selection throughout this radiation. As anoles can evolve very rapidly, it suggests anthropogenic climate change is likely to be less of a conservation threat than other factors, such as habitat loss and invasive species, in this, Lesser Antillean, biodiversity hot spot.

Multilocus phylogeography of the common lizard Zootoca vivipara at the Ibero-Pyrenean suture zone reveals lowland barriers and high-elevation introgression.

BACKGROUND: The geographic distribution of evolutionary lineages and the patterns of gene flow upon secondary contact provide insight into the process of divergence and speciation. We explore the evolutionary history of the common lizard Zootoca vivipara (= Lacerta vivipara) in the Iberian Peninsula and test the role of the Pyrenees and the Cantabrian Mountains in restricting gene flow and driving lineage isolation and divergence. We also assess patterns of introgression among lineages upon secondary contact, and test for the role of high-elevation trans-mountain colonisations in explaining spatial patterns of genetic diversity. We use mtDNA sequence data and genome-wide AFLP loci to reconstruct phylogenetic relationships among lineages, and measure genetic structure. RESULTS: The main genetic split in mtDNA corresponds generally to the French and Spanish sides of the Pyrenees as previously reported, in contrast to genome-wide AFLP data, which show a major division between NW Spain and the rest. Both types of markers support the existence of four distinct and geographically congruent genetic groups, which are consistent with major topographic barriers. Both datasets reveal the presence of three independent contact zones between lineages in the Pyrenean region, one in the Basque lowlands, one in the low-elevation mountains of the western Pyrenees, and one in the French side of the central Pyrenees. The latter shows genetic evidence of a recent, high-altitude trans-Pyrenean incursion from Spain into France. CONCLUSIONS: The distribution and age of major lineages is consistent with a Pleistocene origin and a role for both the Pyrenees and the Cantabrian Mountains in driving isolation and differentiation of Z. vivipara lineages at large geographic scales. However, mountain ranges are not always effective barriers to dispersal, and have not prevented a recent high-elevation trans-Pyrenean incursion that has led to asymmetrical introgression among divergent lineages. Cytonuclear discordance in patterns of genetic structure and introgression at contact zones suggests selection may be involved at various scales. Suture zones are important areas for the study of lineage formation and speciation, and our results show that biogeographic barriers can yield markedly different phylogeographic patterns in different vertebrate and invertebrate taxa.

Optimization of de novo transcriptome assembly from next-generation sequencing data.

Transcriptome analysis has important applications in many biological fields. However, assembling a transcriptome without a known reference remains a challenging task requiring algorithmic improvements. We present two methods for substantially improving transcriptome de novo assembly. The first method relies on the observation that the use of a single k-mer length by current de novo assemblers is suboptimal to assemble transcriptomes where the sequence coverage of transcripts is highly heterogeneous. We present the Multiple-k method in which various k-mer lengths are used for de novo transcriptome assembly. We demonstrate its good performance by assembling de novo a published next-generation transcriptome sequence data set of Aedes aegypti, using the existing genome to check the accuracy of our method. The second method relies on the use of a reference proteome to improve the de novo assembly. We developed the Scaffolding using Translation Mapping (STM) method that uses mapping against the closest available reference proteome for scaffolding contigs that map onto the same protein. In a controlled experiment using simulated data, we show that the STM method considerably improves the assembly, with few errors. We applied these two methods to assemble the transcriptome of the non-model catfish Loricaria gr. cataphracta. Using the Multiple-k and STM methods, the assembly increases in contiguity and in gene identification, showing that our methods clearly improve quality and can be widely used. The new methods were used to assemble successfully the transcripts of the core set of genes regulating tooth development in vertebrates, while classic de novo assembly failed.

Restricted gene flow within and between rapidly diverging Neotropical plant species.

Speciation involves the evolution of traits and genetic differences that contribute to reproductive isolation and the cessation of gene flow, and studying closely related species and divergent populations gives insight into how these phenomena proceed. Here, we document patterns of gene flow within and between two members of a rapid Neotropical species radiation, Costus pulverulentus and Costus scaber (Costaceae). These species co-occur in the tropical rainforest and share pollinators, but are reproductively isolated by a series of prezygotic barriers, some of which show evidence of reinforcement at sympatric sites. Here, we genotype microsatellite markers in plants from eight sites that span the geographical range of both species, including four sympatric sites. We also genotype putative hybrids found at two sympatric sites. We find high levels of genetic isolation among populations within each species and low but detectable levels of introgression between species at sympatric sites. Putative hybrids identified by morphology are consistent with F1 or more advanced hybrids. Our results highlight the effectiveness of prezygotic isolating mechanisms at maintaining species boundaries in young radiations and provide empirical data on levels of gene flow consistent with reinforcement.

Genetic tests for ecological and allopatric speciation in anoles on an island archipelago.

From Darwin's study of the Galapagos and Wallace's study of Indonesia, islands have played an important role in evolutionary investigations, and radiations within archipelagos are readily interpreted as supporting the conventional view of allopatric speciation. Even during the ongoing paradigm shift towards other modes of speciation, island radiations, such as the Lesser Antillean anoles, are thought to exemplify this process. Geological and molecular phylogenetic evidence show that, in this archipelago, Martinique anoles provide several examples of secondary contact of island species. Four precursor island species, with up to 8 mybp divergence, met when their islands coalesced to form the current island of Martinique. Moreover, adjacent anole populations also show marked adaptation to distinct habitat zonation, allowing both allopatric and ecological speciation to be tested in this system. We take advantage of this opportunity of replicated island coalescence and independent ecological adaptation to carry out an extensive population genetic study of hypervariable neutral nuclear markers to show that even after these very substantial periods of spatial isolation these putative allospecies show less reproductive isolation than conspecific populations in adjacent habitats in all three cases of subsequent island coalescence. The degree of genetic interchange shows that while there is always a significant genetic signature of past allopatry, and this may be quite strong if the selection regime allows, there is no case of complete allopatric speciation, in spite of the strong primae facie case for it. Importantly there is greater genetic isolation across the xeric/rainforest ecotone than is associated with any secondary contact. This rejects the development of reproductive isolation in allopatric divergence, but supports the potential for ecological speciation, even though full speciation has not been achieved in this case. It also explains the paucity of anole species in the Lesser Antilles compared to the Greater Antilles.

Quantitative traits and mode of speciation in Martinique anoles.

We investigate extensive quantitative trait variation (dewlap hue, colour pattern, dorsum hue, body proportions and scalation) in the Martinique anole across eight transects representing nascent parapatric ecological speciation, nascent allopatric speciation and allopatric divergence without sufficient genetic structure to suggest speciation. Quantitative trait divergence can be extremely large between adjacent sets of populations, but with one exception that this is associated with difference in habitat rather than past allopatry. Nascent ecological speciation shows the greatest level of quantitative trait divergence across all character sets including those implicated in natural, as well as sexual selection. The sole example of nascent allopatric speciation is associated with fairly strong quantitative trait divergence among most character sets, but not the set most implicated in natural (rather than sexual) selection. The role of sexual selection in ecological speciation is discussed, both in terms of female choice with assortative mating and male-male competition with condition-dependant sexual signals.

Convergent evolution of kin-based sociality in a lizard.

Studies of social birds and mammals have produced extensive theory regarding the formation and dynamics of kin-based social groups in vertebrates. However, comparing kin dynamics in birds and mammals to social reptiles provides the opportunity to identify selective factors that promote independent origins of kin sociality across vertebrates. We combined a 5-year mark-recapture study with a DNA microsatellite analysis of relatedness in a social lizard (Xantusia vigilis) to examine the formation and stability of kin groups. We found that these lizards are highly sedentary and that groups often form through the delayed dispersal of offspring. Groups containing juveniles had higher relatedness than adult-only groups, as juveniles were commonly found in aggregations with at least one parent and/or sibling. Groups containing nuclear family members were more stable than groups of less-related lizards, as predicted by social theory. We conclude that X. vigilis aggregations conform to patterns of kin sociality observed in avian and mammalian systems and represent an example of convergent evolution in social systems. We suggest that kin-based sociality in this and other lizards may be a by-product of viviparity, which can promote delayed juvenile dispersal by allowing prolonged interaction between a neonate and its mother.

The roles of allopatric divergence and natural selection in quantitative trait variation across a secondary contact zone in the lizard Anolis roquet.

Populations of the Caribbean lizard, Anolis roquet, are thought to have experienced long periods of allopatry before recent secondary contact. To elucidate the effects of past allopatry on population divergence in A. roquet, we surveyed parallel transects across a secondary contact zone in northeastern Martinique. We used diagnostic molecular mitochondrial DNA markers to test fine-scale association of mitochondrial DNA lineage and geological region, multivariate statistical techniques to explore quantitative trait pattern, and cline fitting techniques to model trait variation across the zone of secondary contact. We found that lineages were strongly associated with geological regions along both transects, but quantitative trait patterns were remarkably different. Patterns of morphological and mitochondrial DNA variation were consistent with a strong barrier to gene flow on the coast, whereas there were no indications of barriers to gene flow in the transitional forest. Hence, the coastal populations behaved as would be predicted by an allopatric model of divergence in this complex, while those in the transitional forest did not, despite the close proximity of the transects and their shared geological history. Patterns of geographical variation in this species complex, together with environmental data, suggest that on balance, selection regimes on either side of the secondary contact zone in the transitional forest may be more convergent, while those either side of the secondary contact zone on the coast are more divergent. Hence, the evolutionary consequences of allopatry may be strongly influenced by local natural selection regimes.

Microsatellite data show evidence for male-biased dispersal in the Caribbean lizard Anolis roquet.

Dispersal is a key component of an organism's life history and differences in dispersal between sexes appear to be widespread in vertebrates. However, most predictions of sex-biased dispersal have been based on observations of social structure in birds and mammals and more data are needed on other taxa to test whether these predictions apply in other organisms. Caribbean anole lizards are important model organisms in various biological disciplines, including evolutionary biology. However, very little is known about their dispersal strategies despite the importance of dispersal for population structure and dynamics. Here we use nine microsatellite markers to assess signatures of sex-biased dispersal on two spatial sampling scales in Anolis roquet, an anole endemic to the island of Martinique. Significantly higher gene diversity (H(S)) and lower mean assignment value (mAIC) was found in males on the larger spatial sampling scale. Significant heterozygote deficit (F(IS)), lower population differentiation (F(ST)), mAIC and variance of assignment index (vAIC) was found in males on the smaller spatial scale. The observation of male biased dispersal conform with expectations based on the polygynous mating system of Anolis roquet, and contributes to an explanation of the contrasting patterns of genetic structure between maternal and biparental markers that have been reported previously in this, and other anoline, species.

Reconstructing phylogeny from reduced-representation genome sequencing data without assembly or alignment.

Reduced-representation genome sequencing such as RADseq aids the analysis of genomes by reducing the quantity of data, thereby lowering both sequencing costs and computational burdens. RADseq was initially designed for studying genetic variation across genomes at the population level, but has also proved to be suitable for interspecific phylogeny reconstruction. RADseq data pose challenges for standard phylogenomic methods, however, due to incomplete coverage of the genome and large amounts of missing data. Alignment-free methods are both efficient and accurate for phylogenetic reconstructions with whole genomes and are especially practical for nonmodel organisms; nonetheless, alignment-free methods have not been applied with reduced genome sequencing data. Here, we test a full-genome assembly- and alignment-free method, AAF, in application to RADseq data and propose two procedures for reads selection to remove reads from restriction sites that were not found in taxa being compared. We validate these methods using both simulations and real data sets. Reads selection improved the accuracy of phylogenetic construction in every simulated scenario and the two real data sets, making AAF as good or better than a comparable alignment-based method, even though AAF had much lower computational burdens. We also investigated the sources of missing data in RADseq and their effects on phylogeny reconstruction using AAF. The AAF pipeline modified for RADseq or other reduced-representation sequencing data, phyloRAD, is available on github (https://github.com/fanhuan/phyloRAD).

Models of density-dependent genic selection and a new rock-paper-scissors social system.

We describe new ESS models of density regulation driven by genic selection to explain the cyclical dynamics of a social system that exhibits a rock-paper-scissors (RPS) set of three alternative strategies. We tracked changes in morph frequency and fitness of Lacerta vivipara and found conspicuous RPS cycles. Morphs of Uta and Lacerta exhibited parallel survival-performance trade-offs. Frequency cycles in both species of lizards are driven by genic selection. In Lacerta, frequency of each allele in adult cohorts had significant impacts on juvenile recruitment, similar to mutualistic, altruistic, and antagonistic relations of RPS alleles in Uta. We constructed evolutionarily stable strategy (ESS) models in which adults impact juvenile recruitment as a function of self versus nonself color recognition. ESS models suggest that the rapid 4-year RPS cycles exhibited by Lacerta are not possible unless three factors are present: behaviors evolve that discriminate self versus nonself morphs at higher rates than random, self- versus non-self-recognition contributes to density regulation, and context-dependent mate choice evolves in females, which choose sire genotypes to enhance progeny survival. We suggest genic selection coupled to density regulation is widespread and thus fundamental to theories of social system evolution as well as theories of population regulation in diverse animal taxa.

The utility of DNA metabarcoding for studying the response of arthropod diversity and composition to land-use change in the tropics.

Metabarcoding potentially offers a rapid and cheap method of monitoring biodiversity, but real-world applications are few. We investigated its utility in studying patterns of litter arthropod diversity and composition in the tropics. We collected litter arthropods from 35 matched forest-plantation sites across Xishuangbanna, southwestern China. A new primer combination and the MiSeq platform were used to amplify and sequence a wide variety of litter arthropods using simulated and real-world communities. Quality filtered reads were clustered into 3,624 MOTUs at ≥97% similarity and the taxonomy of each MOTU was predicted. We compared diversity and compositional differences between forests and plantations (rubber and tea) for all MOTUs and for eight arthropod groups. We obtained ~100% detection rate after in silico sequencing six mock communities with known arthropod composition. Ordination showed that rubber, tea and forest communities formed distinct clusters. α-diversity declined significantly between forests and adjacent plantations for more arthropod groups in rubber than tea, and diversity of order Orthoptera increased significantly in tea. Turnover was higher in forests than plantations, but patterns differed among groups. Metabarcoding is useful for quantifying diversity patterns of arthropods under different land-uses and the MiSeq platform is effective for arthropod metabarcoding in the tropics.

Development of the uterine shell glands during the preovulatory and early gestation periods in oviparous and viviparous Lacerta vivipara.

The evolutionary process leading to the emergence of viviparity in Squamata consists of lengthening the period of egg retention in utero coupled with marked reduction in the thickness of the eggshell. We used light microscopy and scanning electron microscopy to study uterine structure during the reproductive cycle of oviparous and viviparous females of the reproductively bimodal Lacerta vivipara. We compared the structure of the uterine shell glands, which secrete components of the eggshell, during preovulatory and early gestation phases of the reproductive cycle and also compared histochemistry of the eggshells. The uterine glands of both reproductive forms undergo considerable growth within a period of a few weeks during folliculogenesis and vitellogenesis preceding ovulation. The majority of the proteinaceous fibers of the shell membrane are secreted early in embryonic development and the uterine glands regress shortly thereafter. This supports previous observations indicating that, in Squamata, secretion of the shell membrane occurs very rapidly after ovulation. The most striking differences between reproductive modes were larger uterine glands at late vitellogenesis in oviparous females, 101 microm compared to 60 microm in viviparous females, and greater thickness of the shell membrane during early gestation in oviparous females (52-73 microm) compared to viviparous females (4-8 microm). Our intraspecific comparison supports the conclusions of previous studies that, prior to ovulation, the uterine glandular layer is less developed in viviparous than in oviparous species, and that this is the main factor accounting for differences in the thickness of the shell membrane of the two reproductive forms of squamates.

Variation in eggshell characteristics and in intrauterine egg retention between two oviparous clades of the lizard Lacerta vivipara: insight into the oviparity-viviparity continuum in squamates.

The concept of the oviparity-viviparity continuum refers to the wide range in the length of intrauterine egg retention and, hence, in the stage of embryonic development at oviposition existing in squamates. The evolutionary process underlying this continuum may involve not only a lengthening of egg retention in utero, but also a marked reduction in the thickness of the eggshell. The idea that there may exist a negative correlation between the developmental stage reached by the embryo at oviposition and the eggshell thickness within squamates, although supported by the comparison of oviparous vs. viviparous species, has seldom been evaluated by comparing eggshell thickness of oviparous forms with different lengths of intrauterine egg retention. Eggs of two distinct oviparous clades of the lizard Lacerta vivipara were compared. The eggs laid by females from Slovenian and Italian populations have thicker eggshells, contain embryos on average less developed at the time of oviposition, and require a longer incubation period before hatching than the eggs laid by females from French oviparous populations. Our data and several other examples available from the literature support the idea that the lengthening of intrauterine retention of eggs and the shortening of the subsequent external incubation of eggs are associated with reduction in the thickness of the eggshell, at least in some lineages of oviparous squamates. The current hypotheses that may account for this correlation are presented and a few restrictions and refinements to those hypotheses are discussed. In particular, other changes, such as increased vascularization of the oviduct and of the extraembryonic membranes, may play the same role as the decrease of eggshell thickness in facilitating prolonged intrauterine egg retention in squamates. Future studies should also consider the hypothesis that the length of intrauterine retention might directly depend on the extent of maternal-fetal chemical communication through the eggshell barrier.

Extraembryonic membrane development in a reproductively bimodal lizard, Lacerta (Zootoca) vivipara.

Reproductive mode has been remarkably labile among squamate reptiles and the evolutionary transition from oviparity to viviparity commonly has been accompanied by a shift in the pattern of embryonic nutrition. Structural specializations for placental transfer of nutrients during intrauterine gestation are highly diverse and many features of the extraembryonic membranes of viviparous species differ markedly from those of oviparous species. However, because of a high degree of evolutionary divergence between the species used for comparisons it is likely that the observed differences arose secondarily to the evolution of viviparity. We studied development of the extraembryonic membranes and placentation in the reproductively bimodal lizard Lacerta vivipara because the influence of reproductive mode on the structural/functional relationship between mothers and embryos can best be understood by studying the most recent evolutionary events. Lecithotrophic viviparity has evolved recently within this species and, although populations with different reproductive modes are allopatric, oviparous and viviparous forms interbreed in the laboratory and share many life history characteristics. In contrast to prior comparisons between oviparous and viviparous species, we found no differences in ontogeny or structure of the extraembryonic membranes between populations with different reproductive modes within L. vivipara. However, we did confirm conclusions from previous studies that the tertiary envelope of the egg, the eggshell, is much reduced in the viviparous population. These conclusions support a widely accepted model for the evolution of squamate placentation. We also found support for work published nearly 80 years ago that the pattern of development of the yolk sac of L. vivipara is unusual and that a function of a unique structure of squamate development, the yolk cleft, is hematopoiesis. The structure of the yolk sac splanchnopleure of L. vivipara is inconsistent with a commonly accepted model for amniote yolk sac function and we suggest that a long standing hypothesis that cells from the yolk cleft participate in yolk digestion requires further study.

How does conversion of natural tropical rainforest ecosystems affect soil bacterial and fungal communities in the Nile river watershed of Uganda?

Uganda's forests are globally important for their conservation values but are under pressure from increasing human population and consumption. In this study, we examine how conversion of natural forest affects soil bacterial and fungal communities. Comparisons in paired natural forest and human-converted sites among four locations indicated that natural forest soils consistently had higher pH, organic carbon, nitrogen, and calcium, although variation among sites was large. Despite these differences, no effect on the diversity of dominant taxa for either bacterial or fungal communities was detected, using polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE). Composition of fungal communities did generally appear different in converted sites, but surprisingly, we did not observe a consistent pattern among sites. The spatial distribution of some taxa and community composition was associated with soil pH, organic carbon, phosphorus and sodium, suggesting that changes in soil communities were nuanced and require more robust metagenomic methods to understand the various components of the community. Given the close geographic proximity of the paired sampling sites, the similarity between natural and converted sites might be due to continued dispersal between treatments. Fungal communities showed greater environmental differentiation than bacterial communities, particularly according to soil pH. We detected biotic homogenization in converted ecosystems and substantial contribution of ß-diversity to total diversity, indicating considerable geographic structure in soil biota in these forest communities. Overall, our results suggest that soil microbial communities are relatively resilient to forest conversion and despite a substantial and consistent change in the soil environment, the effects of conversion differed widely among sites. The substantial difference in soil chemistry, with generally lower nutrient quantity in converted sites, does bring into question, how long this resilience will last.