Piwi-interacting RNA (piRNA) expression patterns in pearl oyster (Pinctada fucata) somatic tissues

Piwi-interacting RNAs (piRNAs) belong to a recently discovered class of small non-coding RNAs whose best-understood function is repressing transposable element activity. Most piRNA studies have been conducted on model organisms and little is known about piRNA expression and function in mollusks. We performed high-throughput sequencing of small RNAs extracted from the mantle, adductor muscle, gill, and ovary tissues of the pearl oyster, Pinctada fucata. RNA species with sequences of approximately 30 nt were widely expressed in all tissues. Uridine at the 5′ terminal and protection from β-elimination at the 3′ terminal suggested that these were putative piRNAs. A total of 18.0 million putative piRNAs were assigned to 2.8 million unique piRNAs, and 35,848 piRNA clusters were identified. Mapping to the reference genome showed that 25% of the unique piRNAs mapped to multiple tandem loci on the scaffold. Expression patterns of the piRNA clusters were similar within the somatic tissues, but differed significantly between the somatic and gonadal tissues. These findings suggest that in pearl oysters piRNAs have important and novel functions beyond those in the germ line.