Identification of <Emphasis Type="Italic">Ganoderma</Emphasis>, the causal agent of basal stem rot disease in oil palm using a molecular method

From comparison of the alignments of the internally transcribed spacers (ITS) of ribosomal DNA from Ganoderma associated with oil palm basal stem rot (BSR) and other Ganoderma species, two specific primer pairs were selected to provide a specific DNA amplification of pathogenic Ganoderma in oil palm. Each primer pair produced a single PCR product of about 450 bp (for primer pair IT1–IT2) and 334 bp (for primer pair IT1–IT3) when oil palm Ganoderma DNA was used. No PCR amplification product was observed when other Ganoderma species DNA was used in PCR amplification with these primer pairs. Three specific restriction enzyme sites were identified in the ITS and intergenic spacer (IGS1) regions. The restriction enzymes MluI, SacI and HinfI were used to digest the ITS-PCR product and restriction enzymes TfiI, ScaI and HincII were used to digest the IGS1-PCR product. Of the three restriction enzymes used in each rDNA region, MluI specifically digested the ITS regions, and TfiI specifically digested the IGS1 region of oil palm Ganoderma. Analysis of the published ITS nucleotide sequences of 31 Ganoderma species showed that the MluI restriction site was not present in other Ganoderma species. The use of both specific primers and restriction enzyme analysis can be applied as a standard protocol to identify pathogenic Ganoderma in oil palm. In this study, the use of specific primers and PCR-RFLP analyses of the rDNA gave consistent results for the characterisation of pathogenic Ganoderma, and indicated that Ganoderma strains associated with BSR disease in oil palms belong to a single species.     

Restriction analysis of the amplified ribosomal DNA spacers ITS1 and ITS2 of <Emphasis Type="Italic">Bipolaris sorokiniana</Emphasis> isolates

The purpose of this study was to examine the genotypic variability of Bipolaris sorokiniana by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis using rDNA. Fifty B. sorokiniana isolates from Brazil and other countries, one Bipolaris oryzae and six Drechslera teres isolates were used. The intragenic spacer regions (ITS1 and ITS2) were the regions used for characterization of isolates. The amplification products for both ITS regions, showed two DNA fragments for all isolates. Two B. sorokiniana isolates presented an intraspecific variability showing a third fragment for the ITS1 region. The dendrograms generated with PCR-RFLP data showed intra- and inter-specific groups. The dendrograms showed that most of Brazilian isolates clustered together forming groups between them, and this behavior was repeated with most isolates from other countries. The dendrograms did not enable the separation of B. sorokiniana isolates by their geographic origin or host type. These results suggest the occurrence of gene flow between different populations of the fungus isolated in geographically distant regions and lends cogency to the occurrence of gene flow between species.     

Analysis of nucleotide sequence polymorphism of internal transcribed spacers of ribosomal genes in diploid <Emphasis Type="Italic">Aegilops</Emphasis> (L.) species

The nucleotide sequence of the fragment of the internal transcribed spacer (ITS) of rDNA comprising the full-length ITS1, the gene encoding 5.8S rRNA, and part of the ITS2 sequence was determined in 22 samples of five diploid Aegilops species. The full alignment length of compared sequences was 524 bp. Species-specific substitutions were found in the ITS nucleotide sequence of rDNA of different Aegilops species. Intraspecific differences in ITS structure in diploid Aegilops species were detected for the first time. Polymorphism of the ITS nucleotide sequence within the same sample was revealed, which might be due either to differences between the genomes of individual plants comprising the sample or to the presence of several types of ribosomal genes in the genome of one plant. In general, both interspecific and intraspecific variability of the ITS nucleotide sequences of rDNA is extremely low. In total, 26 variable sites, twelve of which were informative, were identified.__________Translated from Molekulyarnaya Biologiya, Vol. 39, No. 2, 2005, pp. 193–197.Original Russian Text Copyright © 2005 by Goryunova, Chikida, Gori, Kochieva.     

Close sequence identity between ribosomal DNA episomes of the non-pathogenicEntamoeba dispar and pathogenicEntamoeba histolytica

Entamoeba dispar andEntamoeba histolytica are now recognized as two distinct species-the former being nonpathogenic to humans. We had earlier studied the organization of ribosomal RNA genes inE. histolytica. Here we report the analysis of ribosomal RNA genes inE. dispar. The rRNA genes ofE. dispar, like their counterpart inE. histolytica are located on a circular rDNA molecule. From restriction map analysis, the size ofE. dispar rDNA circle was estimated to be 24·4 kb. The size was also confirmed by linearizing the circle withBsaHI, and by limited DNAseI digestion. The restriction map of theE. dispar rDNA circle showed close similarity to EhR1, the rDNA circle ofE. histolytica strain HM-1:IMSS which has two rDNA units per circle. The various families of short tandem repeats found in the upstream and downstream intergenic spacers (IGS) of EhR1 were also present inE. dispar. Partial sequencing of the cloned fragments ofE. dispar rDNA and comparison with EhR1 revealed only 2·6% to 3·8% sequence divergence in the IGS. The region Tr and the adjoiningPvuI repeats in the IGS of EhR1, which are missing in thoseE. histolytica strains that have one rDNA unit per circle, were present in theE. dispar rDNA circle. Such close similarity in the overall organization and sequence of the IGS of rDNAs of two different species is uncommon. In fact the spacer sequences were only slightly more divergent than the 18S rRNA gene sequence which differs by 1·6% in the two species. The most divergent sequence betweenE. histolytica andE. dispar was the internal transcribed spacer, ITS2. Therefore, it was concluded that probes derived from the ITS1 and ITS 2 sequences would be more reliable and reproducible than probes from the IGS regions used earlier for identifying these species.     

Reproductive compatibility and genetic variation between two strains of Aphelinus albipodus (Hymenoptera: Aphelinidae), a parasitoid of the soybean aphid, Aphis glycines (Homoptera: Aphididae)

Aphelinus albipodus Hayat and Fatima is a potential biological control agent of the soybean aphid, Aphis glycines Matsumura, which is a newly introduced soybean pest in the United States. We compared the reproductive compatibility and molecular genetic variation between two geographic strains of A. albipodus. One strain was collected from soybean aphids in Japan and the other recovered from Russian wheat aphid, Diuraphis noxia (Mordvilko), in the western U.S., populations of which were established with parasitoids imported from Eurasia. We present results of crossing experiments between the two strains, genetic differences based on RAPD-PCR markers, rDNA ITS1 and ITS2 gene sequences, and presence of Wolbachia in the two strains using PCR amplification of the wsp gene. We found no reduction in the production of females in reciprocal crosses between strains, but a significant reduction in fecundity when F₁ females stemming from one of the reciprocal crosses were backcrossed to males from either source. The two strains differed by 3.4% in the rDNA ITS1 sequence and by presence/absence of one RAPD-PCR marker from a total of 20 RAPD primers screened, but their rDNA ITS2 sequences were identical. We used restriction enzyme analysis to separate the strains by differential digestion of the ITS1 PCR product. Wolbachia was present in 100% of males and females of both strains of A. albipodus.     

Genetic relationship of Tricholoma matsutake and T. nauseosum from the Northern Hemisphere based on analyses of ribosomal DNA spacer regions

The genetic relationship among Tricholoma matsutake and T. nauseosum strains collected from various parts of the Northern Hemisphere was investigated using sequence analysis of the rDNA ITS region and PCR-RFLP analysis of the rDNA IGS-1 region. ITS sequence similarity between T. matsutake and T. nauseosum ranged between 98.1% and 100%. The strains of T. matsutake from coniferous forests and those from broad-leaved forests showed more than 99.8% similarity in their ITS sequences. Three distinct RFLP types were detected when IGS-1 regions were digested with Cfr13I. RFLP patterns showed no variability among the strains of T. nauseosum and those of T. matsutake from broad-leaved forests. This pattern corresponded to the dominant RFLP type in the Japanese population of T. matsutake. Thus, strains belonging to this RFLP type are widely distributed throughout East Asia and Europe and associated with many tree species of Pinaceae and Fagaceae. The result suggests that T. matsutake in coniferous and broad-leaved forests and T. nauseosum should be treated as the same species genetically.     

Molecular phylogeny in Indian Tinospora species by DNA based molecular markers

The phylogenetic relationships among the three species of Tinospora found in India are poorly understood. Morphology does not fully help to resolve the phylogeny and therefore a fast approach using molecular analysis was explored. Two molecular approaches viz Random Amplified Polymorphic DNA (RAPD) assay and restriction digestion of ITS1-5.8S-ITS2 rDNA (PCR-RFLP) were used to evaluate the genetic similarities between 40 different accessions belonging to three species. Of the 38 random primers used only six generated the polymorphism, while as three out of 11 restriction enzymes used gave polymorphic restriction patterns. The average proportion of polymorphic markers across primers was 95%, however restriction endonucleases showed 92% polymorphism. RAPD alone was found suitable for the species diversions. In contrast PCR- RFLP showed bias in detecting exact species variation. The correlation between the two markers was performed by Jaccard's coefficient of similarity. A significant (r= 0.574) but not very high correlation was obtained. Further to authenticate the results obtained by two markers, sequence analysis of ITS region of ribosomal DNA (ITS1 and ITS2, including 5.8S rDNA) was performed. Three independent clones of each species T. cordifolia, T. malabarica and T. crispa were sequenced. Phylogenetic relationship inferred from ITS sequences is in agreement with RAPD data.     

Nucleotide Sequences of Plasmodium vivaxA-Type rDNA ITS1, 5.8S,and ITS2. Elaboration of Retrospective PCR Diagnostics of Malaria Using Stained Thick Blood Smears

Stages life cycle of the malaria parasite differ in the rate of replication and the structural properties of functionally active A-, S-, and O-type ribosomes. Regions of A-type rDNA including ITS1, 5.8S, and ITS2 from two strains of Plasmodium vivaxwith different incubation periods were amplified and sequenced. No substantial differences in the sequences of two strains were revealed. Phylogenetic analysis of the obtained and homologous sequences of ITS1 rDNA of A, S, and O types of P. vivax; A and S types of P. falciparum; and Cryptosporidium parvum, Eimeria maxima, Toxoplasma gondiias outgroup, by the maximum parsimony method using PAUP 4.0 revealed that divergence of ITS1 might have occurred after speciation and at different rates in individual lineages of the Plasmodiumgenus. Basing on the results of the analysis of orthologous sequences of P. vivaxand P. falciparum, we developed genus- and species-specific primers for PCR diagnostics of malaria, as well as a one-step effective method of DNA isolation from Giemsa–Romanovsky-stained thick blood smears. It was demonstrated that stained preparations could be a reliable source of plasmodial DNA, and the quality of preparations and storage time (10–20 years) did not interfere with the results of PCR analysis.     

Ethiopian soils harbor natural populations of rhizobia that form symbioses with common bean (<Emphasis Type="Italic">Phaseolus vulgaris</Emphasis> L.)

<p>The diversity and taxonomic relationships of 83 bean-nodulating rhizobia indigenous to Ethiopian soils were characterized by PCR-RFLP of the internally transcribed spacer (ITS) region between the 16S and 23S rRNA genes, 16S rRNA gene sequence analysis, multilocus enzyme electrophoresis (MLEE), and amplified fragment-length polymorphism. The isolates fell into 13 distinct genotypes according to PCR-RFLP analysis of the ITS region. Based on MLEE, the majority of these genotypes (70%) was genetically related to the type strain of Rhizobium leguminosarum. However, from analysis of their 16S rRNA genes, the majority was placed with Rhizobium etli. Transfer and recombination of the 16S rRNA gene from presumptively introduced R. etli to local R. leguminosarum is a possible theory to explain these contrasting results. However, it seems unlikely that bean rhizobia originating from the Americas (or Europe) extensively colonized soils of Ethiopia because Rhizobium tropici, Rhizobium gallicum, and Rhizobium giardinii were not detected and only a single ineffective isolate of R. etli that originated from a remote location was identified. Therefore, Ethiopian R. leguminosarum may have acquired the determinants for nodulation of bean from a low number of introduced bean-nodulating rhizobia that either are poor competitors for nodulation of bean or that failed to survive in the Ethiopian environment. Furthermore, it may be concluded from the genetic data presented here that the evidence for separating R. leguminosarum and R. etli into two separate species is inconclusive.     

Identification of Agrobacterium strains by PCR-RFLP analysis of pTi and chromosomal regions

Chromosomes and Ti plasmids of 41 Agrobacterium strains, belonging to biovars 1, 2, 3, and Agrobacterium rubi species were characterized by the restriction fragment length polymorphism of PCR-amplified DNAs. Profiles that were obtained by the analysis of the amplified 16S rDNA confirmed the grouping of the strains according to their species. Higher polymorphism was detected in the intergenic spacer between the 16S rDNA and 23S rDNA genes, allowing efficient discrimination of strains. Identification of most strains was possible, and the genetic relatednesses of Agrobacterium strains could be estimated. The analysis of the plasmid Ti encoded regions between the tmr and nos genes, and the virA and virB₂ genes, allowed fingerprinting of Ti plasmids. Genomic typing by the rapid PCR-RFLP method is thus shown to be useful for an independant identification of strains and of the conjugative Ti plasmids.Abbreviations PCR polymerase chain reaction - RFLP restriction fragment length polymorphism - IGS intergenic spacer Funded by Institut National de la Recherche Agronomique     

Intra-isolate heterogeneity of the ITS region of rDNA in Pythium helicoides

Heterogeneity of the rDNA ITS region in Pythium helicoides and the phylogenetic relationship between P. helicoides and closely related species were investigated. In PCR-RFLP analysis of the rDNA ITS region of six P. helicoides isolates investigated, including the type culture, intraspecific variation was found at the HhaI site. The total length of fragments was longer than before cutting, indicating sequence heterogeneity within isolates. Digestion of the cloned rDNA ITS region derived from seven isolates with HhaI revealed polymorphisms among and within single zoospore isolates, and variability of the region was also present among the clones derived from the same isolate. To test whether the rDNA ITS region of closely related species and other regions in the genome of P. helicoides are also variable, the rDNA ITS region of P. ultimum and the cytochrome oxydase II (cox II) gene encoded in mitochondria were sequenced. P. ultimum had little variation in the rDNA ITS region. The cox II gene sequences of both species revealed only a low intraspecific variability and no intra-isolate variation. In the phylogenic tree based on the rDNA ITS sequences, all clones of P. helicoides formed one large clade that was distinct from the clades comprising morphologically similar species, such as P. oedochilum and P. ostracodes, and was closely related to P. chamaehyphon rather than the other species.     

Molecular divergence of the soil yeasts <Emphasis Type="Italic">Williopsis</Emphasis> sensu stricto

Fifty-three strains of Saturn-spored yeasts were analyzed by means of restriction analysis of the amplified fragment of rDNA comprising the 5.8S rRNA gene and the internal transcribed spacers ITS1 and ITS2. The use of endonucleases HaeIII and MspI enabled clear differentiation of yeast species Williopsis mucosa, W. salicorniae, Zygowilliopsis californica, and Komagataea pratensis and the Williopsis sensu stricto complex. The minisatellite primer M13 was proposed for differentiation between sibling species of Williopsis sensu stricto, which have identical restriction profiles. PCR with primer M13 enabled reidentification of a number of collection strains, species identification of Saturn-spored isolates from the Far East, and detection of three strains affiliated to novel taxa. The latter have unique PCR profiles and differ in the nucleotide sequences of ITS1 and ITS2 fragments of rDNA. Possible variations in the results obtained with different molecular methods are discussed.Translated from Mikrobiologiya, Vol. 73, No. 6, 2004, pp. 768–776.Original Russian Text Copyright © 2004 by Naumova, Gazdiev, Naumov.     

Comparison of the ribosomal RNA genes in four closely relatedCucurbitaceae

Cucurbitaceae are characterized by a high copy number for nuclear ribosomal RNA genes. We have investigated the genomic ribosomal DNA (rDNA) of four closely related species of this family with respect to structure, length heterogeneity, and evolution. InCucumis melo (melon) there are two main length variants of rDNA repeats with 10.7 and 10.55kb.Cucumis sativus (cucumber) shows at least three repeat types with 11.5, 10.5, and 10.2kb.Cucurbita pepo (zucchini) has two different repeat types with 10.0 and 9.3kb. There are also two different repeat types inCucurbita maxima (pumpkin) of about 11.2 and 10.5kb. Restriction enzyme mapping of the genomic rDNA of these four plants and of cloned repeats ofC. sativus shows further heterogeneities which are due to methylation or point mutations. By comparison of the restriction enzyme maps it was possible to trace some evolutionary events in the family ofCucurbitaceae. Some aspects of regulation and function of the middle repetitive rRNA genes (here between 2000 and 10000 copies) are discussed.     

Extensive 5.8S nrDNA polymorphism in <Emphasis Type="Italic">Mammillaria</Emphasis> (Cactaceae) with special reference to the identification of pseudogenic internal transcribed spacer regions

The internal transcribed spacer (ITS) region (ITS1, 5.8S rDNA, ITS2) represents the most widely applied nuclear marker in eukaryotic phylogenetics. Although this region has been assumed to evolve in concert, the number of investigations revealing high degrees of intra-individual polymorphism connected with the presence of pseudogenes has risen. The 5.8S rDNA is the most important diagnostic marker for functionality of the ITS region. In Mammillaria, intra-individual 5.8S rDNA polymorphisms of up to 36% and up to nine different types have been found. Twenty-eight of 30 cloned genomic Mammillaria sequences were identified as putative pseudogenes. For the identification of pseudogenic ITS regions, in addition to formal tests based on substitution rates, we attempted to focus on functional features of the 5.8S rDNA (5.8S motif, secondary structure). The importance of functional data for the identification of pseudogenes is outlined and discussed. The identification of pseudogenes is essential, because they may cause erroneous phylogenies and taxonomic problems. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Allotetraploid origin ofAllium altyncolicum (Alliaceae, Allium sect.Schoenoprasum) as investigated by karyological and molecular markers

The tetraploidAllium altyncolicum (2n = 4x = 32) is considered to be of hybrid origin, because most of its morphological characters are intermediate between those of its putative parents,A. schoenoprasum andA. ledebourianum. In the present work an attempt has been made to ascertain its parentage by several methods: Giemsa C-banding, genomic in situ hybridization (GISH), PCR-RFLP of cpDNA, restriction enzyme mapping of the rDNA, and RAPDs. C-banding and GISH indicates clearly thatA. altyncolicum is a segmental allopolyploid.Allium schoenoprasum andA. ledebourianum are the most likely the parental species and the larger part of the genome ofA. altyncolicum (26 chromosomes) is derived fromA. schoenoprasum. The low genetic divergence between these three species was confirmed by the lack of sequence variation in the ITS sequences of nuclear rRNA genes and of the plastid rbcL-atpB intergenic spacer. Both parental species andA. altyncolicum could be distinguished by RFLP of the rDNA repeats. The geographic origin of the putative parental species was investigated using RAPDs.Dedicated to emer. Univ.-Prof. DrFriedrich Ehrendorfer on the occasion of his 70th birthday     

Phylogenetic analysis of Allium subg. Melanocrommyum infers cryptic species and demands a new sectional classification

Allium subgenus Melanocrommyum (Alliaceae) from Eurasia comprises about 150 mostly diploid species adapted to arid conditions. The group is taxonomically complicated with different and contradictory taxonomic treatments, and was thought to include a considerable number of hybrid species, as the taxa show an admixture of assumed morphological key characters. We studied the phylogeny of the subgenus, covering all existing taxonomic groups and their entire geographic distribution. We analyzed sequences of the nuclear rDNA internal transcribed spacer region (ITS) for multiple individuals of more than 100 species. Phylogenetic analyses of cloned and directly sequenced PCR products confirmed the monophyly of the subgenus, while most sections were either para- or polyphyletic. The splits of the large sections are supported by differences in the anatomy of flower nectaries. ITS data (i) demand a new treatment at sectional level, (ii) do not support the hypotheses of frequent gene flow among species, (iii) indicate that multiple rapid radiations occurred within different monophyletic groups of the subgenus, and (iv) detected separately evolving lineages within three morphologically clearly defined species (cryptic species). In two cases these lineages were close relatives, while in Allium darwasicum they fall in quite different clades in the phylogenetic tree. Fingerprint markers show that this result is not due to ongoing introgression of rDNA (ITS capture) but that genome-wide differences between both lineages exist. Thus, we report one of the rare cases in plants where morphologically indistinguishable diploid species occurring in mixed populations are non-sister cryptic species.     

Ribosomal DNA variation within and among individuals ofLisianthius (Gentianaceae) populations

Restriction endonuclease fragment analysis of nuclear ribosomal DNA (rDNA) was completed on 25 individuals each from seven populations of theLisianthius skinneri (Gentianaceae) species complex in Panama. Seven restriction enzymes were used to determine the amount and type of rDNA variation within and among individuals of the populations. No restriction site variation was seen within populations or individuals although site differences were seen among populations. Spacer length variation within and among individuals of populations was mapped to the internal transcribed spacer (ITS) region between the 18S and 5.8S rRNA genes, a region inLisianthius rDNA that previously was shown to exhibit length differences among populations. This is the first reported case of such variation within and among individuals of populations for the ITS region. Presence or absence of ITS spacer length variation is not correlated with levels of isozymic heterozygosity within populations. No detectable length variation within individuals or populations was seen in the larger intergenic spacer (IGS). Although populations varied with respect to IGS length, all individuals of a given population had a single and equivalent IGS length.     

Two new species of <Emphasis Type="Italic">Trichoderma</Emphasis> from Yunnan,China

Two new species of the fungal genus Trichoderma, Trichoderma compactum and Trichoderma yunnanense, isolated from rhizosphere of tobacco in Yunnan Province, China are described based on morphological characters and phylogenetic analyses of nucleotide sequences. Our DNA sequences included the internal transcribed spacer (ITS) regions of the rDNA cluster (ITS1 and ITS2), and partial sequences of the translation elongation factor 1-alpha (tef1) and a fragment of the gene coding for endochitinase 42 (ech42). The analyses show that T. compactum belongs to the Harzianum clade, and T. yunnanense belongs to the Hamatum clade.     

The sequence of the isoepoxydon dehydrogenase gene of the patulin biosynthetic pathway in <Emphasis Type="Italic">Penicillium</Emphasis> species

Interest in species of the genus Penicillium is related to their ability to produce the mycotoxin patulin and to cause spoilage of fruit products worldwide. The sequence of the isoepoxydon dehydrogenase (idh) gene, a gene in the patulin biosynthetic pathway, was determined for 28 strains representing 12 different Penicillium species known to produce the mycotoxin patulin. Isolates of Penicillium carneum, Penicillium clavigerum, Penicillium concentricum, Penicillium coprobium, Penicillium dipodomyicola, Penicillium expansum, Penicillium gladioli, Penicillium glandicola, Penicillium griseofulvum, Penicillium paneum, Penicillium sclerotigenum and Penicillium vulpinum were compared. Primer pairs for DNA amplification and sequencing were designed from the P. griseofulvum idh gene (GenBank AF006680). The two introns present were removed from the nucleotide sequences, which were translated to produce the IDH sequences of the 12 species for comparison. Phylogenetic relationships among the species were determined from rDNA (ITS1, 5.8 S, ITS2 and partial sequence of 28S rDNA) and from the idh nucleotide sequences minus the two introns. Maximum parsimony analysis showed trees based on rDNA and idh sequences to be congruent. It is anticipated that the genetic information obtained in the present study will aid in the design of probes, specific for patulin biosynthetic pathway genes, to identify the presence of these mycotoxigenic fungi. The U.S. Government's right to retain a non-exclusive, royalty-free license in and to any copyright is acknowledged.     

Morphological and molecular phylogenetic analysis of <Emphasis Type="Italic">Melampsora</Emphasis> species on poplars in China

Many species of Melampsora on Populus have been reported in China, based on morphological characteristics of both uredial and telial states, and on host species, but their morphology and taxonomy are still poorly defined. In this study, 196 specimens representing Melampsora species on poplars and collected from various areas of China were used for morphological observations. The morphological characteristics of urediniospores and teliospores were examined with light and scanning electron microscopy. The specimens could be classified into five groups based on their morphology. For the sequencing of the nuclear large subunit rDNA (D1/D2), 5.8S rDNA and their internal transcribed spacers, ITS1 and ITS2 region, 54 specimens were selected from the specimens used in morphological observations. These specimens were separated into six clades by phylogenetic analyses of the D1/D2 and ITS regions. Correlations among morphological groups and phylogenetic clades based on these results suggest a revision of these species. In particular, no evidence to discriminate specimens of M. acedioides, M. magnusiana, and M. rostrupii was found from either morphological characteristics or sequence analysis.Contribution no. 185 Laboratory of Plant Parasitic Mycology, Institute of Agriculture and Forestry, University of Tsukuba, Japan