Generation of somatic de novo structural variation as a hallmark of cellular senescence in human lung fibroblasts.

Cellular senescence is characterized by replication arrest in response to stress stimuli. Senescent cells accumulate in aging tissues and can trigger organ-specific and possibly systemic dysfunction. Although senescent cell populations are heterogeneous, a key feature is that they exhibit epigenetic changes. Epigenetic changes such as loss of repressive constitutive heterochromatin could lead to subsequent LINE-1 derepression, a phenomenon often described in the context of senescence or somatic evolution. LINE-1 elements decode the retroposition machinery and reverse transcription generates cDNA from autonomous and non-autonomous TEs that can potentially reintegrate into genomes and cause structural variants. Another feature of cellular senescence is mitochondrial dysfunction caused by mitochondrial damage. In combination with impaired mitophagy, which is characteristic of senescent cells, this could lead to cytosolic mtDNA accumulation and, as a genomic consequence, integrations of mtDNA into nuclear DNA (nDNA), resulting in mitochondrial pseudogenes called numts. Thus, both phenomena could cause structural variants in aging genomes that go beyond epigenetic changes. We therefore compared proliferating and senescent IMR-90 cells in terms of somatic de novo numts and integrations of a non-autonomous composite retrotransposons - the so-called SVA elements-that hijack the retropositional machinery of LINE-1. We applied a subtractive and kinetic enrichment technique using proliferating cell DNA as a driver and senescent genomes as a tester for the detection of nuclear flanks of de novo SVA integrations. Coupled with deep sequencing we obtained a genomic readout for SVA retrotransposition possibly linked to cellular senescence in the IMR-90 model. Furthermore, we compared the genomes of proliferative and senescent IMR-90 cells by deep sequencing or after enrichment of nuclear DNA using AluScan technology. A total of 1,695 de novo SVA integrations were detected in senescent IMR-90 cells, of which 333 were unique. Moreover, we identified a total of 81 de novo numts with perfect identity to both mtDNA and nuclear hg38 flanks. In summary, we present evidence for possible age-dependent structural genomic changes by paralogization that go beyond epigenetic modifications. We hypothesize, that the structural variants we observe potentially impact processes associated with replicative aging of IMR-90 cells.

RDA coupled with deep sequencing detects somatic SVA-retrotranspositions and mosaicism in the human brain.

Cells of the developing human brain are affected by the progressive acquisition of genetic and epigenetic alterations that have been reported to contribute to somatic mosaicism in the adult brain and are increasingly considered a possible cause of neurogenetic disorders. A recent work uncovered that the copy-paste transposable element (TE) LINE-1 (L1) is mobilized during brain development, and thus mobile non-autonomous TEs like AluY and SINE-VNTR-Alu (SVA) families can use L1 activity in trans, leading to de novo insertions that may influence the variability of neural cells at genetic and epigenetic levels. In contrast to SNPs and when considering substitutional sequence evolution, the presence or absence of TEs at orthologous loci represents highly informative clade markers that provide insights into the lineage relationships between neural cells and how the nervous system evolves in health and disease. SVAs, as the 'youngest' class of hominoid-specific retrotransposons preferentially found in gene- and GC-rich regions, are thought to differentially co-regulate nearby genes and exhibit a high mobility in the human germline. Therefore, we determined whether this is reflected in the somatic brain and used a subtractive and kinetic enrichment technique called representational difference analysis (RDA) coupled with deep sequencing to compare different brain regions with respect to de novo SINE-VNTR-Alu insertion patterns. As a result, we detected somatic de novo SVA integrations in all human brain regions analyzed, and the majority of de novo insertions can be attributed to lineages of telencephalon and metencephalon, since most of the examined integrations are unique to different brain regions under scrutiny. The SVA positions were used as presence/absence markers, forming informative sites that allowed us to create a maximum parsimony phylogeny of brain regions. Our results largely recapitulated the generally accepted evo-devo patterns and revealed chromosome-wide rates of de novo SVA reintegration targets and preferences for specific genomic regions, e.g., GC- and TE-rich regions as well as close proximity to genes that tend to fall into neural-specific Gene Ontology pathways. We concluded that de novo SVA insertions occur in the germline and somatic brain cells at similar target regions, suggesting that similar retrotransposition modes are effective in the germline and soma.

R-loop landscape in mature human sperm: Regulatory and evolutionary implications.

R-loops are three-stranded nucleic acid structures consisting of an RNA:DNA hybrid and a displaced DNA strand. While R-loops pose a potential threat to genome integrity, they constitute 5% of the human genome. The role of R-loops in transcriptional regulation, DNA replication, and chromatin signature is becoming increasingly clear. R-loops are associated with various histone modifications, suggesting that they may modulate chromatin accessibility. To potentially harness transcription-coupled repair mechanisms in the germline, nearly the entire genome is expressed during the early stages of male gametogenesis in mammals, providing ample opportunity for the formation of a transcriptome-dependent R-loop landscape in male germ cells. In this study, our data demonstrated the presence of R-loops in fully mature human and bonobo sperm heads and their partial correspondence to transcribed regions and chromatin structure, which is massively reorganized from mainly histone to mainly protamine-packed chromatin in mature sperm. The sperm R-loop landscape resembles characteristic patterns of somatic cells. Surprisingly, we detected R-loops in both residual histone and protamine-packed chromatin and localize them to still-active retroposons, ALUs and SINE-VNTR-ALUs (SVAs), the latter has recently arisen in hominoid primates. We detected both evolutionarily conserved and species-specific localizations. Comparing our DNA-RNA immunoprecipitation (DRIP) data with published DNA methylation and histone chromatin immunoprecipitation (ChIP) data, we hypothesize that R-loops epigenetically reduce methylation of SVAs. Strikingly, we observe a strong influence of R-loops on the transcriptomes of zygotes from early developmental stages before zygotic genome activation. Overall, these findings suggest that chromatin accessibility influenced by R-loops may represent a system of inherited gene regulation.

Human sperm heads harbor modified YsRNA as transgenerationally inherited non-coding RNAs.

Most epigenetic information is reprogrammed during gametogenesis and early development. However, some epigenetic information persists and can be inherited, a phenomenon that is common in plants. On the other hand, there are increasing examples of epigenetic inheritance in metazoans, especially for small non-coding RNAs. The presence of regulatory important RNAs in oocytes is undisputed, whereas the corresponding RNA payload in spermatozoa and its regulatory influence in the zygote and early embryogenesis is largely enigmatic. For humans, we herein describe small YRNA fragments (YsRNA) as a paternal contribution to the zygote. First, we trace the biogenesis of these YsRNAs from the source YRNAs with respect to the 5' and 3' modifications. Both the length and modifications make these YsRNAs reminiscent of canonical piRNAs that are not derived from piRNA clusters. Second, from the early stages of spermatogenesis to maturation in the epididymis, we observe distinct YsRNA profile dynamics in the male germline. We detected YsRNAs exclusively in mature sperm heads, the precursor of the male pronucleus in the zygote, suggesting an important role of the epididymis as a site for transmitting and modification of epigenetic information in the form of YsRNA between soma and germline in humans. Since this YsRNA-based epigenetic mechanism is effective across generations, we wondered whether this phenomenon of epigenetic inheritance has an adaptive value. Full-length YRNAs bind to Ro60, an RNA chaperone that additionally binds to non-coding RNAs. We described the profiles of non-coding RNAs bound to Ro60 in the human sperm head and detected specific binding profiles of RNA to Ro60 but no YRNA bound to Ro60. We hypothesize that the sperm head Ro60 system is functional. An adaptive phenotype mediated by the presence of a large amount of YsRNA in the sperm head, and thus as a paternal contribution in the zygote, might be related to an association of YsRNA with YRNA that prevents the adoption of a YRNA secondary structure capable of binding to Ro60. We hypothesize that preventing YRNAs from acting as Ro60-associated gatekeepers for misfolded RNAs in the zygote and early development may enhance RNA chaperoning and, thus, represent the adaptive molecular phenotype.

Protein speciation is likely to increase the chance of proteins to be determined in 2-DE/MS.

Multiple spotting due to protein speciation might increase a protein's chance of being captured in a random selection of 2-DE spots. We tested this expectation in new (PXD015649) and previously published 2-DE/MS data of porcine and human tissues. For comparison, we included bottom-up proteomics studies (BU-LC/MS) of corresponding biological materials. Analyses of altogether ten datasets proposed that amino acid modification fosters multispotting in 2-DE. Thus, the number of 2-DE spots containing a particular protein more tightly associated with a peptide diversity measure accounting for amino acid modification than with an alternative one disregarding it. Furthermore, every 11th amino acid was a post-translational modification candidate site in 2-DE/MS proteins, whereas in BU-LC/MS proteins this was merely the case in every 21st amino acid. Alternative splicing might contribute to multispotting, since genes encoding 2-DE/MS proteins were found to have on average about 0.3 more transcript variants than their counterparts from BU-LC/MS studies. Correspondingly, resolution completeness as estimated from the representation of transcript variant-rich genes was higher in 2-DE/MS than BU-LC/MS datasets. These findings suggest that the ability to resolve proteomes down to protein species can lead to enrichment of multispotting proteins in 2-DE/MS. Low sensitivity of stains and MS instruments appears to enhance this effect.

piRNAclusterDB 2.0: update and expansion of the piRNA cluster database.

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.

Genomics and transcriptomics of epizoic Seisonidea (Rotifera, syn. Syndermata) reveal strain formation and gradual gene loss with growing ties to the host.

BACKGROUND: Seisonidea (also Seisonacea or Seisonidae) is a group of small animals living on marine crustaceans (Nebalia spec.) with only four species described so far. Its monophyletic origin with mostly free-living wheel animals (Monogononta, Bdelloidea) and endoparasitic thorny-headed worms (Acanthocephala) is widely accepted. However, the phylogenetic relationships inside the Rotifera-Acanthocephala clade (Rotifera sensu lato or Syndermata) are subject to ongoing debate, with consequences for our understanding of how genomes and lifestyles might have evolved. To gain new insights, we analyzed first drafts of the genome and transcriptome of the key taxon Seisonidea. RESULTS: Analyses of gDNA-Seq and mRNA-Seq data uncovered two genetically distinct lineages in Seison nebaliae Grube, 1861 off the French Channel coast. Their mitochondrial haplotypes shared only 82% sequence identity despite identical gene order. In the nuclear genome, distinct linages were reflected in different gene compactness, GC content and codon usage. The haploid nuclear genome spans ca. 46 Mb, of which 96% were reconstructed. According to ~ 23,000 SuperTranscripts, gene number in S. nebaliae should be within the range published for other members of Rotifera-Acanthocephala. Consistent with this, numbers of metazoan core orthologues and ANTP-type transcriptional regulatory genes in the S. nebaliae genome assembly were between the corresponding numbers in the other assemblies analyzed. We additionally provide evidence that a basal branching of Seisonidea within Rotifera-Acanthocephala could reflect attraction to the outgroup. Accordingly, rooting via a reconstructed ancestral sequence led to monophyletic Pararotatoria (Seisonidea+Acanthocephala) within Hemirotifera (Bdelloidea+Pararotatoria). CONCLUSION: Matching genome/transcriptome metrics with the above phylogenetic hypothesis suggests that a haploid nuclear genome of about 50 Mb represents the plesiomorphic state for Rotifera-Acanthocephala. Smaller genome size in S. nebaliae probably results from subsequent reduction. In contrast, genome size should have increased independently in monogononts as well as bdelloid and acanthocephalan stem lines. The present data additionally indicate a decrease in gene repertoire from free-living to epizoic and endoparasitic lifestyles. Potentially, this reflects corresponding steps from the root of Rotifera-Acanthocephala via the last common ancestors of Hemirotifera and Pararotatoria to the one of Acanthocephala. Lastly, rooting via a reconstructed ancestral sequence may prove useful in phylogenetic analyses of other deep splits.

Pinpointing the PRDM9-PRDM7 Gene Duplication Event During Primate Divergence.

Studies on the function of PRDM9 in model systems and its evolution during vertebrate divergence shed light on the basic molecular mechanisms of hybrid sterility and its evolutionary consequences. However, information regarding PRDM9-homolog, PRDM7, whose origin is placed in the primate evolutionary tree, as well as information about the fast-evolving DNA-binding zinc finger array of strepsirrhine PRDM9 are scarce. Thus, we aimed to narrow down the date of the duplication event leading to the emergence of PRDM7 during primate evolution by comparing the phylogenetic tree reconstructions of representative primate samples of PRDM orthologs and paralogs. To confirm our PRDM7 paralogization pattern, database-deposited sequences were used to test the presence/absence patterns expected from the paralogization timing. In addition, we extended the existing phylogenetic tree of haplorrhine PRDM9 zinc fingers with their strepsirrhine counterparts. The inclusion of strepsirrhine zinc fingers completes the PRDM9 primate phylogeny. Moreover, the updated phylogeny of PRDM9 zinc fingers showed distinct clusters of strepsirrhine, tarsier, and anthropoid degenerated zinc fingers. Here, we show that PRDM7 emerged on the branch leading to the most recent common ancestor of catarrhines; therefore, its origin is more recent than previously expected. A more detailed character evolutionary study suggests that PRDM7 may have evolved differently in Cercopithecoidea as compared to Hominoidea: it lacks the first four exons in Old World monkeys orthologs and exon 10 in Papionini orthologs. Dating the origin of PRDM7 is essential for further studies investigating why Hominoidea representatives need another putative histone methyltransferase in the testis.

Fertility Relevance Probability Analysis Shortlists Genetic Markers for Male Fertility Impairment.

Impairment of male fertility is one of the major public health issues worldwide. Nevertheless, genetic causes of male sub- and infertility can often only be suspected due to the lack of reliable and easy-to-use routine tests. Yet, the development of a marker panel is complicated by the large quantity of potentially predictive markers. Actually, hundreds or even thousands of genes could have fertility relevance. Thus, a systematic method enabling a selection of the most predictive markers out of the many candidates is required. As a criterion for marker selection, we derived a gene-specific score, which we refer to as fertility relevance probability (FRP). For this purpose, we first categorized 2,753 testis-expressed genes as either candidate markers or non-candidates, according to phenotypes in male knockout mice. In a parallel approach, 2,502 genes were classified as candidate markers or non-candidates based on phenotypes in men. Subsequently, we conducted logistic regression analyses with evolutionary rates of genes (dN/dS), transcription levels in testis relative to other organs, and connectivity of the encoded proteins in a protein-protein interaction network as covariates. In confirmation of the procedure, FRP values showed the expected pattern, thus being overall higher in genes with known relevance for fertility than in their counterparts without corresponding evidence. In addition, higher FRP values corresponded with an increased dysregulation of protein abundance in spermatozoa of 37 men with normal and 38 men with impaired fertility. Present analyses resulted in a ranking of genes according to their probable predictive power as candidate markers for male fertility impairment. Thus, AKAP4, TNP1, DAZL, BRDT, DMRT1, SPO11, ZPBP, HORMAD1, and SMC1B are prime candidates toward a marker panel for male fertility impairment. Additional candidate markers are DDX4, SHCBP1L, CCDC155, ODF1, DMRTB1, ASZ1, BOLL, FKBP6, SLC25A31, PRSS21, and RNF17. FRP inference additionally provides clues for potential new markers, thereunder TEX37 and POU4F2. The results of our logistic regression analyses are freely available at the PreFer Genes website (https://prefer-genes.uni-mainz.de/).

Intra-Protein Coevolution Is Increasingly Functional with Greater Proximity to Fertilization.

Intramolecular coevolution of amino acid sites has repeatedly been studied to improve predictions on protein structure and function. Thereby, the focus was on bacterial proteins with available crystallographic data. However, intramolecular coevolution has not yet been compared between protein sets along a gradient of functional proximity to fertilization. This is especially true for the potential effect of external selective forces on intraprotein coevolution. In this study, we investigated both aspects in equally sized sets of mammalian proteins representing spermatozoa, testis, entire body, and liver. For coevolutionary analyses, we derived the proportion of covarying sites per protein from amino acid alignments of 10 mammalian orthologues each. In confirmation of the validity of our coevolution proxy, we found positive associations with the nonsynonymous or amino acid substitution rate in all protein sets. However, our coevolution proxy negatively correlated with the number of protein interactants (node degree) in male reproductive protein sets alone. In addition, a negative association of our coevolution proxy with protein hydrophobicity was significant in sperm proteins only. Accordingly, the restrictive effect of protein interactants was most pronounced in male reproductive proteins, and the tendency of sperm proteins to form internal structures decreased the more coevolutionary sites they had. Both aspects illustrate that the share of outward and thus functional coevolution increases with greater proximity to fertilization. We found this conclusion confirmed by additional comparisons within sperm proteins. Thus, sperm proteins with high hydrophobicity had the lowest proportions of covarying sites and, according to gene annotations, localized more frequently to internal cellular structures. They should therefore be less exposed to postcopulatory forms of sexual selection. Their counterparts with low hydrophobicity had larger proportions of covarying sites and more often resided at the cell membrane or were secreted. At the cellular level, they are thus closer to externally induced forces of postcopulatory selection which are known for their potential to increase substitution rates. In addition, we show that the intermediary status of the testicular protein set in correlation analyses is probably due to a special combination of reproductive and somatic involvements.

The genome, transcriptome, and proteome of the fish parasite Pomphorhynchus laevis (Acanthocephala).

Thorny-headed worms (Acanthocephala) are endoparasites exploiting Mandibulata (Arthropoda) and Gnathostomata (Vertebrata). Despite their world-wide occurrence and economic relevance as a pest, genome and transcriptome assemblies have not been published before. However, such data might hold clues for a sustainable control of acanthocephalans in animal production. For this reason, we present the first draft of an acanthocephalan nuclear genome, besides the mitochondrial one, using the fish parasite Pomphorhynchus laevis (Palaeacanthocephala) as a model. Additionally, we have assembled and annotated the transcriptome of this species and the proteins encoded. A hybrid assembly of long and short reads resulted in a near-complete P. laevis draft genome of ca. 260 Mb, comprising a large repetitive portion of ca. 63%. Numbers of transcripts and translated proteins (35,683) were within the range of other members of the Rotifera-Acanthocephala clade. Our data additionally demonstrate a significant reorganization of the acanthocephalan gene repertoire. Thus, more than 20% of the usually conserved metazoan genes were lacking in P. laevis. Ontology analysis of the retained genes revealed many connections to the incorporation of carotinoids. These are probably taken up via the surface together with lipids, thus accounting for the orange coloration of P. laevis. Furthermore, we found transcripts and protein sequences to be more derived in P. laevis than in rotifers from Monogononta and Bdelloidea. This was especially the case in genes involved in energy metabolism, which might reflect the acanthocephalan ability to use the scarce oxygen in the host intestine for respiration and simultaneously carry out fermentation. Increased plasticity of the gene repertoire through the integration of foreign DNA into the nuclear genome seems to be another underpinning factor of the evolutionary success of acanthocephalans. In any case, energy-related genes and their proteins may be considered as candidate targets for the acanthocephalan control.

Primate piRNA Cluster Evolution Suggests Limited Relevance of Pseudogenes in piRNA-Mediated Gene Regulation.

PIWI proteins and their guiding Piwi-interacting (pi-) RNAs direct the silencing of target nucleic acids in the animal germline and soma. Although in mammal testes fetal piRNAs are involved in extensive silencing of transposons, pachytene piRNAs have additionally been shown to act in post-transcriptional gene regulation. The bulk of pachytene piRNAs is produced from large genomic loci, named piRNA clusters. Recently, the presence of reversed pseudogenes within piRNA clusters prompted the idea that piRNAs derived from such sequences might direct regulation of their parent genes. Here, we examine primate piRNA clusters and integrated pseudogenes in a comparative approach to gain a deeper understanding about mammalian piRNA cluster evolution and the presumed gene-regulatory role of pseudogene-derived piRNAs. Initially, we provide a broad analysis of the evolutionary relationships of piRNA clusters and their differential activity among six primate species. Subsequently, we show that pseudogenes in reserve orientation relative to piRNA cluster transcription direction generally do not exhibit signs of selection pressure and cause weakly conserved targeting of homologous genes among species, suggesting a lack of functional constraints and thus only a minor significance for gene regulation in most cases. Finally, we report that piRNA-producing loci generally tend to be located in active genomic regions with elevated gene and pseudogene density. Thus, we conclude that the presence of most pseudogenes in piRNA clusters might be regarded as a byproduct of piRNA cluster generation, whereas this does not exclude that some pseudogenes nevertheless play critical roles in individual cases.

Temperature-responsive miRNAs in Drosophila orchestrate adaptation to different ambient temperatures.

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.

Effects of different kinds of essentiality on sequence evolution of human testis proteins.

We asked if essentiality for either fertility or viability differentially affects sequence evolution of human testis proteins. Based on murine knockout data, we classified a set of 965 proteins expressed in human seminiferous tubules into three categories: proteins essential for prepubertal survival ("lethality proteins"), associated with male sub- or infertility ("male sub-/infertility proteins"), and nonessential proteins. In our testis protein dataset, lethality genes evolved significantly slower than nonessential and male sub-/infertility genes, which is in line with other authors' findings. Using tissue specificity, connectivity in the protein-protein interaction (PPI) network, and multifunctionality as proxies for evolutionary constraints, we found that of the three categories, proteins linked to male sub- or infertility are least constrained. Lethality proteins, on the other hand, are characterized by broad expression, many PPI partners, and high multifunctionality, all of which points to strong evolutionary constraints. We conclude that compared with lethality proteins, those linked to male sub- or infertility are nonetheless indispensable, but evolve under more relaxed constraints. Finally, adaptive evolution in response to postmating sexual selection could further accelerate evolutionary rates of male sub- or infertility proteins expressed in human testis. These findings may become useful for in silico detection of human sub-/infertility genes.

The pioneering role of PRDM9 indel mutations in tarsier evolution.

PRDM9 is currently the sole speciation gene found in vertebrates causing hybrid sterility probably due to incompatible alleles. Its role in defining the double strand break loci during the meiotic prophase I is crucial for proper chromosome segregation. Therefore, the rapid turnover of the loci determining zinc finger array seems to be causative for incompatibilities. We here investigated the zinc finger domain-containing exon of PRDM9 in 23 tarsiers. Tarsiers, the most basal extant haplorhine primates, exhibit two frameshifting indels at the 5'-end of the array. The first mutation event interrupts the reading frame and function while the second compensates both. The fixation of this allele variant in tarsiers led to hypothesize that de- and reactivation of the zinc finger domain drove the speciation in early haplorhine or tarsiiform primates. Moreover, the high allelic diversity within Tarsius points to multiple effects of genetic drift reflecting their phylogeographic history since the Miocene.

Stop and Go - Waves of Tarsier Dispersal Mirror the Genesis of Sulawesi Island.

The Indonesian island of Sulawesi harbors a highly endemic and diverse fauna sparking fascination since long before Wallace's contemplation of biogeographical patterns in the region. Allopatric diversification driven by geological or climatic processes has been identified as the main mechanism shaping present faunal distribution on the island. There is both consensus and conflict among range patterns of terrestrial species pointing to the different effects of vicariant events on once co-distributed taxa. Tarsiers, small nocturnal primates with possible evidence of an Eocene fossil record on the Asian mainland, are at present exclusively found in insular Southeast Asia. Sulawesi is hotspot of tarsier diversity, whereby island colonization and subsequent radiation of this old endemic primate lineage remained largely enigmatic. To resolve the phylogeographic history of Sulawesi tarsiers we analyzed an island-wide sample for a set of five approved autosomal phylogenetic markers (ABCA1, ADORA3, AXIN1, RAG1, and TTR) and the paternally inherited SRY gene. We constructed ML and Bayesian phylogenetic trees and estimated divergence times between tarsier populations. We found that their arrival at the Proto-Sulawesi archipelago coincided with initial Miocene tectonic uplift and hypothesize that tarsiers dispersed over the region in distinct waves. Intra-island diversification was spurred by land emergence and a rapid succession of glacial cycles during the Plio-Pleistocene. Some tarsier range boundaries concur with spatial limits in other taxa backing the notion of centers of faunal endemism on Sulawesi. This congruence, however, has partially been superimposed by taxon-specific dispersal patterns.

Piwi proteins and piRNAs in mammalian oocytes and early embryos: From sample to sequence.

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.

piRNAs from Pig Testis Provide Evidence for a Conserved Role of the Piwi Pathway in Post-Transcriptional Gene Regulation in Mammals.

Piwi-interacting (pi-) RNAs guide germline-expressed Piwi proteins in order to suppress the activity of transposable elements (TEs). But notably, the majority of pachytene piRNAs in mammalian testes is not related to TEs. This raises the question of whether the Piwi/piRNA pathway exerts functions beyond TE silencing. Although gene-derived piRNAs were described many times, a possible gene-regulatory function was doubted due to the absence of antisense piRNAs. Here we sequenced and analyzed piRNAs expressed in the adult testis of the pig, as this taxon possesses the full set of mammalian Piwi paralogs while their spermatozoa are marked by an extreme fitness due to selective breeding. We provide an exhaustive characterization of porcine piRNAs and genomic piRNA clusters. Moreover, we reveal that both sense and antisense piRNAs derive from protein-coding genes, while exhibiting features that clearly show that they originate from the Piwi/piRNA-mediated post-transcriptional silencing pathway, commonly referred to as ping-pong cycle. We further show that the majority of identified piRNA clusters in the porcine genome spans exonic sequences of protein-coding genes or pseudogenes, which reveals a mechanism by which primary antisense piRNAs directed against mRNA can be generated. Our data provide evidence that spliced mRNAs, derived from such loci, are not only targeted by piRNAs but are also subject to ping-pong cycle processing. Finally, we demonstrate that homologous genes are targeted and processed by piRNAs in pig, mouse and human. Altogether, this strongly suggests a conserved role for the mammalian Piwi/piRNA pathway in post-transcriptional regulation of protein-coding genes, which did not receive much attention so far.

Piwi proteins and piRNAs in mammalian oocytes and early embryos.

Germ cells of most animals critically depend on piRNAs and Piwi proteins. Surprisingly, piRNAs in mouse oocytes are relatively rare and dispensable. We present compelling evidence for strong Piwi and piRNA expression in oocytes of other mammals. Human fetal oocytes express PIWIL2 and transposon-enriched piRNAs. Oocytes in adult human ovary express PIWIL1 and PIWIL2, whereas those in bovine ovary only express PIWIL1. In human, macaque, and bovine ovaries, we find piRNAs that resemble testis-borne pachytene piRNAs. Isolated bovine follicular oocytes were shown to contain abundant, relatively short piRNAs that preferentially target transposable elements. Using label-free quantitative proteome analysis, we show that these maturing oocytes strongly and specifically express the PIWIL3 protein, alongside other, known piRNA-pathway components. A piRNA pool is still present in early bovine embryos, revealing a potential impact of piRNAs on mammalian embryogenesis. Our results reveal that there are highly dynamic piRNA pathways in mammalian oocytes and early embryos.

Tupaia small RNAs provide insights into function and evolution of RNAi-based transposon defense in mammals.

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.