One- and Two-Locus Population Models With Differential Viability Between Sexes: Parallels Between Haploid Parental Selection and Genomic Imprinting

A model of genomic imprinting with complete inactivation of the imprinted allele is shown to be formally equivalent to the haploid model of parental selection. When single-locus dynamics are considered, an internal equilibrium is possible only if selection acts in the opposite directions in males and females. I study a two-locus version of the latter model, in which maternal and paternal effects are attributed to the single alleles at two different loci. A necessary condition for the allele frequency equilibria to remain on the linkage equilibrium surface is the multiplicative interaction between maternal and paternal fitness parameters. In this case the equilibrium dynamics are independent at both loci and results from the single-locus model apply. When fitness parameters are additive, analytic treatment was not possible but numerical simulations revealed that stable polymorphism characterized by association between loci is possible only in several special cases in which maternal and paternal fitness contributions are precisely balanced. As in the single-locus case, antagonistic selection in males and females is a necessary condition for the maintenance of polymorphism. I also show that the above two-locus results of the parental selection model are very sensitive to the inclusion of weak directional selection on the individual''s own genotypes.PARENTAL genetic effects refer to the influence of the mother''s and father''s genotypes on the phenotypes of their offspring, not attributable just to the transfer of genes. Examples have been documented across a wide range of areas of organism biology; see, for example, Wade (1998) and ​and22 in Rasanen and Kruuk (2007). Parental selection is a more formal concept used in theoretical modeling and concerns situations where the fitness of the offspring depends, besides other factors, on the genotypes of its parent(s) (generalizing from Kirkpatrick and Lande 1989).

TABLE 1

Frequencies of genotypes and fitness parameterizations in model 1

Evolutionary Strata in a Small Mating-Type-Specific Region of the Smut Fungus Microbotryum violaceum

DNA sequence analysis and genetic mapping of loci from mating-type-specific chromosomes of the smut fungus Microbotryum violaceum demonstrated that the nonrecombining mating-type-specific region in this species comprises ∼25% (∼1 Mb) of the chromosome length. Divergence between homologous mating-type-linked genes in this region varies between 0 and 8.6%, resembling the evolutionary strata of vertebrate and plant sex chromosomes.EVOLUTION of mating types or sex-determining systems often involves the suppression of recombination around the primary sex-determining or mating-type-determining locus. In animals and plants, it is often an entire or almost entire chromosome (Y or W in male or female heterogametic species, respectively) that ceases to recombine with its homologous (X or Z) chromosome (Charlesworth and Charlesworth 2000; Charlesworth 2008). Self-incompatibility loci in plants are also thought to be located in regions of suppressed recombination (Charlesworth et al. 2005; Kamau and Charlesworth 2005; Kamau et al. 2007; Li et al. 2007; Yang et al. 2007). Regardless of the phylogenetic position of a species, such nonrecombining regions are known to follow similar evolutionary trajectories. The nonrecombining region on the sex-specific chromosome expands in several steps, forming evolutionary strata—regions of different X/Y (or Z/W) divergence (Lahn and Page 1999; Handley et al. 2004; Sandstedt and Tucker 2004; Nicolas et al. 2005)—and genes in the nonrecombining regions gradually accumulate deleterious mutations that eventually render them dysfunctional (Charlesworth and Charlesworth 2005; Charlesworth 2008).Fungal mating-type systems are very diverse, with the number of mating types varying from two to several hundred (Casselton 2002). Like sex chromosomes in several animals and plants, suppressed recombination has evolved in regions near fungal mating-type loci, including in Ustilago hordei (Lee et al. 1999), Cryptococcus neoformans (Lengeler et al. 2002), and Neurospora tetrasperma (Menkis et al. 2008). These species have two mating types, but no morphologically distinct sexes. The mating-type locus (the region of suppressed recombination) of C. neoformans is small (∼100 kb) compared with known sex chromosomes and contains only ∼20 genes that, unlike many sex chromosomes (Y or W chromosomes), show no obvious signs of genetic degeneration (Lengeler et al. 2002; Fraser et al. 2004). Judging from the divergence between the homologous genes on the two mating-type-specific chromosomes, C. neoformans started to evolve sex chromosomes a long time ago because silent divergence between the two mating types in the most ancient region exceeds 100% (Fraser et al. 2004). Genes in the younger mating-type-specific region are much less diverged between the two sex chromosomes, suggesting that the evolution of the sex locus in C. neoformans might have proceeded through several steps. The nonrecombining region around the mating-type locus of N. tetrasperma is much larger than in C. neoformans (at least 6.6 Mb), and silent divergence between homologous genes on the mating-type-specific chromosomes ranges from zero to 9%, demonstrating that these mating-type-specific chromosomes evolved recently (Menkis et al. 2008).M. violaceum, which causes anther smut disease in Silene latifolia and other species in the family Caryophyllaceae, has two mating types, A1 and A2 (reviewed by Giraud et al. 2008), which are determined by the presence of mating-type-specific chromosomes (hereafter A1 and A2 chromosomes, or sex chromosomes) in the haploid stage of the life cycle (Hood 2002; Hood et al. 2004). The A1 and A2 chromosomes are distinguishable by size in pulsed-field electrophoresis, and it is possible to isolate individual chromosomes electrophoretically (Hood et al. 2004). Random fragments of A1 and A2 chromosomes have previously been isolated from mating-type-specific bands of pulsed-field separated chromosomes of M. violaceum (Hood et al. 2004). These fragments were assumed to be linked to mating type. The same method was used to isolate fragments of non-mating-type-specific chromosomes. On the basis of the analysis of their sequences, (Hood et al. 2004) proposed that mating-type-specific chromosomes in M. violaceum might be degenerate because they contained a lower proportion of protein-coding genes than other chromosomes. However, it was not determined whether the sequences isolated from the mating-type chromosomes originated from the mating-type-specific or from the recombining regions (Hood et al. 2004), and the relative sizes of these regions are not known for these M. violaceum chromosomes. We tested the mating-type specificity of 86 of these fragments and demonstrate that fewer than a quarter of these loci are located in the mating-type-specific region, suggesting that the nonrecombining region on the A1 and A2 chromosomes is quite small, while the rest of the chromosome probably recombines (like pseudoautosomal regions of sex chromosomes) and is therefore not expected to undergo genetic degeneration. Genetic mapping confirms the presence of two pseudoautosomal regions in the M. violaceum mating-type-specific chromosomes.As these chromosomes are mating type specific in the haploid stage of M. violaceum, mating-type-specific loci (or DNA fragments) can be identified by testing whether they are present exclusively in A1 or A2 haploid strains. We therefore prepared haploid A1 and A2 M. violaceum cultures from S. latifolia plants from two geographically remote locations (accessions Sl405 from Sweden and Sl127 from the French Pyrenees). Haploid sporidial cultures were isolated by a standard dilution method (Kaltz and Shykoff 1997; Oudemans and Alexander 1998). Mating types were determined by PCR amplification of each culture with primers designed for A1 and A2 pheromone receptor genes linked to A1 and A2 mating types (Yockteng et al. 2007). The primers were as follows: 5′-TGGCATCCCTCAATGTTTCC-3′ and 5′-CACCTTTTGATGAGAGGCCG-3′ for the A1 pheromone receptor (GenBank accession no. {"type":"entrez-nucleotide","attrs":{"text":"EF584742","term_id":"156254864","term_text":"EF584742"}}EF584742) and 5′-TGACGAGAGCATTCCTACCG-3′ and 5′-GAAGCGGAACTTGCCTTTCT-3′ for the A2 pheromone receptor (GenBank accession no. {"type":"entrez-nucleotide","attrs":{"text":"EF584741","term_id":"156254862","term_text":"EF584741"}}EF584741). Cultures with PCR product amplified only from an A1 or A2 pheromone receptor gene were selected for further use. The mating types of the cultures were verified by conjugating them in all combinations.The GenBank nucleotide database was searched using BLAST for sequences similar to those isolated by Hood et al. (2004). Sequences with similarity to transposable elements (TE) and other repeats were excluded. The resulting set of nonredundant sequences was used to design PCR primers for 98 fragments. Half of these were originally isolated from the A1 and half from the A2 chromosomes and are hereafter called A1-NNN or A2-NNN (where NNN is the locus number; supporting information, Table S1), which does not imply that these loci are A1 or A2 specific, but merely indicates that they were originally isolated from the A1 or A2 chromosomes. Amplification of these regions from new A1 and A2 M. violaceum cultures, independently isolated by ourselves, revealed that only 5 of the 49 loci isolated from the A1 chromosome are indeed A1 specific and only 6 of 49 isolated from the A2 chromosome are A2 specific. All other loci amplified from both A1 and A2 cultures. Figure 1 illustrates some of these results from the Swedish sample (Sl405).Open in a separate windowFigure 1.—Testing of mating-type specificity for loci isolated from A1 and A2 chromosomes. (a) PCR amplifications from haploid cultures from Sl405 using primers designed from six A1-originated loci. Loci in which a PCR product could be amplified only from A1 cultures (boxed) were classified as specific to mating type A1. (b) PCR tests of six A2-originated loci on the same set of haploids as in a. Loci in which a PCR product amplified only from A2 cultures (boxed) were classified as specific to mating type A2. Loci amplified from both A1 and A2 cultures are not mating type specific.The fragments that amplified from both A1 and A2 mating types may be in recombining regions, or they could be present in mating-type-specific regions on both A1 and A2 chromosomes. If they are in recombining regions, the A1- and A2-linked homologs should not be diverged from each other, but if they are in nonrecombining, mating-type-specific regions, the divergence of the A1- and A2-linked homologs should be roughly proportional to the time since recombination stopped in the region. We therefore sequenced and compared PCR fragments amplified from the two mating types of Sl405 or Sl127 cultures (GenBank accession nos. {"type":"entrez-nucleotide","attrs":{"text":"FI855822","term_id":"258612459","term_text":"FI855822"}}FI855822–{"type":"entrez-nucleotide","attrs":{"text":"FI856001","term_id":"258612570","term_text":"FI856001"}}FI856001). Sequencing of PCR products showed that 12 (4 A1 and 8 A2) loci have more than one copy, and they were excluded from further analysis. Sequences of 61 loci were identical between the A1 and A2 strains, and four loci demonstrated low total divergence (0.24–0.61%) between the two mating types (otintseva and D. Filatov, unpublished results). Thus, these loci might be located in the recombining part of the mating-type-specific chromosomes. Ten of 75 loci that amplified in both mating types demonstrated multiple polymorphisms fixed between the mating types rather than between the locations. Given that the strains that we used in the analysis originated from two geographically distant locations, it is highly unlikely that multiple polymorphisms distinguishing the A1 and A2 sequences arose purely by chance; thus, these loci are probably located in the nonrecombining mating-type-specific region of the M. violaceum A1 and A2 chromosomes.

TABLE 1

Loci from mating-type-specific chromosomes of M. violaceum used for PCR analysis and genetic map construction

		With nonzero A1/A2 divergenceb
Loci	Mating type specific	<1%	>1%	With zero A1/A2 divergenceb	Total
A1a	5	2 (1)	3 (3)	35 (3)	45 (7)
A2a	6	2 (0)	7 (7)	26 (3)	41 (10)
Subtotal		4 (1)	10 (10)
Total	11	14 (11)		61 (6)	86 (17)

Open in a separate window^aA1, loci originated from the A1 sex chromosome; A2, loci originated from the A2 sex chromosome.^bThe number of loci used for genetic map construction is in parentheses.To confirm the mating-type-specific or pseudoautosomal locations of the loci with and without A1/A2 divergence, we conducted genetic mapping in a family of 99 individuals, 50 of which were of mating type A1 and 49 of mating type A2. The family was generated by a cross between A1 and A2 M. violaceum strains from S. latifolia accessions Sl405 (Sweden) and Sl127 (France), respectively. The choice of strains from geographically distant locations was motivated by the hope of maximizing the number of DNA sequence differences between them that can be used as molecular genetic markers in segregation analysis. We inoculated S. latifolia seedlings with sporidial cultures of both mating types. For inoculation, petri dishes with 12-day-old seedlings of S. latifolia were flooded with 2.5 ml of inoculum suspension. Inoculum suspension consisted of equal volumes of the A1 and A2 sporidial cultures that were mixed and conjugated overnight at 14° under rotation (Biere and Honders 1996; Van Putten et al. 2003). Seedlings were potted 3 days after inoculation. Two months later, teliospores were collected from the flowers of the infected plant and grown in petri dishes on 3.6% potato dextrose agar medium. Haploid sporidia formed after meiosis were isolated and grown as separate cultures for DNA extraction. The mating types of single sporidia cultures were identified as described above. The loci analyzed in the segregation analysis were sequenced in the two parental haploid strains and in 99 (50 A1 and 49 A2) haploid strains that were generated in the cross. Single nucleotide differences between the parental strains were used as molecular genetic markers for segregation analysis in the progeny. The genetic map was constructed using MAPMAKER/EXP v3.0 (Lincoln et al. 1992) and MapDisto v1.7 (http://mapdisto.free.fr/).The resulting genetic map is shown in Figure 2. As expected, no recombination was observed between the 10 loci with diverged A1- and A2-linked copies. In addition, one marker with no A1/A2 divergence, A2-397, was also completely linked to the loci with significant A1/A2 divergence. This locus either may be very tightly linked to the nonrecombining mating-type-specific region or may have been added to that region more recently than the loci that had already accumulated some divergence between the alleles in the two mating types. The mating-type-specific pheromone receptor locus (Devier et al. 2009) and 11 mating-type-specific loci are also located in this nonrecombining region (Figure 2). Interestingly, the cluster of nonrecombining markers is flanked on both sides with markers that recombine in meiosis, demonstrating that there are pseudoautosomal regions on both ends of the mating-type-specific chromosomes.Open in a separate windowFigure 2.—Genetic map of the mating-type-determining chromosome in M. violaceum. Genetic distance (in centimorgans) and the relative positions of the markers are shown to the left and the right of the chromosome, respectively. The position of the nonrecombining region corresponds to the cluster of linked markers shown on the right of the figure. Total A1/A2 divergence is shown in parentheses. Eleven mating-type-specific markers (for which sequences are available from only one mating type), located in the nonrecombining mating-type-specific region, are not shown.Our results demonstrate that although the loci reported by Hood et al. (2004) were isolated from the A1 and A2 chromosomes, most of these loci are not located in the nonrecombining mating-type-specific regions. In fact, the nonrecombining region might be relatively small: of 86 tested fragments, only 21 appeared to be either mating type specific or linked to the mating-type locus. Assuming that these loci represent a random set of DNA fragments isolated from the A1 and A2 chromosomes, it is possible to estimate the size of the nonrecombining region using the binomial distribution: the nonrecombining region is expected to be 24.4% (95% CI: 16.7–33.6%) of the chromosome length. As the sizes of the A1 and A2 chromosomes are ∼3.4 and 4.2 Mb long (Hood 2002; Hood et al. 2004), the nonrecombining region might be ∼1 Mb long.Interestingly, total A1/A2 divergence for the 11 loci with A1- and A2-linked copies mapped to the nonrecombining region varied from 0% to 8.6% (Figure 2). In addition, 11 loci amplified from only one mating type. These genes could represent degenerated genes, some of which degenerated in A1 strains, and some in A2 strains. Alternatively, they might be highly diverged genes, such that the PCR primers amplify only one allele, and not the other. Variation in divergence may be the result of the stepwise cessation of recombination between the A1 and A2 chromosomes in M. violaceum, resembling the evolutionary strata reported for human, chicken, and white campion sex chromosomes (Lahn and Page 1999; Handley et al. 2004; Bergero et al. 2007). However, only the differences between the most and the least diverged loci are statistically significant (Devier et al. 2009), the M. violaceum mating-type region has at least three strata: one oldest stratum, including the pheromone receptor locus; a younger stratum with ∼5–9% A1/A2 divergence; and the youngest stratum with 1–4% divergence between the two mating types. There may also be an additional very recently evolved stratum containing the locus named A2-397, which is also present in all A1 strains tested, with no fixed differences between the A1 and A2 strains (

Dehalogenation Activities and Distribution of Reductive Dehalogenase Homologous Genes in Marine Subsurface Sediments

Halogenated organic compounds serve as terminal electron acceptors for anaerobic respiration in a diverse range of microorganisms. Here, we report on the widespread distribution and diversity of reductive dehalogenase homologous (rdhA) genes in marine subsurface sediments. A total of 32 putative rdhA phylotypes were detected in sediments from the southeast Pacific off Peru, the eastern equatorial Pacific, the Juan de Fuca Ridge flank off Oregon, and the northwest Pacific off Japan, collected at a maximum depth of 358 m below the seafloor. In addition, significant dehalogenation activity involving 2,4,6-tribromophenol and trichloroethene was observed in sediment slurry from the Nankai Trough Forearc Basin. These results suggest that dehalorespiration is an important energy-yielding pathway in the subseafloor microbial ecosystem.Scientific ocean drilling explorations have revealed that marine subsurface sediments harbor remarkable numbers of microbial cells that account for approximately 1/10 to 1/3 of all living biota on Earth (20, 25, 33). Thermodynamic calculations of pore-water chemistry suggest that subseafloor microbial activities are generally supported by nutrient and energy supplies from the seawater and/or underlying basaltic aquifers (6, 7). Although sulfate, nitrate, Fe(III), Mn(IV), and bicarbonate are known to be potential electron acceptors for anaerobic microbial respiration in marine subsurface sediments (5), the incidence of both the dissimilatory dehalorespiration pathway and microbial activity in halogenated organic substrates remains largely unknown.Previous molecular ecological studies using 16S rRNA gene sequences demonstrated that Chloroflexi is one of the most frequently detected phyla in subseafloor sediments of the Pacific Ocean margins (12-14). Some of the sequences within the Chloroflexi are closely related to sequences in the genus Dehalococcoides, which contains obligatory dehalorespiring bacteria that employ halogenated organic compounds as terminal electron acceptors (21, 29). The frequent detection of Dehalococcoides-related 16S rRNA genes from these environments implies the occurrence of dissimilatory dehalorespiration in marine subsurface sediments.In this study, we detected and phylogenetically analyzed the reductive dehalogenase homologous (rdhA) genes, key functional genes for dehalorespiration pathways, from frozen sediment core samples obtained by Ocean Drilling Program (ODP) Leg 201 (Peru margin and eastern equatorial Pacific) (7, 14); Integrated Ocean Drilling Program (IODP) Expedition 301 (Juan de Fuca Ridge flank) (8, 24); Chikyu Shakedown Expedition CK06-06 (Northwest Pacific off Japan) (20, 23); and IODP Expedition 315 (Nankai Trough Forearc Basin off Japan) (Table ​(Table1).1). DNA was extracted using an ISOIL bead-beating kit (Nippon Gene, Japan) and purified using a MagExtractor DNA fragment purification kit (Toyobo, Japan) according to the manufacturer''s instructions. To increase concentration, DNA was amplified by multiple displacement amplification using the phi29 polymerase supplied with a GenomiPhi kit (GE Healthcare, United Kingdom) (20). Putative rdhA genes were amplified by PCR using Ex Taq polymerase (TaKaRa, Japan) with degenerate primers RRF2 and B1R (17), dehaloF3, dehaloF4, dehaloF5, dehaloR2, dehaloR3, and dehaloR4 (32), and ceRD2S, ceRD2L, and RD7 (26) and the PCR conditions described in those studies. Amplicons of the approximate target size were gel purified and cloned into the pCR2.1 vector (Invitrogen, Japan). Sequence similarity was analyzed using FastGroupII web-based software (34), and sequences with a 95% identity were tentatively assigned to the same phylotype. Amino acid sequences were aligned by ClustalW (31), including known and putative reductive dehalogenase sequences in the genome of Dehalococcoides ethenogenes strain 195 (28), as well as several functionally characterized reductive dehalogenases from other species.

TABLE 1.

Sample locations and results of PCR amplification of rdhA

Spontaneous Quinolone Resistance in the Zoonotic Serovar of Vibrio vulnificus

This work demonstrates that Vibrio vulnificus biotype 2, serovar E, an eel pathogen able to infect humans, can become resistant to quinolone by specific mutations in gyrA (substitution of isoleucine for serine at position 83) and to some fluoroquinolones by additional mutations in parC (substitution of lysine for serine at position 85). Thus, to avoid the selection of resistant strains that are potentially pathogenic for humans, antibiotics other than quinolones must be used to treat vibriosis on farms.Vibrio vulnificus is an aquatic bacterium from warm and tropical ecosystems that causes vibriosis in humans and fish (http://www.cdc.gov/nczved/dfbmd/disease_listing/vibriov_gi.html) (33). The species is heterogeneous and has been subdivided into three biotypes and more than eight serovars (6, 15, 33; our unpublished results). While biotypes 1 and 3 are innocuous for fish, biotype 2 can infect nonimmune fish, mainly eels, by colonizing the gills, invading the bloodstream, and causing death by septicemia (23). The disease is rapidly transmitted through water and can result in significant economic losses to fish farmers. Surviving eels are immune to the disease and can act as carriers, transmitting vibriosis between farms. Interestingly, biotype 2 isolates belonging to serovar E have been isolated from human infections, suggesting that serovar E is zoonotic (2). This serovar is also the most virulent for fish and has been responsible for the closure of several farms due to massive losses of fish. A vaccine, named Vulnivaccine, has been developed from serovar E isolates and has been successfully tested in the field (14). Although the vaccine provides fish with long-term protection from vibriosis, at present its use is restricted to Spain. For this reason, in many fish farms around the world, vibriosis is treated with antibiotics, which are usually added to the food or water.Quinolones are considered the most effective antibiotics against human and fish vibriosis (19, 21, 31). These antibiotics can persist for a long time in the environment (20), which could favor the emergence of resistant strains under selective pressure. In fact, spontaneous resistances to quinolones by chromosomal mutations have been described for some gram-negative bacteria (10, 11, 17, 24, 25, 26). Therefore, improper antibiotic treatment of eel vibriosis or inadequate residue elimination at farms could favor the emergence of human-pathogenic serovar E strains resistant to quinolones by spontaneous mutations. Thus, the main objective of the present work was to find out if the zoonotic serovar of biotype 2 can become quinolone resistant under selective pressure and determine the molecular basis of this resistance.Very few reports on resistance to antibiotics in V. vulnificus have been published; most of them have been performed with biotype 1 isolates. For this reason, the first task of this study was to determine the antibiotic resistance patterns in a wide collection of V. vulnificus strains belonging to the three biotypes that had been isolated worldwide from different sources (see Table S1 in the supplemental material). Isolates were screened for antimicrobial susceptibility to the antibiotics listed in Table S1 in the supplemental material by the agar diffusion disk procedure of Bauer et al. (5), according to the standard guideline (9). The resistance pattern found for each isolate is shown in Table S1 in the supplemental material. Less than 14% of isolates were sensitive to all the antibiotics tested, and more than 65% were resistant to more than one antibiotic, irrespective of their biotypes or serovars. The most frequent resistances were to ampicillin-sulbactam (SAM; 65.6% of the strains) and nitrofurantoin (F; 60.8% of the strains), and the least frequent were to tetracycline (12%) and oxytetracycline (8%). In addition, 15% of the strains were resistant to nalidixic acid (NAL) and oxolinic acid (OA), and 75% of these strains came from fish farms (see Table S1 in the supplemental material). Thus, high percentages of strains of the three biotypes were shown to be resistant to one or more antibiotics, with percentages similar to those found in nonbiotyped environmental V. vulnificus isolates from Asia and North America (4, 27, 34). In those studies, resistance to antibiotics could not be related to human contamination. However, the percentage of quinolone-resistant strains found in our study is higher than that reported in other ones, probably due to the inclusion of fish farm isolates, where the majority of quinolone-resistant strains were concentrated. This fact suggests that quinolone resistance could be related to human contamination due to the improper use of these drugs in therapy against fish diseases, as has been previously suggested (18, 20). Although no specific resistance pattern was associated with particular biotypes or serovars, we found certain differences in resistance distribution, as shown in Table ​Table1.1. In this respect, biotype 3 displayed the narrowest spectrum of resistances and biotype 1 the widest. The latter biotype encompassed the highest number of strains with multiresistance (see Table S1 in the supplemental material). Within biotype 2, there were differences among serovars, with quinolone resistance being restricted to the zoonotic serovar (Table ​(Table11).

TABLE 1.

Percentage of resistant strains distributed by biotypes and serovars

Retrograde Intraflagellar Transport Mutants Identify Complex A Proteins With Multiple Genetic Interactions in Chlamydomonas reinhardtii

The intraflagellar transport machinery is required for the assembly of cilia. It has been investigated by biochemical, genetic, and computational methods that have identified at least 21 proteins that assemble into two subcomplexes. It has been hypothesized that complex A is required for retrograde transport. Temperature-sensitive mutations in FLA15 and FLA17 show defects in retrograde intraflagellar transport (IFT) in Chlamydomonas. We show that IFT144 and IFT139, two complex A proteins, are encoded by FLA15 and FLA17, respectively. The fla15 allele is a missense mutation in a conserved cysteine and the fla17 allele is an in-frame deletion of three exons. The flagellar assembly defect of each mutant is rescued by the respective transgenes. In fla15 and fla17 mutants, bulges form in the distal one-third of the flagella at the permissive temperature and this phenotype is also rescued by the transgenes. These bulges contain the complex B component IFT74/72, but not α-tubulin or p28, a component of an inner dynein arm, which suggests specificity with respect to the proteins that accumulate in these bulges. IFT144 and IFT139 are likely to interact with each other and other proteins on the basis of three distinct genetic tests: (1) Double mutants display synthetic flagellar assembly defects at the permissive temperature, (2) heterozygous diploid strains exhibit second-site noncomplemention, and (3) transgenes confer two-copy suppression. Since these tests show different levels of phenotypic sensitivity, we propose they illustrate different gradations of gene interaction between complex A proteins themselves and with a complex B protein (IFT172).CILIA and flagella are microtubule-based organelles that are found on most mammalian cells. They provide motility to cells and participate in many sensory processes. Defects in or loss of cilia/flagella cause a variety of human diseases that include polycystic kidney disease, retinal degeneration, infertility, obesity, respiratory defects, left–right axis determination, and polydactyly (Fliegauf et al. 2007). Mouse mutants demonstrate that cilia are essential for viability, neural tube closure, and bone development (Eggenschwiler and Anderson 2007; Fliegauf et al. 2007). Cilia and flagella are also present in protists, algae, moss, and some fungi.The assembly and maintenance of cilia and flagella require intraflagellar transport (IFT) (Kozminski et al. 1995). IFT involves the movement of 100- to 200-nm-long protein particles from the basal body located in the cell body to the tip of the flagella using the heterotrimeric kinesin-2 (anterograde movement) (Kozminski et al. 1995) and movement back to the cell body (retrograde movement) using the cytoplasmic dynein complex (Pazour et al. 1999; Porter et al. 1999). IFT particles change their direction of movement as well as their size, speed, and frequency at the ends of the flagella as they switch from anterograde to retrograde movement (Iomini et al. 2001). Biochemical isolation of IFT particles reveals that they are composed of at least 16 proteins and that these particles can be dissociated into two complexes in vitro by changing the salt concentration (Cole et al. 1998; Piperno et al. 1998). Recent genetic and bioinformatics analysis adds at least 7 more proteins to the IFT particle (Follit et al. 2009) (Eggenschwiler and Anderson 2007).

TABLE 1

Proteins and gene names for the intraflagellar transport particles in Chlamydomonas, C. elegans, and mouse

Characterization of a Thermostable Short-Chain Alcohol Dehydrogenase from the Hyperthermophilic Archaeon Thermococcus sibiricus

Short-chain alcohol dehydrogenase, encoded by the gene Tsib_0319 from the hyperthermophilic archaeon Thermococcus sibiricus, was expressed in Escherichia coli, purified and characterized as an NADPH-dependent enantioselective oxidoreductase with broad substrate specificity. The enzyme exhibits extremely high thermophilicity, thermostability, and tolerance to organic solvents and salts.Alcohol dehydrogenases (ADHs; EC 1.1.1.1.) catalyze the interconversion of alcohols to their corresponding aldehydes or ketones by using different redox-mediating cofactors. NAD(P)-dependent ADHs, due to their broad substrate specificity and enantioselectivity, have attracted particular attention as catalysts in industrial processes (5). However, mesophilic ADHs are unstable at high temperatures, sensitive to organic solvents, and often lose activity during immobilization. In this relation, there is a considerable interest in ADHs from extremophilic microorganisms; among them, Archaea are of great interest. The representatives of all groups of NAD(P)-dependent ADHs have been detected in genomes of Archaea (11, 12); however, only a few enzymes have been characterized, and the great majority of them belong to medium-chain (3, 4, 14, 16, 19) or long-chain iron-activated ADHs (1, 8, 9). Up to now, a single short-chain archaeal ADH from Pyrococcus furiosus (10, 18) and only one archaeal aldo-keto reductase also from P. furiosus (11) have been characterized.Thermococcus sibiricus is a hyperthermophilic anaerobic archaeon isolated from a high-temperature oil reservoir capable of growth on complex organic substrates (15). The complete genome sequence of T. sibiricus has been recently determined and annotated (13). Several ADHs are encoded by the T. sibiricus genome, including three short-chain ADHs (Tsib_0319, Tsib_0703, and Tsib_1998) (13). In this report, we describe the cloning and expression of the Tsib_0319 gene from T. sibiricus and the purification and the biochemical characterization of its product, the thermostable short-chain ADH (TsAdh319).The Tsib_0319 gene encodes a protein with a size of 234 amino acids and the calculated molecular mass of 26.2 kDa. TsAdh319 has an 85% degree of sequence identity with short-chain ADH from P. furiosus (AdhA; PF_0074) (18). Besides AdhA, close homologs of TsAdh319 were found among different bacterial ADHs, but not archaeal ADHs. The gene flanked by the XhoI and BamHI sites was PCR amplified using two primers (sense primer, 5′-GTTCTCGAGATGAAGGTTGCTGTGATAACAGGG-3′, and antisense primer, 5′-GCTGGATCCTCAGTATTCTGGTCTCTGGTAGACGG-3′) and cloned into the pET-15b vector. TsAdh319 was overexpressed, with an N-terminal His₆ tag in Escherichia coli Rosetta-gami (DE3) and purified to homogeneity by metallochelating chromatography (Hi-Trap chelating HP column; GE Healthcare) followed by gel filtration on Superdex 200 10/300 GL column (GE Healthcare) equilibrated in 50 mM Tris-HCl (pH 7.5) with 200 mM NaCl. The homogeneity and the correspondence to the calculated molecular mass of 28.7 kDa were verified by SDS-PAGE (7). The molecular mass of native TsAdh319 was 56 to 60 kDa, which confirmed the dimeric structure in solution.The standard ADH activity measurement was made spectrophotometrically at the optimal pH by following either the reduction of NADP (in 50 mM Gly-NaOH buffer; pH 10.5) or the oxidation of NADPH (in 0.1 M sodium phosphate buffer; pH 7.5) at 340 nm at 60°C. The enzyme exhibited a strong preference for NADP(H) and broad substrate specificity (Table ​(Table1).1). The highest oxidation rates were found with pentoses d-arabinose (2.0 U mg⁻¹) and d-xylose (2.46 U mg⁻¹), and the highest reduction rates were found with dimethylglyoxal (5.9 U mg⁻¹) and pyruvaldehyde (2.2 U mg⁻¹). The enzyme did not reduce sugars which were good substrates for the oxidation reaction. The kinetic parameters of TsAdh319 determined for the preferred substrates are shown in Table ​Table2.2. The enantioselectivity of the enzyme was estimated by measuring the conversion rates of 2-butanol enantiomers. TsAdh319 showed an evident preference, >2-fold, for (S)-2-butanol over (RS)-2-butanol. The enzyme stereoselectivity is confirmed by the preferred oxidation of d-arabinose over l-arabinose (Table ​(Table1).1). The fact that TsAdh319 is metal independent was supported by the absence of a significant effect of TsAdh319 preincubation with 10 mM Me²⁺ for 30 min before measuring the activity in the presence of 1 mM Me²⁺ or EDTA (Table ​(Table3).3). TsAdh319 also exhibited a halophilic property, so the enzyme activity increased in the presence of NaCl and KCl and the activation was maintained even at concentration of 4 M and 3 M, respectively (Table ​(Table33).

TABLE 1.

Substrate specificity of TsAdh319

Quantitative Real-Time PCR Assays for Sensitive Detection of Canada Goose-Specific Fecal Pollution in Water Sources

Canada geese (Branta canadensis) are prevalent in North America and may contribute to fecal pollution of water systems where they congregate. This work provides two novel real-time PCR assays (CGOF1-Bac and CGOF2-Bac) allowing for the specific and sensitive detection of Bacteroides 16S rRNA gene markers present within Canada goose feces.The Canada goose (Branta canadensis) is a prevalent waterfowl species in North America. The population density of Canada geese has doubled during the past 15 years, and the population was estimated to be close to 3 million in 2007 (4). Canada geese often congregate within urban settings, likely due to available water sources, predator-free grasslands, and readily available food supplied by humans (6). They are suspected to contribute to pollution of aquatic environments due to the large amounts of fecal matter that can be transported into the water. This can create a public health threat if the fecal droppings contain pathogenic microorganisms (6, 7, 9, 10, 12, 13, 19). Therefore, tracking transient fecal pollution of water due to fecal inputs from waterfowl, such as Canada geese, is of importance for protecting public health.PCR detection of host-specific 16S rRNA gene sequences from Bacteroidales of fecal origin has been described as a promising microbial source-tracking (MST) approach due to its rapidity and high specificity (2, 3). Recently, Lu et al. (15) characterized the fecal microbial community from Canada geese by constructing a 16S rRNA gene sequence database using primers designed to amplify all bacterial 16S rRNA gene sequences. The authors reported that the majority of the 16S rRNA gene sequences obtained were related to Clostridia or Bacilli and to a lesser degree Bacteroidetes, which represent possible targets for host-specific source-tracking assays.The main objective of this study was to identify novel Bacteroidales 16S rRNA gene sequences that are specific to Canada goose feces and design primers and TaqMan fluorescent probes for sensitive and specific quantification of Canada goose fecal contamination in water sources.Primers 32F and 708R from Bernhard and Field (2) were used to construct a Bacteroidales-specific 16S rRNA gene clone library from Canada goose fecal samples (n = 15) collected from grass lawns surrounding Wascana Lake (Regina, SK, Canada) in May 2009 (for a detailed protocol, see File S1 in the supplemental material). Two hundred eighty-eight clones were randomly selected and subjected to DNA sequencing (at the Plant Biotechnology Institute DNA Technologies Unit, Saskatoon, SK, Canada). Representative sequences of each operational taxonomic unit (OTU) were recovered using an approach similar to that described by Mieszkin et al. (16). Sequences that were less than 93% similar to 16S rRNA gene sequences from nontarget host species in GenBank were used in multiple alignments to identify regions of DNA sequence that were putatively goose specific. Subsequently, two TaqMan fluorescent probe sets (targeting markers designated CGOF1-Bac and CGOF2-Bac) were designed using the RealTimeDesign software provided by Biosearch Technologies (http://www.biosearchtech.com/). The newly designed primer and probe set for the CGOF1-Bac assay included CG1F (5′-GTAGGCCGTGTTTTAAGTCAGC-3′) and CG1R (5′-AGTTCCGCCTGCCTTGTCTA-3′) and a TaqMan probe (5′-6-carboxyfluorescein [FAM]-CCGTGCCGTTATACTGAGACACTTGAG-Black Hole Quencher 1 [BHQ-1]-3′), and the CGOF2-Bac assay had primers CG2F (5′-ACTCAGGGATAGCCTTTCGA-3′) and CG2R (5′-ACCGATGAATCTTTCTTTGTCTCC-3′) and a TaqMan probe (5′-FAM-AATACCTGATGCCTTTGTTTCCCTGCA-BHQ-1-3′). Oligonucleotide specificities for the Canada goose-associated Bacteroides 16S rRNA primers were verified through in silico analysis using BLASTN (1) and the probe match program of the Ribosomal Database Project (release 10) (5). Host specificity was further confirmed using DNA extracts from 6 raw human sewage samples from various geographical locations in Saskatchewan and 386 fecal samples originating from 17 different animal species in Saskatchewan, including samples from Canada geese (n = 101) (Table ​(Table1).1). An existing nested PCR assay for detecting Canada goose feces (15) (targeting genetic marker CG-Prev f5) (see Table S1 in the supplemental material) was also tested for specificity using the individual fecal and raw sewage samples (Table ​(Table1).1). All fecal DNA extracts were obtained from 0.25 g of fecal material by using the PowerSoil DNA extraction kit (Mo Bio Inc., Carlsbad, CA) (File S1 in the supplemental material provides details on the sample collection).

TABLE 1.

Specificities of the CGOF1-Bac, CGOF2-Bac, and CG-Prev f5 PCR assays for different species present in Saskatchewan, Canada

Quantitative Fluorescence In Situ Hybridization of Microbial Communities in the Rumens of Cattle Fed Different Diets

At present there is little quantitative information on the identity and composition of bacterial populations in the rumen microbial community. Quantitative fluorescence in situ hybridization using newly designed oligonucleotide probes was applied to identify the microbial populations in liquid and solid fractions of rumen digesta from cows fed barley silage or grass hay diets with or without flaxseed. Bacteroidetes, Firmicutes, and Proteobacteria were abundant in both fractions, constituting 31.8 to 87.3% of the total cell numbers. They belong mainly to the order Bacteroidales (0.1 to 19.2%), hybridizing with probe BAC1080; the families Lachnospiraceae (9.3 to 25.5%) and Ruminococcaceae (5.5 to 23.8%), hybridizing with LAC435 and RUM831, respectively; and the classes Deltaproteobacteria (5.8 to 28.3%) and Gammaproteobacteria (1.2 to 8.2%). All were more abundant in the rumen communities of cows fed diets containing silage (75.2 to 87.3%) than in those of cows fed diets containing hay (31.8 to 49.5%). The addition of flaxseed reduced their abundance in the rumens of cows fed silage-based diets (to 45.2 to 58.7%) but did not change markedly their abundance in the rumens of cows fed hay-based diets (31.8 to 49.5%). Fibrolytic species, including Fibrobacter succinogenes and Ruminococcus spp., and archaeal methanogens accounted for only a small proportion (0.4 to 2.1% and 0.2 to 0.6%, respectively) of total cell numbers. Depending on diet, between 37.0 and 91.6% of microbial cells specifically hybridized with the probes used in this study, allowing them to be identified in situ. The identities of other microbial populations (8.4 to 63.0%) remain unknown.The rumen is an anaerobic ecosystem used by herbivores to convert fibrous plant material into fermentation products that are in turn used as energy by the host. Fibrolytic degradation is accomplished by a complex microbial community which includes specialized fungi, protozoa, and bacteria (14). More than 200 bacterial species (5) have been isolated from rumen, and many of these have been phylogenetically and physiologically characterized. Several of these, including Fibrobacter succinogenes, Ruminococcus albus, and Ruminococcus flavefaciens, have the ability to hydrolyze cellulose in axenic culture (24). Despite the presence of these fibrolytic populations, a large portion of the fiber in low-quality forage diets passes through the rumen undigested. In the rumen, fibrolytic bacteria do not digest plant cell walls in isolation but rather interact with a consortium of bacteria (18). Although culture-dependent studies have improved our understanding of rumen microbiology, the importance of the isolates to the structure and function of the rumen microbial community, with the possible exception of the fibrolytic strains, is still unknown. Expanding our knowledge of the structure and function of the rumen microbial community may provide insights into approaches to improve the efficiency of fiber digestion and biofuel production (14).To provide a high-resolution view of the population structure of the rumen bacterial community, we used quantitative fluorescence in situ hybridization (qFISH) to investigate the composition and distribution of bacterial populations associated with the liquid and solid rumen contents from 12 ruminally cannulated Holstein dairy cows (3 cows were used for each diet) fed (for at least 21 days) grass hay or barley silage diets with or without flaxseed (Table ​(Table1).1). Six new 16S rRNA-targeted FISH probes (Table ​(Table2)2) for not only the fibrolytic groups but also other unclassified bacterial groups in the rumen were designed, using ARB software (17), against the rumen 16S rRNA gene sequences (data not shown) retrieved from the Ribosomal Database Project (RDP) database (6). The new probes target Bacteroidales-related clones (probe BAC1080) (phylum Bacteroidetes), Lachnospiraceae- and Ruminococcaceae-related clones (probes LAC435 and RUM831, respectively) (phylum Firmicutes), Butyrivibrio fibrisolvens-related clones (probe BFI826), and R. albus- and R. flavefaciens-related clones (probes RAL1436 and RFL155, respectively).

TABLE 1.

Composition of diets used in this study

Hierarchical Oligonucleotide Primer Extension as a Time- and Cost-Effective Approach for Quantitative Determination of Bifidobacterium spp. in Infant Feces

The Bifidobacterium spp. present in 10 infant fecal samples (4 from infants with eczema and 6 from healthy infants) were quantified with both hierarchical oligonucleotide primer extension (HOPE) and fluorescence in situ hybridization-flow cytometry. The relative abundances of Bifidobacterium longum and B. catenulatum with respect to the total bifidobacteria had a poor correlation (ρ, <0.600; P value, >0.208), presumably due to differences in primer specificity and the level of hybridization stringency of both methods. In contrast, the relative abundances of organisms of the genus Bifidobacterium against the total amplified 16S rRNA genes and those of B. adolescentis, B. bifidum, and B. breve against the genus Bifidobacterium exhibited a good statistical correlation (ρ, >0.783; P value, <0.066). This good comparability supports HOPE as a method to achieve high-throughput quantitative determination of bacterial targets in a time- and cost-effective manner.The “microflora hygiene” hypothesis states that a lack of exposure to pathogens or certain commensal bacteria in early life may predispose some individuals to allergic disorders (14). However, inconsistent findings on the abundance of health-associated microbes have prevented precise conclusions as to their role in modulating host health. For example, by performing fluorescence in situ hybridization (FISH) on infant feces, Bifidobacterium spp. were found in high abundance in healthy infants (11). In contrast, certain species, like Bifidobacterium pseudocatenulatum, may be more commonly detected in infants with eczema (3). Therefore, to facilitate our understanding of microbial composition and its correlation to human health, it is essential to use a rapid and high-throughput molecular method to determine the abundances of bacterial targets in a large sample size (16). Although FISH-flow cytometry (FISH-FC) is routinely used to quantify the abundances of bacterial targets in feces (9, 11, 17), it does not suffice as a high-throughput method due to the limited range of spectrally distinct fluorophores that are available in the UV spectrum (10, 13). There is a need to develop a high-throughput technique which can complement the existing molecular methods to rapidly evaluate the relative abundance of bacterial targets.A molecular method termed hierarchical oligonucleotide primer extension (HOPE) was developed to rapidly determine the relative abundances of bacterial 16S rRNA genes among total PCR-amplified 16S rRNA genes (19). HOPE uses primers of different lengths that were designed to target bacteria at different phylogenetic levels. The primers anneal to complementary regions of the targeted bacteria and extend with a fluorophore-labeled nucleotide when the bacterial target is present. The extended primers can be differentiated on a genetic analyzer based on primer length and fluorophore color. The relative abundance of the bacterial target against a higher-level primer can then be quantified by calculating the ratio of the peak area of the extended primer with that of a higher-level primer. A subsequent study demonstrated that HOPE can be used for rapid and specific determination of Bacteroides spp. present in feces and wastewaters at different taxonomical levels (5). It also has the versatility to be expanded to include other bacterial groups. This can potentially facilitate the identification and quantification of bacterial populations that modulate the health of an individual at different temporal intervals.This study aimed to demonstrate HOPE as a time- and cost-effective method to quantify the abundances of Bifidobacterium spp. in 10 infant fecal samples (4 from infants with eczema and 6 from healthy infants) that were collected at 1, 3, and 12 months of age. The abundances of the Bifidobacterium spp. as determined, respectively, by HOPE and FISH-FC were also compared to validate the use of HOPE as a quantitative method.To obtain the total PCR-amplified 16S rRNA genes, genomic DNA of the fecal microbiota was extracted based on a previously described protocol (12) prior to 20 cycles of PCR amplification (modified 11F forward primer 5′-GTT YGA TYC TGG CTC AG-3′ and 1492R reverse primer 5′-GGY TAC CTT GTT ACG ACT T-3′) (6, 7). The amplicons were purified, and the concentrations were diluted to 10 ng/μl. For HOPE, a total of 12 primers specifically targeting six Bifidobacterium spp. at different taxonomic levels were designed based on a previously described protocol (5). The specificity of the designed primers was verified in silico against entries in RDP II (2), and the sensitivity of the primers was determined as described previously (19). FISH-FC was performed on the same set of fecal samples based on the protocol described by Lay et al. (9). Table ​Table11 lists the HOPE primers and FISH probes used in this study. A nonparametric Spearson ranked correlation analysis (Minitab) was performed on the abundances of Bifidobacterium spp. as quantified by FISH-FC and HOPE, respectively.

TABLE 1.

HOPE primers and FISH probes included in this study