Rapid Evolution of a Sexual Reproduction Gene in Centric Diatoms of the Genus Thalassiosira

Sexual reproduction is commonly assumed to occur in the vast majority of diatoms due to the intimate association of this process with cell size control. Surprisingly, however, little is known about the impact of sexual events on diatom population dynamics. The Sig1 gene is strongly upregulated during sexual reproduction in the centric diatom Thalassiosira weissflogii and has been hypothesized to encode a protein involved in gamete recognition. In the present study, degenerate PCR primers were designed and used to amplify a portion of Sig1 from three closely related species in the cosmopolitan genus Thalassiosira, Thalassiosira oceanica, Thalassiosira guillardii, and Thalassiosira pseudonana. Identification of Sig1 in these three additional species facilitated development of this gene as a molecular marker for diatom sexual events. Examination of the new sequences indicated that multiple copies of Sig1 are probably present in the genome. Moreover, compared to the housekeeping gene β-tubulin, the Sig1 genes of isolates of T. weissflogii collected from different regions of the Atlantic and Pacific oceans displayed high levels of divergence. The Sig1 genes of the four closely related Thalassiosira species also displayed high levels of sequence divergence compared to the levels observed with a second gene, Fcp, probably explaining why Sig1 could not be amplified from more distantly related species. The high levels of sequence divergence both within and between species suggest that Sig1 is rapidly evolving in a manner reminiscent of the manner observed in other genes that encode gamete recognition proteins. A simple model is presented for Sig1 evolution and the implications of such a rapidly evolving sexual reproduction gene for diatom speciation and population dynamics.     

Identification of a new gene family expressed during the onset of sexual reproduction in the centric diatom Thalassiosira weissflogii.

An intriguing feature of the diatom life cycle is that sexual reproduction and the generation of genetic diversity are coupled to the control of cell size. A PCR-based cDNA subtraction technique was used to identify genes that are expressed as small cells of the centric diatom Thalassiosira weissflogii initiate gametogenesis. Ten genes that are up-regulated during the early stages of sexual reproduction have been identified thus far. Three of the sexually induced genes, Sig1, Sig2, and Sig3, were sequenced to completion and are members of a novel gene family. The three polypeptides encoded by these genes possess different molecular masses and charges but display many features in common: they share five highly conserved domains; they each contain three or more cysteine-rich epithelial growth factor (EGF)-like repeats; and they each display homology to the EGF-like region of the vertebrate extracellular matrix glycoprotein tenascin X. Interestingly, the five conserved domains appear in the same order in each polypeptide but are separated by variable numbers of nonconserved amino acids. SIG1 and SIG2 display putative regulatory domains within the nonconserved regions. A calcium-binding, EF-hand motif is found in SIG1, and an ATP/GTP binding motif is present in SIG2. The striking similarity between the SIG polypeptides and extracellular matrix components commonly involved in cell-cell interactions suggests that the SIG polypeptides may play a role in sperm-egg recognition. The SIG polypeptides are thus important molecular targets for determining when and where sexual reproduction occurs in the field.     

The effect of positive selection on a sexual reproduction gene in Thalassiosira weissflogii (Bacillariophyta): results obtained from maximum-likelihood and parsimony-based methods

Maximum-Likelihood-based and parsimony-based methods were used to test for potential effects of positive selection on the sexually induced gene 1 (Sig1) in Thalassiosira weissflogii. The Sig proteins are thought to play a role in mediating sperm-egg recognition during the sexual reproduction phase. The results obtained from parsimony-based analyses showed that none of the amino acid sites were influenced by positive selection. Maximum-likelihood analyses indicated that positive selection was affecting a maximum of seven and a minimum of four amino acid sites in the polypeptide derived from Sig1. It was concluded that the results obtained from the maximum-likelihood-based method are more reliable than those obtained from the parsimony-based approach. This is apparently the first study that has shown that reproductive proteins in unicellular eukaryotes are influenced by positive selection.     

Identification of a New Gene Family Expressed during the Onset of Sexual Reproduction in the Centric Diatom Thalassiosira weissflogii

An intriguing feature of the diatom life cycle is that sexual reproduction and the generation of genetic diversity are coupled to the control of cell size. A PCR-based cDNA subtraction technique was used to identify genes that are expressed as small cells of the centric diatom Thalassiosira weissflogii initiate gametogenesis. Ten genes that are up-regulated during the early stages of sexual reproduction have been identified thus far. Three of the sexually induced genes, Sig1, Sig2, and Sig3, were sequenced to completion and are members of a novel gene family. The three polypeptides encoded by these genes possess different molecular masses and charges but display many features in common: they share five highly conserved domains; they each contain three or more cysteine-rich epithelial growth factor (EGF)-like repeats; and they each display homology to the EGF-like region of the vertebrate extracellular matrix glycoprotein tenascin X. Interestingly, the five conserved domains appear in the same order in each polypeptide but are separated by variable numbers of nonconserved amino acids. SIG1 and SIG2 display putative regulatory domains within the nonconserved regions. A calcium-binding, EF-hand motif is found in SIG1, and an ATP/GTP binding motif is present in SIG2. The striking similarity between the SIG polypeptides and extracellular matrix components commonly involved in cell-cell interactions suggests that the SIG polypeptides may play a role in sperm-egg recognition. The SIG polypeptides are thus important molecular targets for determining when and where sexual reproduction occurs in the field.     

False-positive selection identified by ML-based methods: examples from the Sig1 gene of the diatom Thalassiosira weissflogii and the tax gene of a human T-cell lymphotropic virus

Sexually induced gene 1 (Sig1) in the centric diatom Thalassiosira weissflogii is considered to encode a gamete recognition protein. Sorhannus (2003) analyzed nucleotide sequences of Sig1 using parsimony analysis and the maximum-likelihood (ML)-based Bayesian method for inferring positive selection at single amino acid sites and reported that positively selected sites were detected by the latter method but not by the former. He then concluded that for this type of study, the ML-based method is more reliable than parsimony analysis. Here we show that his results apparently represent false-positive cases of the ML-based method and that there is no solid evidence that this gene contains positively selected sites. We further demonstrate that in the tax gene of human T-cell lymphotropic virus type I (HTLV-I), all codon sites, including invariable sites, can be inferred as positively selected sites by the ML-based method. These observations indicate that the ML-based method may produce many false-positive sites. One of the main reasons for the occurrence of false positives is that in the ML-based method, codon sites are grouped into several categories, with different nonsynonymous/synonymous rate ratios (omegas), on a purely statistical basis, and positive selection is inferred indirectly by examining whether the average omega for each category is greater than 1. In parsimony analysis, however, the evolutionary change of nucleotides at each codon site is examined. For this reason, parsimony-based methods rarely produce false positives and are safer than ML-based methods for detecting positive selection at individual codon sites, although a large number of sequences are necessary.     

FLOW CYTOMETRIC ANALYSIS OF SPERMATOGENESIS IN THE DIATOM THALASSIOSIRA WEISSFLOGII (BACILLARIOPHYCEAE)1

Flow cytometry was used to detect and quantify sexual differentiation in the centric diatom Thalassiosira weissflogii (Grun.). Size (light scatter), chlorophyll, protein and DNA contents were measured for each cell throughout the process of differentiation. Male gametes were small round cells characterized by one complement of DNA and a lower protein and chlorophyll content than vegetative cells. Male gamete formation was induced by a long period of darkness (2 days) followed by a transfer to continuous light. Up to 30% of the initial cell population produced male gametes which appeared in the culture 14 h after release from darkness. Male gamete production was also detected in exponentially growing cultures in continuous light, but to a much smaller degree.     

Diversity of the cadmium-containing carbonic anhydrase in marine diatoms and natural waters

A recent report of a novel carbonic anhydrase (CDCA1) with Cd as its metal centre in the coastal diatom Thalassiosira weissflogii has led us to search for the occurrence of this Cd enzyme (CDCA) in other marine phytoplankton and in the environment. Using degenerate primers designed from the published sequences from T. weissflogii and a putative sequence in the genome of Thalassiosira pseudonana, we show that CDCA is widespread in diatom species and ubiquitous in the environment. All detected genes share more than 64% amino acid identity with the CDCA of T. pseudonana. Analysis of the amino acid sequence of CDCA shows that the putative Cd binding site resembles that of beta-class carbonic anhydrases (CAs). The prevalence of CAs in diatoms that presumably contain Cd at their active site probably reflects the very low concentration of Zn in the marine environment and the difficulty in acquiring inorganic carbon for photosynthesis. The cdca primers developed in this study should be useful for detecting cdca genes in the field, and studying the conditions under which they are expressed.     

The evolution of sperm competition genes: The effect of mating system on levels of genetic variation within and between species

It is widely established that proteins involved in reproduction diverge between species more quickly than other proteins. For male sperm proteins, rapid divergence is believed to be caused by postcopulatory sexual selection and/or sexual conflict. Here, we derive the expected levels of gene diversity within populations and divergence between them for male sperm protein genes evolving by postcopulatory, prezygotic fertility competition, i.e. the function imputed for some sperm and seminal fluid genes. We find that, at the mutation‐selection equilibrium, both gene diversity within species and divergence between them are elevated relative to genes with similar selection coefficients expressed by both sexes. We show that their expected level of diversity is a function of the harmonic mean number of mates per female, which affects the strength of fertility selection stemming from male–male sperm competition. Our predictions provide a null hypothesis for distinguishing between other selective hypotheses accounting for the rapid evolution of male reproductive genes.     

Chemotaxis of Marinobacter adhaerens and Its Impact on Attachment to the Diatom Thalassiosira weissflogii

Alga-bacterium interactions are crucial for aggregate formation and carbon cycling in aquatic systems. To understand the initiation of these interactions, we investigated bacterial chemotaxis within a bilateral model system. Marinobacter adhaerens HP15 has been demonstrated to attach to the diatom Thalassiosira weissflogii and induce transparent exopolymeric particle and aggregate formation. M. adhaerens possesses one polar flagellum and is highly motile. Bacterial cells were attracted to diatom cells, as demonstrated by addition of diatom cell homogenate or diatom culture supernatant to soft agar, suggesting that chemotaxis might be important for the interaction of M. adhaerens with diatoms. Three distinct chemotaxis-associated gene clusters were identified in the genome sequence of M. adhaerens, with the clusters showing significant sequence similarities to those of Pseudomonas aeruginosa PAO1. Mutations in the genes cheA, cheB, chpA, and chpB, which encode histidine kinases and methylesterases and which are putatively involved in either flagellum-associated chemotaxis or pilus-mediated twitching motility, were generated and mutants with the mutations were phenotypically analyzed. ΔcheA and ΔcheB mutants were found to be swimming deficient, and all four mutants were impaired in biofilm formation on abiotic surfaces. Comparison of the HP15 wild type and its chemotaxis mutants in cocultures with the diatom revealed that the fraction of bacteria attaching to the diatom decreased significantly for mutants in comparison to that for the wild type. Our results highlight the importance of M. adhaerens chemotaxis in initiation of its interaction with the diatom. In-depth knowledge of these basic processes in interspecies interactions is pivotal to obtain a systematic understanding of organic matter flux and nutrient cycling in marine ecosystems.     

Sexual selection at the protein level drives the extraordinary divergence of sex-related genes during sympatric speciation

An increasing number of molecular studies are indicating that, in a wide variety of species, genes directly related to fertilization evolve at extraordinarily high rates. We try to gain insight into the dynamics of this rapid evolution and its underlying mechanisms by means of a simple theoretical model. In the model, sexual selection and sympatric speciation act together in order to drive rapid divergence of gamete recognition proteins. In this process, intraspecific competition for fertilizations enlarges male gamete protein variation by means of evolutionary branching, which initiates sympatric speciation. In addition, avoidance of competition for fertilizations between the incipient species drives the rapid evolution of gamete recognition proteins. This mechanism can account for both strong stabilizing selection on gamete recognition proteins within species and rapid divergence between species. Moreover, it can explain the empirical finding that the rate of divergence of fertilization genes is not constant, but highest between closely related species.     

GAMETE INTERACTIONS AND GENETIC DIFFERENTIATION AMONG THREE SYMPATRIC POLYCHAETES

The evolution of gamete incompatibility between free-spawning marine invertebrate species has been explained by three hypotheses: (1) independent divergence at gamete recognition loci; (2) selection against hybrids; and (3) a process of sexual selection involving polymorphic gamete recognition loci (Metz and Palumbi 1996). The first two hypotheses predict that gamete incompatibility appears only after gene flow has been halted for other reasons and the third that gamete incompatibility appears simultanously with blocks to gene flow. Here I show that gametes of three sympatric polychaetes in the genus Arctonoe are compatible in all crosses, over a broad range of gamete concentrations and contact times. Although at least some hybrid crosses produce fertile adults, allozyme and mitochondrial DNA sequence data indicate that the three species do not regularly exchange genes. These data are consistent with predictions of the first two hypotheses for the evolution of gamete incompatibility, but allow rejection of the third hypothesis. Gametes of the three species are compatible despite estimated divergence times of 1–3 M.Y.B.P.; in several other marine invertebrates, divergence times of the same magnitude are associated with asymmetric or complete gamete incompatibility. It appears likely that segregation of symbiotic adults on their respective host species restricts mating opportunities, and thus gene flow, among Arctonoe species.     

Complete genome sequence of Marinobacter adhaerens type strain (HP15), a diatom-interacting marine microorganism

Marinobacter adhaerens HP15 is the type strain of a newly identified marine species, which is phylogenetically related to M. flavimaris, M. algicola, and M. aquaeolei. It is of special interest for research on marine aggregate formation because it showed specific attachment to diatom cells. In vitro it led to exopolymer formation and aggregation of these algal cells to form marine snow particles. M. adhaerens HP15 is a free-living, motile, rod-shaped, Gram-negative gammaproteobacterium, which was originally isolated from marine particles sampled in the German Wadden Sea. M. adhaerens HP15 grows heterotrophically on various media, is easy to access genetically, and serves as a model organism to investigate the cellular and molecular interactions with the diatom Thalassiosira weissflogii. Here we describe the complete and annotated genome sequence of M. adhaerens HP15 as well as some details on flagella-associated genes. M. adhaerens HP15 possesses three replicons; the chromosome comprises 4,422,725 bp and codes for 4,180 protein-coding genes, 51 tRNAs and three rRNA operons, while the two circular plasmids are ~187 kb and ~42 kb in size and contain 178 and 52 protein-coding genes, respectively.     

Light and dark control of the cell cycle in two marine phytoplankton species

The effect of light and dark on growth, DNA replication and cell division of two marine phytoplankters Thalassiosira weissflogii (a diatom) and Hymenomonas carterae (a coccolithophorid) was investigated using flow cytometry. The two species displayed very differing behavior. When transferred from light to prolonged darkness, all coccolithophorid cells were arrested at the beginning of the G1 stage of the cell cycle. When shifted back into light, they resumed cycling at a rate slightly slower than prior to arrest. In contrast, diatom cells were arrested either in the G1 or G2 stage of the cell cycle in the dark. Upon re-exposure to light, cells which had been dark-arrested in G1 resumed cycling at the same rate as prior to arrest, while cells arrested in G2 cycled much more slowly. These results suggest that in both species, light control of cell cycle progression is effective only over a restricted part of the cell cycle, as has been hypothesized by Spudich & Sager (J cell biol 83 (1980) 136) [38] for Chlamydomonas. In the coccolithophorid there is a single light-dependent segment located at the beginning of G1, whereas the diatom appears to have two such segments, one in G1 and the other in G2, corresponding to two different light requiring processes.     

The active site structure of Thalassiosira weissflogii carbonic anhydrase 1

X-ray absorption spectroscopy at the Zn K-edge indicates that the active site of the marine diatom Thalassiosira weissflogii carbonic anhydrase is strikingly similar to that of mammalian alpha-carbonic anhydrase enzymes. The zinc has three histidine ligands and a single water at 1.98 A. This is quite different from the beta-carbonic anhydrases of higher plants in which zinc is coordinated by two cysteine thiolates, one histidine, and a water molecule. The diatom carbonic anhydrase shows no significant sequence similarity with other carbonic anhydrases and may represent an example of convergent evolution at the molecular level.     

Regulation of carbonic anhydrase expression by zinc, cobalt, and carbon dioxide in the marine diatom Thalassiosira weissflogii

TWCA1 is the major Zn-requiring isoform of carbonic anhydrase (CA) in the marine diatom Thalassiosira weissflogii. We have examined the roles that trace metals and CO(2) play in the regulation of TWCA1 expression over ranges of concentrations that bracket those encountered in the marine environment. Both steady-state levels of TWCA1 and the kinetics of induction were measured by western analysis. TWCA1 levels correlated well with cellular CA activity levels. TWCA1 was induced at a low CO(2) concentration but the level of induction, as determined by western analysis, was dependent on the availability of Zn. Co effectively substituted for Zn in regulating TWCA1 expression and promoting TWCA1 activity. Upon shift from low to high CO(2), the concentration of TWCA1 decreased. The expression of TWCA1 is diel cycle regulated, and cellular TWCA1 decreased during the dark phase. These results provide the basis for studying the expression of CA in field populations and, taken together with previous radiolabeling studies, provide strong evidence of in vivo metal substitution of Co for Zn in a CA. Our data also support the conclusion that TWCA1 plays a central role in carbon acquisition in T. weissflogii.     

Extraordinary sequence divergence at Tsga8, an X-linked gene involved in mouse spermiogenesis

The X chromosome plays an important role in both adaptive evolution and speciation. We used a molecular evolutionary screen of X-linked genes potentially involved in reproductive isolation in mice to identify putative targets of recurrent positive selection. We then sequenced five very rapidly evolving genes within and between several closely related species of mice in the genus Mus. All five genes were involved in male reproduction and four of the genes showed evidence of recurrent positive selection. The most remarkable evolutionary patterns were found at Testis-specific gene a8 (Tsga8), a spermatogenesis-specific gene expressed during postmeiotic chromatin condensation and nuclear transformation. Tsga8 was characterized by extremely high levels of insertion-deletion variation of an alanine-rich repetitive motif in natural populations of Mus domesticus and M. musculus, differing in length from the reference mouse genome by up to 89 amino acids (27% of the total protein length). This population-level variation was coupled with striking divergence in protein sequence and length between closely related mouse species. Although no clear orthologs had previously been described for Tsga8 in other mammalian species, we have identified a highly divergent hypothetical gene on the rat X chromosome that shares clear orthology with the 5' and 3' ends of Tsga8. Further inspection of this ortholog verified that it is expressed in rat testis and shares remarkable similarity with mouse Tsga8 across several general features of the protein sequence despite no conservation of nucleotide sequence across over 60% of the rat-coding domain. Overall, Tsga8 appears to be one of the most rapidly evolving genes to have been described in rodents. We discuss the potential evolutionary causes and functional implications of this extraordinary divergence and the possible contribution of Tsga8 and the other four genes we examined to reproductive isolation in mice.