Search

Background Planarians reliably regenerate all body parts after injury, including a fully functional head and central nervous system. But until now, the expression dynamics and functional role of miRNAs and other small RNAs during the process of head regeneration are not well understood. Furthermore, little is known about the evolutionary conservation of the relevant small RNAs pathways, rendering it difficult to assess whether insights from planarians will apply to other taxa. Results In this study, we applied high throughput sequencing to identify miRNAs, tRNA fragments and piRNAs that are dynamically expressed during head regeneration in Dugesia japonica. We further show that knockdown of selected small RNAs, including three novel Dugesia-specific miRNAs, during head regeneration induces severe defects including abnormally small-sized eyes, cyclopia and complete absence of eyes. Conclusions Our findings suggest that a complex pool of small RNAs takes part in the process of head regeneration in Dugesia japonica and provide novel insights into global small RNA expression profiles and expression changes in response to head amputation. Our study reveals the evolutionary conserved role of miR-124 and brings further promising candidate small RNAs into play that might unveil new avenues for inducing restorative programs in non-regenerative organisms via small RNA mimics based therapies.

Ciliopathies often comprise retinal degeneration since the photoreceptors outer segment is an adapted primary cilium. CEP162 is a distal end centriolar protein required for proper transition zone assembly during ciliogenesis and whose loss causes ciliopathy in zebrafish. CEP162 has so far not been implicated in human disease. Here, we identified a homozygous CEP162 frameshift variant, c.1935dupA (p.(E646R*5)), in retinitis pigmentosa patients from two unrelated Moroccan families, likely representing a founder allele. We found that even though mRNA levels were reduced, the truncated CEP162-E646R*5 protein was expressed and localized to the mitotic spindle during mitosis, but not at the basal body of the cilium. In CEP162 knockdown cells, expression of the truncated CEP162-E646R*5 protein is unable to restore ciliation indicating its loss of function at the cilium. In patient fibroblasts, cilia overcome the absence of CEP162 from the primary cilium by delaying ciliogenesis through the persistence of CP110 at the mother centriole. The patient fibroblasts are ultimately able to extend some abnormally long cilia that are missing key transition zone components. Defective transition zone formation likely disproportionately affects the long-living ciliary outer segment of photoreceptors resulting in retinal dystrophy. CEP162 is expressed in human retina, and we show that wild-type CEP162, but not truncated CEP162-E646R*5, specifically localizes to the distal end of centrioles of mouse photoreceptor cilia. Together, our genetic, cell-based, and in vivo modeling establish that CEP162 deficiency causes retinal ciliopathy in humans. SIGNIFICANCE STATEMENTCiliopathies often comprise retinal degeneration since the photoreceptors outer segment is an adapted primary cilium. CEP162 is a basal body protein required for proper transition zone assembly during ciliogenesis and has so far not been implicated in human disease. Using genetic, cell-based, and in vivo studies, we show that a biallelic mutation of CEP162 causes late-onset human retinal degeneration. This mutation specifically hinders CEP162 function at the cilia, leading to impaired transition zone assembly and delayed formation of dysmorphic cilia. Studies in mouse photoreceptors support that absence of CEP162 could lead to defective outer segments. The discovery of novel ciliopathy genes, such as CEP162, advances our insight into cell-specific functions of the primary cilium.

Fragments of mature tRNAs have long been considered as mere degradation products without physiological function. However, recent reports show that tRNA-derived small RNAs (tsRNAs) play prominent roles in diverse cellular processes across a wide spectrum of species. Contrasting the situation in other small RNA pathways the mechanisms behind these effects appear more diverse, more complex, and are generally less well understood. In addition, surprisingly little is known about the expression profiles of tsRNAs across different tissues and species. Here, we provide an initial overview of tsRNA expression in different species and tissues, revealing very high levels of 5′ tRNA halves (5′ tRHs) particularly in the primate hippocampus. We further modulated the regulation capacity of selected 5′ tRHs in human cells by transfecting synthetic tsRNA mimics (“overexpression”) or antisense-RNAs (“inhibition”) and identified differentially expressed transcripts based on RNA-seq. We then used a novel k-mer mapping approach to dissect the underlying targeting rules, suggesting that 5′ tRHs silence genes in a sequence-specific manner, while the most efficient target sites align to the mid-region of the 5′ tRH and are located within the CDS or 3′ UTR of the target. This amends previous observations that tsRNAs guide Argonaute proteins to silence their targets via a miRNA-like 5′ seed match and suggests a yet unknown mechanism of regulation. Finally, our data suggest that some 5′ tRHs that are also able to sequence-specifically stabilize mRNAs as up-regulated mRNAs are also significantly enriched for 5′ tRH target sites.

PIWI-interacting RNAs (piRNAs) and their partnering PIWI proteins defend the animal germline against transposable elements and play a crucial role in fertility. Numerous studies in the past have uncovered many additional functions of the piRNA pathway, including gene regulation, anti-viral defense, and somatic transposon repression. Further, comparative analyses across phylogenetic groups showed that the PIWI/piRNA system evolves rapidly and exhibits great evolutionary plasticity. However, the presence of so-called piRNA clusters as the major source of piRNAs is common to nearly all metazoan species. These genomic piRNA-producing loci are highly divergent across taxa and critically influence piRNA populations in different evolutionary lineages. We launched the initial version of the piRNA cluster database to facilitate research on regulation and evolution of piRNA-producing loci across tissues und species. In recent years the amount of small RNA sequencing data that was generated and the abundance of species that were studied has grown rapidly. To keep up with this recent progress, we have released a major update for the piRNA cluster database (https://www.smallrnagroup.uni-mainz.de/piRNAclusterDB), expanding it from 12 to a total of 51 species with hundreds of new datasets, and revised its overall structure to enable easy navigation through this large amount of data.

PIWI proteins and PIWI-interacting RNAs (piRNAs) suppress transposon activity in animals, thus protecting their genomes from detrimental insertion mutagenesis. Here, we reveal that PIWI genes and piRNAs are ubiquitously expressed in mollusks, similar to the situation in arthropods. We describe lineage-specific adaptations of transposon composition in piRNA clusters in the great pond snail and the pacific oyster, likely reflecting differential transposon activity in gastropods and bivalves. We further show that different piRNA clusters with unique transposon composition are dynamically expressed during oyster development. Finally, bioinformatics analyses suggest that different populations of piRNAs presumably bound to different PIWI paralogs participate in homotypic and heterotypic ping-pong amplification loops in a tissue- and sex-specific manner. Together with recent findings from other animal species, our results support the idea that somatic piRNA expression represents the ancestral state in metazoans.

Temperature has a major impact on gene expression in ectotherms. But until recently, it was not clear in which way, if any, small non-coding RNAs such as miRNAs or piRNAs contribute to thermosensitive gene regulation. We have recently shown that temperature-responsive miRNAs in Drosophila drive adaptation to different ambient temperatures on the transcriptome level. Moreover, we demonstrated that higher temperatures lead to a more efficient piRNA-dependent transposon silencing, possibly due to heat-induced unfolding of RNA secondary structures. In this commentary, we will dwell upon particular interesting aspects connected to our findings, hoping that our point of view may encourage other scientists to address some of the questions raised here. We will particularly focus on aspects related to climate-dependent transposon propagation in evolution and putative transgenerational epigenetic effects of altered small RNA transcriptomes. We further briefly indicate how temperature-responsive miRNAs may confound the interpretation of data obtained from experiments comprising heat-shock treatment which is a widely used technique not only in Drosophila genetics.

The majority of Drosophila genes are expressed in a temperature-dependent manner, but the way in which small RNAs may contribute to this effect is completely unknown as we currently lack an idea of how small RNA transcriptomes change as a function of temperature. Applying high-throughput sequencing techniques complemented by quantitative real-time PCR experiments, we demonstrate that altered ambient temperature induces drastic but reversible changes in sequence composition and total abundance of both miRNA and piRNA populations. Further, mRNA sequencing reveals that the expression of miRNAs and their predicted target transcripts correlates inversely, suggesting that temperature-responsive miRNAs drive adaptation to different ambient temperatures on the transcriptome level. Finally, we demonstrate that shifts in temperature affect both primary and secondary piRNA pools, and the observed aberrations are consistent with altered expression levels of the involved Piwi-pathway factors. We further reason that enhanced ping–pong processing at 29°C is driven by dissolved RNA secondary structures at higher temperatures, uncovering target sites that are not accessible at low temperatures. Together, our results show that small RNAs are an important part of epigenetic regulatory mechanisms that ensure homeostasis and adaptation under fluctuating environmental conditions.

Piwi-interacting RNAs (piRNAs) are 26–30-nucleotide germ line-specific small non-coding RNAs that have evolutionarily conserved function in mobile genetic element (transposons) silencing and maintenance of genome integrity. Drosophila Hsp70/90-organizing protein homolog (Hop), a co-chaperone, interacts with piRNA-binding protein Piwi and mediates silencing of phenotypic variations. However, it is not known whether Hop has a direct role in piRNA biogenesis and transposon silencing. Here, we show that knockdown of Hop in the germ line nurse cells (GLKD) of Drosophila ovaries leads to activation of transposons. Hop GLKD females can lay eggs at the same rate as wild-type counterparts, but the eggs do not hatch into larvae. Hop GLKD leads to the accumulation of γ-H2Av foci in the germ line, indicating increased DNA damage in the ovary. We also show that Hop GLKD-induced transposon up-regulation is due to inefficient piRNA biogenesis. Based on these results, we conclude that Hop is a critical component of the piRNA pathway and that it maintains genome integrity by silencing transposons.

Fragments of mature tRNAs have long been considered as mere degradation products without physiological function. However, recent reports show that tRNA fragments (tRFs) play prominent roles in diverse cellular processes across a wide spectrum of species. Contrasting the situation in other small RNA pathways the mechanisms behind these effects appear more diverse, more complex and are generally less well understood. In addition, surprisingly little is known about the expression profiles of tRFs across different tissues and species. Here, we provide an initial overview of tRF expression in different species and tissues, revealing very high tRF-levels particularly in the primate hippocampus. We further modulated the regulation capacity of selected tRFs in human cells by transfecting synthetic tRF mimics (“overexpression”) or antisense-RNAs (“inhibition”) and identified differentially expressed transcripts based on RNAseq. We then used a novel k-mer mapping approach to dissect the underlying targeting rules, demonstrating that 5’ tRNA halves (5’ tRHs) silence genes in a sequence-specific manner, while the most efficient target sites align to the mid-region of the 5’ tRH and are located within the CDS or 3’ UTR of the target. This amends previous observations that tRFs guide Argonaut proteins to silence their targets via a miRNA-like 5’ seed match and suggests a yet unknown mechanism of regulation. Finally, our data suggests that some 5’ tRHs are also able to sequence-specifically stabilize mRNAs as upregulated mRNAs are also significantly enriched for 5’ tRH target sites.

Codon composition, GC-content and local RNA secondary structures can have a profound effect on gene expression and mutations affecting these parameters, even though they do not alter the protein sequence, are not neutral in terms of selection. Although evidence exists that in some cases selection favors more stable RNA secondary structures, we currently lack a concrete idea of how many genes are affected within a species, and if this is a universal phenomenon in nature.We searched for signs of structural selection in a global manner, analyzing a set of one million coding sequences from 73 species representing all domains of life, as well as viruses, by means of our newly developed software PACKEIS. We show that codon composition and amino acid identity are main determinants of RNA secondary structure. In addition, we show that the arrangement of synonymous codons within coding sequences is non-random, yielding extremely high, but also extremely low secondary structures significantly more often than expected by chance.Together, we demonstrate that selection for high and low secondary structure is a widespread phenomenon. Our results provide another line of evidence that synonymous mutations are less neutral than commonly thought, which is of importance for many evolutionary models.

A monophyletic origin of endoparasitic thorny-headed worms (Acanthocephala) and wheel-animals (Rotifera) is widely accepted. However, the phylogeny inside the clade, be it called Syndermata or Rotifera, has lacked validation by mitochondrial (mt) data. Herein, we present the first mt genome of the key taxon Seison and report conflicting results of phylogenetic analyses: while mt sequence-based topologies showed monophyletic Lemniscea (Bdelloidea + Acanthocephala), gene order analyses supported monophyly of Pararotatoria (Seisonidea + Acanthocephala) and Hemirotifera (Bdelloidea + Pararotatoria). Sequence-based analyses obviously suffered from substitution saturation, compositional bias, and branch length heterogeneity; however, we observed no compromising effects in gene order analyses. Moreover, gene order-based topologies were robust to changes in coding (genes vs. gene pairs, two-state vs. multistate, aligned vs. non-aligned), tree reconstruction methods, and the treatment of the two monogonont mt genomes. Thus, mt gene order verifies seisonids as sister to acanthocephalans within monophyletic Hemirotifera, while deviating results of sequence-based analyses reflect artificial signal. This conclusion implies that the complex life cycle of extant acanthocephalans evolved from a free-living state, as retained by most monogononts and bdelloids, via an epizoic state with a simple life cycle, as shown by seisonids. Hence, Acanthocephala represent a rare example where ancestral transitional stages have counterparts amongst the closest relatives.

Piwi proteins and their guiding small RNAs, termed Piwi-interacting (pi-) RNAs, are essential for silencing of transposons in the germline of animals. A substantial fraction of piRNAs originates from genomic loci termed piRNA clusters and sequences encoded in these piRNA clusters determine putative targets for the Piwi/piRNA system. In the past decade, studies of piRNA transcriptomes in different species revealed additional roles for piRNAs beyond transposon silencing, reflecting the astonishing plasticity of the Piwi/piRNA system along different phylogenetic branches. Moreover, piRNA transcriptomes can change drastically during development and vary across different tissues. Since piRNA clusters crucially shape piRNA profiles, analysis of these loci is imperative for a thorough understanding of functional and evolutionary aspects of the piRNA pathway. But despite the ever-growing amount of available piRNA sequence data, we know little about the factors that determine differential regulation of piRNA clusters, nor the evolutionary events that cause their gain or loss. In order to facilitate addressing these subjects, we established a user-friendly piRNA cluster database (http://www.smallrnagroup-mainz.de/piRNAclusterDB.html) that provides comprehensive data on piRNA clusters in multiple species, tissues and developmental stages based on small RNA sequence data deposited at NCBI's Sequence Read Archive (SRA).

Throughout the domains of life, transposon activity represents a serious threat to genome integrity and evolution has realized different molecular mechanisms that aim to inhibit the transposition of mobile DNA. Small noncoding RNAs that function as guides for Argonaute effector proteins represent a key feature of so-called RNA interference (RNAi) pathways and specialized RNAi pathways exist to repress transposon activity on the transcriptional and posttranscriptional level. Transposon transcription can be diminished by targeted DNA methylation or chromatin remodeling via repressive Histone modifications. Posttranscriptional transposon silencing bases on degradation of transposon transcripts to prevent either reverse transcription followed by genomic reintegration or translation into proteins that mediate the transposition process. In plants, Argonaute-like proteins guided by short interfering RNAs (siRNAs) are essential for transposon repression on the epigenetic and posttranscriptional level. In the germline of animals, these tasks are often assumed by a second subclass of Argonaute proteins referred to as Piwi-like proteins, which bind a distinct class of small noncoding RNAs named piwi-interacting RNAs (piRNAs). Though the principals of RNAi pathways are essentially the same in all eukaryotic organisms, remarkable differences can be observed even in closely related species reflecting the astonishing plasticity and diversity of these pathways.

PIWI proteins and a specific class of small non-coding RNAs, termed Piwi interacting RNAs (piRNAs), suppress transposon activity in animals on the transcriptional and post-transcriptional level, thus protecting genomes from detrimental insertion mutagenesis. While in vertebrates the PIWI/piRNA system appears to be restricted to the germline, somatic expression of piRNAs directed against transposons is widespread in arthropods, likely representing the ancestral state for this phylum. Here, we show that somatic expression of PIWI genes and piRNAs directed against transposons is conserved in mollusks, suggesting that somatic PIWI/piRNA expression was already realized in an early bilaterian ancestor. We further describe lineage specific adaptations regarding transposon composition of piRNA clusters and show that different piRNA clusters are dynamically expressed during oyster development. Finally, bioinformatics analyses suggest that different populations of piRNAs participate in the ping-pong amplification loop in a tissue specific manner.

Argonaute proteins comprising Piwi-like and Argonaute-like proteins and their guiding small RNAs combat mobile DNA on the transcriptional and post-transcriptional level. While Piwi-like proteins and associated piRNAs are generally restricted to the germline, Argonaute-like proteins and siRNAs have been linked with transposon control in the germline as well as in the soma. Intriguingly, evolution has realized distinct Argonaute subfunctionalization patterns in different species but our knowledge about mammalian RNA interference pathways relies mainly on findings from the mouse model. However, mice differ from other mammals by absence of functional Piwil3 and expression of an oocyte-specific Dicer isoform. Thus, studies beyond the mouse model are required for a thorough understanding of function and evolution of mammalian RNA interference pathways. We high-throughput sequenced small RNAs from the male Tupaia belangeri germline, which represents a close outgroup to primates, hence phylogenetically links mice with humans. We identified transposon-derived piRNAs as well as siRNAs clearly contrasting the separation of piRNA- and siRNA-pathways into male and female germline as seen in mice. Genome-wide analysis of tree shrew transposons reveal that putative siRNAs map to transposon sites that form foldback secondary structures thus representing suitable Dicer substrates. In contrast piRNAs target transposon sites that remain accessible. With this we provide a basic mechanistic explanation how secondary structure of transposon transcripts influences piRNA- and siRNA-pathway utilization. Finally, our analyses of tree shrew piRNA clusters indicate A-Myb and the testis-expressed transcription factor RFX4 to be involved in the transcriptional regulation of mammalian piRNA clusters.

The transcription repressor FOXP2 is a crucial player in nervous system evolution and development of humans and songbirds. In order to provide an additional insight into its functional role we compared target gene expression levels between human neuroblastoma cells (SH-SY5Y) stably overexpressing either human FOXP2 cDNA or its orthologues from the common chimpanzee, Rhesus monkey, and marmoset, respectively. Subsequent RNA-seq led to identification of 27 genes with differential regulation under the control of human FOXP2, which were previously reported to have FOXP2-driven and/or songbird song-related expression regulation. Importantly, RT-qPCR and Western blotting indicated differential regulation of additional 13 new target genes in response to overexpression of human FOXP2. These genes may be directly regulated by FOXP2 considering numerous matches of established FOXP2 binding motifs as well as publicly available FOXP2-ChIP-seq reads within their putative promoters. Ontology analysis of the new and reproduced targets, along with their interactors in a network, revealed an enrichment of terms relating to cellular signaling and communication, metabolism and catabolism, cellular migration and differentiation, and expression regulation. Notably, terms including the words “neuron” or “axonogenesis” were also enriched. Complementary literature screening uncovered many connections to human developmental (autism spectrum disease, schizophrenia, Down syndrome, agenesis of corpus callosum, trismus-pseudocamptodactyly, ankyloglossia, facial dysmorphology) and neurodegenerative diseases and disorders (Alzheimer’s, Parkinson’s, and Huntington’s diseases, Lewy body dementia, amyotrophic lateral sclerosis). Links to deafness and dyslexia were detected, too. Such relations existed for single proteins (e.g., DCDC2, NURR1, PHOX2B, MYH8, MYH13) and groups of proteins which conjointly function in mRNA processing, ribosomal recruitment, cell-cell adhesion (e.g., CDH4), cytoskeleton organization, neuro-inflammation, and processing of amyloid precursor protein. Conspicuously, many links pointed to an involvement of the FOXP2-driven network in JAK/STAT signaling and the regulation of the ezrin-radixin-moesin complex. Altogether, the applied phylogenetic perspective substantiated FOXP2’s importance for nervous system development, maintenance, and functioning. However, the study also disclosed new regulatory pathways that might prove to be useful for understanding the molecular background of the aforementioned developmental disorders and neurodegenerative diseases.

BackgroundNext generation sequencing is a key technique in small RNA biology research that has led to the discovery of functionally different classes of small non-coding RNAs in the past years. However, reliable annotation of the extensive amounts of small non-coding RNA data produced by high-throughput sequencing is time-consuming and requires robust bioinformatics expertise. Moreover, existing tools have a number of shortcomings including a lack of sensitivity under certain conditions, limited number of supported species or detectable sub-classes of small RNAs.ResultsHere we introduce unitas, an out-of-the-box ready software for complete annotation of small RNA sequence datasets, supporting the wide range of species for which non-coding RNA reference sequences are available in the Ensembl databases (currently more than 800). unitas combines high quality annotation and numerous analysis features in a user-friendly manner. A complete annotation can be started with one simple shell command, making unitas particularly useful for researchers not having access to a bioinformatics facility. Noteworthy, the algorithms implemented in unitas are on par or even outperform comparable existing tools for small RNA annotation that base on available ncRNA sequence information.Conclusionsunitas brings together annotation and analysis features that hitherto required the installation of numerous different bioinformatics tools which can pose a challenge for the non-expert user. With this, unitas overcomes the problem of read normalization. Moreover, the high quality of sequence annotation and analysis, paired with the ease of use, make unitas a valuable tool for researchers in all fields connected to small RNA biology.

PIWI proteins and their associated piRNAs have been the focus of intensive research in the past decade; therefore, their participation in the maintenance of genomic integrity during spermatogenesis has been well established. Recent studies have suggested important roles for the PIWI/piRNA system outside of gametogenesis, based on the presence of piRNAs and PIWI proteins in several somatic tissues, cancers, and the early embryo. Here, we investigated the small RNA complement present in bovine gonads, gametes, and embryos through next-generation sequencing. A distinct piRNA population was present in the testis as expected. However, we also found a large population of slightly shorter, 24–27 nt piRNA-like RNA (pilRNAs) in pools of oocytes and zygotes. These oocyte and embryo pilRNAs exhibited many of the canonical characteristics of piRNAs including a 1U bias, the presence of a ‘ping-pong’ signature, genomic clustering, and transposable element targeting. Some of the major transposons targeted by oocyte and zygote pilRNA were from the LINE RTE and ERV1 classes. We also identified pools of pilRNA potentially derived from, or targeted at, specific mRNA sequences. We compared the frequency of these gene-associated pilRNAs to the fold change in the expression of respective mRNAs from two previously reported transcriptome datasets. We observed significant negative correlations between the number of pilRNAs targeting mRNAs, and their fold change in expression between the 4–8 cell and 8–16 cell stages. Together, these results represent one of the first characterizations of the PIWI/piRNA pathway in the translational bovine model, and in the novel context of embryogenesis.

article i nfo Human Fetuin-A (FETUA, also termed AHSG) is a serum protein composed of two cystatin-like domains D1 and D2 of together 235 amino acids (aa) and an unrelated domain D3 of 114 aa. Though the protein plays a role in diverse physiological and pathological processes, comparably little is known about sequence evolution of FETUA. We therefore analyzed its molecular evolution on the basis of coding sequences of 16 primate species. Ratios of non-synonymous to synonymous substitution rates (dn/ds= ω) suggest that a previously reported acceleration of sequence evolution of exon 7, which encodes domain D3, is driven by positive selection instead of neutral evolution. Irrespective of the maximum likelihood model used, a remarkable 40% of codon sites is estimated to be positively selected in domain D3, with an average ω of 4.3. Analyses of GC/AT- and AT/GC-substitutions reveal that GC-biased gene conversion (gBGC) cannot explain the observed pattern. We thus reject the extended null model of molecular evolution (neutral evolution and gBGC) and conclude that positive Darwinian selection contributes to the sequence evolution of FETUA domain D3. Remarkably, domain D3 of human FETUA is characterized by a significant accumulation of proteolytic cleavage sites. Though speculative at present, the accelerated evolution of domain D3 might be due to the evolution of differential cleavage site patterns across primates reflecting different paths to achieve one goal, the decomposition of FETUA under pathological conditions. Present analyses use a PERL-script ("Substitution Bias Indicator", SBI) for the extraction of the numbers of GC/AT- and AT/GC-substitutions from supplemental results files of PAML analyses. The script is available upon request and the presented approach of testing the extended null hypothesis of molecular evolution can be transferred to every other protein coding locus.

Video abstract A video abstract by the authors of this paper is available. video-abstract8870.mov, The PCR-amplification of unknown homologous or paralogous genes generally relies on PCR primers predicted from multi sequence alignments. But increasing sequence divergence can induce the need to use degenerate primers which entails the problem of testing the characteristics, unwanted interactions and potential mispriming of degenerate primers. Here I introduce easyPAC, a new software for the prediction of degenerate primers from multi sequence alignments or single consensus sequences. As a major innovation, easyPAC allows to apply all customary primer test procedures to degenerate primer sequences including fast mapping to reference files. Thus, easyPAC simplifies and expedites the designing of specific degenerate primers enormously. Degenerate primers suggested by easyPAC were used in PCR amplification with subsequent de novo sequencing of TDRD1 exon 11 homologs from several representatives of the haplorrhine primate phylogeny. The results demonstrate the efficient performance of the suggested primers and therefore show that easyPAC can advance upcoming comparative genetic studies.

To assess the relative impact of functional constraint and post-mating sexual selection on sequence evolution of reproductive proteins, we examined 169 primate sperm proteins. In order to recognize potential genome-wide trends, we additionally analysed a sample of altogether 318 non-reproductive (brain and postsynaptic) proteins. Based on cDNAs of eight primate species (Anthropoidea), we observed that pre-mating sperm proteins engaged in sperm composition and assembly show significantly lower incidence of site-specific positive selection and overall lower non-synonymous to synonymous substitution rates ( d N / d S ) across sites as compared with post-mating sperm proteins involved in capacitation, hyperactivation, acrosome reaction and fertilization. Moreover, database screening revealed overall more intracellular protein interaction partners in pre-mating than in post-mating sperm proteins. Finally, post-mating sperm proteins evolved at significantly higher evolutionary rates than pre-mating sperm and non-reproductive proteins on the branches to multi-male breeding species, while no such increase was observed on the branches to unimale and monogamous species. We conclude that less protein–protein interactions of post-mating sperm proteins account for lowered functional constraint, allowing for stronger impact of post-mating sexual selection, while the opposite holds true for pre-mating sperm proteins. This pattern is particularly strong in multi-male breeding species showing high female promiscuity.

Background Several mechanistic models aim to explain the diversification of the multitude of endemic species on Madagascar. The island's biogeographic history probably offered numerous opportunities for secondary contact and subsequent hybridization. Existing diversification models do not consider a possible role of these processes. One key question for a better understanding of their potential importance is how they are influenced by different environmental settings. Here, we characterized a contact zone between two species of mouse lemurs, Microcebus griseorufus and M. murinus, in dry spiny bush and mesic gallery forest that border each other sharply without intermediate habitats between them. We performed population genetic analyses based on mtDNA sequences and nine nuclear microsatellites and compared the results to a known hybrid zone of the same species in a nearby wide gradient from dry spiny bush over transitional forest to humid littoral forest. Results In the spiny-gallery system, Microcebus griseorufus is restricted to the spiny bush; Microcebus murinus occurs in gallery forest and locally invades the dryer habitat of its congener. We found evidence for bidirectional introgressive hybridization, which is closely linked to increased spatial overlap within the spiny bush. Within 159 individuals, we observed 18 hybrids with mitochondrial haplotypes of both species. Analyses of simulated microsatellite data indicate that we identified hybrids with great accuracy and that we probably underestimated their true number. We discuss short-term climatic fluctuations as potential trigger for the dynamic of invasion and subsequent hybridization. In the gradient hybrid zone in turn, long-term aridification could have favored unidirectional nuclear introgression from Microcebus griseorufus into M. murinus in transitional forest. Conclusions Madagascar's southeastern transitional zone harbors two very different hybrid zones of mouse lemurs in different environmental settings. This sheds light on the multitude of opportunities for the formation of hybrid zones and indicates an important influence of environmental factors on secondary contact and hybridization. Our findings suggest that hybridization could enhance the adaptability of mouse lemurs without necessarily leading to a loss of distinctiveness. They point to a potential role of hybridization in Madagascar's diversification history that requires further investigation.

The role of the Piwi/piRNA pathway during mammalian oogenesis has remained enigmatic thus far, especially since experiments with Piwi knockout mice did not reveal any phenotypic defects in female individuals. This is in striking contrast with results obtained from other species including flies and zebrafish. In mouse oocytes, however, only low levels of piRNAs are found and they are not required for their function. We recently demonstrated dynamic expression of PIWIL1, PIWIL2, and PIWIL3 during mammalian oogenesis and early embryogenesis. In addition, small RNA analysis of human, crab-eating macaque and cattle revealed that piRNAs are also expressed in the female germline and closely resemble piRNAs from testis. Here, we thoroughly describe the experimental and computational methods that we applied for the generation, processing and analyses of next generation sequencing (NGS) data associated with our study on Piwi proteins and piRNAs in mammalian oocytes and embryos (Roovers et al., 2015). The complete sequence data is available at NCBI's Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) under the accession GSE64942.