Does stress mess with rodents' heads? Influence of habitat amount and genetic factors in mandible fluctuating asymmetry in South American water rats (Nectomys squamipes, Sigmodontinae) from Brazilian Atlantic rainforest remnants

Loss of developmental stability can lead to deviations from bilateral symmetry (i.e. Fluctuating Asymmetry - FA), and is thought to be caused by environmental and genetic factors associated with habitat loss and stress. Therefore, levels of FA might be a valuable tool to monitor wild populations if FA serves as an indicator of exposure to stress due to impacts of habitat loss and fragmentation. In studies examining FA and habitat fragmentation, FA levels are often explained by loss of genetic variation, though few studies have addressed FA's use as indicator of environmental impact. Here, we investigated whether habitat loss, genetic variation, and/or inbreeding affect the developmental instability in Brazilian Atlantic forest populations of a Neotropical water rat (Nectomys squamipes). We sampled individuals from eight sites within Atlantic forest remnants with different amounts of available forest habitat and assessed FA levels with geometric morphometric techniques using adult mandibles. We used observed heterozygosity (Ho) and inbreeding coefficient (Fis), from seven microsatellite markers, as a proxy of genetic variation at individual and population levels. Populations were not significantly different for shape or size FA levels. Furthermore, interindividual variation in both shape and size FA levels and interpopulational differences in size FA levels were best explained by chance. However, habitat amount was negatively associated with both interpopulational variance and average shape FA levels. This association was stronger in populations living in areas with <28% of forest cover, which presented higher variance and higher average FA, suggesting that Nectomys squamipes might have a tolerance threshold to small availability of habitat. Our work is one of the first to use FA to address environmental stress caused by habitat loss in small mammal populations from a Neotropical biome. We suggest that shape FA might serve as a conservation tool to monitor human impact on natural animal populations.

Influence of female cuticular hydrocarbon (CHC) profile on male courtship behavior in two hybridizing field crickets Gryllus firmus and Gryllus pennsylvanicus

Background

The hybridizing field crickets, Gryllus firmus and Gryllus pennsylvanicus have several barriers that prevent gene flow between species. The behavioral pre-zygotic mating barrier, where males court conspecifics more intensely than heterospecifics, is important because by acting earlier in the life cycle it has the potential to prevent a larger fraction of hybridization. The mechanism behind such male mate preference is unknown. Here we investigate if the female cuticular hydrocarbon (CHC) profile could be the signal behind male courtship.

Results

While males of the two species display nearly identical CHC profiles, females have different, albeit overlapping profiles and some females (between 15 and 45%) of both species display a male-like profile distinct from profiles of typical females. We classified CHC females profile into three categories: G. firmus-like (F; including mainly G. firmus females), G. pennsylvanicus-like (P; including mainly G. pennsylvanicus females), and male-like (ML; including females of both species). Gryllus firmus males courted ML and F females more often and faster than they courted P females (p < 0.05). Gryllus pennsylvanicus males were slower to court than G. firmus males, but courted ML females more often (p < 0.05) than their own conspecific P females (no difference between P and F). Both males courted heterospecific ML females more often than other heterospecific females (p < 0.05, significant only for G. firmus males).

Conclusions

Our results suggest that male mate preference is at least partially informed by female CHC profile and that ML females elicit high courtship behavior in both species. Since ML females exist in both species and are preferred over other heterospecific females, it is likely that this female type is responsible for most hybrid offspring production.

Mapping reduced introgression loci to the X chromosome of the hybridizing field crickets, Gryllus firmus and G. pennsylvanicus

The genomic architecture of barriers to gene exchange during the speciation process is poorly understood. The genomic islands model suggests that loci associated with barriers to gene exchange prevent introgression of nearby genomic regions via linkage disequilibrium. But few analyses of the actual genomic location of non-introgressing loci in closely related species exist. In a previous study Maroja et al. showed that in the hybridizing field crickets, Gryllus firmus and G. pennsylvanicus, 50 non-introgressing loci are localized on two autosomal regions and the X chromosome, but they were not able to map the loci along the X chromosome because they used a male informative cross. Here, we localize the introgressing and non-introgressing loci on the X chromosome, and reveal that all X-linked non-introgressing loci are restricted to a 50-cM region with 10 of these loci mapped to a single location. We discuss the implications of this finding to speciation.

Rapid sexual and genomic isolation in sympatric Drosophila without reproductive character displacement

The rapid evolution of sexual isolation in sympatry has long been associated with reinforcement (i.e., selection to avoid maladaptive hybridization). However, there are many species pairs in sympatry that have evolved rapid sexual isolation without known costs to hybridization. A major unresolved question is what evolutionary processes are involved in driving rapid speciation in such cases. Here, we focus on one such system; the Drosophila athabasca species complex, which is composed of three partially sympatric and interfertile semispecies: WN, EA, and EB. To study speciation in this species complex, we assayed sexual and genomic isolation within and between these semispecies in both sympatric and allopatric populations. First, we found no evidence of reproductive character displacement (RCD) in sympatric zones compared to distant allopatry. Instead, semispecies were virtually completely sexually isolated from each other across their entire ranges. Moreover, using spatial approaches and coalescent demographic simulations, we detected either zero or only weak heterospecific gene flow in sympatry. In contrast, within each semispecies we found only random mating and little population genetic structure, except between highly geographically distant populations. Finally, we determined that speciation in this system is at least an order of magnitude older than previously assumed, with WN diverging first, around 200K years ago, and EA and EB diverging 100K years ago. In total, these results suggest that these semispecies should be given full species status and we adopt new nomenclature: WN-D. athabasca, EA-D. mahican, and EB-D. lenape. While the lack of RCD in sympatry and interfertility do not support reinforcement, we discuss what additional evidence is needed to further decipher the mechanisms that caused rapid speciation in this species complex.

The wavy Mutation Maps to the Inositol 1,4,5-Trisphosphate 3-Kinase 2 (IP3K2) Gene of Drosophila and Interacts with IP3R to Affect Wing Development

Inositol 1,4,5-trisphosphate (IP3) regulates a host of biological processes from egg activation to cell death. When IP3-specific receptors (IP3Rs) bind to IP3, they release calcium from the ER into the cytoplasm, triggering a variety of cell type- and developmental stage-specific responses. Alternatively, inositol polyphosphate kinases can phosphorylate IP3; this limits IP3R activation by reducing IP3 levels, and also generates new signaling molecules altogether. These divergent pathways draw from the same IP3 pool yet cause very different cellular responses. Therefore, controlling the relative rates of IP3R activation vs. phosphorylation of IP3 is essential for proper cell functioning. Establishing a model system that sensitively reports the net output of IP3 signaling is crucial for identifying the controlling genes. Here we report that mutant alleles of wavy (wy), a classic locus of the fruit fly Drosophila melanogaster, map to IP3 3-kinase 2 (IP3K2), a member of the inositol polyphosphate kinase gene family. Mutations in wy disrupt wing structure in a highly specific pattern. RNAi experiments using GAL4 and GAL80(ts) indicated that IP3K2 function is required in the wing discs of early pupae for normal wing development. Gradations in the severity of the wy phenotype provide high-resolution readouts of IP3K2 function and of overall IP3 signaling, giving this system strong potential as a model for further study of the IP3 signaling network. In proof of concept, a dominant modifier screen revealed that mutations in IP3R strongly suppress the wy phenotype, suggesting that the wy phenotype results from reduced IP4 levels, and/or excessive IP3R signaling.

Genes with Restricted Introgression in a Field Cricket (Gryllus firmus/Gryllus pennsylvanicus) Hybrid Zone Are Concentrated on the X Chromosome and a Single Autosome

Characterizing the extent of genomic differentiation between recently diverged lineages provides an important context for understanding the early stages of speciation. When such lineages form discrete hybrid zones, patterns of differential introgression allow direct estimates of which genome regions are likely involved in speciation and local adaptation. Here we use a backcross experimental design to construct a genetic linkage map for the field crickets Gryllus firmus and Gryllus pennsylvanicus, which interact in a well-characterized hybrid zone in eastern North America. We demonstrate that loci with major allele frequency differences between allopatric populations are not randomly distributed across the genome. Instead, most are either X-linked or map to a few small autosomal regions. Furthermore, the subset of those highly differentiated markers that exhibit restricted introgression across the cricket hybrid zone are also concentrated on the X chromosome (39 of 50 loci) and in a single 7-cM region of one autosome. Although the accumulation on the sex chromosome of genes responsible for postzygotic barriers is a well-known phenomenon, less attention has been given to the genomic distribution of genes responsible for prezygotic barriers. We discuss the implications of our results for speciation, both in the context of the role of sex chromosomes and also with respect to the likely causes of heterogeneous genomic divergence. Although we do not yet have direct evidence for the accumulation of ecological, behavioral, or fertilization prezygotic barrier genes on the X chromosome, faster-X evolution could make these barriers more likely to be X-linked.

Barriers to gene exchange in hybridizing field crickets: the role of male courtship effort and cuticular hydrocarbons

Background

Pre-zygotic barriers often involve some form of sexual selection, usually interpreted as female choice, as females are typically the choosier sex. However, males typically show some mate preferences, which are increasingly reported. Here we document previously uncharacterized male courtship behavior (effort and song) and cuticular hydrocarbon (CHC) profiles in the hybridizing crickets Gryllus firmus and G. pennsylvanicus. These two species exhibit multiple barriers to gene exchange that act throughout their life history, including a behavioral barrier that results in increased time to mate in heterospecific pairs.

Results

We demonstrated that male mate choice (as courtship effort allocation) plays a more important role in the prezygotic behavioral barrier than previously recognized. In gryllids females ultimately decide whether or not to mate, yet we found males were selective by regulating courtship effort intensity toward the preferred (conspecific) females. Females were also selective by mating with more intensely courting males, which happened to be conspecifics. We report no differences in courtship song between the two species and suggest that the mechanism that allows males to act differentially towards conspecific and heterospecific females is the cuticular hydrocarbon (CHC) composition. CHC profiles differed between males and females of both species, and there were clear differences in CHC composition between female G. firmus and G. pennsylvanicus but not between the males of each species.

Conclusion

Although many barriers to gene exchange are known in this system, the mechanism behind the mate recognition leading to reduced heterospecific mating remains unknown. The CHC profiles might be the phenotypic cue that allow males to identify conspecifics and thus to adjust their courtship intensity accordingly, leading to differential mating between species.

Development and characterization of 10 microsatellite markers in Sagina nodosa (Caryophyllaceae)

PREMISE OF THE STUDY

We developed 10 novel microsatellite loci for Sagina nodosa, a diploid perennial arctic-alpine herb. To our knowledge, these are the first microsatellite loci for a Sagina species. •

METHODS AND RESULTS

We performed a low-coverage 454 next-generation sequencing of enriched genomic fragments derived from one individual to generate a massive library of contigs containing potential polymorphic microsatellites. We present data for 10 novel polymorphic microsatellite loci containing di-, tri-, tetra-, and hexanucleotide repeats with two to nine alleles per locus assessed in 29 individuals. •

CONCLUSIONS

These polymorphic microsatellite loci in S. nodosa will provide insights on the population structure and life history of S. nodosa in Isle Royale and other North American populations.

High-throughput microsatellite marker development for the distylous herb Primula mistassinica (Primulaceae)

PREMISE OF THE STUDY

Twelve microsatellite markers were developed for Primula mistassinica, a distylous, diploid arctic-alpine plant. The markers will be used to investigate the landscape genetics of a disjunct population on Isle Royale, Michigan, and the phylogeographic patterns of the species. •

METHODS AND RESULTS

We used Roche/454 high-throughput technology to sequence microsatellite-enriched regions in the P. mistassinica genome. We developed 12 polymorphic microsatellite primer sets. These loci contained di-, tri-, and tetranucleotide repeats with two to nine alleles per locus when assessed in 23 individuals. •

CONCLUSIONS

Understanding the historical movements of P. mistassinica will provide insight to the survival prospects of current Arctic plant populations, which face the pressures of global, anthropogenic climate change.

Characterisation and expression of microRNAs in developing wings of the neotropical butterfly Heliconius melpomene

BACKGROUND

Heliconius butterflies are an excellent system for studies of adaptive convergent and divergent phenotypic traits. Wing colour patterns are used as signals to both predators and potential mates and are inherited in a Mendelian manner. The underlying genetic mechanisms of pattern formation have been studied for many years and shed light on broad issues, such as the repeatability of evolution. In Heliconius melpomene, the yellow hindwing bar is controlled by the HmYb locus. MicroRNAs (miRNAs) are important post-transcriptional regulators of gene expression that have key roles in many biological processes, including development. miRNAs could act as regulators of genes involved in wing development, patterning and pigmentation. For this reason we characterised miRNAs in developing butterfly wings and examined differences in their expression between colour pattern races.

RESULTS

We sequenced small RNA libraries from two colour pattern races and detected 142 Heliconius miRNAs with homology to others found in miRBase. Several highly abundant miRNAs were differentially represented in the libraries between colour pattern races. These candidates were tested further using Northern blots, showing that differences in expression were primarily due to developmental stage rather than colour pattern. Assembly of sequenced reads to the HmYb region identified hme-miR-193 and hme-miR-2788; located 2380 bp apart in an intergenic region. These two miRNAs are expressed in wings and show an upregulation between 24 and 72 hours post-pupation, indicating a potential role in butterfly wing development. A search for miRNAs in all available H. melpomene BAC sequences (~2.5 Mb) did not reveal any other miRNAs and no novel miRNAs were predicted.

CONCLUSIONS

Here we describe the first butterfly miRNAs and characterise their expression in developing wings. Some show differences in expression across developing pupal stages and may have important functions in butterfly wing development. Two miRNAs were located in the HmYb region and were expressed in developing pupal wings. Future work will examine the expression of these miRNAs in different colour pattern races and identify miRNA targets among wing patterning genes.

Genealogical discordance and patterns of introgression and selection across a cricket hybrid zone

In recently diverged species, ancestral polymorphism and introgression can cause incongruence between gene and species trees. In the face of hybridization, few genomic regions may exhibit reciprocal monophyly, and these regions, usually evolving rapidly under selection, may be important for the maintenance of species boundaries. In animals with internal fertilization, genes encoding seminal protein are candidate barrier genes. Recently diverged hybridizing species such as the field crickets Gryllus firmus and G. pennsylvanicus, offer excellent opportunities to investigate the origins of barriers to gene exchange. These recently diverged species form a well-characterized hybrid zone, and share ancestral polymorphisms across the genome. We analyzed DNA sequence divergence for seminal protein loci, housekeeping loci, and mtDNA, using a combination of analytical approaches and extensive sampling across both species and the hybrid zone. We report discordant genealogical patterns and differential introgression rates across the genome. The most dramatic outliers, showing near-zero introgression and more structured species trees, are also the only two seminal protein loci under selection. These are candidate barrier genes with possible reproductive functions. We also use genealogical data to examine the demographic history of the field crickets and the current structure of the hybrid zone.

Searching for candidate speciation genes using a proteomic approach: seminal proteins in field crickets

In many animals, male seminal proteins influence gamete interactions and fertilization ability and are probably involved in barriers to gene flow between diverging lineages. Here we use a proteomic approach to identify seminal proteins that are transferred to females during copulation and that may be involved in fertilization barriers between two hybridizing field crickets (Gryllus firmus and Gryllus pennsylvanicus). Analyses of patterns of divergence suggest that much of the field cricket genome has remained undifferentiated following the evolution of reproductive isolation. By contrast, seminal protein genes are highly differentiated. Tests of selection reveal that positive selection is likely to be responsible for patterns of differentiation. Together, our observations suggest that some of the loci encoding seminal proteins may indeed play a role in fertilization barriers in field crickets.

Phylogeography of spruce beetles (Dendroctonus rufipennis Kirby) (Curculionidae: Scolytinae) in North America

Tree-feeding insects that are widespread in north temperate regions are excellent models for studying how past glaciations have impacted differentiation and speciation. We used mitochondrial DNA (mtDNA) sequences and allele frequencies at nine microsatellite loci to examine genetic population structure across the current range of the spruce beetle (Dendroctonus rufipennis), an economically important insect in North America. Two major haplotype groups occur across northern North America, from Newfoundland to Alaska, on white spruce (Picea glauca), and a third distinctive haplotype group occurs throughout the Rocky Mountains on Engelmann spruce (Picea engelmannii). The two mtDNA lineages found in northern populations are 3-4% divergent from each other and from the lineages found in the Rocky Mountains. Analyses of microsatellite data also suggest the existence of major population groupings associated with different geographical regions. In the Pacific Northwest, concordant contact zones for genetically distinct populations of spruce beetles and their principal hosts appear to reflect recent secondary contact. Although we could detect no evidence of historical mtDNA gene flow between allopatric population groups, patterns of variation in the Pacific Northwest suggest recent hybridization and introgression. Together with the pollen record for spruce, they also suggest that beetles have spread from at least three glacial refugia. A minimum estimate of divergence time between the Rocky Mountain and northern populations was 1.7 Myr (million years), presumably reflecting the combined effects of isolation during multiple glacial cycles.

Small mammal populations of an agroecosystem in the Atlantic Forest domain, southeastern Brazil

This study reports 2 years of the population dynamics and reproduction of a small mammal community using the removal method. The study was conducted in a rural area of the Atlantic Forest, in Sumidouro, Rio de Janeiro State, Brazil. The population sizes, age structure and reproduction were studied for the four most common species in the study area. The overall diversity was 1.67 and ranged between 0.8 to 1.67. The species richness was 13 considering the whole study. The most abundant species were the rodents Nectomys squamipes (n = 133), Akodon cursor (n = 74), Oligoryzomys nigripes (n = 25) and the marsupials Didelphis aurita (n = 58) and Philander frenatus (n = 50). Seven other rodents were captured once: Necromys lasiurus, Akodon montensis, Sooretamys angouya, Oecomys catherine, Oxymycterus judex, Euryzygomatomys spinosus and Trinomys iheringi. There were higher peaks for diversity and species richness during the winter (dry) months, probably due to higher food availability. The marsupials had a seasonal reproduction with highest population sizes at the end of the rainy seasons. Nectomys squamipes reproduced mostly during rainy periods. Akodon cursor reproduced predominantly in the winter with the highest population peaks occurring during this season. The analysis of the population dynamics of the rodent species indicated that no species behaved as an agricultural pest, probably due to the heterogeneous landscape of high rotativity of vegetable cultivation. Rodent populations were more susceptible to the removal procedure than marsupial ones.

Identification and comparative analysis of accessory gland proteins in Orthoptera

Accessory reproductive gland proteins (Acps) in Drosophila evolve quickly and appear to play an important role in ensuring the fertilization success of males. Moreover, Acps are thought to be involved in establishing barriers to fertilization between closely related species. While accessory glands are known to occur in the males of many insect groups, the proteins that are passed on to females by males during mating have not been well characterized outside of Drosophila. To gain a better understanding of these proteins, we characterized ESTs from the accessory glands of two cricket species, Allonemobius fasciatus and Gryllus firmus. Using an expressed sequence tag (EST) approach, followed by bioinformatic and evolutionary analyses, we found that many proteins are secreted and, therefore, available for transfer to the female during mating. Further, we found that most ESTs are novel, showing little sequence similarity between taxa. Evolutionary analyses suggest that cricket proteins are subject to diversifying selection and indicate that Allonemobius is much less polymorphic than Gryllus. Despite rapid nucleotide sequence divergence, there appears to be functional conservation of protein classes among Drosophila and cricket taxa.

Molecular evolution of seminal proteins in field crickets

In sexually reproducing organisms, male ejaculates are complex traits that are potentially subject to many different selection pressures. Recent experimental evidence supports the hypothesis that postmating sexual selection, and particularly sexual conflict, may play a key role in the evolution of the proteinaceous components of ejaculates. However, this evidence is based almost entirely on the study of Drosophila, a species with a mating system characterized by a high cost of mating for females. In this paper, we broaden our understanding of the role of selection on the evolution of seminal proteins by characterizing these proteins in field crickets, a group of insects in which females appear to benefit from mating multiply. We have used an experimental protocol that can be applied to other organisms for which complete genome sequences are not yet available. By combining an evolutionary expressed sequence tag screen of the male accessory gland in 2 focal species (Gryllus firmus and Gryllus pennsylvanicus) with a bioinformatics approach, we have been able to identify as many as 30 seminal proteins. Evolutionary analyses among 5 species of the genus Gryllus suggest that seminal protein genes evolve more rapidly than genes encoding proteins that are not involved with reproduction. The rates of synonymous substitution (dS) are similar in genes encoding seminal proteins and genes encoding "housekeeping" proteins. For the same comparison, the rate of fixation of nonsynonymous substitutions (dN) is 3 times higher in genes encoding seminal proteins, suggesting that the divergence of seminal proteins in field crickets has been accelerated by positive Darwinian selection. In spite of the contrasting characteristics of the Drosophila and Gryllus mating systems, the mean selection parameter omega and the proportion of loci estimated to be affected by positive selection are very similar.

Nectomys squamipes microsatellites and homologous loci in sigmodontine rodents

Three monomorphic and four highly polymorphic microsatellites of Nectomys squamipes were isolated and characterised in a sample of 141 specimens from eight different Brazilian localities. These seven microsatellites and four others previously described in this species were tested in seven other nonfocus sigmodontine species. At least three loci were successfully amplified in every species, but none was amplified in all species. All sequenced products in nonfocus species showed (GT)(n) motifs as in N. squamipes. Several loci were amplified in Nectomys rattus and Oligoryzomys nigripes, while absence of PCR products was observed more frequently in Oxymycterus dasythricus and Akodon cursor. Two of three monomorphic loci in N. squamipes were polymorphic in other species.

Small non-flying mammals from conserved and altered areas of Atlantic forest and Cerrado: comments on their potential use for monitoring environment

Two Atlantic Forests and two Cerrado areas in Brazil were sampled for non-flying small mammal fauna. In each biome one area with altered and another with almost unaltered vegetation (national parks), were chosen to investigate these fauna. Species richness of Atlantic Forest and Cerrado was comparable in the conserved as well as in the altered areas. Data suggested that species could be divided into different ecological categories according to distribution, use of altered and/or relatively unaltered vegetation and habitat specificity. Within these ecological categories some species are appropriate indicators for monitoring environmental quality and degradation. Useful guidelines for wildlife management planning, including selecting areas for conservation units and their better boundary delimitation can ensue.

The parasitism of Schistosoma mansoni (Digenea-Trematoda) in a naturally infected population of water rats, Nectomys squamipes (Rodentia-Sigmodontinae) in Brazil

Schistosomiasis is a health problem in Brazil and the role of rodents in maintaining the schistosome life-cycle requires further clarification. The influence of Schistosoma mansoni on a population of Nectomys squamipes was studied by capture-recapture (1st phase, from June 1991 to November 1995) and removal (2nd phase, from April 1997 to March 1999) studies at Sumidouro, Rio de Janeiro, Brazil. During both phases coproscopic examinations were performed. At the 2nd phase the rodents were perfused and worms were counted. The population dynamics of parasites was studied. During the 1st phase, female reproductive parameters, longevity, recruitment and survivorship rates and migration patterns were studied in relation to schistosome prevalence. Water contamination (source of miracidia), abundance intermediate host and rodent migration were related to prevalence. The N. squamipes population was not obviously influenced by the infection, as shown by the high number of reproductive infected females, high longevity of infected individuals and the absence of a relationship between recruitment or survivorship rates and the intensity of schistosome infection. The data indicate that N. squamipes can increase transmission of S. mansoni in endemic areas and carry it to non-infected areas. Furthermore, this rodent can be used as an indicator of a transmission focus.