Genome Sequence and Analysis of Nicotiana benthamiana, the Model Plant for Interactions between Organisms.

Nicotiana benthamiana is widely used as a model plant for dicotyledonous angiosperms. In fact, the strains used in research are highly susceptible to a wide range of viruses. Accordingly, these strains are subject to plant pathology and plant-microbe interactions. In terms of plant-plant interactions, N. benthamiana is one of the plants that exhibit grafting affinity with plants from different families. Thus, N. benthamiana is a good model for plant biology and has been the subject of genome sequencing analyses for many years. However, N. benthamiana has a complex allopolyploid genome, and its previous reference genome is fragmented into 141,000 scaffolds. As a result, molecular genetic analysis is difficult to perform. To improve this effort, de novo whole-genome assembly was performed in N. benthamiana with Hifi reads, and 1,668 contigs were generated with a total length of 3.1 Gb. The 21 longest scaffolds, regarded as pseudomolecules, contained a 2.8-Gb sequence, occupying 95.6% of the assembled genome. A total of 57,583 high-confidence gene sequences were predicted. Based on a comparison of the genome structures between N. benthamiana and N. tabacum, N. benthamiana was found to have more complex chromosomal rearrangements, reflecting the age of interspecific hybridization. To verify the accuracy of the annotations, the cell wall modification genes involved in grafting were analyzed, which revealed not only the previously indeterminate untranslated region, intron and open reading frame sequences but also the genomic locations of their family genes. Owing to improved genome assembly and annotation, N. benthamiana would increasingly be more widely accessible.

A chromosome-level genome sequence of Chrysanthemum seticuspe, a model species for hexaploid cultivated chrysanthemum.

Chrysanthemums are one of the most industrially important cut flowers worldwide. However, their segmental allopolyploidy and self-incompatibility have prevented the application of genetic analysis and modern breeding strategies. We thus developed a model strain, Gojo-0 (Chrysanthemum seticuspe), which is a diploid and self-compatible pure line. Here, we present the 3.05 Gb chromosome-level reference genome sequence, which covered 97% of the C. seticuspe genome. The genome contained more than 80% interspersed repeats, of which retrotransposons accounted for 72%. We identified recent segmental duplication and retrotransposon expansion in C. seticuspe, contributing to arelatively large genome size. Furthermore, we identified a retrotransposon family, SbdRT, which was enriched in gene-dense genome regions and had experienced a very recent transposition burst. We also demonstrated that the chromosome-level genome sequence facilitates positional cloning in C. seticuspe. The genome sequence obtained here can greatly contribute as a reference for chrysanthemum in front-line breeding including genome editing.

Root-knot nematode genetic diversity associated with host compatibility to sweetpotato cultivars.

Plant parasitic root-knot nematodes (RKN) such as Meloidogyne incognita cause significant crop losses worldwide. Although RKN are polyphagous, with wide host ranges, races with differing host compatibilities have evolved. Associations between genotype and infection phenotype in M. incognita have not yet been discovered. In this study, 48 M. incognita isolates were collected from geographically diverse fields in Japan and their genomes sequenced. The isolates exhibited various infection compatibilities to five sweetpotato (SP) cultivars and were assigned to SP races. Genome-wide association analysis identified 743 SNPs affecting gene coding sequences, a large number of which (575) were located on a single 1 Mb region. To examine how this polymorphic region evolved, nucleotide diversity (Pi) was scanned at the whole genome scale. The SNP-rich 1 Mb region exhibited high Pi values and was clearly associated with the SP races. SP1 and 2 races showed high Pi values in this region whereas the Pi values of SP3, 4, and 6 were low. Principal component analysis of isolates from this study and globally collected isolates showed selective divergence in this 1 Mb region. Our results suggest for the first time that the host could be a key determining factor stimulating the genomic divergence of M. incognita.

The persimmon genome reveals clues to the evolution of a lineage-specific sex determination system in plants.

Most angiosperms bear hermaphroditic flowers, but a few species have evolved outcrossing strategies, such as dioecy, the presence of separate male and female individuals. We previously investigated the mechanisms underlying dioecy in diploid persimmon (D. lotus) and found that male flowers are specified by repression of the autosomal gene MeGI by its paralog, the Y-encoded pseudo-gene OGI. This mechanism is thought to be lineage-specific, but its evolutionary path remains unknown. Here, we developed a full draft of the diploid persimmon genome (D. lotus), which revealed a lineage-specific whole-genome duplication event and provided information on the architecture of the Y chromosome. We also identified three paralogs, MeGI, OGI and newly identified Sister of MeGI (SiMeGI). Evolutionary analysis suggested that MeGI underwent adaptive evolution after the whole-genome duplication event. Transformation of tobacco plants with MeGI and SiMeGI revealed that MeGI specifically acquired a new function as a repressor of male organ development, while SiMeGI presumably maintained the original function. Later, a segmental duplication event spawned MeGI's regulator OGI on the Y-chromosome, completing the path leading to dioecy, and probably initiating the formation of the Y-chromosome. These findings exemplify how duplication events can provide flexible genetic material available to help respond to varying environments and provide interesting parallels for our understanding of the mechanisms underlying the transition into dieocy in plants.

Genome-wide association study overcomes the genome complexity in autohexaploid chrysanthemum and tags SNP markers onto the flower color genes.

The use of DNA markers has revolutionized selection in crop breeding by linkage mapping and QTL analysis, but major problems still remain for polyploid species where marker-assisted selection lags behind the situation in diploids because of its high genome complexity. To overcome the complex genetic mode in the polyploids, we investigated the development of a strategy of genome-wide association study (GWAS) using single-dose SNPs, which simplify the segregation patterns associated polyploids, with respect to the development of DNA markers. In addition, we employed biparental populations for the GWAS, wherein the SNP allele frequency could be predicted. The research investigated whether the method could be used to effectively develop DNA markers for petal color in autohexaploid chrysanthemum (Chrysanthemum morifolium; 2n = 6x = 54). The causal gene for this trait is already-known CmCCD4a encoding a dioxygenase which cleaves carotenoids in petals. We selected 9,219 single-dose SNPs, out of total 52,489 SNPs identified by dd-RAD-Seq, showing simplex (1 × 0) and double-simplex (1 × 1) inheritance pattern according to alternative allele frequency with respect to the SNP loci in the F1 population. GWAS, using these single-dose SNPs, discovered highly reproducible SNP markers tightly linked to the causal genes. This is the first report of a straightforward GWAS-based marker developing system for use in autohexaploid species.

Current status in whole genome sequencing and analysis of Ipomoea spp.

The recent advances of next-generation sequencing have made it possible to construct reference genome sequences in divergent species. However, de novo assembly at the chromosome level remains challenging in polyploid species, due to the existence of more than two pairs of homoeologous chromosomes in one nucleus. Cultivated sweet potato (Ipomoea batatas (L.) Lam) is a hexaploid species with 90 chromosomes (2n = 6X = 90). Although the origin of sweet potato is also still under discussion, diploid relative species, I. trifida and I. triloba have been considered as one of the most possible progenitors. In this manuscript, we review the recent results and activities of whole-genome sequencing in the genus Ipomoea series Batatas, I. trifida, I. triloba and sweet potato (I. batatas). Most of the results of genome assembly suggest that the genomes of sweet potato consist of two pairs and four pairs of subgenomes, i.e., B1B1B2B2B2B2. The results also revealed the relation between sweet potato and other Ipomoea species. Together with the development of bioinformatics approaches, the large-scale publicly available genome and transcript sequence resources and international genome sequencing streams are expected to promote the genome sequence dissection in sweet potato.

Establishment of a genome-wide and quantitative protocol for assessment of transcriptional activity at human retrotransposon L1 antisense promoters.

Long interspersed element 1 (L1) retrotransposon sequences are widespread in the human genome, occupying ~500,000 locations. The majority of L1s have lost their retrotransposition capability, although a significant population of human L1s maintains bidirectional transcriptional activity from the internal promoter. While the sense promoter drives transcription of the entire L1 mRNA and leads to L1 retrotransposition, the antisense promoter (ASP) transcribes L1-gene chimeric RNAs that include neighboring exon sequences. Activation mechanisms and functional impacts of L1ASP transcription are thought to vary at every L1ASP location. To explore the locus-specific regulation and function of L1ASP transcription, quantitative methodology is necessary for identifying the genomic positions of highly active L1ASPs on a genome-wide scale. Here, we employed deep-sequencing techniques and built a 3' RACE-based experimental and bioinformatics protocol, named the L1 antisense transcriptome protocol (LATRAP). In LATRAP, the PCR primer and the read mapping scheme were designed to reduce false positives and negatives, which may have been included as hits in previous cloning studies. LATRAP was here applied to the A549 human lung cancer cell line, and 313 L1ASP loci were detected to have transcriptional activity but differed in the number of mapped reads by four orders of magnitude. This indicates that transcriptional activities of the individual L1ASPs can vary greatly and that only a small population of L1ASP loci is active within individual nuclei. LATRAP is the first experimental method for ranking L1ASPs according to their transcriptional activity and will thus open a new avenue to unveiling the locus-specific biology of L1ASPs.

Discrimination of candidate subgenome-specific loci by linkage map construction with an S1 population of octoploid strawberry (Fragaria × ananassa).

BACKGROUND: The strawberry, Fragaria × ananassa, is an allo-octoploid (2n = 8x = 56) and outcrossing species. Although it is the most widely consumed berry crop in the world, its complex genome structure has hindered its genetic and genomic analysis, and thus discrimination of subgenome-specific loci among the homoeologous chromosomes is needed. In the present study, we identified candidate subgenome-specific single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) loci, and constructed a linkage map using an S1 mapping population of the cultivar 'Reikou' with an IStraw90 Axiom® SNP array and previously published SSR markers. RESULTS: The 'Reikou' linkage map consisted of 11,574 loci (11,002 SNPs and 572 SSR loci) spanning 2816.5 cM of 31 linkage groups. The 11,574 loci were located on 4738 unique positions (bin) on the linkage map. Of the mapped loci, 8999 (8588 SNPs and 411 SSR loci) showed a 1:2:1 segregation ratio of AA:AB:BB allele, which suggested the possibility of deriving loci from candidate subgenome-specific sequences. In addition, 2575 loci (2414 SNPs and 161 SSR loci) showed a 3:1 segregation of AB:BB allele, indicating they were derived from homoeologous genomic sequences. Comparative analysis of the homoeologous linkage groups revealed differences in genome structure among the subgenomes. CONCLUSIONS: Our results suggest that candidate subgenome-specific loci are randomly located across the genomes, and that there are small- to large-scale structural variations among the subgenomes. The mapped SNPs and SSR loci on the linkage map are expected to be seed points for the construction of pseudomolecules in the octoploid strawberry.

A high-density SNP genetic map consisting of a complete set of homologous groups in autohexaploid sweetpotato (Ipomoea batatas).

Sweetpotato (Ipomoea batatas) is an autohexaploid species with 90 chromosomes (2n = 6x = 90) and a basic chromosome number of 15, and is therefore regarded as one of the most challenging species for high-density genetic map construction. Here, we used single nucleotide polymorphisms (SNPs) identified by double-digest restriction site-associated DNA sequencing based on next-generation sequencing technology to construct a map for sweetpotato. We then aligned the sequence reads onto the reference genome sequence of I. trifida, a likely diploid ancestor of sweetpotato, to detect SNPs. In addition, to simplify analysis of the complex genetic mode of autohexaploidy, we used an S1 mapping population derived from self-pollination of a single parent. As a result, 28,087 double-simplex SNPs showing a Mendelian segregation ratio in the S1 progeny could be mapped onto 96 linkage groups (LGs), covering a total distance of 33,020.4 cM. Based on the positions of the SNPs on the I. trifida genome, the LGs were classified into 15 groups, each with roughly six LGs and six small extra groups. The molecular genetic techniques used in this study are applicable to high-density mapping of other polyploid plant species, including important crops.

Target Amplicon Sequencing for Genotyping Genome-Wide Single Nucleotide Polymorphisms Identified by Whole-Genome Resequencing in Peanut.

Genome-wide genotyping data regarding breeding materials are essential resources for improving breeding efficiency, especially in plants with complex genomes with a high degree of polyploidy. Several current breeding efforts in cultivated peanut ( L.), which has a tetraploid genome, are devoted to developing high oleic acid cultivars. Genetic maps for such breeding programs have been developed using simple-sequence repeat (SSR) markers, the use of which requires time-consuming electrophoretic analyses. Next-generation sequencing (NGS) technology can overcome this technical hurdle. Initially, we attempted double-digest restriction site-associated DNA sequencing on peanut breeding materials used to develop high oleic acid cultivars. However, this method was not effective because few single nucleotide polymorphism (SNPs) were available because of low genetic diversity of the lines. The genome sequences of the probable diploid ancestors of cultivated peanut, Krapov. & W. C. Greg. and Krapov. & W. C. Greg., are available. Therefore, we next employed whole-genome resequencing analysis to obtain genome-wide SNP data. In this analysis, we observed large biases in the numbers and genomic positions of interspecific and intraspecific SNPs. For genome-wide genotyping, we selected a subset of SNPs covering the peanut genome as the targets of amplicon sequencing analysis. Using this technique, genome-wide genotypes of the breeding materials were easily and rapidly determined. The SNP information and analytic methods developed in this study should accelerate genetics, genomics, and breeding in peanut.

Clustering of CARMA1 through SH3-GUK domain interactions is required for its activation of NF-κB signalling.

CARMA1-mediated NF-κB activation controls lymphocyte activation through antigen receptors and survival of some malignant lymphomas. CARMA1 clusters are formed on physiological receptor-mediated activation or by its oncogenic mutation in activated B-cell-diffuse large B-cell lymphomas (ABC-DLBCLs) with constitutive NF-κB activation. However, regulatory mechanisms and relevance of CARMA1 clusters in the NF-κB pathway are unclear. Here we show that SH3 and GUK domain interactions of CARMA1 link CARMA1 clustering to signal activation. SH3 and GUK domains of CARMA1 interact by either intra- or intermolecular mechanisms, which are required for activation-induced assembly of CARMA1. Disruption of these interactions abolishes the formation of CARMA1 microclusters at the immunological synapse, CARMA-regulated signal activation following antigen receptor stimulation as well as spontaneous CARMA1 clustering and NF-κB activation by the oncogenic CARMA1 mutation in ABC-DLBCLs. Thus, the SH3-GUK interactions that regulate CARMA1 cluster formations are promising therapeutic targets for ABC-DLBCLs.

A reference genetic linkage map of apomictic Hieracium species based on expressed markers derived from developing ovule transcripts.

BACKGROUND AND AIMS: Apomixis in plants generates clonal progeny with a maternal genotype through asexual seed formation. Hieracium subgenus Pilosella (Asteraceae) contains polyploid, highly heterozygous apomictic and sexual species. Within apomictic Hieracium, dominant genetic loci independently regulate the qualitative developmental components of apomixis. In H. praealtum, LOSS OF APOMEIOSIS (LOA) enables formation of embryo sacs without meiosis and LOSS OF PARTHENOGENESIS (LOP) enables fertilization-independent seed formation. A locus required for fertilization-independent endosperm formation (AutE) has been identified in H. piloselloides. Additional quantitative loci appear to influence the penetrance of the qualitative loci, although the controlling genes remain unknown. This study aimed to develop the first genetic linkage maps for sexual and apomictic Hieracium species using simple sequence repeat (SSR) markers derived from expressed transcripts within the developing ovaries. METHODS: RNA from microdissected Hieracium ovule cell types and ovaries was sequenced and SSRs were identified. Two different F1 mapping populations were created to overcome difficulties associated with genome complexity and asexual reproduction. SSR markers were analysed within each mapping population to generate draft linkage maps for apomictic and sexual Hieracium species. KEY RESULTS: A collection of 14 684 Hieracium expressed SSR markers were developed and linkage maps were constructed for Hieracium species using a subset of the SSR markers. Both the LOA and LOP loci were successfully assigned to linkage groups; however, AutE could not be mapped using the current populations. Comparisons with lettuce (Lactuca sativa) revealed partial macrosynteny between the two Asteraceae species. CONCLUSIONS: A collection of SSR markers and draft linkage maps were developed for two apomictic and one sexual Hieracium species. These maps will support cloning of controlling genes at LOA and LOP loci in Hieracium and should also assist with identification of quantitative loci that affect the expressivity of apomixis. Future work will focus on mapping AutE using alternative populations.

Draft genome sequence of eggplant (Solanum melongena L.): the representative solanum species indigenous to the old world.

Unlike other important Solanaceae crops such as tomato, potato, chili pepper, and tobacco, all of which originated in South America and are cultivated worldwide, eggplant (Solanum melongena L.) is indigenous to the Old World and in this respect it is phylogenetically unique. To broaden our knowledge of the genomic nature of solanaceous plants further, we dissected the eggplant genome and built a draft genome dataset with 33,873 scaffolds termed SME_r2.5.1 that covers 833.1 Mb, ca. 74% of the eggplant genome. Approximately 90% of the gene space was estimated to be covered by SME_r2.5.1 and 85,446 genes were predicted in the genome. Clustering analysis of the predicted genes of eggplant along with the genes of three other solanaceous plants as well as Arabidopsis thaliana revealed that, of the 35,000 clusters generated, 4,018 were exclusively composed of eggplant genes that would perhaps confer eggplant-specific traits. Between eggplant and tomato, 16,573 pairs of genes were deduced to be orthologous, and 9,489 eggplant scaffolds could be mapped onto the tomato genome. Furthermore, 56 conserved synteny blocks were identified between the two species. The detailed comparative analysis of the eggplant and tomato genomes will facilitate our understanding of the genomic architecture of solanaceous plants, which will contribute to cultivation and further utilization of these crops.

The rice endosperm ADP-glucose pyrophosphorylase large subunit is essential for optimal catalysis and allosteric regulation of the heterotetrameric enzyme.

Although an alternative pathway has been suggested, the prevailing view is that starch synthesis in cereal endosperm is controlled by the activity of the cytosolic isoform of ADPglucose pyrophosphorylase (AGPase). In rice, the cytosolic AGPase isoform is encoded by the OsAGPS2b and OsAGPL2 genes, which code for the small (S2b) and large (L2) subunits of the heterotetrameric enzyme, respectively. In this study, we isolated several allelic missense and nonsense OsAGPL2 mutants by N-methyl-N-nitrosourea (MNU) treatment of fertilized egg cells and by TILLING (Targeting Induced Local Lesions in Genomes). Interestingly, seeds from three of the missense mutants (two containing T139I and A171V) were severely shriveled and had seed weight and starch content comparable with the shriveled seeds from OsAGPL2 null mutants. Results from kinetic analysis of the purified recombinant enzymes revealed that the catalytic and allosteric regulatory properties of these mutant enzymes were significantly impaired. The missense heterotetramer enzymes and the S2b homotetramer had lower specific (catalytic) activities and affinities for the activator 3-phosphoglycerate (3-PGA). The missense heterotetramer enzymes showed more sensitivity to inhibition by the inhibitor inorganic phosphate (Pi) than the wild-type AGPase, while the S2b homotetramer was profoundly tolerant to Pi inhibition. Thus, our results provide definitive evidence that starch biosynthesis during rice endosperm development is controlled predominantly by the catalytic activity of the cytoplasmic AGPase and its allosteric regulation by the effectors. Moreover, our results show that the L2 subunit is essential for both catalysis and allosteric regulatory properties of the heterotetramer enzyme.

Dissection of the octoploid strawberry genome by deep sequencing of the genomes of Fragaria species.

Cultivated strawberry (Fragaria x ananassa) is octoploid and shows allogamous behaviour. The present study aims at dissecting this octoploid genome through comparison with its wild relatives, F. iinumae, F. nipponica, F. nubicola, and F. orientalis by de novo whole-genome sequencing on an Illumina and Roche 454 platforms. The total length of the assembled Illumina genome sequences obtained was 698 Mb for F. x ananassa, and ∼200 Mb each for the four wild species. Subsequently, a virtual reference genome termed FANhybrid_r1.2 was constructed by integrating the sequences of the four homoeologous subgenomes of F. x ananassa, from which heterozygous regions in the Roche 454 and Illumina genome sequences were eliminated. The total length of FANhybrid_r1.2 thus created was 173.2 Mb with the N50 length of 5137 bp. The Illumina-assembled genome sequences of F. x ananassa and the four wild species were then mapped onto the reference genome, along with the previously published F. vesca genome sequence to establish the subgenomic structure of F. x ananassa. The strategy adopted in this study has turned out to be successful in dissecting the genome of octoploid F. x ananassa and appears promising when applied to the analysis of other polyploid plant species.

The genome of the mesopolyploid crop species Brassica rapa.

We report the annotation and analysis of the draft genome sequence of Brassica rapa accession Chiifu-401-42, a Chinese cabbage. We modeled 41,174 protein coding genes in the B. rapa genome, which has undergone genome triplication. We used Arabidopsis thaliana as an outgroup for investigating the consequences of genome triplication, such as structural and functional evolution. The extent of gene loss (fractionation) among triplicated genome segments varies, with one of the three copies consistently retaining a disproportionately large fraction of the genes expected to have been present in its ancestor. Variation in the number of members of gene families present in the genome may contribute to the remarkable morphological plasticity of Brassica species. The B. rapa genome sequence provides an important resource for studying the evolution of polyploid genomes and underpins the genetic improvement of Brassica oil and vegetable crops.

Sequence analysis of the genome of an oil-bearing tree, Jatropha curcas L.

The whole genome of Jatropha curcas was sequenced, using a combination of the conventional Sanger method and new-generation multiplex sequencing methods. Total length of the non-redundant sequences thus obtained was 285 858 490 bp consisting of 120 586 contigs and 29 831 singlets. They accounted for ~95% of the gene-containing regions with the average G + C content was 34.3%. A total of 40 929 complete and partial structures of protein encoding genes have been deduced. Comparison with genes of other plant species indicated that 1529 (4%) of the putative protein-encoding genes are specific to the Euphorbiaceae family. A high degree of microsynteny was observed with the genome of castor bean and, to a lesser extent, with those of soybean and Arabidopsis thaliana. In parallel with genome sequencing, cDNAs derived from leaf and callus tissues were subjected to pyrosequencing, and a total of 21 225 unigene data have been generated. Polymorphism analysis using microsatellite markers developed from the genomic sequence data obtained was performed with 12 J. curcas lines collected from various parts of the world to estimate their genetic diversity. The genomic sequence and accompanying information presented here are expected to serve as valuable resources for the acceleration of fundamental and applied research with J. curcas, especially in the fields of environment-related research such as biofuel production. Further information on the genomic sequences and DNA markers is available at http://www.kazusa.or.jp/jatropha/.

A protein secretion system linked to bacteroidete gliding motility and pathogenesis.

Porphyromonas gingivalis secretes strong proteases called gingipains that are implicated in periodontal pathogenesis. Protein secretion systems common to other Gram-negative bacteria are lacking in P. gingivalis, but several proteins, including PorT, have been linked to gingipain secretion. Comparative genome analysis and genetic experiments revealed 11 additional proteins involved in gingipain secretion. Six of these (PorK, PorL, PorM, PorN, PorW, and Sov) were similar in sequence to Flavobacterium johnsoniae gliding motility proteins, and two others (PorX and PorY) were putative two-component system regulatory proteins. Real-time RT-PCR analysis revealed that porK, porL, porM, porN, porP, porT, and sov were down-regulated in P. gingivalis porX and porY mutants. Disruption of the F. johnsoniae porT ortholog resulted in defects in motility, chitinase secretion, and translocation of a gliding motility protein, SprB adhesin, to the cell surface, providing a link between a unique protein translocation system and a motility apparatus in members of the Bacteroidetes phylum.

Role of disulfide bonds in goose-type lysozyme.

The role of the two disulfide bonds (Cys4-Cys60 and Cys18-Cys29) in the activity and stability of goose-type (G-type) lysozyme was investigated using ostrich egg-white lysozyme as a model. Each of the two disulfide bonds was deleted separately or simultaneously by substituting both Cys residues with either Ser or Ala. No remarkable differences in secondary structure or catalytic activity were observed between the wild-type and mutant proteins. However, thermal and guanidine hydrochloride unfolding experiments revealed that the stabilities of mutants lacking one or both of the disulfide bonds were significantly decreased relative to those of the wild-type. The destabilization energies of mutant proteins agreed well with those predicted from entropic effects in the denatured state. The effects of deleting each disulfide bond on protein stability were found to be approximately additive, indicating that the individual disulfide bonds contribute to the stability of G-type lysozyme in an independent manner. Under reducing conditions, the thermal stability of the wild-type was decreased to a level nearly equivalent to that of a Cys-free mutant (C4S/C18S/C29S/C60S) in which all Cys residues were replaced by Ser. Moreover, the optimum temperature of the catalytic activity for the Cys-free mutant was downshifted by about 20 degrees C as compared with that of the wild-type. These results indicate that the formation of the two disulfide bonds is not essential for the correct folding into the catalytically active conformation, but is crucial for the structural stability of G-type lysozyme.

Genome-wide analysis of Chlamydophila pneumoniae gene expression at the late stage of infection.

Chlamydophila pneumoniae, an obligate intracellular eubacterium, changes its form from a vegetative reticulate body into an infectious elementary body during the late stage of its infection cycle. Comprehension of the molecular events in the morphological change is important to understand the switching mechanism between acute and chronic infection, which is deemed to relate to the pathogenesis of atherosclerosis. Herein, we have attempted to screen genes expressed in the late stage with a genome-wide DNA microarray, resulting in nomination of 17 genes as the late-stage genes. Fourteen of the 17 genes and six other genes predicted as late-stage genes were confirmed to be up-regulated in the late stage with a quantitative reverse transcriptase-polymerase chain reaction. These 20 late-stage genes were classified into two groups by clustering analysis: 'drastically induced' and 'moderately induced' genes. Out of eight drastically induced genes, four contain sigma(28) promoter-like sequences and the other four contain an upstream common sequence. It suggests that besides sigma(28), there are certain up-regulatory mechanisms at the late stage, which may be involved in the chlamydial morphological change and thus pathogenesis.