Differences in Mitochondrial Genome Organization of Cryptococcus Neoformans Strains

The organization of the mitochondrial genomes in two strains belonging in different varieties of Cryptococcus neoformans was analysed. Physical maps of the mtDNA of the IFM5844 (var. neoformans) and IFO410 (var. grubii) strains were constructed by using EcoRI and EcoRV restriction enzymes; functional maps were constructed by hybridization, cloning and sequencing. Most of the genes important in the mitochondrial function (ND1, ND2, ND3, ND4, ND4L, ND5, ND6, ATP6, ATP9, COX1, COX2 and COB) and protein synthesis (SsrRNA and LsrRNA) were localized. We did not find any differences between the strains in the order of these genes. However, they differed significantly in the sizes of the mtDNAs: 32.6 kb for IFM5844, and 24.1 kb for IFO410. This can be attributed to two large regions of the mtDNA. In these regions, differences were found in the numbers of introns in COX1 (no intron in var. grubii, 5 introns in var. neoformans), COB (1 intron in var. grubii, 2 introns in var. neoformans), LsrRNA (no intron in var. grubii, 2 introns in var. neoformans), and ND5 (no intron in var. grubii, 1 intron in var. neoformans) genes. In several introns of the COB and COX1 genes LAGLIDADG motifs were found. Differences were also observed in the nucleotide sequences of some genes and in the sizes and sequences of intergenic regions. The nucleotide sequences of the genes of the IFM and IFO strains were compared with those of the H-99 and JEC 21 strains from the database. Surprisingly high similarities were found between the strains belonging in var. grubii (IFO 410 and H-99) and var. neoformans (IFM 5844 and JEC 21).     

Chromosomal Rearrangements between Serotype A and D Strains in Cryptococcus neoformans

Cryptococcus neoformans is a major human pathogenic fungus that can cause meningoencephalitis in immunocompromised hosts. It contains two divergent varieties, var. grubii (serotype A) and var. neoformans (serotype D), as well as hybrids (serotype AD) between these two varieties. In this study, we investigated the extent of chromosomal rearrangements between the two varieties, estimated the effects of chromosomal rearrangements on recombination frequencies, and surveyed the potential polymorphisms of the rearrangements among natural strains of the three serotypes. Through the analyses of two sequenced genomes from strains H99 (representing var. grubii) and JEC21 (representing var. neoformans), we revealed a total of 32 unambiguous chromosome rearrangements, including five translocations, nine simple inversions, and 18 complex rearrangements. Our analyses identified that overall, rearranged regions had recombination frequencies about half of those around syntenic regions. Using a direct PCR screening strategy, we examined the potential polymorphisms of 11 rearrangements among 64 natural C. neoformans strains from five countries. We found no polymorphism within var. neoformans and very limited polymorphism within var. grubii. However, strains of serotype AD showed significant polymorphism, consistent with their hybrid origins coupled with differential loss of heterozygosity. We discuss the implications of these results on the genome structure, ecology, and evolution of C. neoformans.     

Test of intron predictions reveals novel splice sites, alternatively spliced mRNAs and new introns in meiotically regulated genes of yeast

Correct identification of all introns is necessary to discern the protein-coding potential of a eukaryotic genome. The existence of most of the spliceosomal introns predicted in the genome of Saccharomyces cerevisiae remains unsupported by molecular evidence. We tested the intron predictions for 87 introns predicted to be present in non-ribosomal protein genes, more than a third of all known or suspected introns in the yeast genome. Evidence supporting 61 of these predictions was obtained, 20 predicted intron sequences were not spliced and six predictions identified an intron-containing region but failed to specify the correct splice sites, yielding a successful prediction rate of <80%. Alternative splicing has not been previously described for this organism, and we identified two genes (YKL186C/MTR2 and YML034W) which encode alternatively spliced mRNAs; YKL186C/MTR2 produces at least five different spliced mRNAs. One gene (YGR225W/SPO70) has an intron whose removal is activated during meiosis under control of the MER1 gene. We found eight new introns, suggesting that numerous introns still remain to be discovered. The results show that correct prediction of introns remains a significant barrier to understanding the structure, function and coding capacity of eukaryotic genomes, even in a supposedly simple system like yeast.     

The plastid aldolase gene fromChlamydomonas reinhardtii: Intron/exon organization,evolution, and promoter structure

Genomic clones encoding the plastidic fructose- 1,6-bisphosphate aldolase ofChlamydomonas reinhardtii were isolated and sequenced. The gene contains three introns which are located within the coding sequence for the mature protein. No introns are located within or near the sequence encoding the transit-peptide, in contrast to the genes for plastidic aldolases of higher plants. Neither the number nor the positions of the three introns of theC. reinhardtii aldolase gene are conserved in the plastidic or cytosolic aldolase genes of higher plants and animals. The 5′ border sequences of introns in the aldolase gene ofC. reinhardtii exhibit the conserved plant consensus sequence. The 3′ acceptor splice sites for introns 1 and 3 show much less similarity to the eukaryotic consensus sequences than do those of intron 2. The plastidic aldolase gene has two tandemly repeated CAAT box motifs in the promoter region. Genomic Southern blots indicate that the gene is encoded by a single locus in theC. reinhardtii genome.     

The complete nucleotide sequence and RNA editing content of the mitochondrial genome of rapeseed (Brassica napus L.): comparative analysis of the mitochondrial genomes of rapeseed and Arabidopsis thaliana

The entire mitochondrial genome of rapeseed (Brassica napus L.) was sequenced and compared with that of Arabidopsis thaliana. The 221 853 bp genome contains 34 protein-coding genes, three rRNA genes and 17 tRNA genes. This gene content is almost identical to that of Arabidopsis. However the rps14 gene, which is a pseudo-gene in Arabidopsis, is intact in rapeseed. On the other hand, five tRNA genes are missing in rapeseed compared to Arabidopsis, although the set of mitochondrially encoded tRNA species is identical in the two Cruciferae. RNA editing events were systematically investigated on the basis of the sequence of the rapeseed mitochondrial genome. A total of 427 C to U conversions were identified in ORFs, which is nearly identical to the number in Arabidopsis (441 sites). The gene sequences and intron structures are mostly conserved (more than 99% similarity for protein-coding regions); however, only 358 editing sites (83% of total editings) are shared by rapeseed and Arabidopsis. Non-coding regions are mostly divergent between the two plants. One-third (about 78.7 kb) and two-thirds (about 223.8 kb) of the rapeseed and Arabidopsis mitochondrial genomes, respectively, cannot be aligned with each other and most of these regions do not show any homology to sequences registered in the DNA databases. The results of the comparative analysis between the rapeseed and Arabidopsis mitochondrial genomes suggest that higher plant mitochondria are extremely conservative with respect to coding sequences and somewhat conservative with respect to RNA editing, but that non-coding parts of plant mitochondrial DNA are extraordinarily dynamic with respect to structural changes, sequence acquisition and/or sequence loss.     

Introns of Entamoeba histolytica and Entamoeba dispar

The genome of Entamoeba histolytica is considered to possess very few intervening sequences (introns), as only 5 intron-containing genes from this protozoan parasite have been reported so far. However, while sequencing a number of genomic contigs as well as three independent genes coding for ribosomal protein L27a, we have identified 9 additional intron-containing genes of E. histolytica and the closely related species Entamoeba dispar,indicating that introns are more common in these organisms than previously suggested. The various amoeba introns are relatively short comprising between 46 and 115 nucleotides only and have a higher AT-content compared to thecorresponding exon sequences. In contrast to higher eukaryotes, amoeba introns do not contain a well-conserved branch point consensus, and have extended donor and acceptor splice sites of the sequences GTTTGTT and TAG, respectively. Consistent with the close phylogenetic relationship of E. histolytica and E. dispar, the position and length of introns is conserved between the two species but the degree of sequence identity is reduced compared to orthologous coding regions.     

The six genes of the Rubisco small subunit multigene family from Mesembryanthemum crystallinum,a facultative CAM plant

The nucleotide sequences of the entire gene family, comprising six genes, that encodes the Rubisco small subunit (rbcS) multigene family in Mesembryanthemum crystallinum (common ice plant), were determined. Five of the genes are arranged in a tandem array spanning 20 kb, while the sixth gene is not closely linked to this array. The mature small subunit coding regions are highly conserved and encode four distinct polypeptides of equal lengths with up to five amino acid differences distinguishing individual genes. The transit peptide coding regions are more divergent in both amino acid sequence and length, encoding five distinct peptide sequences that range from 55 to 61 amino acids in length. Each of the genes has two introns located at conserved sites within the mature peptide-coding regions. The first introns are diverse in sequence and length ranging from 122 by to 1092 bp. Five of the six second introns are highly conserved in sequence and length. Two genes, rbcS-4 and rbcS-5, are identical at the nucleotide level starting from 121 by upstream of the ATG initiation codon to 9 by downstream of the stop codon including the sequences of both introns, indicating recent gene duplication and/or gene conversion. Functionally important regulatory elements identified in rbcS promoters of other species are absent from the upstream regions of all but one of the ice plant rbcS genes. Relative expression levels were determined for the rbcS genes and indicate that they are differentially expressed in leaves.     

Understanding the Differences between Genome Sequences of Escherichia coli B Strains REL606 and BL21(DE3) and Comparison of the E. coli B and K-12 Genomes

Each difference between the genome sequences of Escherichia coli B strains REL606 and BL21(DE3) can be interpreted in light of known laboratory manipulations plus a gene conversion between ribosomal RNA operons. Two treatments with 1-methyl-3-nitro-1-nitrosoguanidine in the REL606 lineage produced at least 93 single-base-pair mutations (∼ 90% GC-to-AT transitions) and 3 single-base-pair GC deletions. Two UV treatments in the BL21(DE3) lineage produced only 4 single-base-pair mutations but 16 large deletions. P1 transductions from K-12 into the two B lineages produced 317 single-base-pair differences and 9 insertions or deletions, reflecting differences between B DNA in BL21(DE3) and integrated restriction fragments of K-12 DNA inherited by REL606. Two sites showed selective enrichment of spontaneous mutations. No unselected spontaneous single-base-pair mutations were evident. The genome sequences revealed that a progenitor of REL606 had been misidentified, explaining initially perplexing differences. Limited sequencing of other B strains defined characteristic properties of B and allowed assembly of the inferred genome of the ancestral B of Delbrück and Luria. Comparison of the B and K-12 genomes shows that more than half of the 3793 proteins of their basic genomes are predicted to be identical, although ∼ 310 appear to be functional in either B or K-12 but not in both. The ancestral basic genome appears to have had ∼ 4039 coding sequences occupying ∼ 4.0 Mbp. Repeated horizontal transfer from diverged Escherichia coli genomes and homologous recombination may explain the observed variable distribution of single-base-pair differences. Fifteen sites are occupied by phage-related elements, but only six by comparable elements at the same site. More than 50 sites are occupied by IS elements in both B and K, 16 in common, and likely founding IS elements are identified. A signature of widespread cryptic phage P4-type mobile elements was identified. Complex deletions (dense clusters of small deletions and substitutions) apparently removed nonessential genes from ∼ 30 sites in the basic genomes.     

Proteogenomic analysis of pathogenic yeast Cryptococcus neoformans using high resolution mass spectrometry

Background

Cryptococcus neoformans, a basidiomycetous fungus of universal occurrence, is a significant opportunistic human pathogen causing meningitis. Owing to an increase in the number of immunosuppressed individuals along with emergence of drug-resistant strains, C. neoformans is gaining importance as a pathogen. Although, whole genome sequencing of three varieties of C. neoformans has been completed recently, no global proteomic studies have yet been reported.

Results

We performed a comprehensive proteomic analysis of C. neoformans var. grubii (Serotype A), which is the most virulent variety, in order to provide protein-level evidence for computationally predicted gene models and to refine the existing annotations. We confirmed the protein-coding potential of 3,674 genes from a total of 6,980 predicted protein-coding genes. We also identified 4 novel genes and corrected 104 predicted gene models. In addition, our studies led to the correction of translational start site, splice junctions and reading frame used for translation in a number of proteins. Finally, we validated a subset of our novel findings by RT-PCR and sequencing.

Conclusions

Proteogenomic investigation described here facilitated the validation and refinement of computationally derived gene models in the intron-rich genome of C. neoformans, an important fungal pathogen in humans.     

Cyclic AMP-dependent protein kinase catalytic subunits have divergent roles in virulence factor production in two varieties of the fungal pathogen Cryptococcus neoformans

Our earlier findings established that cyclic AMP-dependent protein kinase functions in a signaling cascade that regulates mating and virulence of Cryptococcus neoformans var. grubii (serotype A). Mutants lacking the serotype A protein kinase A (PKA) catalytic subunit Pka1 are unable to mate, fail to produce melanin or capsule, and are avirulent in animal models, whereas mutants lacking the PKA regulatory subunit Pkr1 overproduce capsule and are hypervirulent. Because other mutations have been observed to confer different phenotypes in two diverged varieties of C. neoformans (grubii variety [serotype A] and neoformans variety [serotype D]), we analyzed the functions of the PKA genes in the serotype D neoformans variety. Surprisingly, the Pka1 catalytic subunit was not required for mating, haploid fruiting, or melanin or capsule production of serotype D strains. Here we identify a second PKA catalytic subunit gene, PKA2, that is present in both serotype A and D strains of C. neoformans. The divergent Pka2 catalytic subunit was found to regulate mating, haploid fruiting, and virulence factor production in serotype D strains. In contrast, Pka2 has no role in mating, melanin production, or capsule formation in serotype A strains. Our studies illustrate how different components of signaling pathways can be co-opted and functionally specialized during the evolution of related but distinct varieties or subspecies of a human fungal pathogen.     

Microsatellite Typing of Clinical and Environmental Cryptococcus neoformans var. grubii Isolates from Cuba Shows Multiple Genetic Lineages

Background

Human cryptococcal infections have been associated with bird droppings as a likely source of infection. Studies toward the local and global epidemiology of Cryptococcus spp. have been hampered by the lack of rapid, discriminatory, and exchangeable molecular typing methods.

Methodology/Principal Findings

We selected nine microsatellite markers for high-resolution fingerprinting from the genome of C. neoformans var. grubii. This panel of markers was applied to a collection of clinical (n = 122) and environmental (n = 68; from pigeon guano) C. neoformans var. grubii isolates from Cuba. All markers proved to be polymorphic. The average number of alleles per marker was 9 (range 5–51). A total of 104 genotypes could be distinguished. The discriminatory power of this panel of markers was 0.993. Multiple clusters of related genotypes could be discriminated that differed in only one or two microsatellite markers. These clusters were assigned as microsatellite complexes. The majority of environmental isolates (>70%) fell into 1 microsatellite complex containing only few clinical isolates (49 environmental versus 2 clinical). Clinical isolates were segregated over multiple microsatellite complexes.

Conclusions/Significance

A large genotypic variation exists in C. neoformans var. grubii. The genotypic segregation between clinical and environmental isolates from pigeon guano suggests additional source(s) of human cryptococcal infections. The selected panel of microsatellite markers is an excellent tool to study the epidemiology of C. neoformans var. grubii.     

Genome-wide discovery of DNA polymorphism in Brassica rapa

Single nucleotide polymorphisms (SNPs) and/or insertion/deletions (InDels) are frequent sequence variations in the plant genome, which can be developed as molecular markers for genetic studies on crop improvement. The ongoing Brassica rapa genome sequencing project has generated vast amounts of sequence data useful in genetic research. Here, we report a genome-wide survey of DNA polymorphisms in the B. rapa genome based on the 557 bacterial artificial clone sequences of B. rapa ssp. pekinensis cv. Chiifu. We identified and characterized 21,311 SNPs and 6,753 InDels in the gene space of the B. rapa genome by re-sequencing 1,398 sequence-tagged sites (STSs) in eight genotypes. Comparison of our findings with a B. rapa genetic linkage map confirmed that STS loci were distributed randomly over the B. rapa whole genome. In the 1.4 Mb of aligned sequences, mean nucleotide polymorphism and diversity were θ = 0.00890 and π = 0.00917, respectively. Additionally, the nucleotide diversity in introns was almost three times greater than that in exons, and the frequency of observed InDel was almost 17 times higher in introns than in exons. Information regarding SNPs/InDels obtained here will provide an important resource for genetic studies and breeding programs of B. rapa.     

Choosing a Benchtop Sequencing Machine to Characterise Helicobacter pylori Genomes

The fully annotated genome sequence of the European strain, 26695 was first published in 1997 and, in 1999, it was directly compared to the USA isolate J99, promoting two standard laboratory isolates for Helicobacter pylori (H. pylori) research. With the genomic scaffolds available from these important genomes and the advent of benchtop high-throughput sequencing technology, a bacterial genome can now be sequenced within a few days. We sequenced and analysed strains J99 and 26695 using the benchtop-sequencing machines Ion Torrent PGM and the Illumina MiSeq Nextera and Nextera XT methodologies. Using publically available algorithms, we analysed the raw data and interrogated both genomes by mapping the data and by de novo assembly. We compared the accuracy of the coding sequence assemblies to the originally published sequences. With the Ion Torrent PGM, we found an inherently high-error rate in the raw sequence data. Using the Illumina MiSeq, we found significantly more non-covered nucleotides when using the less expensive Illumina Nextera XT compared with the Illumina Nextera library creation method. We found the most accurate de novo assemblies using the Nextera technology, however, extracting an accurate multi-locus sequence type was inconsistent compared to the Ion Torrent PGM. We found the cagPAI failed to assemble onto a single contig in all technologies but was more accurate using the Nextera. Our results indicate the Illumina MiSeq Nextera method is the most accurate for de novo whole genome sequencing of H. pylori.     

Isolation and Characterization of Cryptococcus neoformans Varieties Recovered from Natural Sources in Bogotá, Colombia, and Study of Ecological Conditions in the Area

Cryptococcus neoformans, the etiological agent of cryptococcosis, has been associated with avian droppings and certain trees in different countries, including Colombia. C neoformans environmental isolates were obtained in urban areas in Bogotá, Colombia, and the strains recovered were phenotypically characterized. Attempts to determine the ecological conditions (micro- and macroclimatic) possibly related to their habitat were also undertaken. Four hundred and eighty samples from bark, soil around trunk bases, and detritus inside hollows of 32 trees were collected in three urban areas during a 5-month period, as well as 89 avian droppings samples from different places. Of plant samples, 6.7% collected from nine tree species yielded C. neoformans var. gattii, serotype B strains in 99% of the cases, and C. neoformans var. grubii, serotype A in 1%. The yeast was more frequently recovered from bark than from soil or detritus inside hollows, and from trees with hollows or rotted wood rather than from trees in which birds nest. C. neoformans was present with higher frequency and density in the rainy season than in the dry season; we found that slightly higher temperature and humidity values of the microhabitat, as compared to those of the environment, favored fungal occurrence, but the phenological state of the tree did not. Of dropping samples, 7.9% yielded C. neoformans strains, all of them C. neoformans var. grubii, serotype A. The yeast was obtained more frequently from dry droppings than from moist ones, but neither the sunlight exposure nor the site of collection of samples was correlated with this occurrence. Population density was significantly higher in droppings than in tree samples. Under laboratory conditions, isolates of different serotype showed similar capsular sizes. Water content and pH ranges were wide and did not show any significant difference between positive and negative samples.     

Uniparental Mitochondrial Transmission in Sexual Crosses in Cryptococcus neoformans

Restriction fragment length polymorphism (RFLP) in the large ribosomal RNA region of the mitochondrial DNA (mtDNA) was developed as a genetic marker for investigating mitochondrial transmission in sexual crosses of the human pathogenic basidiomycetous yeast Cryptococcus neoformans. Strain JEC20 of C. neoformans var. neoformans (mat a) was mated with six strains of C. neoformans var. grubii (mat alpha). Successful mating was indicated by the formation of hyphae and basidiospores. These basidiospores were examined for mtDNA RFLP genotypes. All 570 basidiospores examined from the six crosses showed the mtDNA genotype of strain JEC20. The failure to recover the C. neoformans var. grubii mtDNA in any cross indicates that the C. neoformans var. grubii mtDNA is either selectively eliminated in the newly formed dikaryon or selectively excluded in the immediate dikaryotic hyphae of the newly formed dikaryon. Received: 20 September 1999 / Accepted: 12 November 1999