Soil fungal communities differ between shaded and sun-intensive coffee plantations in El Salvador.

Coffea arabica is a highly traded commodity worldwide, and its plantations are habitat to a wide range of organisms. Coffee farmers are shifting away from traditional shade coffee farms in favor of sun-intensive, higher yield farms, which can impact local biodiversity. Using plant-associated microorganisms in biofertilizers, particularly fungi collected from local forests, to increase crop yields has gained traction among coffee producers. However, the taxonomic and spatial distribution of many fungi in coffee soil, nearby forests and biofertilizers is unknown. We collected soil samples from a sun coffee system, shade coffee system, and nearby forest from Izalco, Sonsonate, El Salvador. At each coffee system, we collected soil from the surface (upper) and 10 cm below the surface (lower), and from the coffee plant drip line (drip line) and the walkway between two plants (walkway). Forest soils were collected from the surface only. We used ITS metabarcoding to characterize fungal communities in soil and in the biofertilizer (applied in both coffee systems), and assigned fungal taxa to functional guilds using FUNGuild. In the sun and shade coffee systems, we found that drip line soil had higher richness in pathotrophs, symbiotrophs, and saprotrophs than walkway soil, suggesting that fungi select for microhabitats closer to coffee plants. Upper and lower soil depths did not differ in fungal richness or composition, which may reflect the shallow root system of Coffea arabica. Soil from shade, sun, and forest had similar numbers of fungal taxa, but differed dramatically in community composition, indicating that local habitat differences drive fungal species sorting among systems. Yet, some fungal taxa were shared among systems, including seven fungal taxa present in the biofertilizer. Understanding the distribution of coffee soil mycobiomes can be used to inform sustainable, ecologically friendly farming practices and identify candidate plant-growth promoting fungi for future studies.

Phylogeography and connectivity of molluscan parasites: Perkinsus spp. in Panama and beyond.

Panama is a major hub for commercial shipping between two oceans, making it an ideal location to examine parasite biogeography, potential invasions, and the spread of infectious agents. Our goals were to (i) characterise the diversity and genetic connectivity of Perkinsus spp. haplotypes across the Panamanian Isthmus and (ii) combine these data with sequences from around the world to evaluate the current phylogeography and genetic connectivity of these widespread molluscan parasites. We collected 752 bivalves from 12 locations along the coast of Panama including locations around the Bocas del Toro archipelago and the Caribbean and Pacific entrances to the Panama Canal, from December 2012 to February 2013. We used molecular genetic methods to screen for Perkinsus spp. and obtained internal transcribed spacer region (ITS) ribosomal DNA (rDNA) sequences for all positive samples. Our sequence data were used to evaluate regional haplotype diversity and distribution across both coasts of Panama, and were then combined with publicly available sequences to create global haplotype networks. We found 26 ITS haplotypes from four Perkinsus spp. (1-12 haplotypes per species) in Panama. Perkinsus beihaiensis haplotypes had the highest genetic diversity, were the most regionally widespread, and were associated with the greatest number of hosts. On a global scale, network analyses demonstrated that some haplotypes found in Panama were cosmopolitan (Perkinsus chesapeaki, Perkinsus marinus), while others were more geographically restricted (Perkinsus olseni, P. beihaiensis), indicating different levels of genetic connectivity and dispersal. We found some Perkinsus haplotypes were shared across the Isthmus of Panama and several regions around the world, including across ocean basins. We also found that haplotype diversity is currently underestimated and directly related to the number of sequences. Nevertheless, our results demonstrate long-range dispersal and global connectivity for many haplotypes, suggesting that dispersal through shipping probably contributes to these biogeographical patterns.

Protistan Biogeography: A Snapshot Across a Major Shipping Corridor Spanning Two Oceans.

Deciphering patterns of protistan taxa is a crucial step for understanding anthropogenic and environmental impacts on biogeography. We characterized and compared protistan communities from environmental samples collected along a major shipping corridor, the Panama Canal, and the Bocas del Toro archipelago. We used metabarcoding with high throughput sequencing (HTS) with the V4 hypervariable region of the ribosomal gene complex (rDNA). We detected many protistan taxa, including a variety of parasitic and toxic taxa. There were 1,296 OTUs shared across all three regions, with an additional 342-1,526 OTUs occurring across two or more regions, suggesting some mixing within the Caribbean and across the Isthmus. In general, this mixing did not impact community similarity, which was primarily distinct across regions. When OTUs identified as gregarines were analyzed separately, most samples grouped by region and communities were distinct across the Canal. Shipping traffic through the Panama Canal could move some taxa across regions; however, different environmental conditions in the two oceans may limit their establishment. Overall our results suggest that contemporary protistan biogeographic patterns are likely caused by a complex combination of factors, including anthropogenic dispersal and environmental tolerance.

Richness and distribution of tropical oyster parasites in two oceans.

Parasites can exert strong effects on population to ecosystem level processes, but data on parasites are limited for many global regions, especially tropical marine systems. Characterizing parasite diversity and distributions are the first steps towards understanding the potential impacts of parasites. The Panama Canal serves as an interesting location to examine tropical parasite diversity and distribution, as it is a conduit between two oceans and a hub for international trade. We examined metazoan and protistan parasites associated with ten oyster species collected from both Panamanian coasts, including the Panama Canal and Bocas del Toro. We found multiple metazoan taxa (pea crabs, Stylochus spp., Urastoma cyrinae). Our molecular screening for protistan parasites detected four species of Perkinsus (Perkinsus marinus, Perkinsus chesapeaki, Perkinsus olseni, Perkinsus beihaiensis) and several haplosporidians, including two genera (Minchinia, Haplosporidium). Species richness was higher for the protistan parasites than for the metazoans, with haplosporidian richness being higher than Perkinsus richness. Perkinsus species were the most frequently detected and most geographically widespread among parasite groups. Parasite richness and overlap differed between regions, locations and oyster hosts. These results have important implications for tropical parasite richness and the dispersal of parasites due to shipping associated with the Panama Canal.

Conservation and divergence in the frog immunome: pyrosequencing and de novo assembly of immune tissue transcriptomes.

BACKGROUND: Frogs are a diverse group of vertebrates for which limited genomic resources are available. Natural frog populations face a multitude of threats, including habitat degradation, infectious disease, and environmental change. Characterizing the functional genomics of anuran tissues in general - and the immune system in particular - will enhance our knowledge of genetic and epigenetic responses to environmental threats and inform conservation and recovery efforts. RESULTS: To increase the number of species with genomic datasets and characterize gene expression in immune-related tissues, we sequenced the transcriptomes of three tissues from two frogs (Espadarana prosoblepon and Lithobates yavapaiensis) on the Roche 454 GS FLX platform. Our sequencing produced 8881 E. prosoblepon and 5428 L. yavapaiensis annotated gene products after de novo assembly and Gene Ontology classification. Transcripts of the innate and acquired immune system were expressed in all three tissues. Inflammatory response and acquired immunity transcripts were significantly more diverged between E. prosoblepon and L. yavapaiensis compared to innate immunity and immune system development transcripts. Immune-related transcripts did not show an overall elevated rate of functional evolution, with the exception of glycosyl proteases, which include lysozymes, central bacterial and fungal-killing enzymes of the innate immune system. CONCLUSIONS: The three frog transcriptomes provide more than 600 Mbp of new genomic data, and will serve as a valuable framework for future comparative studies of non-model anurans. Additionally, we show that immune gene divergence varies by functional group and that transcriptome studies can be useful in comparing rates of evolutionary change across gene families.

Cryptic population structuring and the role of the Isthmus of Tehuantepec as a gene flow barrier in the critically endangered Central American River Turtle.

The critically endangered Central American River Turtle (Dermatemys mawii) is the only remaining member of the Dermatemydidae family, yet little is known about its population structuring. In a previous study of mitochondrial (mt) DNA in the species, three main lineages were described. One lineage (Central) was dominant across most of the range, while two other lineages were restricted to Papaloapan (PAP; isolated by the Isthmus of Tehuantepec and the Sierra de Santa Marta) or the south-eastern part of the range (1D). Here we provide data from seven polymorphic microsatellite loci and the R35 intron to re-evaluate these findings using DNA from the nuclear genome. Based on a slightly expanded data set of a total of 253 samples from the same localities, we find that mtDNA and nuclear DNA markers yield a highly congruent picture of the evolutionary history and population structuring of D. mawii. While resolution provided by the R35 intron (sequenced for a subset of the samples) was very limited, the microsatellite data revealed pronounced population structuring. Within the Grijalva-Usumacinta drainage basin, however, many populations separated by more than 300 kilometers showed signals of high gene flow. Across the entire range, neither mitochondrial nor nuclear DNA show a significant isolation-by-distance pattern, but both genomes highlight that the D. mawii population in the Papaloapan basin is genetically distinctive. Further, both marker systems detect unique genomic signals in four individuals with mtDNA clade 1D sampled on the southeast edge of the Grijalva-Usumacinta basin. These individuals may represent a separate cryptic taxon that is likely impacted by recent admixture.

Whole genome SNP discovery and analysis of genetic diversity in Turkey (Meleagris gallopavo).

BACKGROUND: The turkey (Meleagris gallopavo) is an important agricultural species and the second largest contributor to the world's poultry meat production. Genetic improvement is attributed largely to selective breeding programs that rely on highly heritable phenotypic traits, such as body size and breast muscle development. Commercial breeding with small effective population sizes and epistasis can result in loss of genetic diversity, which in turn can lead to reduced individual fitness and reduced response to selection. The presence of genomic diversity in domestic livestock species therefore, is of great importance and a prerequisite for rapid and accurate genetic improvement of selected breeds in various environments, as well as to facilitate rapid adaptation to potential changes in breeding goals. Genomic selection requires a large number of genetic markers such as e.g. single nucleotide polymorphisms (SNPs) the most abundant source of genetic variation within the genome. RESULTS: Alignment of next generation sequencing data of 32 individual turkeys from different populations was used for the discovery of 5.49 million SNPs, which subsequently were used for the analysis of genetic diversity among the different populations. All of the commercial lines branched from a single node relative to the heritage varieties and the South Mexican turkey population. Heterozygosity of all individuals from the different turkey populations ranged from 0.17-2.73 SNPs/Kb, while heterozygosity of populations ranged from 0.73-1.64 SNPs/Kb. The average frequency of heterozygous SNPs in individual turkeys was 1.07 SNPs/Kb. Five genomic regions with very low nucleotide variation were identified in domestic turkeys that showed state of fixation towards alleles different than wild alleles. CONCLUSION: The turkey genome is much less diverse with a relatively low frequency of heterozygous SNPs as compared to other livestock species like chicken and pig. The whole genome SNP discovery study in turkey resulted in the detection of 5.49 million putative SNPs compared to the reference genome. All commercial lines appear to share a common origin. Presence of different alleles/haplotypes in the SM population highlights that specific haplotypes have been selected in the modern domesticated turkey.

Mitochondrial and nuclear DNA sequences reveal recent divergence in morphologically indistinguishable petrels.

Often during the process of divergence, genetic markers will only gradually obtain the signal of isolation. Studies of recently diverged taxa utilizing both mitochondrial and nuclear data sets may therefore yield gene trees with differing levels of phylogenetic signal as a result of differences in coalescence times. However, several factors can lead to this same pattern, and it is important to distinguish between them to gain a better understanding of the process of divergence and the factors driving it. Here, we employ three nuclear intron loci in addition to the mitochondrial Cytochrome b gene to investigate the magnitude and timing of divergence between two endangered and nearly indistinguishable petrel taxa: the Galapagos (GAPE) and Hawaiian (HAPE) petrels (Pterodroma phaeopygia and P. sandwichensis). Phylogenetic analyses indicated reciprocal monophyly between these two taxa for the mitochondrial data set, but trees derived from the nuclear introns were unresolved. Coalescent analyses revealed effectively no migration between GAPE and HAPE over the last 100,000 generations and that they diverged relatively recently, approximately 550,000 years ago, coincident with a time of intense ecological change in both the Galapagos and Hawaiian archipelagoes. This indicates that recent divergence and incomplete lineage sorting are causing the difference in the strength of the phylogenetic signal of each data set, instead of insufficient variability or ongoing male-biased dispersal. Further coalescent analyses show that gene flow is low even between islands within each archipelago suggesting that divergence may be continuing at a local scale. Accurately identifying recently isolated taxa is becoming increasingly important as many clearly recognizable species are already threatened by extinction.

Long-term isolation of a highly mobile seabird on the Galapagos.

The Galapagos Islands are renowned for their high degree of endemism. Marine taxa inhabiting the archipelago might be expected to be an exception, because of their utilization of pelagic habitats-the dispersal barrier for terrestrial taxa-as foraging grounds. Magnificent frigatebirds (Fregata magnificens) have a highly vagile lifestyle and wide geographical distribution around the South and Central American coasts. Given the potentially high levels of gene flow among populations, the species provides a good test of the effectiveness of the Galapagos ecosystem in isolating populations of highly dispersive marine species. We studied patterns of genetic (mitochondrial DNA, microsatellites and nuclear introns) and morphological variation across the distribution of magnificent frigatebirds. Concordant with predictions from life-history traits, we found signatures of extensive gene flow over most of the range, even across the Isthmus of Panama, which is a major barrier to gene flow in other tropical seabirds. In contrast, individuals from the Galapagos were strongly differentiated from all conspecifics, and have probably been isolated for several hundred thousand years. Our finding is a powerful testimony to the evolutionary uniqueness of the taxa inhabiting the Galapagos archipelago and its associated marine ecosystems.

Mid-Pleistocene divergence of Cuban and North American ivory-billed woodpeckers.

We used ancient DNA analysis of seven museum specimens of the endangered North American ivory-billed woodpecker (Campephilus principalis) and three specimens of the species from Cuba to document their degree of differentiation and their relationships to other Campephilus woodpeckers. Analysis of these mtDNA sequences reveals that the Cuban and North American ivory bills, along with the imperial woodpecker (Campephilus imperialis) of Mexico, are a monophyletic group and are roughly equidistant genetically, suggesting each lineage may be a separate species. Application of both internal and external rate calibrations indicates that the three lineages split more than one million years ago, in the Mid-Pleistocene. We thus can exclude the hypothesis that Native Americans introduced North American ivory-billed woodpeckers to Cuba. Our sequences of all three woodpeckers also provide an important DNA barcoding resource for identification of non-invasive samples or remains of these critically endangered and charismatic woodpeckers.

Vocal dialects, sex-biased dispersal, and microsatellite population structure in the parrot Amazona auropalliata.

Geographic variation in microsatellite allele frequencies was assessed at nine sites in two regional vocal dialects of the parrot Amazona auropalliata (yellow-naped amazon) to test for correspondence between dialects and population structure. There was no relationship between the genetic distances between individuals and their dialect membership. High rates of gene flow were estimated between vocal dialects based on genetic differentiation. In addition, 5.5% of pairs of individuals compared across the dialect boundary were estimated to be related at the level of half siblings, indicating that dispersal is ongoing. The number of effective migrants per generation between dialects estimated with the microsatellite data was roughly one-seventh the number estimated with mitochondrial control region sequence data from the same individuals, suggesting that gene flow may be female-biased. Together, these results suggest that the observed mosaic pattern of geographic variation in vocalizations is maintained by learning of local call types by immigrant birds after dispersal. We found no evidence that ongoing habitat fragmentation has contributed to cryptic population structure.

Male dominance, paternity, and relatedness in the Jamaican fruit-eating bat (Artibeus jamaicensis).

We analysed variation at 14 nuclear microsatellite loci to assess the genetic structure, relatedness, and paternity of polygynous Jamaican fruit-eating bats. A total of 84 adults captured in two caves exhibited little genetic differentiation between caves (FST = 0.008). Average relatedness among adult females in 10 harem groups was very low (R = 0.014 +/- 0.011), providing no evidence of harem structure. Dominant and subordinate males shared paternity in large groups, while dominant and satellite males shared paternity in smaller groups. However, our results suggest that male rank influences paternity. Dominant males fathered 69% of 40 offspring, followed by satellite (22%) and subordinate males (9%). Overall adult male bats are not closely related, however, in large harem groups we found that subordinate and dominant males exhibited relatedness values consistent with a father-offspring relationship. Because dominant and subordinate males also sired all the pups in large groups, we propose that their association provides inclusive fitness to them.