Analysis of the basal chordate Botryllus schlosseri reveals a set of genes associated with fertility
© Rodriguez et al.; licensee BioMed Central. 2014
Received: 10 July 2014
Accepted: 20 November 2014
Published: 26 December 2014
Gonad differentiation is an essential function for all sexually reproducing species, and many aspects of these developmental processes are highly conserved among the metazoa. The colonial ascidian, Botryllus schlosseri is a chordate model organism which offers two unique traits that can be utilized to characterize the genes underlying germline development: a colonial life history and variable fertility. These properties allow individual genotypes to be isolated at different stages of fertility and gene expression can be characterized comprehensively.
Here we characterized the transcriptome of both fertile and infertile colonies throughout blastogenesis (asexual development) using differential expression analysis. We identified genes (as few as 7 and as many as 647) regulating fertility in Botryllus at each stage of blastogenesis. Several of these genes appear to drive gonad maturation, as they are expressed by follicle cells surrounding both testis and oocyte precursors. Spatial and temporal expression of differentially expressed genes was analyzed by in situ hybridization, confirming expression in developing gonads.
We have identified several genes expressed in developing and mature gonads in B. schlosseri. Analysis of genes upregulated in fertile animals suggests a high level of conservation of the mechanisms regulating fertility between basal chordates and vertebrates.
KeywordsAscidian Tunicate Fertility Infertility Gonad formation Germline
B. schlosseri belongs to a subset of ascidians that are colonial, and grow not by increasing in size, but by a lifelong, recurring asexual budding process during which entire bodies are regenerated de novo every week, resulting in an expanding colony of genetically identical individuals, called blastozooids (hereafter shortened to zooids). Zooids have the same body plan as the oozooid, and arrange themselves into star-shaped structures called systems (Figure 1A). A single tadpole larva can give rise to a colony consisting of thousands of systems with tens of thousands of zooids. The zooids are connected by a common, extracorporeal vasculature that ramifies throughout the colony. However, while linked by a common vasculature, the zooids and buds are independent, and zooids or systems can be separated from the colony and continue to grow (called subcloning). The ability to collect pieces of the same genotype at different time points is a powerful characteristic of B. schlosseri as an experimental model for these studies .
While the entire colony can have a lifespan ranging from nine months to several years, the zooids themselves are transient. After completing a two-week developmental program, zooids have a defined lifespan of one week as an adult. During that week, new zooids are regenerating in a process called blastogenesis. This process is coordinated throughout the colony and arranged spatially such that a colony consists of three co-existing asexual generations within each system (Figure 1A and B) [3–7]. The center of the colony is occupied by the zooids, which are actively feeding and, when fertile, sexually reproducing. They are joined peripherally by primary buds, which are completing development. In turn, these are connected to secondary buds, which are in the initial stages of development. Neither primary nor secondary buds feed or sexually reproduce.
The zooid has a lifespan of one week, then dies in a process called takeover. During takeover, all zooids undergo simultaneous apoptosis and are removed by phagocytic cells in the blood. The primary bud then migrates into the newly vacated region of the colony, opens its siphons and becomes a zooid, the secondary bud becomes the primary bud, and a new secondary bud begins developing from a region of the primary bud called the peribranchial epithelium (Figure 1B,C).
Blastogenesis is synchronized throughout a colony and takes 14 days under laboratory conditions. This process can be divided into distinct visual stages, each lasting one day and demarcated by major morphogenetic changes (illustrated in Figure 1B). As outlined in Figure 1 (panel B), a new generation starts as a secondary bud, first visible as a thickening of the peribranchial epithelium of a primary bud (Stage A1), which evaginates and forms a closed vesicle (Stage B2). Next, a series of epithelial invaginations and protrusions appear as development proceeds (Stage C1). After 7 days (Stage D), a takeover event occurs, the secondary bud transitions to a primary bud, and over the next 7 days completes organogenesis, including initiating new secondary bud(s). At day 14, the zooid of the preceding generation undergoes takeover and is removed, the bud migrates to the new zooid position, and opens its siphons.
Ascidians are hermaphrodites, and adult animals can develop both testes and eggs. Following metamorphosis, colonies undergo an average of eight to twelve developmental cycles prior to the development of gametes (sexual maturity) [8, 9]. In addition, natural populations show seasonal fertility, and in the lab cycle in and out of reproductive (fertile) and non-reproductive (infertile) states [10–13]. In turn, sexuality is also plastic, and colonies can be either hermaphroditic, or male only (with no oocytes visible) (Figure 1D-E).
When a colony is fertile, development of the gametes is synchronized with somatic development [10–13]. Germline stem cells (GSCs) seed the secondary bud (Figure 1B) and begin to proliferate. However, development of oocytes and testis are not linked. Testes develop in situ, and are visible in the secondary bud at stage B2, and completing development approximately 10 days later, when that secondary bud becomes a zooid [10–13]. In contrast, oocytes appear to take several cycles to mature (described in detail, below), and oocytes at different stages of development (pre-vitellogenic and vitellogenic) will seed the secondary bud along with the GSCs.
B. schlosseri is ovoviviparous and mature eggs are about 300 μm in diameter, and are surrounded by a series of cell layers. These include a layer of outer follicle cells (OFC) a layer of inner follicle cells (IFC) and test cells (TC) . Each egg resembles a small ovary consisting of an oocyte, the envelopes (OFC, IFC and TC) and its own oviduct. Oocyte precursors first appear in secondary buds, ripen in primary buds, and ovulate when the primary buds are transitioning to mature zooids during takeover [14–16]. Oocyte development takes place in the primary bud from stages A1 to C2 and can be described as 5 stages: At stage 1, oocyte precursors are 10 μm in diameter, with a high nuclear to cytoplasmic ratio. At stage 2, the oocyte diameter increases up to 50 μm. At stage 3 oocyte diameter is approximately 80 μ. At stage 4, the oocyte diameter is 120 μm, and vitellogenesis, the process of secretion of nutrients into the oocyte by follicle cells, is initiated. The diameter of stage 5 oocytes is 220 μm, vitellogenesis is almost completed, with yolk granules becoming larger. At the end of the blastogenic cycle one to four oocytes mature and are ovulated into the peribranchial space during takeover (stage D). Usually a number (approximately 6 per primary bud) of oocytes remain arrested at stage 2 (in a previtellogenic stage) which can move from the zooid into the blood stream and reach a new primary bud where they can potentially mature in successive blastogenic generations .
To better understand the mechanisms of gonad differentiation and fertility, we used mRNA-seq to identify genes that are differentially expressed between fertile and infertile colonies. Using differential expression analysis combined with in situ hybridization, we have identified a number of genes associated with fertility in Botryllus, particularly genes expressed by follicle cells - surrounding both testis and oocyte precursors - that may drive maturation. This study is the first step to comprehensively characterize the mechanisms of gonad maturation and fertility in Botryllus schlosseri.
Results and discussion
The recent publication of a draft genome sequence for Botryllus schlosseri has identified homologs of vertebrate genes known for their function in lymphoid-mediated immunity as well as eye, ear and heart function . Although this is an important step towards establishing Botryllus as a model organism, an analysis of gene expression dynamics throughout the blastogenic cycle is an important tool to further strengthen this non-conventional model system. One of the limitations of a newly established model organism is that we can infer the function of most proteins in Botryllus only based on comparison to homologs in other species or analysis of conserved protein domains.
Number of differentially expressed EST’s (10% FDR) at each stage of the blastogenic cycle and number of those EST’s with human homologs
Stage of blastogenesis
Number of differetially expressed genes (10%FDR)
Number of human homologs
Differentially expressed genes upregulated in fertile versus infertile colonies and their human and/or ascidian homologs
DE at Stage(s) of Blastogenesis
EST 5-11 contig name
Homo sapiens or Ascidian homologue
Function or Expression
A1*, A2, B1, C2 and D
Regulates gene transcription through chromatin remodeling, facilitating expression of fertility-related genes.
A1, B1*, B2, C1, C2 and D
Involved in secretion of nutrients into the oocyte by follicle cells
A1*, B1, and B2
F-box only protein 24 isoform 3
Substrate-recognition component of the SCF (SKP1-CUL1-F-box protein)-type E3 ubiquitin ligase complex
A1 and B1*
Erythrocyte band 7 integral membrane protein isoform a
Expressed on red blood cells.
A1 and B1*
Testis-specific serine/threonine-protein kinase 2
Expressed in both mouse and human sperm.
Microfibrillar-associated protein 4
Involved in photoprotection of the skin.
A1*, B1 and B2
Kelch-like protein 10
Encodes a component of intercellular bridges in Drosophila egg chambers.
A1, A2, B1*, and B2
Histone H2B type 1-L
Required to condense chromatin in sperm, removes heterochromatin marks, facilitating gene expression
A1, A2, B1* and B2
Creatine kinase S-type mitochondrial
Is expressed in cardiac and striated muscle.
A1*, B1 and D
Mammalian tolloid (BMP-1) is a proteinase involved in ovarian tissue remodeling
A1*, A2, B1 and B2
Testis-specific serine/threonine-protein kinase 1
Testicular germ cell-specific expression, possible role at and after the meiotic phase of spermatogenesis
A1, B1, B2*, C2 and D
Low density lipoprotein-related protein 1
Mediates endocytosis of cholesterol-rich LDL. In mice, LDLR is expressed in germ cells
B1, B2 and C2*
Zona Pellucida sperm binding protein-1
Involved in oocyte development, protection, sperm binding.
A1, A2, B1* and D
Expressed in the oolema of oocytes in hamsters and humans, as well as in sperm following the acrosomal reaction
A1 and B1*
Ubiquitin-conjugating enzyme E2 R2
Ubiquitin-protein ligase complex, mediates ubiquitination and proteosomal degradation of target proteins during spermatogenesis
A1 and B1*
NAD+-specific isocitrate dehydrogenase beta precursor
Possible role in oxidation of isocitrate to alpha-ketoglutarate in the citric acid cycle.
A1* and B1
Outer dense fiber of sperm tails 3-like 2
Structural protein surrounding the axoneme in both the middle piece and principal piece of the sperm tail.
A1* and B1
CD81 antigen (Tetraspanin-8)
Expressed by granulosa cells surrounding the oocyte in mice. Reduced fertility in CD81-/- null mice, possible role in acrosomal reaction
A1*, B1, B2
Mediates cleavage of proteoglycans during the release of the oocyte in mammals
ADAM metallopeptidase domain 12
ADAMs play roles in spermatogenesis and sperm function, potentially by effecting maturation of sperm and their adhesion and migration
A1, A2*, and B1
Involved in cell cycle.
A1 and B1*
Expressed in human testes,pecifically in spermatogonia, spermatocytes, spermatids, Leydig and Sertoli cells
A1, B1*, B2
Fibrous sheath-interacting protein 2
Expressed in late spermatocyte development.
A1*, B1 and D
Facilitates the gliding of the developing gonad along epithelial basement membranes and germline cellularization.
A1*, A2 and B1
Testis-specific serine/threonine-protein kinase 6
Essential for spem production and function
A1*, B1 and D
Expressed in sensory epithelia of the inner ear.t Has also been classified as a testis-selective cancer/testis (CT) gene
A1, A2, B1*, C2 and D
Mouse homozygous-null embryos showed folding of the egg cylinder, leakage of maternal red blood cells into the yolk sac cavity.
A1*, A2 and D
ATP-binding cassette sub-family B member 5 isoform 1
In Drosophila, ABC transporters efflux prenylated peptides out of somatic gonadal precursors to serve as chemoattractants for migrating germ cells
Retinol dehydrogenase 12
In Botryllus schlosseri retinoic acid signaling is involved in gonad formation
A1, B1* and D
Tubulin beta-4B chain
Tubulin is the major constituent of microtubules; post-translational modification bu monoglycalation specifi incorporation into axonemes.
A1 and B1*
Ciliary dynein heavy chain 9
Possible role in sperm development or motility
A1*, A2, B1 and D
Cyclic nucleotide-gated olfactory channel
Expressed in the flagellum of mature sperm, involved in sperm movement
Involved in cell cycle regulation
The number of differentially expressed genes varied greatly between stages of the blastogenic cycle, from a minimum of 7 at stage C1 to a maximum of 647 at stage B1 (Table 1).
In Botryllus spermatogenesis begins at stage A1 in the testis of the zooid and continues until stage B2 when mature sperm is released through the oral siphon out into the water . Consistent with this timing, we found that based on homology, a significant portion of the genes upregulated in stages A1 and A2 are functionally related to sperm formation (Additional file 3: Table S1, Additional file 4: Table S2, Additional file 14: Table S12, Additional file 15: Table S13, Additional file 16: Table S14, Additional file 17: Table S15, Additional file 18: Table S16, Additional file 19: Table S17, Additional file 20: Table S18, Additional file 21: Table S19). Since most of these homologs are involved in spermatogenesis in vertebrates (references in Table 2), these data suggest some degree of conservation of the mechanisms governing spermatogenesis between ascidians and vertebrates.
To investigate the spatial and temporal expression pattern of DE genes, we performed quantitative RT-PCR (qPCR) and fluorescent in situ hybridization (FISH) for genes that were selected from the most highly DE genes with human homologs: p-selectin (slep), otoancorin (otoa), tetraspanin-8 (tspan8), low-density lipoprotein receptor (ldlr), testis-specific serine/threonine kinase 1 (tsk1) and 2 (tsk2) and genes with known roles in fertility in other organisms: estradiol 17 β dehydrogenase-8 (hsd17β8), zona pellucida sperm binding protein 1 (zp1), vitellogenin-1 (vtg1), tolloid-like-1 (tll1). All of these genes were analyzed in fertile colonies at each stage of the blastogenic cycle by qPCR (Additional file 22: Figure S2).
Testis-specific serine/threonine protein kinases phosphorylate myelin basic protein and histones in vitro. The testicular germ cell-specific expression of the mouse homologs suggests that these genes play an important role at and after the meiotic phase of spermatogenesis . In Botryllus, Testis-specific serine/threonine protein kinases, tsk-1, tsk-2 and tsk-6 are upregulated in fertile animals based on our DE analysis (tsk-1: 203-fold upregulated at stage A1; tsk-2: 462-fold upregulated at stage B1; tsk-6: 57-fold upregulated at stage A1). qPCR analysis showed that both tsk-1 and tsk-2 are highly upregulated between stages A2 and B2, and in situ hybridization indicates that the mRNA is expressed on developing testes (Figure 2) consistent with a crucial role of these genes in spermatogenesis.
In mammals, Otoancorin is specific to sensory epithelia of the inner ear but in humans it has also been classified as a testis-selective cancer/testis gene . Based on our DE analysis, the Botryllus homolog of otoancorin is upregulated 54-fold in fertile animals at stage A1. By in situ hybridization, otoancorin shows strong expression from stage A1 through B1 in the testes of the primary bud (pb) (Figure 2A), indicating that in Botryllus, otoancorin might be involved in spermatogenesis. By qPCR, its expression level is highest in A1, and decreases as the blastogenic cycle continues. Interestingly, testes of the zooid also show signal for otoancorin mRNA, but expression is restricted to the peripheral cells of the testes (Figure 2B).
Tetraspanin/CD9-like is an oocyte factor required for sperm-oocyte fusion . Consistent with these findings in mammals, in situ hybridization shows expression of the Botryllus homolog of tetraspanin-8 on stage 2 oocytes. The ADAM-integrin-tetraspanin complex, known to constitute a network of membrane microdomains called the tetraspanin web, is potentially involved in the migration of prespermatogonia from the center to the periphery of the testicular cords and in the reinitiation of mitotic activity during the initial wave of spermatogenesis . By qPCR, expression of tetraspanin-8 in Botryllus is highest at stage A1, coinciding with the initiation of spermatogenesis. Furthermore, our DE analysis revealed a 125-fold upregulation of tetraspanin-8 in fertile animals at stage A1. In situ hybridization showed cells expressing tetraspanin-8 mRNA in the testes of the primary bud, and suggests that the role of Tetraspanin proteins in spermatogenesis is likely to be conserved between Botryllus and vertebrates.
Vitellogenins are lipophosphoglycoproteins that are produced under female hormonal control in large central organs (fat body in insects; liver in higher animals) and are transported in the circulation to the female gonads . In Botryllus, a vitellogenin homolog is 927-fold upregulated in fertile animals at stage B1. Two novel isoforms of vitellogenin, both of which possess vWF-D and CT domains but not a lipovitellin or phosvitin domain, are expressed in the gonad of the ascidian Halocynthia roretzi. In situ hybridization revealed that mRNAs of these proteins are specifically expressed in oocytes and test cells, accessory cells in the perivitelline space of ascidian eggs. Immunohistochemistry showed that these proteins are localized around the surface of test cells in immature oocytes . In Botryllus, vitellogenin expression is highest at stage B2 (Additional file 22: Figure S2), and, consistent with the findings in Halocynthia, in situ hybridization shows expression on maturing eggs of the primary bud (Figure 2).
Interestingly, vitellogenin mRNA is also expressed by follicle cells on the testes of both the primary bud and the zooid (Figure 2). These cells are restricted to the periphery of the testes and seem to surround them, in contrast to the mRNA’s of otoancorin, testis-specific serine/threonine-protein kinase 1 and tetraspanin-8 (Figure 2), which are expressed by cells in the inner compartment of the testes.
To date, all molecularly characterized vitellogenin receptors belong to the low-density lipoprotein receptor (ldlr) supergene family. Ldlr receptors localized in coated pits on the surface of growth-competent oocytes are able to accumulate in the yolk high concentrations of vitellogenin and other ligands they recognize . A Botryllus schlosseri homolog of ldlr is upregulated 199-fold in fertile animals at stage B2, and qPCR analysis shows that expression of ldlr in Botryllus is highest at stage A2 (Additional file 22: Figure S2). By in situ hybridization, we found that ldlr is expressed in follicle cells of maturing eggs in the primary bud (Figure 2), suggesting that it might function as a vitellogenin-receptor during oocyte maturation.
In Halocynthia roretzi, a 120-kDa transmembrane protein with a zona pellucida domain and a 13 EGF-like repeats (HrVC120), has been shown to be the receptor in the vitelline coat that binds to sperm . Ciona intestinalis has several genes related to HrVC120 that are exclusively expressed in developing oocytes but not in eggs . We found a gene containing a Zona Pellucida domain, zp-1, that is highly upregulated in fertile animals (197-fold change at stage C2) and highly expressed at stage A1 (Additional file 22: Figure S2). By in situ hybridization, we found that zp-1 is expressed by the follicle cells that surround the egg in the primary bud. To our knowledge this is the first report of a gene with a Zona Pellucida domain in a colonial ascidian. This finding reflects the close relationship of these animals with vertebrates and the fact that this domain has been maintained in deuterostomes suggests that it is indispensable for sperm-egg interaction [27, 33].
Mammalian Tolloid (Bmp-1) is a proteinase involved in ovulation. During ovarian tissue remodeling, it contributes to the maturation of procollagen molecules and the deposition of collagen fibrils . In Botryllus, a tolloid-like protein 1 homolog is expressed on follicle cells of maturing oocytes of the primary bud and is upregulated 261-fold in fertile animals at stage A1 (Figure 2).
A Botryllus schlosseri homolog of estradiol 17 β dehydrogenase-8 (hsd17β8) is highly upregulated in fertile animals (5.4-fold at stage B1), and as expected, given its involvement in estrogen metabolism , its mRNA localizes to oocytes at stage 2 (Figure 2).
Measurements of p-selectin positive cells, mean and standard deviation (SD) in the primary bud throughout the blastogenic cycle
Range of diameters (μm)
Unexpectedly, we found that stage 2 oocytes were positive for the mRNA of testis-specific serine/threonine-protein kinase 2 gene at stage B1 and B2 of the blastogenic cycle. Moreover, a strong nucleolus like signal was observed in about 30-40% of positive cells. This finding could indicate a previously unknown function of this gene in oogenesis.
The set of fertility-related genes we have identified represents the first reported molecular markers for a variety of Botryllus schlosseri anatomical structures. Specifically, we found transcripts expressed by the testes (both developing and mature testes), by follicle cells surrounding both maturing oocytes and testes, by oocytes at stage 2 and by the siphons (oral and excurrent). These markers will undoubtedly serve as a valuable resource for other researchers examining fertility in Botryllus.
To obtain a broader perspective and explore the biological roles of the genes identified in our DE analysis, we performed a gene ontology (GO) analysis on the human homologs of differentially expressed genes (Table 1) at each stage of the blastogenic cycle (Human genes were used as a proxy in our GO analysis because Botryllus ESTs lack GO annotation) [39–41]. Specifically, annotation based on biological processes provided significant insight into functional changes related to fertility throughout the blastogenic cycle, and are discussed here. Other GO annotations (All annotations and molecular function) are included in Additional file 14: Table S12, Additional file 15: Table S13, Additional file 16: Table S14, Additional file 17: Table S15, Additional file 18: Table S16, Additional file 19: Table S17, Additional file 20: Table S18, Additional file 21: Table S19, Additional file 23: Table S20, Additional file 24: Table S21, Additional file 25: Table S22, Additional file 26: Table S23, Additional file 27: Table S24, Additional file 28: Table S25, Additional file 29: Table S26, Additional file 30: Table S27, Additional file 31: Table S28, Additional file 32: Table S29, Additional file 33: Table S30, Additional file 34: Table S31, Additional file 35: Table S32, Additional file 36: Table S33, Additional file 37: Table S34, Additional file 38: Table S35, Additional file 39: Table S36, Additional file 40: Table S37 and Additional file 41: Table S38. Furthermore, we conducted Inter Pro analysis to find conserved protein domains (Additional file 17: Table S15, Additional file 21: Table S19, Additional file 26: Table S23, Additional file 30: Table S27, Additional file 37: Table S34, and Additional file 41: Table S38).
At stages A1 and A2, 592 (290 with human homologs) and 136 (70 with human homologs) genes are upregulated in fertile animals, respectively. GO analysis revealed that genes related to male meiosis, spermatogenesis, cilium, proteolysis, microtubule based movement, meiosis, spermatid development, cilliary or flagellar motility, female gamete generation and male meiosis are upregulated (Additional file 14: Table S12, Additional file 15: Table S13, Additional file 16: Table S14, Additional file 18: Table S16, Additional file 19: Table S17 and Additional file 24: Table S18), correlating with the fact that spermatogenesis begins in this timeframe. Inter Pro analysis revealed that these genes contain domains such as EGF, LDLR, BMP, Selectin and Meiotic recombination Spo11, among others (Additional file 21: Table S19 and Additional file 26: Table S23). These results are congruent with previous non-molecular observations indicating that gametogenesis starts during these stages.
At stages B1 and B2, 647 and 125 genes are DE with 304 and 59 human homologs respectively. Genes related to male meiosis, spermatogenesis, spermatid development, microtubule-based movement, cell differentiation, motility, filipodium assembly, flagellar motility, cell adhesion, female gamete generation and apoliprotein binding are upregulated in fertile animals (Additional file 23: Table S20, Additional file 24: Table S21, Additional file 25: Table S22, Additional file 27, Table S24, Additional file 28: Table S25 and Additional file 29: Table S26). Inter Pro analysis revealed that these genes contain domains such as EGF, Kringle, CUB, Speract/Scavenger Kelch, Calcium binding, ADAM, and Peptidase S1/S6/M12A, BMP1/Tolloid, Zona pellucida sperm-binding protein and LDLR among others (Additional file 26: Table S23 and Additional file 30: Table S27). These results reflect the complex biological processes occurring during these stages, which also have the highest number of DE genes, and are consistent with the timing of sperm maturation and release.
At stage C1 and C2, 7 and 61 genes are DE with 2 and 32 human homologs respectively. The predominant GO terms at these stages included oocyte maturation, oocyte growth, lipoprotein metabolic process, regulation of cholesterol transport, positive regulation of estrogen receptor, extracellular matrix disassembly, and tissue remodeling, among others (Additional file 31: Table S28, Additional file 32: Table S29, Additional file 34: Table S31, Additional file 35: Table S32 and Additional file 36: Table S33). The Inter Pro analysis for this stage revealed domains such as: Apple-like, LDLR, Kringe, Aquaporin, Tetraspanin, EGF, Zona Pellucida sperm-binding protein, Peptidase S1/S6/M12A and FY rich, among others (Additional file 33: Table S30 and Additional file 37: Table S34). At these stages maturing oocytes rapidly increase in size and are preparing to move into the branchial basket of the zooid.
Finally at stage D 221 genes are DE (with 122 human homologs). GO analysis revealed the presence of annotations including cell adhesion, carbohydrate metabolism, regulation of developmental process and response to corticosteroids stimulus, proteolysis, lipoprotein transport, cholesterol transport, estrogen receptor signaling pathway, oocyte growth (Additional file 38: Table S35, Additional file 39: Table S36 and Additional file 40: Table S37) and the Inter Pro analysis of the genes upregulated at this stage revealed the following protein domains: EGF, Notch, CUB, LDLR, Selectin, and Peptidases S1/S1A/S6, among others (Additional file 41: Table S38). At this stage mature oocytes begin to transition from the peribranchial cavity into the branchial basket of the new zooid.
In summary, our GO analysis has revealed a number of genes with potential involvement in spermatogenesis at stages A1-B2, correlating with what has been described for Botryllus schlosseri. The Botryllus homologs of these genes are promising candidates for future studies on testes maturation and sperm production. Homologs of genes with reported roles in lipid and steroid metabolism are upregulated at stages B1-D, representing key processes involved in oocyte maturation in other species that may be conserved in Botryllus[42–44].
This manuscript represents the first step in the comprehensive characterization of gonad differentiation and fertility in Botryllus schlosseri. We have analyzed the transcriptomes of both infertile and fertile colonies and characterized the spatial and temporal expression of highly expressed genes involved in fertility. Our findings suggest a high level of conservation of the mechanisms regulating fertility between basal chordates and vertebrates. Some of the genes characterized in this study were found to have unexpected spatial localization, such as testis serine/threonine kinase 2, which is expressed by egg precursors. Another example is vitellogenin, which localizes to both follicle cells surrounding the egg and equivalent cells surrounding the testes.
Interestingly, known germline markers such as vasa, nanos and piwi were not significantly differentially expressed between fertile and infertile colonies, most likely because both juveniles and infertile colonies have immature germ cells expressing these genes. We have identified specific markers for both developing gonads (and mature testes) throughout the blastogenic cycle in zooids, primary buds and secondary buds. Some of these markers are now being used to investigate the phenomenon of induced juvenile fertility in Botryllus schlosseri (Di Maio et al. unpublished).
In future studies, we will investigate differential gene expression between male and hermaphroditic colonies, which will allow us to characterize genes used exclusively in male or female gamete development and identify genes governing the process by which Botryllus schlosseri assumes a given reproductive state.
All mariculture procedures have been described previously (Boyd et al., 1986). Briefly, we collected pregnant Botryllus schlosseri colonies from the Santa Barbara Marina, California. Individual colonies were tied to 3 × 5-cm glass slides and placed 5 cm opposite another glass slide (called the settlement slide) in a slide rack. The slide rack was placed into an aquarium, and within a few days the tadpoles hatched, swam to the settlement slide, and metamorphosed into the adult body plan (oozooid). Over 95% of the tadpoles hatch and immediately swim to and settle on the juxtaposed slide . Animals are reared in 5 L tanks supplemented with food in suspension daily, and food is not limiting . Only isolated genotypes reared in our mariculture facility were used for this study.
To explore only the genes involved in gonad formation and fertility, mRNA seq analysis was performed at every stage of the blastogenic cycle (A1, A2, B1, B2, C1, C2 and D, each stage as described by bud morphogenesis based on in vivo and histological features ) on a total of 3 fertile samples (SB801, SB841 and Mix genotypes SB802a, SB804, SB722 ) and 3 infertile samples (SB802d, SB-BB-001 and Mix genotypes SB831, SB645, SB714, SB841, SB842 ).
RT-PCR, cloning, and quantitative PCR
Total RNA was isolated from whole colonies using Nucleo-Spin II columns from Macherey Nagel, and mRNA was isolated using the NEB Magnetic Bead Isolation kit. cDNA was made using Superscript II from Invitrogen. For RT-PCR, we used Clontech Advantage following the manufacturers′ recommendations for PCR conditions.
PCR products were isolated using Qiagen columns, then cloned into the Promega Easy-T vector and transformed using competent NEB5α from NEB. Single colonies were picked and inserts amplified using colony PCR and Clontech Advantage. PCR products were prepared for sequencing then sent for Sanger sequencing at the UC Berkeley Sequencing Facility.
Quantitative PCR (QPCR) analysis was done as described previously . Briefly, Q-PCR was carried out using a LightCycler 480 II (Roche) and LightCycler DNA Master SYBR Green I detection (Roche, 12015099001). The thermocycling profile was as follows: 5 minutes at 95°C, 45 cycles of 95°C for 10 sec, 50-60°C for 10 sec, and 72°C for 10 sec. All gene expression data was normalized to elongation factor 1-α (EF1- α) as a reference house keeping gene and reported as relative expression using the 2-ΔΔCt method.
Three biological and technical replicates were used for each gene.
Primers sequence for QPCR
Primers sequence 5′-3′
Total RNA was isolated from whole colonies using the Nucleo-Spin RNA II kit (MN). Libraries were prepared and sequenced at the USC Epigenome Center using kits from Illumina following the manufacturers’ instructions. Paired-end libraries were generated for each sample and sequenced with an Illumina Hi-Seq 2000. The data supporting the results of this article have been deposited in NCBI’s Gene Expression Omnibus [48, 49] and are accessible through GEO Series accession number GSE62112 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE62112). After RNASeq each paired end (PE) library was first checked for quality control (QC) using the software FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ ) with an average score of 28 across all bases (Illumina 1.5 encoding). Using FASTQ Trimmer we removed the adapters of the sequences by removing the first 12 bases of each read. All reads pairs passing QC were mapped using BOWTIE (2012 version 0.12.7)  to the public EST database Bot_asmb assembly (04.05.2011, A. Gracey) (http://octopus.obs-vlfr.fr/public/botryllus/blast_botryllus.php). Bot_asmb assembly (04.05.2011, A. Gracey), consisting of 50,107 contigs, is an assembly of two cDNA libraries generated from a mixture of Botryllus schlosseri genotypes at each stage of the blastogenic cycle and representing both infertile and fertile states. The libraries used to construct Bot_asmb assembly (04.05.2011, A. Gracey) consist of: (1) 10,000 arrayed clones sequenced from 5’ to 3’ ends using the Sanger method (200,000 sequences total with an average length of 950 base pairs) and (2) 454 sequencing of the same RNA with an average length of 400 base pairs. The number of reads mapping to each EST were obtained with Sam2Counts (Samtools version 0.1.18) and differential expression analysis was performed with DESeq 1.10.1 using triplicates for the analysis and a false discovery rate of 10% . Fragments per kilobase per million (FPKM) were calculated using RSEM 1.2.3 for each paired-end library for each genotype at each stage of the blastogenic cycle for both infertile and fertile colonies. To identify putative homologs we blasted (BLASTX) to the non-redundant protein databases of Homo sapiens Mus musculus, Drosophila melanogaster, Ciona intestinalis and Danio reirio with an E-value cutoff of 1 × 10-4.
All bioinformatics analysis was performed on a Desktop computer running Ubuntu 12.10 and equipped with 16 GB of RAM.
Gene Ontology and InterPro analysis
The list of human homologs of differentially expressed Botryllus contigs was submitted to GeneCodis 3, an on-line modular enrichment tool [39–41]. For this analysis, the following annotations were selected: GO Biological Process, GO Molecular Function and InterPro Motifs. The statistical parameters for these analyses were as follows: First, for co-occurrence analyses of annotations, a minimum support of 3 genes was required. Second, the statistical method to compute p-values was the hypergeometric test. Finally, to correct p-values for multiple hypothesis testing, FDR estimation was utilized. Of particular interest to this study were processes and pathways with a known role in fertility and gametogenesis.
Fluorescent in situhybridization (FISH)
Fragments of Botryllus schlosseri Estradiol 17-β-dehydrogenase 8 (hsd17b8), P-Selectin (selp), Otoancorin (otoa), Testis-specific serine/threonine-protein kinase 2 (tsk2), Vitellogenin (vtg1), Tetraspanin 8 (tspan8), Zona-pellucida sperm binding protein 1 (zp1), Low-Density Lipoprotein receptor (ldlr), Testis-specific serine/threonine-protein kinase 1 (tsk1), Tolliod-like 1 (tll1), ranging in size from 200 bp – 1 kb were amplified with Advantage cDNA polymerase (Clontech) and cloned into the pGEM-T Easy vector (Promega). SP6 or T7 RNA polymerase (Roche) were used to generate antisense RNA probes labeled with digoxigenin (Roche) or DNP (Perkin Elmer).
Single systems of fertile or infertile Botryllus were fixed with 4% formaldehyde in 0.1 M MOPS pH 7.5, 0.5 M NaCl for 3 hours and then transferred to methanol. To eliminate pigmentation, animals were bleached in 6% H2O2/methanol and then stored in methanol at -20°C. After rehydration, animals were permeabilized with proteinase K (10 μg/ml for 30 min) and then post-fixed with 4% formaldehyde. Prehybridization was carried out for 4 hours at 65°C in hybridization buffer (65% formamide, 5X SSC, 1X Denhardt's solution, 0.1% Tween-20, 5 mg/ml torula yeast RNA, 50 μg/ml heparin), followed by hybridization with DIG- and DNP-labeled probes in hybridization buffer overnight at 65°C. After washing off unbound probes, DIG-labeled probe was detected with an anti-DIG HRP-conjugated antibody (Roche) and TSA Plus fluorescent tyramide conjugate (Perkin Elmer). For double labeling, the HRP-conjugated anti-DIG antibody was then inactivated with 2% H2O2 in PBS with 0.1% Triton X-100 and DNP-labeled probes were detected with anti-DNP HRP-conjugated antibody (Perkin Elmer) and TSA Plus fluorescent tyramide conjugate (Perkin Elmer). After washing in PBS, specimens were flat-mounted with Vectashield (Vector Labs) and imaged on an Olympus Fluoview 1000 Spectral confocal microscope equipped with a 40X oil immersion objective. All FISH experiments were performed in triplicates. Only consistent reproducible staining throughout the whole-mount colonies was reported.
Primer sequences for in situ hybridization probes
Primers sequence used for each gene 5′-3′
Measurement of stage 2 oocytes
After FISH (above) the diameter of selp positive cells were measured using FIJI image processing software. Each primary bud (average of five buds per experiment) with positive cells (average of ten positive oocytes stage 2) were measured of total of three colonies at each stage of the blastogenic cycle. Standard deviation was calculated using Microsoft Office Excel 2011.
The authors would like to acknowledge the Grant Funding NIH R01-AG037699 and the NIR-MCDB Microscopy Facility and the Spectral Laser Scanning Confocal supported by the Office of The Director, National Institutes of Health of the NIH under Award # S10OD010610. We thank Dr. Susannah Kassmer for critical reading of the manuscript and helpful discussions. We thank Amy N. Marten for experiments on the early stages of this project. We would like to thank Mike Caun for his expert care of the De Tomaso laboratory mariculture facility.
- Berrill NJ: The developmental cycle of Botrylloides. Q J Microsc Sci. 1947, 88 (Pt 4): 393-407.PubMedGoogle Scholar
- Chadwick-Furman NE, Weissman IL: Life history plasticity in chimaeras of the colonial ascidian Botryllus schlosseri. Proc Biol Sci. 1995, 262 (1364): 157-162. 10.1098/rspb.1995.0190.PubMedView ArticleGoogle Scholar
- Manni L, Zaniolo G, Cima F, Burighel P, Ballarin L: Botryllus schlosseri: a model ascidian for the study of asexual reproduction. Dev Dyn. 2007, 236 (2): 335-352. 10.1002/dvdy.21037.PubMedView ArticleGoogle Scholar
- Tiozzo S, Christiaen L, Deyts C, Manni L, Joly JS, Burighel P: Embryonic versus blastogenetic development in the compound ascidian Botryllus schlosseri: insights from Pitx expression patterns. Dev Dyn. 2005, 232 (2): 468-478. 10.1002/dvdy.20250.PubMedView ArticleGoogle Scholar
- Chadwick-Furman NE, Weissman IL: Life histories and senescence of Botryllus schlosseri (Chordata, Ascidiacea) in Monterey Bay. Biol Bull. 1995, 189 (1): 36-41. 10.2307/1542199.PubMedView ArticleGoogle Scholar
- Stoner DS, Weissman IL: Somatic and germ cell parasitism in a colonial ascidian: possible role for a highly polymorphic allorecognition system. Proc Natl Acad Sci U S A. 1996, 93 (26): 15254-15259. 10.1073/pnas.93.26.15254.PubMed CentralPubMedView ArticleGoogle Scholar
- Tiozzo S, De Tomaso AW: Functional analysis of Pitx during asexual regeneration in a basal chordate. Evol Dev. 2009, 11 (2): 152-162. 10.1111/j.1525-142X.2009.00316.x.PubMedView ArticleGoogle Scholar
- Sabbadin A: Self- and cross-fertilization in the compound ascidian Botryllus schlosseri. Dev Biol. 1971, 24 (3): 379-391. 10.1016/0012-1606(71)90087-X.PubMedView ArticleGoogle Scholar
- Zaniolo G, Sabbadin A, Resola C: Dynamics of the colonial cycle in the ascidian, Botryllus schlosseri. The fate of isolated buds. Acta Embryol Exp (Palermo). 1976, 2: 205-213.Google Scholar
- Johnson SL, Yund PO: Remarkable longevity of dilute sperm in a free-spawning colonial ascidian. Biol Bull. 2004, 206 (3): 144-151. 10.2307/1543638.PubMedView ArticleGoogle Scholar
- Mukai H: Comparative studies on the structure of reproductive organs of four botryllid ascidians. J Morphol. 1977, 152 (3): 363-380. 10.1002/jmor.1051520307.PubMedView ArticleGoogle Scholar
- Mukai H, Watanabe H: Studies on the formation of germ cells in a compound ascidian Botryllus primigenus Oka. J Morphol. 1976, 148 (3): 377-362.PubMedView ArticleGoogle Scholar
- Sabbadin AaZ G: Sexual Differentiation and Germ Cell Transfer in the Colonial Ascidian Botryllus schlosseri. J Exp Zool. 1979, 207: 279-301.Google Scholar
- Manni L, Zaniolo G, Burighel P: Egg Envelop Cytodifferentiation in the Colonial Ascidian Botryllus schlosseri (Tunicata). Acta Zoologica. 1993, 74 (2): 103-113. 10.1111/j.1463-6395.1993.tb01226.x.View ArticleGoogle Scholar
- Manni L, Zaniolo G, Burighel P: An unusual membrane system in the oocyte of the ascidian Botryllus schlosseri. Tissue Cell. 1994, 26 (3): 403-412. 10.1016/0040-8166(94)90023-X.PubMedView ArticleGoogle Scholar
- Manni LZ G, Burighel P: Ultrastructural Study of Oogenesis in the Compound Ascidian Botryllus schlosseri (Tunicata). Acta Zoologica. 1994, 75 (2): 101-113. 10.1111/j.1463-6395.1994.tb01115.x.View ArticleGoogle Scholar
- Voskoboynik A, Neff NF, Sahoo D, Newman AM, Pushkarev D, Koh W, Passarelli B, Fan HC, Mantalas GL, Palmeri KJ, Ishizuka KJ, Gissi C, Griggio F, Ben-Shlomo R, Corey DM, Penland L, White RA, Weissman IL, Quake SR: The genome sequence of the colonial chordate, Botryllus schlosseri. ELife. 2013, 2: e00569-PubMedView ArticleGoogle Scholar
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11 (10): R106-10.1186/gb-2010-11-10-r106.PubMed CentralPubMedView ArticleGoogle Scholar
- Brown FD, Tiozzo S, Roux MM, Ishizuka K, Swalla BJ, De Tomaso AW: Early lineage specification of long-lived germline precursors in the colonial ascidian Botryllus schlosseri. Development. 2009, 136 (20): 3485-3494. 10.1242/dev.037754.PubMed CentralPubMedView ArticleGoogle Scholar
- Kawamura K, Sunanaga T: Role of Vasa, Piwi, and Myc-expressing coelomic cells in gonad regeneration of the colonial tunicate, Botryllus primigenus. Mech Dev. 2011, 128 (7–10): 457-470.PubMedView ArticleGoogle Scholar
- Rinkevich Y, Voskoboynik A, Rosner A, Rabinowitz C, Paz G, Oren M, Douek J, Alfassi G, Moiseeva E, Ishizuka KJ, Palmeri KJ, Weissman IL, Rinkevich B: Repeated, long-term cycling of putative stem cells between niches in a basal chordate. Dev Cell. 2013, 24 (1): 76-88. 10.1016/j.devcel.2012.11.010.PubMedView ArticleGoogle Scholar
- Rosner A, Moiseeva E, Rabinowitz C, Rinkevich B: Germ lineage properties in the urochordate Botryllus schlosseri - from markers to temporal niches. Dev Biol. 2013, 384 (2): 356-374. 10.1016/j.ydbio.2013.10.002.PubMedView ArticleGoogle Scholar
- Sunanaga T, Inubushi H, Kawamura K: Piwi-expressing hemoblasts serve as germline stem cells during postembryonic germ cell specification in colonial ascidian, Botryllus primigenus. Dev Growth Differ. 2010, 52 (7): 603-614. 10.1111/j.1440-169X.2010.01196.x.PubMedView ArticleGoogle Scholar
- Bucko-Justyna M, Lipinski L, Burgering BM, Trzeciak L: Characterization of testis-specific serine-threonine kinase 3 and its activation by phosphoinositide-dependent kinase-1-dependent signalling. FEBS J. 2005, 272 (24): 6310-6323. 10.1111/j.1742-4658.2005.05018.x.PubMedView ArticleGoogle Scholar
- Toshima J, Koji T, Mizuno K: Stage-specific expression of testis-specific protein kinase 1 (TESK1) in rat spermatogenic cells. Biochem Biophys Res Commun. 1998, 249 (1): 107-112. 10.1006/bbrc.1998.9099.PubMedView ArticleGoogle Scholar
- Hofmann O, Caballero OL, Stevenson BJ, Chen YT, Cohen T, Chua R, Maher CA, Panji S, Schaefer U, Kruger A, Lehvaslaiho M, Carninci P, Hayashizaki Y, Jongeneel CV, Simpson AJ, Old LJ, Hide W: Genome-wide analysis of cancer/testis gene expression. Proc Natl Acad Sci U S A. 2008, 105 (51): 20422-20427. 10.1073/pnas.0810777105.PubMed CentralPubMedView ArticleGoogle Scholar
- Sutovsky P: Sperm proteasome and fertilization. Reproduction. 2011, 142 (1): 1-14. 10.1530/REP-11-0041.PubMedView ArticleGoogle Scholar
- Tres LL, Kierszenbaum AL: The ADAM-integrin-tetraspanin complex in fetal and postnatal testicular cords. Birth Defects Res C Embryo Today. 2005, 75 (2): 130-141. 10.1002/bdrc.20041.PubMedView ArticleGoogle Scholar
- Schneider WJ: Vitellogenin receptors: oocyte-specific members of the low-density lipoprotein receptor supergene family. Int Rev Cytol. 1996, 166: 103-137.PubMedView ArticleGoogle Scholar
- Akasaka M, Kato KH, Kitajima K, Sawada H: Identification of novel isoforms of vitellogenin expressed in ascidian eggs. J Exp Zool B Mol Dev Evol. 2013, 320 (2): 118-128. 10.1002/jez.b.22488.PubMedView ArticleGoogle Scholar
- Sawada H, Sakai N, Abe Y, Tanaka E, Takahashi Y, Fujino J, Kodama E, Takizawa S, Yokosawa H: Extracellular ubiquitination and proteasome-mediated degradation of the ascidian sperm receptor. Proc Natl Acad Sci U S A. 2002, 99 (3): 1223-1228. 10.1073/pnas.032389499.PubMed CentralPubMedView ArticleGoogle Scholar
- Kurn U, Sommer F, Bosch TC, Khalturin K: In the urochordate Ciona intestinalis zona pellucida domain proteins vary among individuals. Dev Comp Immunol. 2007, 31 (12): 1242-1254. 10.1016/j.dci.2007.03.011.PubMedView ArticleGoogle Scholar
- Mengerink KJ, Vacquier VD: Glycobiology of sperm-egg interactions in deuterostomes. Glycobiology. 2001, 11 (4): 37R-43R. 10.1093/glycob/11.4.37R.PubMedView ArticleGoogle Scholar
- Ohnishi J, Ohnishi E, Shibuya H, Takahashi T: Functions for proteinases in the ovulatory process. Biochim Biophys Acta. 2005, 1751 (1): 95-109. 10.1016/j.bbapap.2005.05.002.PubMedView ArticleGoogle Scholar
- Napoli JL: 17beta-Hydroxysteroid dehydrogenase type 9 and other short-chain dehydrogenases/reductases that catalyze retinoid, 17beta- and 3alpha-hydroxysteroid metabolism. Mol Cell Endocrinol. 2001, 171 (1–2): 103-109.PubMedView ArticleGoogle Scholar
- Fusi FM, Montesano M, Bernocchi N, Panzeri C, Ferrara F, Villa A, Bronson RA: P-selectin is expressed on the oolemma of human and hamster oocytes following sperm adhesion and is also detected on the equatorial region of acrosome-reacted human spermatozoa. Mol Hum Reprod. 1996, 2 (5): 341-347. 10.1093/molehr/2.5.341.PubMedView ArticleGoogle Scholar
- Geng JG, Raub TJ, Baker CA, Sawada GA, Ma L, Elhammer AP: Expression of a P-selectin ligand in zona pellucida of porcine oocytes and P-selectin on acrosomal membrane of porcine sperm cells. Potential implications for their involvement in sperm-egg interactions. J Cell Biol. 1997, 137 (3): 743-754. 10.1083/jcb.137.3.743.PubMed CentralPubMedView ArticleGoogle Scholar
- Oehninger S, Patankar M, Seppala M, Clark GF: Involvement of selectin-like carbohydrate binding specificity in human gamete interaction. Andrologia. 1998, 30 (4–5): 269-274.PubMedGoogle Scholar
- Carmona-Saez P, Chagoyen M, Tirado F, Carazo JM, Pascual-Montano A: GENECODIS: a web-based tool for finding significant concurrent annotations in gene lists. Genome Biol. 2007, 8 (1): R3-10.1186/gb-2007-8-1-r3.PubMed CentralPubMedView ArticleGoogle Scholar
- Nogales-Cadenas R, Carmona-Saez P, Vazquez M, Vicente C, Yang X, Tirado F, Carazo JM, Pascual-Montano A: GeneCodis: interpreting gene lists through enrichment analysis and integration of diverse biological information. Nucleic Acids Res. 2009, 37 (Web Server issue): W317-W322.PubMed CentralPubMedView ArticleGoogle Scholar
- Tabas-Madrid D, Nogales-Cadenas R, Pascual-Montano A: GeneCodis3: a non-redundant and modular enrichment analysis tool for functional genomics. Nucleic Acids Res. 2012, 40 (Web Server issue): W478-W483.PubMed CentralPubMedView ArticleGoogle Scholar
- Babin PJ: Apolipoproteins and the association of egg yolk proteins with plasma high density lipoproteins after ovulation and follicular atresia in the rainbow trout (Salmo gairdneri). J Biol Chem. 1987, 262 (9): 4290-4296.PubMedGoogle Scholar
- Fujimoto VY, Kane JP, Ishida BY, Bloom MS, Browne RW: High-density lipoprotein metabolism and the human embryo. Hum Reprod Update. 2010, 16 (1): 20-38. 10.1093/humupd/dmp029.PubMedView ArticleGoogle Scholar
- Walzem RL, Hansen RJ, Williams DL, Hamilton RL: Estrogen induction of VLDLy assembly in egg-laying hens. J Nutr. 1999, 129 (2S Suppl): 467S-472S.PubMedGoogle Scholar
- Boyd HC, Brown SK, Harpi JA, Weissmanz IL: Growth and Sexual Maturation of Laboratory Cultured Monterey Botryllus Schlosseri. Biol Bull. 1986, 170: 91-109. 10.2307/1541383.View ArticleGoogle Scholar
- Lauzon RJ, Rinkevich B, Patton CW, Weissman IL: A morphological study of nonrandom senescence in a colonial urochordate. Biol Bull. 2000, 198 (3): 367-378. 10.2307/1542692.PubMedView ArticleGoogle Scholar
- McKitrick TR, Muscat CC, Pierce JD, Bhattacharya D, De Tomaso AW: Allorecognition in a basal chordate consists of independent activating and inhibitory pathways. Immunity. 2011, 34 (4): 616-626. 10.1016/j.immuni.2011.01.019.PubMedView ArticleGoogle Scholar
- Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M, Yefanov A, Lee H, Zhang N, Robertson CL, Serova N, Davis S, Soboleva A: NCBI GEO: archive for functional genomics data sets–update. Nucleic Acids Res. 2013, 41 (Database issue): D991-D995.PubMed CentralPubMedView ArticleGoogle Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30 (1): 207-210. 10.1093/nar/30.1.207.PubMed CentralPubMedView ArticleGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10 (3): R25-10.1186/gb-2009-10-3-r25.PubMed CentralPubMedView ArticleGoogle Scholar
- Lapatsina L, Brand J, Poole K, Daumke O, Lewin GR: Stomatin-domain proteins. Eur J Cell Biol. 2012, 91 (4): 240-245. 10.1016/j.ejcb.2011.01.018.PubMedView ArticleGoogle Scholar
- Yang MH, Su JZ, Liu J, Li YH: Cloning and expression of Tssk1 & Tssk2 in mice and the presence & localization of them in mature sperm. Dongwuxue Yanjiu. 2012, 33 (4): 381-388.PubMedGoogle Scholar
- Kasamatsu S, Hachiya A, Fujimura T, Sriwiriyanont P, Haketa K, Visscher MO, Kitzmiller WJ, Bello A, Kitahara T, Kobinger GP, et al: Essential role of microfibrillar-associated protein 4 in human cutaneous homeostasis and in its photoprotection. Sci Rep. 2011, 1: 164-PubMed CentralPubMedView ArticleGoogle Scholar
- Xue F, Cooley L: kelch encodes a component of intercellular bridges in Drosophila egg chambers. Cell. 1993, 72 (5): 681-693. 10.1016/0092-8674(93)90397-9.PubMedView ArticleGoogle Scholar
- Carrell DT: Epigenetics of the male gamete. Fertil Steril. 2012, 97 (2): 267-274. 10.1016/j.fertnstert.2011.12.036.PubMedView ArticleGoogle Scholar
- Richard I, Devaud C, Cherif D, Cohen D, Beckmann JS: The gene for creatine kinase, mitochondrial 2 (sarcomeric; CKMT2), maps to chromosome 5q13.3. Genomics. 1993, 18 (1): 134-136. 10.1006/geno.1993.1437.PubMedView ArticleGoogle Scholar
- Toshima J, Nakagawara K, Mori M, Noda T, Mizuno K: Structural organization and chromosomal localization of the mouse tesk1 (testis-specific protein kinase 1) gene. Gene. 1998, 206 (2): 237-245. 10.1016/S0378-1119(97)00591-X.PubMedView ArticleGoogle Scholar
- Yamaguchi S, Yamane T, Takahashi-Niki K, Kato I, Niki T, Goldberg MS, Shen J, Ishimoto K, Doi T, Iguchi-Ariga SM, et al: Transcriptional activation of low-density lipoprotein receptor gene by DJ-1 and effect of DJ-1 on cholesterol homeostasis. PLoS One. 2012, 7 (5): e38144-10.1371/journal.pone.0038144.PubMed CentralPubMedView ArticleGoogle Scholar
- Herlyn H, Zischler H: The molecular evolution of sperm zonadhesin. Int J Dev Biol. 2008, 52 (5–6): 781-790.PubMedView ArticleGoogle Scholar
- Wang S, Zheng H, Esaki Y, Kelly F, Yan W: Cullin3 is a KLHL10-interacting protein preferentially expressed during late spermiogenesis. Biol Reprod. 2006, 74 (1): 102-108.PubMedView ArticleGoogle Scholar
- Hartong DT, Dange M, McGee TL, Berson EL, Dryja TP, Colman RF: Insights from retinitis pigmentosa into the roles of isocitrate dehydrogenases in the Krebs cycle. Nat Genet. 2008, 40 (10): 1230-1234. 10.1038/ng.223.PubMed CentralPubMedView ArticleGoogle Scholar
- Govindaraju A, Dogan S, Rodriguez-Osorio N, Grant K, Kaya A, Memili E: Delivering value from sperm proteomics for fertility. Cell Tissue Res. 2012, 349 (3): 783-793. 10.1007/s00441-012-1452-2.PubMedView ArticleGoogle Scholar
- Tanigawa M, Miyamoto K, Kobayashi S, Sato M, Akutsu H, Okabe M, Mekada E, Sakakibara K, Miyado M, Umezawa A, et al: Possible involvement of CD81 in acrosome reaction of sperm in mice. Mol Reprod Dev. 2008, 75 (1): 150-155. 10.1002/mrd.20709.PubMedView ArticleGoogle Scholar
- Edwards DR, Handsley MM, Pennington CJ: The ADAM metalloproteinases. Mol Aspects Med. 2008, 29 (5): 258-289. 10.1016/j.mam.2008.08.001.PubMedView ArticleGoogle Scholar
- Zhong JL, Poghosyan Z, Pennington CJ, Scott X, Handsley MM, Warn A, Gavrilovic J, Honert K, Kruger A, Span PN, et al: Distinct functions of natural ADAM-15 cytoplasmic domain variants in human mammary carcinoma. Mol Cancer Res. 2008, 6 (3): 383-394. 10.1158/1541-7786.MCR-07-2028.PubMedView ArticleGoogle Scholar
- Costache V, McDougall A, Dumollard R: Cell cycle arrest and activation of development in marine invertebrate deuterostomes. Biochem Biophys Res Commun. 2014Google Scholar
- Lum L, Blobel CP: Evidence for distinct serine protease activities with a potential role in processing the sperm protein fertilin. Dev Biol. 1997, 191 (1): 131-145. 10.1006/dbio.1997.8609.PubMedView ArticleGoogle Scholar
- Brown PR, Miki K, Harper DB, Eddy EM: A-kinase anchoring protein 4 binding proteins in the fibrous sheath of the sperm flagellum. Biol Reprod. 2003, 68 (6): 2241-2248.PubMedView ArticleGoogle Scholar
- Vogel BE, Hedgecock EM: Hemicentin, a conserved extracellular member of the immunoglobulin superfamily, organizes epithelial and other cell attachments into oriented line-shaped junctions. Development. 2001, 128 (6): 883-894.PubMedGoogle Scholar
- Spiridonov NA, Wong L, Zerfas PM, Starost MF, Pack SD, Paweletz CP, Johnson GR: Identification and characterization of SSTK, a serine/threonine protein kinase essential for male fertility. Mol Cell Biol. 2005, 25 (10): 4250-4261. 10.1128/MCB.25.10.4250-4261.2005.PubMed CentralPubMedView ArticleGoogle Scholar
- Kurahashi H, Taniguchi M, Meno C, Taniguchi Y, Takeda S, Horie M, Otani H, Toda T: Basement membrane fragility underlies embryonic lethality in fukutin-null mice. Neurobiol Dis. 2005, 19 (1–2): 208-217.PubMedView ArticleGoogle Scholar
- Ricardo S, Lehmann R: An ABC transporter controls export of a Drosophila germ cell attractant. Science. 2009, 323 (5916): 943-946. 10.1126/science.1166239.PubMed CentralPubMedView ArticleGoogle Scholar
- Konno A, Setou M, Ikegami K: Ciliary and flagellar structure and function–their regulations by posttranslational modifications of axonemal tubulin. Int Rev Cell Mol Biol. 2012, 294: 133-170.PubMedView ArticleGoogle Scholar
- Milisav I, Jones MH, Affara NA: Characterization of a novel human dynein-related gene that is specifically expressed in testis. Mamm Genome. 1996, 7 (9): 667-672. 10.1007/s003359900202.PubMedView ArticleGoogle Scholar
- Wiesner B, Weiner J, Middendorff R, Hagen V, Kaupp UB, Weyand I: Cyclic nucleotide-gated channels on the flagellum control Ca2+ entry into sperm. J Cell Biol. 1998, 142 (2): 473-484. 10.1083/jcb.142.2.473.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.