Genomic and phenomic analysis of island ant community assembly.

Island biodiversity has long fascinated biologists as it typically presents tractable systems for unpicking the eco-evolutionary processes driving community assembly. In general, two recurring themes are of central theoretical interest. First, immigration, diversification, and extinction typically depend on island geographical properties (e.g., area, isolation, and age). Second, predictable ecological and evolutionary trajectories readily occur after colonization, such as the evolution of adaptive trait syndromes, trends toward specialization, adaptive radiation, and eventual ecological decline. Hypotheses such as the taxon cycle draw on several of these themes to posit particular constraints on colonization and subsequent eco-evolutionary dynamics. However, it has been challenging to examine these integrated dynamics with traditional methods. Here, we combine phylogenomics, population genomics and phenomics, to unravel community assembly dynamics among Pheidole (Hymenoptera, Formicidae) ants in the isolated Fijian archipelago. We uphold basic island biogeographic predictions that isolated islands accumulate diversity primarily through in situ evolution rather than dispersal, and population genomic support for taxon cycle predictions that endemic species have decreased dispersal ability and demography relative to regionally widespread taxa. However, rather than trending toward island syndromes, ecomorphological diversification in Fiji was intense, filling much of the genus-level global morphospace. Furthermore, while most endemic species exhibit demographic decline and reduced dispersal, we show that the archipelago is not an evolutionary dead-end. Rather, several endemic species show signatures of population and range expansion, including a successful colonization to the Cook islands. These results shed light on the processes shaping island biotas and refine our understanding of island biogeographic theory.

Conserved community structure and simultaneous divergence events in the fig wasps associated with Ficus benjamina in Australia and China.

BACKGROUND: Localised patterns of species diversity can be influenced by many factors, including regional species pools, biogeographic features and interspecific interactions. Despite recognition of these issues, we still know surprisingly little about how invertebrate biodiversity is structured across geographic scales. In particular, there have been few studies of how insect communities vary geographically while using the same plant host. We compared the composition (species, genera) and functional structure (guilds) of the chalcid wasp communities associated with the widespread fig tree, Ficus benjamina, towards the northern (Hainan province, China) and southern (Queensland, Australia) edges of its natural range. Sequence data were generated for nuclear and mtDNA markers and used to delimit species, and Bayesian divergence analyses were used to test patterns of community cohesion through evolutionary time. RESULTS: Both communities host at least 14 fig wasp species, but no species are shared across continents. Community composition is similar at the genus level, with six genera shared although some differ in species diversity between China and Australia; a further three genera occur in only China or Australia. Community functional structure remains very similar in terms of numbers of species in each ecological guild despite community composition differing a little (genera) or a lot (species), depending on taxonomic level. Bayesian clustering analyses favour a single community divergence event across continents over multiple events for different ecological guilds. Molecular dating estimates of lineage splits between nearest inter-continental species pairs are broadly consistent with a scenario of synchronous community divergence from a shared "ancestral community". CONCLUSIONS: Fig wasp community structure and genus-level composition are largely conserved in a wide geographic comparison between China and Australia. Moreover, dating analyses suggest that the functional community structure has remained stable for long periods during historic range expansions. This suggests that ecological interactions between species may play a persistent role in shaping these communities, in contrast to findings in some comparable temperate systems.

Phylogenetic analysis of RAD-seq data: examining the influence of gene genealogy conflict on analysis of concatenated data.

One of the major issues in phylogenetic analysis is that gene genealogies from different gene regions may not reflect the true species tree or history of speciation. This has led to considerable debate about whether concatenation of loci is the best approach for phylogenetic analysis. The application of Next-generation sequencing techniques such as RAD-seq generates thousands of relatively short sequence reads from across the genomes of the sampled taxa. These data sets are typically concatenated for phylogenetic analysis leading to data sets that contain millions of base pairs per taxon. The influence of gene region conflict among so many loci in determining the phylogenetic relationships among taxa is unclear. We simulated RAD-seq data by sampling 100 and 500 base pairs from alignments of over 6000 coding regions that each produce one of three highly supported alternative phylogenies of seven species of Drosophila. We conducted phylogenetic analyses on different sets of these regions to vary the sampling of loci with alternative gene trees to examine the effect on detecting the species tree. Irrespective of sequence length sampled per region and which subset of regions was used, phylogenetic analyses of the concatenated data always recovered the species tree. The results suggest that concatenated alignments of Next-generation data that consist of many short sequences are robust to gene tree/species tree conflict when the goal is to determine the phylogenetic relationships among taxa.

Molecular species delimitation of a symbiotic fig-pollinating wasp species complex reveals extreme deviation from reciprocal partner specificity.

BACKGROUND: Symbiotic relationships have contributed to major evolutionary innovations, the maintenance of fundamental ecosystem functions, and the generation and maintenance of biodiversity. However, the exact nature of host/symbiont associations, which has important consequences for their dynamics, is often poorly known due to limited understanding of symbiont taxonomy and species diversity. Among classical symbioses, figs and their pollinating wasps constitute a highly diverse keystone resource in tropical forest and savannah environments. Historically, they were considered to exemplify extreme reciprocal partner specificity (one-to-one host-symbiont species relationships), but recent work has revealed several more complex cases. However, there is a striking lack of studies with the specific aims of assessing symbiont diversity and how this varies across the geographic range of the host. RESULTS: Here, we use molecular methods to investigate cryptic diversity in the pollinating wasps of a widespread Australian fig species. Standard barcoding genes and methods were not conclusive, but incorporation of phylogenetic analyses and a recently developed nuclear barcoding gene (ITS2), gave strong support for five pollinator species. Each pollinator species was most common in a different geographic region, emphasising the importance of wide geographic sampling to uncover diversity, and the scope for divergence in coevolutionary trajectories across the host plant range. In addition, most regions had multiple coexisting pollinators, raising the question of how they coexist in apparently similar or identical resource niches. CONCLUSION: Our study offers a striking example of extreme deviation from reciprocal partner specificity over the full geographical range of a fig-wasp system. It also suggests that superficially identical species may be able to co-exist in a mutualistic setting albeit at different frequencies in relation to their fig host's range. We show that comprehensive sampling and molecular taxonomic techniques may be required to uncover the true structure of cryptic biodiversity underpinning intimate ecological interactions.

Enhanced and asymmetric signatures of hybridization at climatic margins: Evidence from closely related dioecious fig species.

Hybridization plays a significant role in biological evolution. However, it is not clear whether ecological contingency differentially influences likelihood of hybridization, particularly at ecological margins where parental species may exhibit reduced fitnesses. Moreover, it is unknown whether future ecosystem change will increase the prevalence of hybridization. Ficus heterostyla and F. squamosa are closely related species co-distributed from southern Thailand to southwest China where hybridization, yielding viable seeds, has been documented. As a robust test of ecological factors driving hybridization, we investigated spatial hybridization signatures based on nuclear microsatellites from extensive population sampling across a widespread contact range. Both species showed high population differentiation and strong patterns of isolation by distance. Admixture estimates exposed asymmetric interspecific gene flow. Signatures of hybridization increase significantly towards higher latitude zones, peaking at the northern climatic margins. Geographic variation in reproductive phenology combined with ecologically challenging marginal habitats may promote this phenomenon. Our work is a first systematic evaluation of such patterns in a comprehensive, latitudinally-based clinal context, and indicates that tendency to hybridize appears strongly influenced by environmental conditions. Moreover, that future climate change scenarios will likely alter and possibly augment cases of hybridization at ecosystem scales.

Maintenance of specificity in sympatric host-specific fig/wasp pollination mutualisms.

Background: Fig/wasp pollination mutualisms are extreme examples of species-specific plant-insect symbioses, but incomplete specificity occurs, with potentially important evolutionary consequences. Why pollinators enter alternative hosts, and the fates of pollinators and the figs they enter, are unknown. Methods: We studied the pollinating fig wasp, Ceratosolen emarginatus, which concurrently interacts with its typical host Ficus auriculata and the locally sympatric alternative host F. hainanensis, recording frequencies of the wasp in figs of the alternative hosts. We measured ovipositor lengths of pollinators and style lengths in female and male figs in the two host species. Volatile organic compounds (VOCs) emitted by receptive figs of each species were identified using GC-MS. We tested the attraction of wasps to floral scents in choice experiments, and detected electrophysiologically active compounds by GC-EAD. We introduced C. emarginatus foundresses into figs of both species to reveal the consequences of entry into the alternative host. Results: C. emarginatus entered a low proportion of figs of the alternative host, and produced offspring in a small proportion of them. Despite differences in the VOC profiles of the two fig species, they included shared semiochemicals. Although C. emarginatus females prefer receptive figs of F. auriculata, they are also attracted to those of F. hainanensis. C. emarginatus that entered male figs of F. hainanensis produced offspring, as their ovipositors were long enough to reach the bottom of the style; however, broods were larger and offspring smaller than in the typical host. Female figs of F. hainanensis failed to produce seeds when visited by C. emarginatus. These findings advance our current understanding of how these species-specific mutualisms usually remain stable and the conditions that allow their diversification.

riceExplorer: Uncovering the Hidden Potential of a National Genomic Resource Against a Global Database.

Agricultural crop breeding programs, particularly at the national level, typically consist of a core panel of elite breeding cultivars alongside a number of local landrace varieties (or other endemic cultivars) that provide additional sources of phenotypic and genomic variation or contribute as experimental materials (e.g., in GWAS studies). Three issues commonly arise. First, focusing primarily on core development accessions may mean that the potential contributions of landraces or other secondary accessions may be overlooked. Second, elite cultivars may accumulate deleterious alleles away from nontarget loci due to the strong effects of artificial selection. Finally, a tendency to focus solely on SNP-based methods may cause incomplete or erroneous identification of functional variants. In practice, integration of local breeding programs with findings from global database projects may be challenging. First, local GWAS experiments may only indicate useful functional variants according to the diversity of the experimental panel, while other potentially useful loci-identifiable at a global level-may remain undiscovered. Second, large-scale experiments such as GWAS may prove prohibitively costly or logistically challenging for some agencies. Here, we present a fully automated bioinformatics pipeline (riceExplorer) that can easily integrate local breeding program sequence data with international database resources, without relying on any phenotypic experimental procedure. It identifies associated functional haplotypes that may prove more robust in determining the genotypic determinants of desirable crop phenotypes. In brief, riceExplorer evaluates a global crop database (IRRI 3000 Rice Genomes) to identify haplotypes that are associated with extreme phenotypic variation at the global level and recorded in the database. It then examines which potentially useful variants are present in the local crop panel, before distinguishing between those that are already incorporated into the elite breeding accessions and those only found among secondary varieties (e.g., landraces). Results highlight the effectiveness of our pipeline, identifying potentially useful functional haplotypes across the genome that are absent from elite cultivars and found among landraces and other secondary varieties in our breeding program. riceExplorer can automatically conduct a full genome analysis and produces annotated graphical output of chromosomal maps, potential global diversity sources, and summary tables.

Predicting distributions of Wolbachia strains through host ecological contact-Who's manipulating whom?

Reproductive isolation in response to divergent selection is often mediated via third-party interactions. Under these conditions, speciation is inextricably linked to ecological context. We present a novel framework for understanding arthropod speciation as mediated by Wolbachia, a microbial endosymbiont capable of causing host cytoplasmic incompatibility (CI). We predict that sympatric host sister-species harbor paraphyletic Wolbachia strains that provide CI, while well-defined congeners in ecological contact and recently diverged noninteracting congeners are uninfected due to Wolbachia redundancy. We argue that Wolbachia provides an adaptive advantage when coupled with reduced hybrid fitness, facilitating assortative mating between co-occurring divergent phenotypes-the contact contingency hypothesis. To test this, we applied a predictive algorithm to empirical pollinating fig wasp data, achieving up to 91.60% accuracy. We further postulate that observed temporal decay of Wolbachia incidence results from adaptive host purging-adaptive decay hypothesis-but implementation failed to predict systematic patterns. We then account for post-zygotic offspring mortality during CI mating, modeling fitness clines across developmental resources-the fecundity trade-off hypothesis. This model regularly favored CI despite fecundity losses. We demonstrate that a rules-based algorithm accurately predicts Wolbachia infection status. This has implications among other systems where closely related sympatric species encounter adaptive disadvantage through hybridization.

Identification of Bacterial Blight Resistance Loci in Rice (Oryza sativa L.) against Diverse Xoo Thai Strains by Genome-Wide Association Study.

Bacterial leaf blight (BLB) is a serious disease affecting global rice agriculture caused by Xanthomonas oryzae pv. oryzae (Xoo). Most resistant rice lines are dependent on single genes that are vulnerable to resistance breakdown caused by pathogen mutation. Here we describe a genome-wide association study of 222 predominantly Thai rice accessions assayed by phenotypic screening against 20 Xoo isolates. Loci corresponding to BLB resistance were detected using >142,000 SNPs. We identified 147 genes according to employed significance thresholds across chromosomes 1-6, 8, 9 and 11. Moreover, 127 of identified genes are located on chromosomal regions outside estimated Linkage Disequilibrium influences of known resistance genes, potentially indicating novel BLB resistance markers. However, significantly associated SNPs only occurred across a maximum of six Xoo isolates indicating that the development of broad-spectrum Xoo strain varieties may prove challenging. Analyses indicated a range of gene functions likely underpinning BLB resistance. In accordance with previous studies of accession panels focusing on indica varieties, our germplasm displays large numbers of SNPs associated with resistance. Despite encouraging data suggesting that many loci contribute to resistance, our findings corroborate previous inferences that multi-strain resistant varieties may not be easily realised in breeding programs without resorting to multi-locus strategies.

Development of 19 universal microsatellite loci for three closely related Ficus species (Moraceae) by high-throughput sequencing.

The genus Ficus is striking for its species diversity, ecological significance, and its often species-specific relationship with coevolved pollinating fig wasps, which has long fascinated biologists. The three closely related and generally co-distributed dioecious species Ficus hispida, F. heterostyla and F. squamosa provide an ideal system for the study of speciation, hybridization (caused by pollinator sharing) and comparative phylogeography to infer historical biogeography. We aimed to develop microsatellite markers for these allied species to facilitate the outlined study investigations. A DNA library was constructed from one F. heterostyla sample, and 19 microsatellite loci were developed based on high-throughput sequencing. These markers showed relatively high polymorphism in all three fig species. The mean number of alleles per locus was 3.594-5.286, and the mean observed and expected heterozygosity were 0.469-0.546 and 0.467-0.528, respectively. Principal coordinate, STRUCTURE and AMOVA analyses revealed different degrees of genetic differentiation within species, and, despite some observed genetic admixture, indicated the presence of clear boundaries between different species. In summary, we successfully developed universal microsatellite markers for three closely related Ficus species. These markers will be of great value for investigating patterns of biodiversity among the species in this model system for coevolutionary studies.

Genome-Wide Association Analysis Identifies Resistance Loci for Bacterial Leaf Streak Resistance in Rice (Oryza sativa L.).

Bacterial leaf streak (BLS) caused by Xanthomonas oryzae pv. oryzicola (Xoc) is one of the most devastating diseases in rice production areas, especially in humid tropical and subtropical zones throughout Asia and worldwide. A genome-wide association study (GWAS) analysis conducted on a collection of 236 diverse rice accessions, mainly indica varieties, identified 12 quantitative trait loci (QTLs) on chromosomes 1, 2, 3, 4, 5, 8, 9 and 11, conferring resistance to five representative isolates of Thai Xoc. Of these, five QTLs conferred resistance to more than one Xoc isolates. Two QTLs, qBLS5.1 and qBLS2.3, were considered promising QTLs for broad-spectrum resistance to BLS. The xa5 gene was proposed as a potential candidate gene for qBLS5.1 and three genes, encoding pectinesterase inhibitor (OsPEI), eukaryotic zinc-binding protein (OsRAR1), and NDP epimerase function, were proposed as candidate genes for qBLS2.3. Results from this study provide an insight into the potential QTLs and candidate genes for BLS resistance in rice. The recessive xa5 gene is suggested as a potential candidate for strong influence on broad-spectrum resistance and as a focal target in rice breeding programs for BLS resistance.

Partial begging: an empirical model for the early evolution of offspring signalling.

Species where, from birth, the offspring feed themselves in addition to begging for food from the parents can be described as 'partially begging'. Such species provide a unique opportunity to examine the evolution of offspring begging from non-signalling offspring foraging strategies. We used the partially begging burying beetle Nicrophorus vespilloides to test specific hypotheses concerning the coexistence of begging and self-feeding. We first tested whether the cessation of larval begging coincided with an increase in the efficiency of self-feeding. As predicted, begging ceased when the efficiency of self-feeding reached the point where the larvae grew just as well without as with access to food provided by the parent. We next tested whether the transition to nutritional independence was under parental or offspring control. The parent did not change its behaviour towards the larvae over time, while the larvae changed their behaviour by reducing the time spent begging in the presence of the parent. Food allocation during the transition to nutritional independence was therefore under offspring control. Our results on partial begging provide a starting point for new theoretical models for the origin of begging. We suggest that these should be constructed as scramble-competition models because the offspring control food allocation.