'Next generation' biogeography: towards understanding the drivers of species diversification and persistence

The drivers of species diversification and persistence are of great interest to current biogeography, especially in those global biodiversity hotspots' harbouring most of Earth's animal and plant life. Classical multispecies biogeographical work has yielded fascinating insights into broad-scale patterns of diversification, and DNA-based intraspecific phylogeographical studies have started to complement this picture at much finer temporal and spatial scales. The advent of novel next-generation sequencing (NGS) technologies provides the opportunity to greatly scale up the numbers of individuals, populations and species sampled, potentially merging intraspecific and interspecific approaches to biogeographical inference. Here, we outline these prospects and issues by using the example of an undisputed hotspot, the Cape of southern Africa. We outline the current state of knowledge on the biogeography of species diversification within the Cape, review the literature for phylogeographical evidence of its likely drivers and mechanisms, and suggest possible ways forward based on NGS approaches. We demonstrate the potential of these methods and current bioinformatic issues with the help of restriction-site-associated DNA (RAD) sequencing data for three highly divergent species of the Restionaceae, an important plant radiation in the Cape. A thorough understanding of the mechanisms that facilitate species diversification and persistence in spatially structured, species-rich environments will require the adoption of novel genomic and bioinformatic tools in biogeographical studies.

Genetic variability among populations of Lutzomyia (Psathyromyia) shannoni (Dyar 1929) (Diptera: Psychodidae: Phlebotominae) in Colombia

Polyacrylamide gel electrophoresis was used to elucidate genetic variation at 13 isozyme loci among forest populations of Lutzomyia shannoni from three widely separated locations in Colombia: Palambí (Nariño Department), Cimitarra (Santander Department) and Chinácota (Norte de Santander Department). These samples were compared with a laboratory colony originating from the Magdalena Valley in Central Colombia. The mean heterozygosity ranged from 16 to 22%, with 2.1 to 2.6 alleles detected per locus. Nei's genetic distances among populations were low, ranging from 0.011 to 0.049. The estimated number of migrants (Nm=3.8) based on Wright's F-Statistic, F ST, indicated low levels of gene flow among Lu. shannoni forest populations. This low level of migration indicates that the spread of stomatitis virus occurs via infected host, not by infected insect. In the colony sample of 79 individuals, the Gpi locus was homozygotic (0.62/0.62) in all females and heterozygotic (0.62/0.72) in all males. Although this phenomenon is probably a consequence of colonization, it indicates that Gpi is linked to a sex determining locus.

Genetic diversity in Brazilian populations of Aedes albopictus

Random amplified polymorphic DNA (RAPD) analysis technique was undertaken in Aedes albopictus populations from three states in Brazil, Rio de Janeiro (RJ), Minas Gerais (MG) and Pernambuco (PE), to estimate the level of genetic variability and levels of genetic exchange between populations. Allele and genotype frequencies were measured on 47 RAPD loci. Average observed heterozigosity (Ho) ranged from 0.282 in MG to 0.355 in Casa Forte (PE) population. Genetic distances estimates indicated that RJ and MG were more genetically similar than populations from PE. Genetic variation observed in local Brazilian populations was attributed to genetic drift associated with restricted gene flow in recently established populations.

Low relatedness among cooperatively breeding workers of the greenhead ant Rhytidoponera metallica

The greenhead ant Rhytidoponera metallica has long been recognized as posing a potential challenge to kin selection theory because it has large queenless colonies where apparently many of the morphological workers are mated and reproducing. However this species has never been studied genetically and important elements of its breeding system and kin structure remain uncertain. We used microsatellite markers to measure the relatedness among nestmates unravel the fine-scale population genetic structure and infer the breeding system of R. metallica. The genetic relatedness among worker nestmates is very low but significantly greater than zero (r = 0.082 +/- 0.015) which demonstrates that nests contain many distantly related breeders. The inbreeding coefficient is very close to and not significantly different from zero indicating random mating and lack of microgeographic genetic differentiation. On average. closely located nests are not more similar genetically than distant nests which is surprising as new colonies form by budding and female dispersal is restricted. Lack of inbreeding and absence of population viscosity indicates high gene flow mediated by males. Overall the genetic pattern detected in R. metallica suggests that a high number of moderately related workers mate with unrelated males from distant nests. This breeding system results in the lowest relatedness among nestmates reported for social insect species where breeders and helpers are not morphologically differentiated. [References: 69]

Phenotypic divergence but not genetic distance predicts assortative mating among species of a cichlid fish radiation.

The hypothesis of ecological divergence giving rise to premating isolation in the face of gene flow is controversial. However, this may be an important mechanism to explain the rapid multiplication of species during adaptive radiation following the colonization of a new environment when geographical barriers to gene flow are largely absent but underutilized niche space is abundant. Using cichlid fish, we tested the prediction of ecological speciation that the strength of premating isolation among species is predicted by phenotypic rather than genetic distance. We conducted mate choice experiments between three closely related, sympatric species of a recent radiation in Lake Mweru (Zambia/DRC) that differ in habitat use and phenotype, and a distantly related population from Lake Bangweulu that resembles one of the species in Lake Mweru. We found significant assortative mating among all closely related, sympatric species that differed phenotypically, but none between the distantly related allopatric populations of more similar phenotype. Phenotypic distance between species was a good predictor of the strength of premating isolation, suggesting that assortative mating can evolve rapidly in association with ecological divergence during adaptive radiation. Our data also reveals that distantly related allopatric populations that have not diverged phenotypically, may hybridize when coming into secondary contact, e.g. upon river capture because of diversion of drainage systems.

Isoenzymatic analysis of four Anopheles (Kerteszia) bellator Dyar & Knab (Diptera: Culicidae) populations

Anopheles bellator is a small silvatic bromelia-breeding mosquito and is a primary human malaria vector species in Southern Brazil. The bromelia-breeding habitat of the species should accompany the Atlantic forest coastal distribution, where bromeliads are abundant. Nonetheless, records on An. bellator collections show a gap in the species geographical distribution. An. bellator has been recorded in Southern Brazil and in the Brazilian states of Bahia and Paraíba. It appers again in the island of Trinidad, in Trinidad and Tobago. The aim of this work was to measure gene flow between different populations of An. bellator collected in the northern and southern extremes of the geographic distribution of this species. Mosquitoes were captured in forest borders in Santa Catarina, São Paulo, and Bahia states in Brazil and in the island of Trinidad in Republic of Trinidad and Tobago. Genetic distances varied between 0.076 and 0.680, based on enzymatic profiles from 11 distinct isoenzymes. Results indicate the existence of low-level gene flow between Brazilian populations of An. bellator, and a gene flow was even lower between the Brazilian and the Trinidad populations. This finding lead us to hypothesize that An. bellator did not spread along the coast, but reached northeastern areas through inland routes.

Isoenzimatic analysis of four Anopheles (Kerteszia) cruzii (Diptera: Culicidae) populations of Brazil

Anopheles cruzii is a small sylvatic mosquito and primary human Plasmodium vector in Southern Brazil. The distribution of this bromeliad-breeding mosquito follows the Atlantic forest coastal distribution, where bromeliads are abundant. Morphological, genetic, and molecular polymorphisms among different populations have been reported and it has recently been suggested that An. cruzii is a complex of cryptic species. The aim of this work is to analyze the gene flow between different populations of An. cruzii collected in four localities within the geographic distribution range of the species, and to examine if An. cruzii is a complex of cryptic species. The genetic distances show that populations of the states of Santa Catarina, São Paulo, and Rio de Janeiro are genetically closer (0.032 to 0.083) than populations of Bahia (0.364 to 0.853) based on profiles from 10 distinct isoenzyme loci. The Fst was lower (0.077) when the Bahia population was excluded than when it was included (0.300) in the analyses. The inferred number of migrants per generation was 2.99 individuals among populations from the states of Santa Catarina, São Paulo, and Rio de Janeiro and 0.58 migrants per generation among all populations. Results suggest that An. cruzii is a complex of species and that the specimes of state of Bahia can be considered as belonging to a species that is distinct from other three closely-related populations studied.

Genetic lineages in the yellow fever mosquito Aedes (Stegomyia) aegypti (Diptera: Culicidae) from Peru

The yellow fever mosquito Aedes aegypti was introduced in Peru in 1852 and was considered to be eradicated in 1958. In 2001, Ae. aegypti had been recorded in 15 out of 24 Peruvian Departments. Peru has great ecological differences between the east and west sides of Andes. Because of this, we consider that Ae. aegypti populations of both east and west sides can have a genetically distinct population structure. In this study we examined genetic variability and genealogical relationships among three Ae. aegypti Peruvian populations: Lima, Piura (west Andes), and Iquitos (east Andes) using a fragment of the ND4 gene of the mitochondrial genome. Three haplotypes were detected among 55 samples. Lima and Iquitos showed the same haplotype (Haplotype I), whereas Piura has two haplotypes (Haplotype II and III). Haplotype II is four mutational steps apart from Haplotype I, while Haplotype III is 13 mutational steps apart from Haplotype I in the network. The analysis of molecular variation showed that mostly of the detected genetic variation occurs at interpopulational level. The significant value phist suggests that Piura population is structured in relation to Lima and Iquitos populations and the gene flow of the ND4 is restricted in Piura when compared to Lima and Iquitos. Genetic relationship between haplotype I and haplotype II suggests introduction of the same mtDNA lineage into those localities. However the existence of a genetically distant haplotype III also suggests introduction of at least two Ae. aegypti lineages in Peru.

Patterns of phenotypic variation reveal substantial differentiation in sexual dimorphism of three Psammodromus (Squamata, Lacertidae) species

The Spanish sand racer (Psammodromus hispanicus) has been recently split into three distinct species: P. hispanicus, P. edwardsianus, and P. occidentalis. Some morphological differences have been reported but there is as yet no description allowing unambiguous identification of the three species. Here, we describe differentiation in body measurements, scalation traits, and colour traits as well as in the degree of sexual dimorphism. Our results show that P. edwardsianus can be easily distinguished by the presence of a supralabial scale below the subocular scale, which is absent in the other two species. Psammodromus hispanicus and P. occidentalis can be distinguished by the number of femoral pores, throat scales and ocelli, and the relative width of the anal scale. The degree of sexual size dimorphism and sexual colour dimorphism substantially differs among species, suggesting that different scenarios of sexual and natural selection may exist for each species. Moreover, sexually selected traits (nuptial colouration, ocelli, and femoral pores) significantly differ among species, suggesting that visual and chemical communication may also differ among species. Such differences could prevent reproduction and gene flow at secondary contact zones, potentially reinforcing isolation and speciation within this group of lizards.

Evolution of life history traits in a natural population of the shrew Crocidura russula

Abstract Life history traits encompass all the decisions concerning fitness an individual is faced with during his life. The study of these traits is crucial to understand the factors shaping the biology of living organisms. Up until now, most of the information on the evolution of life history traits comes from laboratory studies. While these studies are interesting to test the effect of specific parameters, their conclusions are difficult to extrapolate to natural populations. Investigating the evolution of life history traits in natural populations is of great interest. This may be tricky because it requires information on reproduction, survival and morphology of individuals. Mark-recapture methods allow most of this information to be obtained. However, when direct observations of a species are not possible due to its ecology, indirect methods must be used to infer lifetime reproductive success. In this case, molecular markers are particularly helpful in assessing the genetic relationships between individuals and allow the construction of a pedigree. This thesis focuses on a natural population of a small insectivorous mammal, the greater white-toothed shrew, Crocidura russula. Because of its hidden lifestyle, the two complementary techniques mentioned above were combined to gather information on this population. The data were used to explore diverse aspects of evolutionary biology. We demonstrated that the high genetic variance displayed by the species was not maintained by its mating system because this shrew was less monogamous than previously thought. The large genetic diversity was most likely promoted by gene flow from the neighborhood. Dispersal was thus a central topic in this thesis. We showed that dispersal was not driven by inbreeding avoidance. In addition, we did not find any inbreeding depression in the population. Dispersal was promoted by a high number of vacant territories in the population for both sexes, meaning that territory acquisition played an important role in driving dispersal. Moreover, dispersal propensity was shown to have a genetic basis and, once achieved, to have no effect on individual fitness. Body mass was found to be a life history trait strongly influenced by sexual and viability selection in both sexes. Larger individuals had higher access to reproduction through territory acquisition and defense than lighter ones. By contrast, intermediate size individuals were favored by viability selection presumably because of ecological constraints and metabolic costs. Finally, we demonstrated that the majority of the life history traits in our shrew population has the potential to evolve because they maintained substantial amounts of additive genetic variance. Nonetheless, life history traits had no significant heritability due to their high level of nonadditive or environmental variance. Résumé Les traits d'histoire de vie comprennent toutes les décisions auxquelles un individu est confronté au cours de sa vie et qui concernent sa valeur adaptative. L'étude de ces traits est cruciale pour comprendre les facteurs qui façonnent la biologie des êtres vivants. Jusqu'à ce jour, la majorité des informations sur l'évolution des traits d'histoire de vie provient d'études réalisées en laboratoire. Alors que ces études sont intéressantes pour tester l'effet de paramètres spécifiques, leurs conclusions sont difficilement extrapolables aux populations naturelles. Il est particulièrement intéressant d'étudier l'évolution des traits d'histoire de vie dans des populations naturelles. Toutefois, ces études peuvent se révéler difficiles parce qu'elles requièrent des informations sur la reproduction, la survie et la morphologie des individus. Des méthodes de marquage-recapture permettent d'obtenir ces informations. Cependant, lorsque l'écologie de l'espèce rend les obervations directes impossibles, des méthodes indirectes doivent être utilisées pour obtenir le succès reproducteur des individus. Dans ce cas, les marqueurs moléculaires sont particulièrement utiles pour évaluer les relations génétiques entre individus et permettre la construction d'un pedigree. Cette thèse porte sur une population naturelle d'un petit mammifère insectivore, la musaraigne musette, Crocidura russula. Parce que cette espèce présente un mode de vie souterrain, les deux techniques complémentaires mentionnées ci-dessus ont été combinées pour acquérir les informations nécessaires. Les données ont été utilisées pour explorer divers aspects de biologie evolutive. Nous avons montré que la grande quantité de variance génétique trouvée chez cette espèce n'est pas maintenue par son système d'appariement. Celle-ci s'est en effet avérée être moins monogame que ce qui était admis jusqu'ici. Sa grande diversité génétique est plutôt entretenue par le flux de gènes provenant du voisinage. La dispersion a donc été un sujet phare dans cette thèse. Nous avons montré qu'elle n'est pas provoquée par un évitement de la consanguinité et nous n'avons pas trouvé de dépression de consanguité dans notre population. L'acquisition d'un territoire joue par contre un rôle important dans la dispersion. En outre, la dispersion possède une base génétique chez cette espèce. De plus, une fois qu'ils ont dispersé, les individus n'ont pas une valeur adaptative differente d'individus philopatriques. Le poids s'est avéré être un trait d'histoire de vie fortement influencé par la sélection sexuelle et de viabilité chez les deux sexes. Les gros individus ont accès à la reproduction parce qu'ils acquièrent et défendent un territoire plus facilement que les plus légers. Au contraire, les individus de taille intermédiaire sont favorisés par la sélection de viabilité, certainement à cause de contraintes écologiques et de coûts métaboliques. Finalement, nous avons montré que la majorité des traits d'histoire de vie dans notre population a le potentiel d'évoluer parce qu'elle maintient des quantités considérables de variance génétique additive. Néanmoins, l'héritabilité de ces traits d'histoire de vie n'est pas significative à cause de la grande quantité de variance non-additive ou environmentale associée à ces traits.

Higher genetic variation estimated by microsatellites compared to isoenzyme markers in Aedes aegypti from Rio de Janeiro

Aedes aegypti populations from five districts in Rio de Janeiro were analyzed using five microsatellites and six isoenzyme markers, to assess the amount of variation and patterns of gene flow at local levels. Microsatellite loci were polymorphic enough to detect genetic differentiation of populations collected at small geographic scales (e.g. within a city). Ae. aegypti populations were highly differentiated as well in the city center as in the outskirt. Thus, dengue virus propagation by mosquitoes could be as efficient in the urban area as in the outskirt of Rio de Janeiro, the main entry point of dengue in Brazil.

A population genetics study of Anopheles darlingi (Diptera: Culicidae) from Colombia based on random amplified polymorphic DNA-polymerase chain reaction and amplified fragment lenght polymorphism markers

The genetic variation and population structure of three populations of Anopheles darlingi from Colombia were studied using random amplified polymorphic markers (RAPDs) and amplified fragment length polymorphism markers (AFLPs). Six RAPD primers produced 46 polymorphic fragments, while two AFLP primer combinations produced 197 polymorphic fragments from 71 DNA samples. Both of the evaluated genetic markers showed the presence of gene flow, suggesting that Colombian An. darlingi populations are in panmixia. Average genetic diversity, estimated from observed heterozygosity, was 0.374 (RAPD) and 0.309 (AFLP). RAPD and AFLP markers showed little evidence of geographic separation between eastern and western populations; however, the F ST values showed high gene flow between the two western populations (RAPD: F ST = 0.029; Nm: 8.5; AFLP: F ST = 0.051; Nm: 4.7). According to molecular variance analysis (AMOVA), the genetic distance between populations was significant (RAPD:phiST = 0.084; AFLP:phiST = 0.229, P < 0.001). The F ST distances and AMOVAs using AFLP loci support the differentiation of the Guyana biogeographic province population from those of the Chocó-Magdalena. In this last region, Chocó and Córdoba populations showed the highest genetic flow.

Population genetic structure of the major malaria vector Anopheles darlingi (Diptera: Culicidae) from the Brazilian Amazon, using microsatellite markers

The population genetic structure of Anopheles darlingi, the major human malaria vector in the Neotropics, was examined using seven microsatellite loci from nine localities in central and western Amazonian Brazil. High levels of genetic variability were detected (5-25 alleles per locus; H E = 0.519-0.949). There was deviation from Hardy-Weinberg Equilibrium for 59.79% of the tests due to heterozygote deficits, while the analysis of linkage disequilibrium was significant for only two of 189 (1.05%) tests, most likely caused by null alleles. Genetic differentiation (F ST = 0.001-0.095; Nm = 4.7-363.8) indicates that gene flow is extensive among locations < 152 km apart (with two exceptions) and reduced, but not absent, at a larger geographic scale. Genetic and geographic distances were significantly correlated (R² = 0.893, P < 0.0002), supporting the isolation by distance (IBD) model. The overall estimate of Ne was 202.4 individuals under the linkage disequilibrium model, and 8 under the heterozygote excess model. Analysis of molecular variance showed that nearly all variation (~ 94%) was within sample locations. The UPGMA phenogram clustered the samples geographically, with one branch including 5/6 of the state of Amazonas localities and the other branch the Acre, Rondônia, and remaining Amazonas localities. Taken together, these data suggest little genetic structure for An. darlingi from central and western Amazonian Brazil. These findings also imply that the IBD model explains nearly all of the differentiation detected. In practical terms, populations of An. darlingi at distances < 152 km should respond similarly to vector control measures, because of high gene flow.

A preliminary analysis of the population genetics and molecular phylogenetics of Onchocerca volvulus (Nematoda: Filarioidea) using nuclear ribosomal second internal transcribed spacer sequences

Nuclear internal transcribed spacer 2 (ITS2) rDNA sequences were used for a molecular phylogenetics analysis of five Onchocerca species. The sister species of the human parasite O. volvulus was found to be the cattle parasite O. ochengi and not O. gibsoni, contrary to chromosomal evidence. The genetic differentiation of two African populations (representing the two African strains) and a Brazilian population of O. volvulus was also studied. Phylogenetic and network reconstruction did not show any clustering of ITS2 alleles on geographic or strain grounds. Furthermore, population genetics tests showed no indication of population differentiation but suggested gene flow among the three populations.

Population structure of the malaria vector Anopheles darlingi in Rondônia, Brazilian Amazon, based on mitochondrial DNA

Anopheles darlingi is the most important Brazilian malaria vector, with a widespread distribution in the Amazon forest. Effective strategies for vector control could be better developed through knowledge of its genetic structure and gene flow among populations, to assess the vector diversity and competence in transmitting Plasmodium. The aim of this study was to assess the genetic diversity of An. darlingi collected at four locations in Porto Velho, by sequencing a fragment of the ND4 mitochondrial gene. From 218 individual mosquitoes, we obtained 20 different haplotypes with a diversity index of 0.756, equivalent to that found in other neotropical anophelines. The analysis did not demonstrate significant population structure. However, haplotype diversity within some populations seems to be over-represented, suggesting the presence of sub-populations, but the presence of highly represented haplotypes complicates this analysis. There was no clear correlation among genetic and geographical distance and there were differences in relation to seasonality, which is important for malarial epidemiology.