Arabidopsis LON2 Is Necessary for Peroxisomal Function and Sustained Matrix Protein Import

Relatively little is known about the small subset of peroxisomal proteins with predicted protease activity. Here, we report that the peroxisomal LON2 (At5g47040) protease facilitates matrix protein import into Arabidopsis (Arabidopsis thaliana) peroxisomes. We identified T-DNA insertion alleles disrupted in five of the nine confirmed or predicted peroxisomal proteases and found only two—lon2 and deg15, a mutant defective in the previously described PTS2-processing protease (DEG15/At1g28320)—with phenotypes suggestive of peroxisome metabolism defects. Both lon2 and deg15 mutants were mildly resistant to the inhibitory effects of indole-3-butyric acid (IBA) on root elongation, but only lon2 mutants were resistant to the stimulatory effects of IBA on lateral root production or displayed Suc dependence during seedling growth. lon2 mutants displayed defects in removing the type 2 peroxisome targeting signal (PTS2) from peroxisomal malate dehydrogenase and reduced accumulation of 3-ketoacyl-CoA thiolase, another PTS2-containing protein; both defects were not apparent upon germination but appeared in 5- to 8-d-old seedlings. In lon2 cotyledon cells, matrix proteins were localized to peroxisomes in 4-d-old seedlings but mislocalized to the cytosol in 8-d-old seedlings. Moreover, a PTS2-GFP reporter sorted to peroxisomes in lon2 root tip cells but was largely cytosolic in more mature root cells. Our results indicate that LON2 is needed for sustained matrix protein import into peroxisomes. The delayed onset of matrix protein sorting defects may account for the relatively weak Suc dependence following germination, moderate IBA-resistant primary root elongation, and severe defects in IBA-induced lateral root formation observed in lon2 mutants.Peroxisomes are single-membrane-bound organelles found in most eukaryotes. Peroxin (PEX) proteins are necessary for various aspects of peroxisome biogenesis, including matrix protein import (for review, see Distel et al., 1996; Schrader and Fahimi, 2008). Most matrix proteins are imported into peroxisomes from the cytosol using one of two targeting signals, a C-terminal type 1 peroxisome-targeting signal (PTS1) or a cleavable N-terminal type 2 peroxisome-targeting signal (PTS2) (Reumann, 2004). PTS1- and PTS2-containing proteins are bound in the cytosol by soluble matrix protein receptors, escorted to the peroxisome membrane docking complex, and translocated into the peroxisome matrix (for review, see Platta and Erdmann, 2007). Once in the peroxisome, many matrix proteins participate in metabolic pathways, such as β-oxidation, hydrogen peroxide decomposition, and photorespiration (for review, see Gabaldon et al., 2006; Poirier et al., 2006).In addition to metabolic enzymes, several proteases are found in the peroxisome matrix. Only one protease, DEG15/Tysnd1, has a well-defined role in peroxisome biology. The rat Tysnd1 protease removes the targeting signal after PTS2-containing proteins enter the peroxisome and also processes certain PTS1-containing β-oxidation enzymes (Kurochkin et al., 2007). Similarly, the Arabidopsis (Arabidopsis thaliana) Tysnd1 homolog DEG15 (At1g28320) is a peroxisomal Ser protease that removes PTS2 targeting signals (Helm et al., 2007; Schuhmann et al., 2008).In contrast with DEG15, little is known about the other eight Arabidopsis proteins that are annotated as proteases in the AraPerox database of putative peroxisomal proteins (Reumann et al., 2004; Carter et al., 2004; Shimaoka et al., 2004), which, in combination with the minor PTS found in both of these predicted proteases (Reumann, 2004), suggests that these enzymes may not be peroxisomal. Along with DEG15, only two of the predicted peroxisomal proteases, an M16 metalloprotease (At2g41790), which we have named PXM16 for peroxisomal M16 protease, and a Lon-related protease (At5g47040/LON2; Ostersetzer et al., 2007), are found in the proteome of peroxisomes purified from Arabidopsis suspension cells (Eubel et al., 2008). DEG15 and LON2 also have been validated as peroxisomally targeted using GFP fusions (Ostersetzer et al., 2007; Schuhmann et al., 2008).

Table I.

Putative Arabidopsis proteases predicted or demonstrated to be peroxisomal

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

Adaptive Divergence in Experimental Populations of Pseudomonas fluorescens. IV. Genetic Constraints Guide Evolutionary Trajectories in a Parallel Adaptive Radiation

The capacity for phenotypic evolution is dependent upon complex webs of functional interactions that connect genotype and phenotype. Wrinkly spreader (WS) genotypes arise repeatedly during the course of a model Pseudomonas adaptive radiation. Previous work showed that the evolution of WS variation was explained in part by spontaneous mutations in wspF, a component of the Wsp-signaling module, but also drew attention to the existence of unknown mutational causes. Here, we identify two new mutational pathways (Aws and Mws) that allow realization of the WS phenotype: in common with the Wsp module these pathways contain a di-guanylate cyclase-encoding gene subject to negative regulation. Together, mutations in the Wsp, Aws, and Mws regulatory modules account for the spectrum of WS phenotype-generating mutations found among a collection of 26 spontaneously arising WS genotypes obtained from independent adaptive radiations. Despite a large number of potential mutational pathways, the repeated discovery of mutations in a small number of loci (parallel evolution) prompted the construction of an ancestral genotype devoid of known (Wsp, Aws, and Mws) regulatory modules to see whether the types derived from this genotype could converge upon the WS phenotype via a novel route. Such types—with equivalent fitness effects—did emerge, although they took significantly longer to do so. Together our data provide an explanation for why WS evolution follows a limited number of mutational pathways and show how genetic architecture can bias the molecular variation presented to selection.UNDERSTANDING—and importantly, predicting—phenotypic evolution requires knowledge of the factors that affect the translation of mutation into phenotypic variation—the raw material of adaptive evolution. While much is known about mutation rate (e.g., Drake et al. 1998; Hudson et al. 2002), knowledge of the processes affecting the translation of DNA sequence variation into phenotypic variation is minimal.Advances in knowledge on at least two fronts suggest that progress in understanding the rules governing the generation of phenotypic variation is possible (Stern and Orgogozo 2009). The first stems from increased awareness of the genetic architecture underlying specific adaptive phenotypes and recognition of the fact that the capacity for evolutionary change is likely to be constrained by this architecture (Schlichting and Murren 2004; Hansen 2006). The second is the growing number of reports of parallel evolution (e.g., Pigeon et al. 1997; ffrench-Constant et al. 1998; Allender et al. 2003; Colosimo et al. 2004; Zhong et al. 2004; Boughman et al. 2005; Shindo et al. 2005; Kronforst et al. 2006; Woods et al. 2006; Zhang 2006; Bantinaki et al. 2007; McGregor et al. 2007; Ostrowski et al. 2008)—that is, the independent evolution of similar or identical features in two or more lineages—which suggests the possibility that evolution may follow a limited number of pathways (Schluter 1996). Indeed, giving substance to this idea are studies that show that mutations underlying parallel phenotypic evolution are nonrandomly distributed and typically clustered in homologous genes (Stern and Orgogozo 2008).While the nonrandom distribution of mutations during parallel genetic evolution may reflect constraints due to genetic architecture, some have argued that the primary cause is strong selection (e.g., Wichman et al. 1999; Woods et al. 2006). A means of disentangling the roles of population processes (selection) from genetic architecture is necessary for progress (Maynard Smith et al. 1985; Brakefield 2006); also necessary is insight into precisely how genetic architecture might bias the production of mutations presented to selection.Despite their relative simplicity, microbial populations offer opportunities to advance knowledge. The wrinkly spreader (WS) morphotype is one of many different niche specialist genotypes that emerge when experimental populations of Pseudomonas fluorescens are propagated in spatially structured microcosms (Rainey and Travisano 1998). Previous studies defined, via gene inactivation, the essential phenotypic and genetic traits that define a single WS genotype known as LSWS (Spiers et al. 2002, 2003) (Figure 1). LSWS differs from the ancestral SM genotype by a single nonsynonymous nucleotide change in wspF. Functionally (see Figure 2), WspF is a methyl esterase and negative regulator of the WspR di-guanylate cyclase (DGC) (Goymer et al. 2006) that is responsible for the biosynthesis of c-di-GMP (Malone et al. 2007), the allosteric activator of cellulose synthesis enzymes (Ross et al. 1987). The net effect of the wspF mutation is to promote physiological changes that lead to the formation of a microbial mat at the air–liquid interface of static broth microcosms (Rainey and Rainey 2003).Open in a separate windowFigure 1.—Outline of experimental strategy for elucidation of WS-generating mutations and their subsequent identity and distribution among a collection of independently evolved, spontaneously arising WS genotypes. The strategy involves, first, the genetic analysis of a specific WS genotype (e.g., LSWS) to identify the causal mutation, and second, a survey of DNA sequence variation at specific loci known to harbor causal mutations among a collection of spontaneously arising WS genotypes. For example, suppressor analysis of LSWS using a transposon to inactivate genes necessary for expression of the wrinkly morphology delivered a large number of candidate genes (top left) (Spiers et al. 2002). Genetic and functional analysis of these candidate genes (e.g., Goymer et al. 2006) led eventually to the identity of the spontaneous mutation (in wspF) responsible for the evolution of LSWS from the ancestral SM genotype (Bantinaki et al. 2007). Subsequent analysis of the wspF sequence among 26 independent WS genotypes (bottom) showed that 50% harbored spontaneous mutations (of different kinds; see Open in a separate windowFigure 2.—Network diagram of DGC-encoding pathways underpinning the evolution of the WS phenotype and their regulation. Overproduction of c-di-GMP results in overproduction of cellulose and other adhesive factors that determine the WS phenotype. The ancestral SBW25 genome contains 39 putative DGCs, each in principle capable of synthesizing the production of c-di-GMP, and yet WS genotypes arise most commonly as a consequence of mutations in just three DGC-containing pathways: Wsp, Aws, and Mws. In each instance, the causal mutations are most commonly in the negative regulatory component: wspF, awsX, and the phosphodiesterase domain of mwsR (see text).To determine whether spontaneous mutations in wspF are a common cause of the WS phenotype, the nucleotide sequence of this gene was obtained from a collection of 26 spontaneously arising WS genotypes (WS_A-Z) taken from 26 independent adaptive radiations, each founded by the same ancestral SM genotype (Figure 1): 13 contained mutations in wspF (Bantinaki et al. 2007). The existence of additional mutational pathways to WS provided the initial motivation for this study.

TABLE 1

Mutational causes of WS

Efficient Transposition of Tol2 in the Mouse Germline

Insertional mutagenesis screens play an integral part in the annotating of functional data for all sequenced genes in the postgenomic era. Chemical mutagenesis screens are highly efficient but identifying the causative gene can be a laborious task. Other mutagenesis platforms, such as transposable elements, have been successfully applied for insertional mutagenesis screens in both the mouse and rat. However, relatively low transposition efficiency has hampered their use as a high-throughput forward genetic mutagenesis screen. Here we report the first evidence of germline activity in the mouse using a naturally active DNA transposon derived from the medaka fish called Tol2, as an alternative system for high-throughput forward genetic mutagenesis screening tool.THE Tol2 (transposable element of Oryzias latipes, number 2) element belongs to the hAT family (hobo of Drosophilia, Activator of maize and Tam3 of snapdragon) of transposons and was the first known autonomously active vertebrate type II transposable element (Koga et al. 1996; Kawakami et al. 1998). Unlike other DNA-type transposons like Sleeping Beauty (SB) (Ivics et al. 1997) or piggyBac (PB) (Fraser et al. 1996), Tol2 does not exhibit any known strong site specificity for integration nor does it exhibit any significant overexpression inhibition activity (Kawakami and Noda 2004; Balciunas et al. 2006) as seen in SB (Geurts et al. 2003). Recently, Tol2 was shown to effectively carry large DNA cargo of up to 10 kb in human and mouse cells without affecting its transposition efficiency (Balciunas et al. 2006). To date, Tol2 has also been demonstrated to transpose efficiently in zebrafish, frog, chicken, mouse cells, and human cells (Kawakami et al. 2000, 2004; Koga et al. 2003; Kawakami and Noda 2004; Balciunas et al. 2006; Hamlet et al. 2006; Sato et al. 2007).Germline mutagenesis using the SB transposon system has been demonstrated in both the mouse (Dupuy et al. 2001; Horie et al. 2001) and rat (Kitada et al. 2007; Lu et al. 2007). In addition, PB germline mutagenesis in mice has also been demonstrated (Ding et al. 2005; Wu et al. 2007). However, the relatively low germline transposition efficiency of both transposon systems reported so far has hampered their use in a high-throughput forward genetic mutagenesis screen (Keng et al. 2005; Kitada et al. 2007).

TABLE 1

Germline transposition frequency in various transposon systems

Evidence for a New Avian Paramyxovirus Serotype 10 Detected in Rockhopper Penguins from the Falkland Islands

The biological, serological, and genomic characterization of a paramyxovirus recently isolated from rockhopper penguins (Eudyptes chrysocome) suggested that this virus represented a new avian paramyxovirus (APMV) group, APMV10. This penguin virus resembled other APMVs by electron microscopy; however, its viral hemagglutination (HA) activity was not inhibited by antisera against any of the nine defined APMV serotypes. In addition, antiserum generated against this penguin virus did not inhibit the HA of representative viruses of the other APMV serotypes. Sequence data produced using random priming methods revealed a genomic structure typical of APMV. Phylogenetic evaluation of coding regions revealed that amino acid sequences of all six proteins were most closely related to APMV2 and APMV8. The calculation of evolutionary distances among proteins and distances at the nucleotide level confirmed that APMV2, APMV8, and the penguin virus all were sufficiently divergent from each other to be considered different serotypes. We propose that this isolate, named APMV10/penguin/Falkland Islands/324/2007, be the prototype virus for APMV10. Because of the known problems associated with serology, such as antiserum cross-reactivity and one-way immunogenicity, in addition to the reliance on the immune response to a single protein, the hemagglutinin-neuraminidase, as the sole base for viral classification, we suggest the need for new classification guidelines that incorporate genome sequence comparisons.Viruses from the Paramyxoviridae family have caused disease in humans and animals for centuries. Over the last 40 years, many paramyxoviruses isolated from animals and people have been newly described (16, 17, 22, 29, 31, 32, 36, 42, 44, 46, 49, 58, 59, 62-64). Viruses from this family are pleomorphic, enveloped, single-stranded, nonsegmented, negative-sense RNA viruses that demonstrate serological cross-reactivity with other paramyxoviruses related to them (30, 46). The subfamily Paramyxovirinae is divided into five genera: Respirovirus, Morbillivirus, Rubulavirus, Henipavirus, and Avulavirus (30). The Avulavirus genus contains nine distinct avian paramyxovirus (APMV) serotypes (Table ​(Table1),1), and information on the discovery of each has been reported elsewhere (4, 6, 7, 9, 12, 34, 41, 50, 51, 60, 68).

TABLE 1.

Characteristics of prototype viruses APMV1 to APMV9 and the penguin virus

Comparative Analysis of Myxococcus Predation on Soil Bacteria

Predator-prey relationships among prokaryotes have received little attention but are likely to be important determinants of the composition, structure, and dynamics of microbial communities. Many species of the soil-dwelling myxobacteria are predators of other microbes, but their predation range is poorly characterized. To better understand the predatory capabilities of myxobacteria in nature, we analyzed the predation performance of numerous Myxococcus isolates across 12 diverse species of bacteria. All predator isolates could utilize most potential prey species to effectively fuel colony expansion, although one species hindered predator swarming relative to a control treatment with no growth substrate. Predator strains varied significantly in their relative performance across prey types, but most variation in predatory performance was determined by prey type, with Gram-negative prey species supporting more Myxococcus growth than Gram-positive species. There was evidence for specialized predator performance in some predator-prey combinations. Such specialization may reduce resource competition among sympatric strains in natural habitats. The broad prey range of the Myxococcus genus coupled with its ubiquity in the soil suggests that myxobacteria are likely to have very important ecological and evolutionary effects on many species of soil prokaryotes.Predation plays a major role in shaping both the ecology and evolution of biological communities. The population and evolutionary dynamics of predators and their prey are often tightly coupled and can greatly influence the dynamics of other organisms as well (1). Predation has been invoked as a major cause of diversity in ecosystems (11, 12). For example, predators may mediate coexistence between superior and inferior competitors (2, 13), and differential trajectories of predator-prey coevolution can lead to divergence between separate populations (70).Predation has been investigated extensively in higher organisms but relatively little among prokaryotes. Predation between prokaryotes is one of the most ancient forms of predation (27), and it has been proposed that this process may have been the origin of eukaryotic cells (16). Prokaryotes are key players in primary biomass production (44) and global nutrient cycling (22), and predation of some prokaryotes by others is likely to significantly affect these processes. Most studies of predatory prokaryotes have focused on Bdellovibrionaceae species (e.g., see references 51, 55, and 67). These small deltaproteobacteria prey on other Gram-negative cells, using flagella to swim rapidly until they collide with a prey cell. After collision, the predator cells then enter the periplasmic space of the prey cell, consume the host cell from within, elongate, and divide into new cells that are released upon host cell lysis (41). Although often described as predatory, the Bdellovibrionaceae may also be considered to be parasitic, as they typically depend (apart from host-independent strains that have been observed [60]) on the infection and death of their host for their reproduction (47).In this study, we examined predation among the myxobacteria, which are also deltaproteobacteria but constitute a monophyletic clade divergent from the Bdellovibrionaceae (17). Myxobacteria are found in most terrestrial soils and in many aquatic environments as well (17, 53, 74). Many myxobacteria, including the model species Myxococcus xanthus, exhibit several complex social traits, including fruiting body formation and spore formation (14, 18, 34, 62, 71), cooperative swarming with two motility systems (64, 87), and group (or “wolf pack”) predation on both bacteria and fungi (4, 5, 8, 9, 15, 50). Using representatives of the genus Myxococcus, we tested for both intra- and interspecific variation in myxobacterial predatory performance across a broad range of prey types. Moreover, we examined whether prey vary substantially in the degree to which they support predatory growth by the myxobacteria and whether patterns of variation in predator performance are constant or variable across prey environments. The latter outcome may reflect adaptive specialization and help to maintain diversity in natural populations (57, 59).Although closely related to the Bdellovibrionaceae (both are deltaproteobacteria), myxobacteria employ a highly divergent mode of predation. Myxobacteria use gliding motility (64) to search the soil matrix for prey and produce a wide range of antibiotics and lytic compounds that kill and decompose prey cells and break down complex polymers, thereby releasing substrates for growth (66). Myxobacterial predation is cooperative both in its “searching” component (6, 31, 82; for details on cooperative swarming, see reference 64) and in its “handling” component (10, 29, 31, 32), in which secreted enzymes turn prey cells into consumable growth substrates (56, 83). There is evidence that M. xanthus employs chemotaxis-like genes in its attack on prey cells (5) and that predation is stimulated by close contact with prey cells (48).Recent studies have revealed great genetic and phenotypic diversity within natural populations of M. xanthus, on both global (79) and local (down to centimeter) scales (78). Phenotypic diversity includes variation in social compatibility (24, 81), the density and nutrient thresholds triggering development (33, 38), developmental timing (38), motility rates and patterns (80), and secondary metabolite production (40). Although natural populations are spatially structured and both genetic diversity and population differentiation decrease with spatial scale (79), substantial genetic diversity is present even among centimeter-scale isolates (78). No study has yet systematically investigated quantitative natural variation in myxobacterial predation phenotypes across a large number of predator genotypes.Given the previous discovery of large variation in all examined phenotypes, even among genetically extremely similar strains, we anticipated extensive predatory variation as well. Using a phylogenetically broad range of prey, we compared and contrasted the predatory performance of 16 natural M. xanthus isolates, sampled from global to local scales, as well as the commonly studied laboratory reference strain DK1622 and representatives of three additional Myxococcus species: M. flavescens (86), M. macrosporus (42), and M. virescens (63) (Table ​(Table1).1). In particular, we measured myxobacterial swarm expansion rates on prey lawns spread on buffered agar (31, 50) and on control plates with no nutrients or with prehydrolyzed growth substrate.

TABLE 1.

List of myxobacteria used, with geographical origin

Natural Infection of Burkholderia pseudomallei in an Imported Pigtail Macaque (Macaca nemestrina) and Management of the Exposed Colony

Identification of the select agent Burkholderia pseudomallei in macaques imported into the United States is rare. A purpose-bred, 4.5-y-old pigtail macaque (Macaca nemestrina) imported from Southeast Asia was received from a commercial vendor at our facility in March 2012. After the initial acclimation period of 5 to 7 d, physical examination of the macaque revealed a subcutaneous abscess that surrounded the right stifle joint. The wound was treated and resolved over 3 mo. In August 2012, 2 mo after the stifle joint wound resolved, the macaque exhibited neurologic clinical signs. Postmortem microbiologic analysis revealed that the macaque was infected with B. pseudomallei. This case report describes the clinical evaluation of a B. pseudomallei-infected macaque, management and care of the potentially exposed colony of animals, and protocols established for the animal care staff that worked with the infected macaque and potentially exposed colony. This article also provides relevant information on addressing matters related to regulatory issues and risk management of potentially exposed animals and animal care staff.Abbreviations: CDC, Centers for Disease Control and Prevention; IHA, indirect hemagglutination assay; PEP, postexposure prophylacticBurkholderia pseudomallei, formerly known as Pseudomonas pseudomallei, is a gram-negative, aerobic, bipolar, motile, rod-shaped bacterium. B. pseudomallei infections (melioidosis) can be severe and even fatal in both humans and animals. This environmental saprophyte is endemic to Southeast Asia and northern Australia, but it has also been found in other tropical and subtropical areas of the world.^7,22,32,42 The bacterium is usually found in soil and water in endemic areas and is transmitted to humans and animals primarily through percutaneous inoculation, ingestion, or inhalation of a contaminated source.^{8, 22,28,32,42} Human-to-human, animal-to-animal, and animal-to-human spread are rare.^8,32 In December 2012, the National Select Agent Registry designated B. pseudomallei as a Tier 1 overlap select agent.³⁹ Organisms classified as Tier 1 agents present the highest risk of deliberate misuse, with the most significant potential for mass casualties or devastating effects to the economy, critical infrastructure, or public confidence. Select agents with this status have the potential to pose a severe threat to human and animal health or safety or the ability to be used as a biologic weapon.³⁹Melioidosis in humans can be challenging to diagnose and treat because the organism can remain latent for years and is resistant to many antibiotics.^12,37,41 B. pseudomallei can survive in phagocytic cells, a phenomenon that may be associated with latent infections.^19,38 The incubation period in naturally infected animals ranges from 1 d to many years, but symptoms typically appear 2 to 4 wk after exposure.^13,17,35,38 Disease generally presents in 1 of 2 forms: localized infection or septicemia.²² Multiple methods are used to diagnose melioidosis, including immunofluorescence, serology, and PCR analysis, but isolation of the bacteria from blood, urine, sputum, throat swabs, abscesses, skin, or tissue lesions remains the ‘gold standard.’^9,22,40,42 The prognosis varies based on presentation, time to diagnosis, initiation of appropriate antimicrobial treatment, and underlying comorbidities.^7,28,42 Currently, there is no licensed vaccine to prevent melioidosis.There are several published reports of naturally occurring melioidosis in a variety of nonhuman primates (NHP; 2,10,13,17,25,30,31,35 The first reported case of melioidosis in monkeys was recorded in 1932, and the first published case in a macaque species was in 1966.³⁰ In the United States, there have only been 7 documented cases of NHP with B. pseudomallei infection.^2,13,17 All of these cases occurred prior to the classification of B. pseudomallei as a select agent. Clinical signs in NHP range from subclinical or subacute illness to acute septicemia, localized infection, and chronic infection. NHP with melioidosis can be asymptomatic or exhibit clinical signs such as anorexia, wasting, purulent drainage, subcutaneous abscesses, and other soft tissue lesions. Lymphadenitis, lameness, osteomyelitis, paralysis and other CNS signs have also been reported.^{2,7,10,22,28,32} In comparison, human''s clinical signs range from abscesses, skin ulceration, fever, headache, joint pain, and muscle tenderness to abdominal pain, anorexia, respiratory distress, seizures, and septicemia.^7,9,21,22

Table 1.

Summary of reported cases of naturally occurring Burkholderia pseudomalleiinfections in nonhuman primates

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

TrichOME: A Comparative Omics Database for Plant Trichomes

Plant secretory trichomes have a unique capacity for chemical synthesis and secretion and have been described as biofactories for the production of natural products. However, until recently, most trichome-specific metabolic pathways and genes involved in various trichome developmental stages have remained unknown. Furthermore, only a very limited amount of plant trichome genomics information is available in scattered databases. We present an integrated “omics” database, TrichOME, to facilitate the study of plant trichomes. The database hosts a large volume of functional omics data, including expressed sequence tag/unigene sequences, microarray hybridizations from both trichome and control tissues, mass spectrometry-based trichome metabolite profiles, and trichome-related genes curated from published literature. The expressed sequence tag/unigene sequences have been annotated based upon sequence similarity with popular databases (e.g. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and Transporter Classification Database). The unigenes, metabolites, curated genes, and probe sets have been mapped against each other to enable comparative analysis. The database also integrates bioinformatics tools with a focus on the mining of trichome-specific genes in unigenes and microarray-based gene expression profiles. TrichOME is a valuable and unique resource for plant trichome research, since the genes and metabolites expressed in trichomes are often underrepresented in regular non-tissue-targeted cDNA libraries. TrichOME is freely available at http://www.planttrichome.org/.Plant trichomes are epidermal tissues located on the surfaces of leaves, petals, stems, petioles, peduncles, and seed coats depending on species. By virtue of their physical properties (size, density), trichome hairs can directly serve to protect buds of plants from insect damage, reduce leaf temperature, increase light reflectance, prevent loss of water, and reduce leaf abrasion (Wagner, 1991; Wagner et al., 2004).Although the morphology of trichomes varies greatly, they can be generally classified into two types: simple trichomes (STs) and glandular secreting trichomes (GSTs; Wagner et al., 2004). STs of Arabidopsis (Arabidopsis thaliana) have been chosen as models for studying cell fate and differentiation (Wagner, 1991; Breuer et al., 2009; Marks et al., 2009). In Arabidopsis, STs on leaves consist of a unicellular structure with a stalk and three to four branches (Fig. 1B). Although the STs are referred to as “nonglandular” (presumably nonsecreting), expression of genes involved in anthocyanin, flavonoid, and glucosinolate pathways can nevertheless be detected in STs, indicating the roles of STs in the biosynthesis of secondary compounds and defense (Wang et al., 2002; Jakoby et al., 2008). GSTs are found on about one-third of vascular plants. GSTs have a multicellular structure with a stalk terminating in a glandular head (Fig. 1, A and C–G). GSTs are initiated from a single protodermal cell that undergoes vertical enlargement and multiple divisions to give rise to fully developed trichomes. GSTs often produce and accumulate terpenoid and phenylpropanoid oils (Wagner et al., 2004). However, alkaloids, the third major class of plant secondary compounds, are not common in GST exudates (Laue et al., 2000). The amount of exudates produced by GSTs may reach 30% of mature leaf dry weight, as found in certain Australian desert plants (Dell and McComb, 1978). Plant GSTs can impact pathogen defense, pest resistance, pollinator attraction, and water retention based on the phytochemicals they secrete.Open in a separate windowFigure 1.Representative scanning electron microscopy images of trichomes on plants. A, Erect glandular trichome on the stem of M. sativa. B, Nonglandular trichome on a rosette leaf of Arabidopsis. C, Procumbent trichome on the petiole of M. truncatula. D, Field of glandular trichomes on a female bract of Cannabis sativa. E, Glandular trichomes on a bract of H. lupulus. F, Nonglandular trichome on a leaf of M. truncatula. G, Types VI (small arrow) and I (large arrow) trichomes on a leaf of S. lycopersicum. All the scanning electron microscopy images were generated as described previously (Ahlstrand, 1996; Esch et al., 2004) using an Emitech Technologies (www.emitech.co.uk) K1150 cyropreparation system and a Hitachi High Technologies (www.hitachi-hhta.com) S3500N scanning electron microscope. Bars = 100 μm.GSTs on the aerial organ surfaces have a unique capacity for synthesis and secretion of chemicals (largely plant secondary metabolites), and they have been described as “chemical factories” for the production of high-value natural products (Mahmoud and Croteau, 2002; Wagner et al., 2004; Schilmiller et al., 2008). Secondary metabolites play important roles in protecting the plant against insect predation and other biotic challenges (Peter and Shanower, 1998), and they are potential sources for pharmaceutical and nutraceutical product development. For example, the trichome-borne artemisinin from Artemisia annua is still the most effective drug against malaria, and the early steps of its biosynthetic pathway have been extensively studied (Duke et al., 1994; Arsenault et al., 2008). Recently, the mechanisms by which plant glandular trichomes make, transport, store, and secrete a great variety of unique compounds, especially terpenoids and flavoniods, have received extended research interest because of the potential use of these compounds in pharmaceutical and nutraceutical applications. Seminal studies have reported the assignment of gene functions to specific metabolic pathways in glandular trichomes of several plant species, including mint (Mentha × piperita; Alonso et al., 1992; Rajaonarivony et al., 1992; Lange et al., 2000), basil (Ocimum basilicum; Gang et al., 2002; Iijima et al., 2004; Xie et al., 2008), Artemisia (Teoh et al., 2006; Zhang et al., 2008), tomato (Solanum lycopersicum; Fridman et al., 2005; Besser et al., 2009; Schilmiller et al., 2009), and hop (Humulus lupulus; Nagel et al., 2008; Wang et al., 2008).

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 (No. of sites analyzedWithin A1

---

Within A2

---

Fixed differences between A1 and A2A1/A2 divergence (%)Loci^aS^b totalSπ (%)^cSπ (%)^cA1/A2 divergence <1%A1-23645630020.4410.44A1-0456544000040.61A2-568413220.4820.4800.24A2-411480210.210010.31A1/A2 divergence >1%A1-2176679000091.35A1-12856990010.1881.49A1-199618130010.16122.02A2-4223449000092.62A2-516470140000142.98A2-404508200030.59173.64A2-4355062220.3920.39183.95A2-4734572310.2210.22214.81A2-4573031710.3300165.54A2-5755034750.9930.59398.55Open in a separate window^aA1, loci originated from the A1 sex chromosome; A2, loci originated from the A2 sex chromosome.^bS, number of polymorphic sites.^cπ (%), average number of differences per 100 nucleotides.

TABLE 3

P-values for the 2 × 2 G-tests for significance of differences in A1/A2 divergence between the loci in the nonrecombining region

Variation in Genomic Recombination Rates Among Heterogeneous Stock Mice

We used a large panel of pedigreed, genetically admixed house mice to study patterns of recombination rate variation in a leading mammalian model system. We found considerable inter-individual differences in genomic recombination rates and documented a significant heritable component to this variation. These findings point to clear variation in recombination rate among common laboratory strains, a result that carries important implications for genetic analysis in the house mouse.THE rate of recombination—the amount of crossing over per unit DNA—is a key parameter governing the fidelity of meiosis. Recombination rates that are too high or too low frequently give rise to aneuploid gametes or prematurely arrest the meiotic cell cycle (Hassold and Hunt 2001). As a consequence, recombination rates should experience strong selective pressures to lie within the range defined by the demands of meiosis (Coop and Przeworski 2007). Nonetheless, classical genetic studies in Drosophila (Chinnici 1971; Kidwell 1972; Brooks and Marks 1986), crickets (Shaw 1972), flour beetles (Dewees 1975), and lima beans (Allard 1963) have shown that considerable inter-individual variation for recombination rate is present within populations. Recent studies examining the transmission of haplotypes in human pedigrees have corroborated these findings (Broman et al. 1998; Kong et al. 2002; Coop et al. 2008).Here, we use a large panel of heterogeneous stock (HS) mice to study variation in genomic recombination rates in a genetic model system. These mice are genetically admixed, derived from an initial generation of pseudorandom mating among eight common inbred laboratory strains (DBA/2J, C3H/HeJ, AKR/J, A/J, BALB/cJ, CBA/J, C57BL/6J, and LP/J), followed by >50 generations of pseudorandom mating in subsequent hybrid cohorts (Mott et al. 2000; Demarest et al. 2001). The familial relationships among animals in recent generations were tracked to organize the mice into pedigrees. In total, this HS panel includes ∼2300 animals comprising 85 families, 8 of which span multiple generations. The remainder consists of nuclear families (sibships) that range from 1 to 34 sibs, with an average of 9.6 sibs (Valdar et al. 2006) (Mott et al. 2000; Demarest et al. 2001; Shifman et al. 2006).

TABLE 1

Heterogeneous stock mouse pedigrees

Immunomodulation by Mesenchymal Stem Cells in Veterinary Species

Mesenchymal stem cells (MSC) are adult-derived multipotent stem cells that have been derived from almost every tissue. They are classically defined as spindle-shaped, plastic-adherent cells capable of adipogenic, chondrogenic, and osteogenic differentiation. This capacity for trilineage differentiation has been the foundation for research into the use of MSC to regenerate damaged tissues. Recent studies have shown that MSC interact with cells of the immune system and modulate their function. Although many of the details underlying the mechanisms by which MSC modulate the immune system have been defined for human and rodent (mouse and rat) MSC, much less is known about MSC from other veterinary species. This knowledge gap is particularly important because the clinical use of MSC in veterinary medicine is increasing and far exceeds the use of MSC in human medicine. It is crucial to determine how MSC modulate the immune system for each animal species as well as for MSC derived from any given tissue source. A comparative approach provides a unique translational opportunity to bring novel cell-based therapies to the veterinary market as well as enhance the utility of animal models for human disorders. The current review covers what is currently known about MSC and their immunomodulatory functions in veterinary species, excluding laboratory rodents.Abbreviations: AT, adipose tissue; BM, Bone marrow; CB, umbilical cord blood; CT, umbilical cord tissue; DC, dendritic cell; IDO, indoleamine 2;3-dioxygenase; MSC, mesenchymal stem cells; PGE₂, prostaglandin E2; VEGF, vascular endothelial growth factorMesenchymal stem cells (MSC, alternatively known as mesenchymal stromal cells) were first reported in the literature in 1968.³⁹ MSC are thought to be of pericyte origin (cells that line the vasculature)^21,22 and typically are isolated from highly vascular tissues. In humans and mice, MSC have been isolated from fat, placental tissues (placenta, Wharton jelly, umbilical cord, umbilical cord blood), hair follicles, tendon, synovial membrane, periodontal ligament, and every major organ (brain, spleen, liver, kidney, lung, bone marrow, muscle, thymus, pancreas, skin).^23,121 For most current clinical applications, MSC are isolated from adipose tissue (AT), bone marrow (BM), umbilical cord blood (CB), and umbilical cord tissue (CT; 11,87,99 Clinical trials in human medicine focus on the use of MSC both for their antiinflammatory properties (graft-versus-host disease, irritable bowel syndrome) and their ability to aid in tissue and bone regeneration in combination with growth factors and bone scaffolds (clinicaltrials.gov).¹³¹ For tissue regeneration, the abilities of MSC to differentiate and to secrete mediators and interact with cells of the immune system likely contribute to tissue healing (Figure 1). The current review will not address the specific use of MSC for orthopedic applications and tissue regeneration, although the topic is covered widely in current literature for both human and veterinary medicine.^57,62,90

Table 1.

Tissues from which MSC have been isolated

Complexation and Toxicity of Copper in Higher Plants. II. Different Mechanisms for Copper versus Cadmium Detoxification in the Copper-Sensitive Cadmium/Zinc Hyperaccumulator Thlaspi caerulescens (Ganges Ecotype)

The cadmium/zinc hyperaccumulator Thlaspi caerulescens is sensitive toward copper (Cu) toxicity, which is a problem for phytoremediation of soils with mixed contamination. Cu levels in T. caerulescens grown with 10 μm Cu²⁺ remained in the nonaccumulator range (<50 ppm), and most individuals were as sensitive toward Cu as the related nonaccumulator Thlaspi fendleri. Obviously, hyperaccumulation and metal resistance are highly metal specific. Cu-induced inhibition of photosynthesis followed the “sun reaction” type of damage, with inhibition of the photosystem II reaction center charge separation and the water-splitting complex. A few individuals of T. caerulescens were more Cu resistant. Compared with Cu-sensitive individuals, they recovered faster from inhibition, at least partially by enhanced repair of chlorophyll-protein complexes but not by exclusion, since the content of Cu in their shoots was increased by about 25%. Extended x-ray absorption fine structure (EXAFS) measurements on frozen-hydrated leaf samples revealed that a large proportion of Cu in T. caerulescens is bound by sulfur ligands. This is in contrast to the known binding environment of cadmium and zinc in the same species, which is dominated by oxygen ligands. Clearly, hyperaccumulators detoxify hyperaccumulated metals differently compared with nonaccumulated metals. Furthermore, strong features in the Cu-EXAFS spectra ascribed to metal-metal contributions were found, in particular in the Cu-resistant specimens. Some of these features may be due to Cu binding to metallothioneins, but a larger proportion seems to result from biomineralization, most likely Cu(II) oxalate and Cu(II) oxides. Additional contributions in the EXAFS spectra indicate complexation of Cu(II) by the nonproteogenic amino acid nicotianamine, which has a very high affinity for Cu(II) as further characterized here.Many heavy metals are well known to be essential microelements for plants, but elevated concentrations of these metals cause toxicity (for review, see Prasad and Hagemeyer, 1999; Küpper and Kroneck, 2005), as explained in more detail in our companion article (Küpper et al., 2009) on copper (Cu) complexation and toxicity in Crassula helmsii. Plants developed a number of strategies to resist this toxicity, including active efflux, sequestration, and binding of heavy metals inside the cells by strong ligands. Among the zinc (Zn) and cadmium (Cd) hyperaccumulators (Brooks, 1998; Lombi et al., 2000), the best known species is Thlaspi caerulescens, which has been proposed as a hyperaccumulator model species by several authors (Assunção et al., 2003; Peer et al., 2003, 2006). An enhanced uptake of metals into the root symplasm was found in T. caerulescens compared with the related nonaccumulator Thlaspi arvense (Lasat et al., 1996, 1998), and a reduced sequestration into the root vacuoles was associated with the higher root-to-shoot translocation efficiency of T. caerulescens (Shen et al., 1997; Lasat et al., 1998). This is likely related to elevated expression of xylem-loading transporters in the roots (Papoyan and Kochian, 2004; Weber et al., 2004). One of these, the Cd/Zn-pumping P_1b-type ATPase TcHMA4, was recently purified as a protein, which revealed posttranslational processing, and its biochemical characterization showed Cd and Zn transport affinity in the submicromolar range (Parameswaran et al., 2007). Studies of cellular metal compartmentation have shown that in most hyperaccumulators, the metal is sequestered preferentially into compartments where it does no harm to the metabolism. These are, in most cases studied so far, the epidermal vacuoles (Küpper et al., 1999, 2001; Frey et al., 2000; Bidwell et al., 2004; Broadhurst et al., 2004), where concentrations of several hundred mmol L⁻¹ can be reached in the large metal storage cells (Küpper et al., 1999, 2001). The latter showed that hyperaccumulation must be mediated by active pumping of the heavy metals into their storage sites, which was shown to be achieved by an extremely increased expression of metal transport proteins in leaves of hyperaccumulators compared with nonaccumulators (Pence et al., 2000; Assunção et al., 2001; Becher et al., 2004; Papoyan and Kochian, 2004; Küpper et al., 2007b). Strong sulfur (S) ligands like phytochelatins were shown not to be relevant for Cd detoxification in the Cd hyperaccumulator T. caerulescens. Phytochelatin levels are lower in this plant than in the related nonaccumulator T. arvense (Ebbs et al., 2002), inhibition of phytochelatin synthase in hyperaccumulators does not affect their Cd resistance (Schat et al., 2002), and direct measurements of the Cd ligands by EXAFS (Küpper et al., 2004). Thus, the main detoxification strategy in hyperaccumulators is clearly not binding to strong ligands but sequestration of the hyperaccumulated heavy metals. However, the nonproteogenic amino acid nicotianamine (NA) seems to play an important role in metal homeostasis of plants. According to several studies, it binds iron, Zn, and Cu, mainly for long-distance transport in the vascular bundle (Stephan and Scholz, 1993; Pich et al., 1994; Schmidke and Stephan, 1995; Stephan et al., 1996, Pich and Scholz, 1996; von Wiren et al., 1999; Liao et al., 2000), and NA synthase has been shown to be highly overexpressed in hyperaccumulators compared with nonaccumulator plants (Becher et al., 2004; Weber et al., 2004; van de Mortel et al., 2006, 2008).

Table I.

Explanation of technical terms

A Targeted Multilocus Genotyping Assay for Lineage,Serogroup, and Epidemic Clone Typing of Listeria monocytogenes

A 30-probe assay was developed for simultaneous classification of Listeria monocytogenes isolates by lineage (I to IV), major serogroup (4b, 1/2b, 1/2a, and 1/2c), and epidemic clone (EC) type (ECI, ECIa, ECII, and ECIII). The assay was designed to facilitate rapid strain characterization and the integration of subtype data into risk-based inspection programs.Listeria monocytogenes is a facultative intracellular pathogen that can cause serious invasive illness (listeriosis) in humans and other animals. L. monocytogenes is responsible for over 25% of food-borne-disease-related deaths attributable to known pathogens and is a leading cause of food recalls due to microbial adulteration (12, 21). However, not all L. monocytogenes subtypes contribute equally to human illness, and substantial differences in the ecologies and virulence attributes of different L. monocytogenes subtypes have been identified (9, 13, 14, 23, 24, 33, 35, 36). Among the four major evolutionary lineages of L. monocytogenes, only lineages I and II are commonly isolated from contaminated food and human listeriosis patients (19, 27, 29, 33). Lineage I strains are overrepresented among human listeriosis isolates, particularly those associated with epidemic outbreaks, whereas lineage II strains are overrepresented in foods and the environment (13, 14, 24). Lineage III strains account for approximately 1% of human listeriosis cases but are common among animal listeriosis isolates and appear to be a host-adapted group that is poorly adapted to food-processing environments (6, 34-36). The ecological and virulence attributes of lineage IV are poorly understood, as this lineage is rare and was only recently described based on a small number of strains (19, 26, 29, 33).L. monocytogenes is differentiated into 13 serotypes; however, four major serogroups (4b, 1/2b, 1/2a, and 1/2c) from within lineages I and II account for more than 98% of human and food isolates (16, 31). Serogroups refer to evolutionary complexes typified by a predominant serotype but which include very rare serotypes that represent minor evolutionary variants (7, 9, 33). Phylogenetic analyses have indicated that rare serotypes may have evolved recently, or even multiple times, from one of the major serotypes (9), and numerous molecular methods fail to discriminate minor serotypes as independent groups (1, 4, 7, 9, 18, 22, 33, 38, 39). Serotyping is one of the most common methods for L. monocytogenes subtyping, and serogroup classifications are a useful component of strain characterization because ecotype divisions appear largely congruent with serogroup distinctions (16, 34). Serogroup 4b strains are of particular public health concern because contamination with these strains appears to increase the probability that a ready-to-eat (RTE) food will be implicated in listeriosis (16, 28). Serogroup 4b strains account for approximately 40% of sporadic listeriosis and also are responsible for the majority of listeriosis outbreaks despite being relatively rare contaminants of food products (9, 13, 17, 30, 34). In addition, serogroup 4b strains are associated with more severe clinical presentations and higher mortality rates than other serogroups (11, 16, 20, 31, 34). Serogroups 1/2a and 1/2b are overrepresented among food isolates but also contribute significantly to human listeriosis, whereas serogroup 1/2c rarely causes human illness and may pose a lower risk of listeriosis for humans (16). Serogroup-specific differences in association with human listeriosis are consistent with the prevalence of virulence-attenuating mutations in inlA within these serogroups (32, 34); however, a number of additional factors likely contribute to these differences.Four previously described epidemic clones (ECs; ECI, ECIa, ECII, and ECIII) of L. monocytogenes have been implicated in numerous listeriosis outbreaks and have contributed significantly to sporadic illness (15, 34). ECI, ECIa, and ECII are distinct groups within serogroup 4b that were each responsible for repeated outbreaks of listeriosis in the United States and Europe. ECIII is a lineage II clone of serotype 1/2a that persisted in the same processing facility for more than a decade prior to causing a multistate outbreak linked to contaminated turkey (15, 25). While there has been speculation that epidemic clones possess unique adaptations that explain their frequent involvement in listeriosis outbreaks (9, 34, 37), it is not clear that epidemic clones are more virulent than other strains with the same serotype. However, contamination of RTE food with EC strains would be cause for increased concern due to the previous involvement of these clones in major outbreaks of listeriosis (16).As a result of the L. monocytogenes subtype-specific differences in ecology, virulence, and association with human illness, molecular subtyping technologies have the potential to inform assessments of relative risk and to improve risk-based inspection programs. The objective of the present study was to develop a single assay for rapid and accurate classification of L. monocytogenes isolates by lineage, major serogroup, and epidemic clone in order to facilitate strain characterization and the integration of subtype data into inspection programs that are based on assessment of relative risk.A database of more than 5.3 Mb of comparative DNA sequences from 238 L. monocytogenes isolates (9, 33-35) was scanned for single nucleotide polymorphisms that could be used to differentiate lineages, major serogroups, and epidemic clones via a targeted multilocus genotyping (TMLGT) approach. The acronym TMLGT is used to distinguish this approach from previously published multilocus genotyping (MLGT) assays that were lineage specific and designed for haplotype discrimination (9, 33). To provide for simultaneous interrogation of the selected polymorphisms via TMLGT, six genomic regions (Table ​(Table1)1) were coamplified in a multiplex PCR. While the previous MLGT assays were based on three lineage-specific multiplexes and required prior identification of lineage identity, TMLGT was designed to target variation across all of the lineages simultaneously and is based on a unique set of amplicons. PCR was performed in 50-μl volumes with 1× High Fidelity PCR buffer (Invitrogen Life Technologies), 2 mM MgSO₄, 100 μM deoxynucleoside triphosphate (dNTP), 300 nM primer, 1.5 U Platinum Taq high-fidelity DNA polymerase (Invitrogen Life Technologies), and 100 ng of genomic DNA. PCR consisted of an initial denaturation of 90 s at 96°C, followed by 40 cycles of 30 s at 94°C, 30 s at 50°C, and 90 s at 68°C. Amplification products were purified using Montage PCR cleanup filter plates (Millipore) and served as a template for allele-specific primer extension (ASPE) reactions utilizing subtype-specific probes.

TABLE 1.

Primers used in multiplex amplification for the TMLGT assay