Search

Chemolithoautotrophic bacterium Ralstonia eutropha is a versatile host for production of various useful compounds including polyhydroxyalkanoates (PHAs) under both heterotrophic and autotrophic conditions. In this bacterium, Calvin–Benson–Bassham (CBB) cycle is functional even under heterotrophic conditions on sugars and reutilizes CO2 emitted through sugar metabolisms into PHA, leading to increase in yield of the storage polyester. This study focused on isopropanol production from glucose by engineered strains of R. eutropha. The isopropanol-producing strains were constructed by introduction of codon-optimized genes of acetoacetate decarboxylase (adc) and primary-secondary alcohol dehydrogenase (adh) from clostridia into glucose-utilizing and PHA-negative (ΔphaC1) strain of R. eutropha. Several genetic modifications showed that high expression of the isopropanol synthesis genes by using a strong synthetic promoter and deletion of NAD+-dependent (S)-3-hydroxybutyryl-CoA dehydrogenase genes (paaH1 and had) in addition to NADPH-dependent acetoacetyl-CoA reductase genes (phaB1 and phaB3) were effective for improving isopropanol production with low by-production of acetone. Isopropanol titer of 4.13 g/L was achieved by two-stage cultivation of the strain IP-007/pBj5c2-adh-adc, corresponding to overall yield of 0.6 mol mol-glucose−1. The fixation of sugar-derived CO2 during isopropanol synthesis was evaluated by 13C-labelling of the isopropanol produced from [1–13C]-glucose. The 13C-abundance in isopropanol synthesized by the engineered strain was significantly increased up to 4.8%, demonstrating actual reassimilation of CO2 emitted from glucose moiety by decarboxylation and potential contribution towards increase in the carbon yield of isopropanol on glucose.

Advances in genomics have expedited the improvement of several agriculturally important crops but similar efforts in wheat (Triticum spp.) have been more challenging. This is largely owing to the size and complexity of the wheat genome1, and the lack of genome-assembly data for multiple wheat lines2,3. Here we generated ten chromosome pseudomolecule and five scaffold assemblies of hexaploid wheat to explore the genomic diversity among wheat lines from global breeding programs. Comparative analysis revealed extensive structural rearrangements, introgressions from wild relatives and differences in gene content resulting from complex breeding histories aimed at improving adaptation to diverse environments, grain yield and quality, and resistance to stresses4,5. We provide examples outlining the utility of these genomes, including a detailed multi-genome-derived nucleotide-binding leucine-rich repeat protein repertoire involved in disease resistance and the characterization of Sm16, a gene associated with insect resistance. These genome assemblies will provide a basis for functional gene discovery and breeding to deliver the next generation of modern wheat cultivars., Comparison of multiple genome assemblies from wheat reveals extensive diversity that results from the complex breeding history of wheat and provides a basis for further potential improvements to this important food crop.

Excessive amounts of metal ions in soil are toxic for most plant species, yet metal can also facilitate plant survival by elemental defense against herbivores and pathogens. Zinc and cadmium hyperaccumulation in Arabidopsis halleri is known to be effective for the defense against natural enemies. The allotetraploid species A. kamchatica, derived from A. halleri and a non‐hyperaccumulator A. lyrata, has a lower hyperaccumulation level of zinc than A. halleri, but its significance for elemental defense remains unknown. In this study, we evaluated the accumulation levels of zinc and cadmium in the allotetraploid compared with its diploid progenitors, and evaluated the contribution of metal treatments to anti‐herbivore resistance under field conditions. The accumulation level of zinc in A. kamchatica was intermediate between the progenitors, but that of cadmium was lower than in both diploid progenitors. The elemental defense of A. kamchatica and A. halleri was supported by a field experiment comparing the herbivory level between a control group and metal‐supplemented plants. Moreover, the effect of elemental defense was lower in A. kamchatica than in the hyperaccumulator progenitor A. halleri, which is consistent with the metal accumulation level. This result reveals that the allotetraploid plant inherited its hyperaccumulating ability from one progenitor as an advantageous trait but at an intermediate level.

Bread wheat is a major crop that has long been the focus of basic and breeding research. Assembly of its genome has been difficult because of its large size and allohexaploid nature (AABBDD genome). Following the first reported assembly of the genome of the experimental strain Chinese Spring (CS), the 10+ Wheat Genomes Project was launched to produce multiple assemblies of worldwide modern cultivars. The only Asian cultivar in the project is Norin 61, a representative Japanese cultivar adapted to grow across a broad latitudinal range, mostly characterized by a wet climate and a short growing season. Here, we characterize the key aspects of its chromosome-scale genome assembly spanning 15 Gb with a raw scaffold N50 of 22 Mb. Analysis of the repetitive elements identified chromosomal regions unique to Norin 61 that encompass a tandem array of the pathogenesis-related 13 family. We report novel copy-number variations in the B homeolog of the florigen gene FT1/VRN3, pseudogenization of its D homeolog and the association of its A homeologous alleles with the spring/winter growth habit. Furthermore, the Norin 61 genome carries typical East Asian functional variants different from CS, ranging from a single nucleotide to multi-Mb scale. Examples of such variation are the Fhb1 locus, which confers Fusarium head-blight resistance, Ppd-D1a, which confers early flowering, Glu-D1f for Asian noodle quality and Rht-D1b, which introduced semi-dwarfism during the green revolution. The adoption of Norin 61 as a reference assembly for functional and evolutionary studies will enable comprehensive characterization of the underexploited Asian bread wheat diversity.

Wild Japanese macaques (Macaca fuscata Blyth) living in the highland and lowland areas of Yakushima are known to have different diets, with highland individuals consuming more leaves. We aim to clarify whether and how these differences in diet are also reflected by gut microbial composition and fermentation ability. Therefore, we conduct an in vitro fermentation assay using fresh feces from macaques as inoculum and dry leaf powder of Eurya japonica Thunb. as a substrate. Fermentation activity was higher for feces collected in the highland, as evidenced by higher gas and butyric acid production and lower pH. Genetic analysis indicated separation of highland and lowland in terms of both community structure and function of the gut microbiota. Comparison of feces and suspension after fermentation indicated that the community structure changed during fermentation, and the change was larger for lowland samples. Analysis of the 16S rRNA V3-V4 barcoding region of the gut microbiota showed that community structure was clearly clustered between the two areas. Furthermore, metagenomic analysis indicated separation by gene and pathway abundance patterns. Two pathways (glycogen biosynthesis I and D-galacturonate degradation I) were enriched in lowland samples, possibly related to the fruit-eating lifestyle in the lowland. Overall, we demonstrated that the more leaf-eating highland Japanese macaques harbor gut microbiota with higher leaf fermentation ability compared with the more fruit-eating lowland ones. Broad, non-specific taxonomic and functional gut microbiome differences suggest that this pattern may be driven by a complex interplay between many taxa and pathways rather than single functional traits.

Interspecific gene flow by hybridization may weaken species barriers and adaptive divergence, but can also initiate reinforcement of reproductive isolation trough natural and sexual selection. The extent of interspecific gene flow and its consequences for the initiation and maintenance of species barriers in natural systems remain poorly understood, however. To assess genome-wide patterns of gene flow between the two closely related European dung fly species Sepsis cynipsea and Sepsis neocynipsea (Diptera: Sepsidae), we tested for historical gene flow with the aid of ABBA-BABA test using whole-genome resequencing data from pooled DNA of male specimens originating from natural and laboratory populations. We contrasted genome-wide variation in DNA sequence differences between samples from sympatric populations of the two species in France and Switzerland with that of interspecific differences between pairs of samples involving allopatric populations from Estonia and Italy. In the French Cevennes, we detected a relative excess of DNA sequence identity, suggesting interspecific gene flow in sympatry. In contrast, at two sites in Switzerland, we observed a relative depletion of DNA sequence identity compatible with reinforcement of species boundaries in sympatry. Our results suggest that the species boundaries between S. cynipsea and S. neocynipsea in Europe depend on the eco-geographic context.

Background: Whole genome duplication (WGD) events are common in the evolutionary history of many living organisms. For decades, researchers have been trying to understand the genetic and epigenetic impact of WGD and its underlying molecular mechanisms. Particular attention was given to allopolyploid study systems, species resulting from an hybridization event accompanied by WGD. Investigating the mechanisms behind the survival of a newly formed allopolyploid highlighted the key role of DNA methylation. With the improvement of high-throughput methods, such as whole genome bisulfite sequencing (WGBS), an opportunity opened to further understand the role of DNA methylation at a larger scale and higher resolution. However, only a few studies have applied WGBS to allopolyploids, which might be due to lack of genomic resources combined with a burdensome data analysis process. To overcome these problems, we developed the Automated Reproducible Polyploid EpiGenetic GuIdance workflOw (ARPEGGIO): the first workflow for the analysis of epigenetic data in polyploids. This workflow analyzes WGBS data from allopolyploid species via the genome assemblies of the allopolyploid’s parent species. ARPEGGIO utilizes an updated read classification algorithm (EAGLE-RC), to tackle the challenge of sequence similarity amongst parental genomes. ARPEGGIO offers automation, but more importantly, a complete set of analyses including spot checks starting from raw WGBS data: quality checks, trimming, alignment, methylation extraction, statistical analyses and downstream analyses. A full run of ARPEGGIO outputs a list of genes showing differential methylation. ARPEGGIO was made simple to set up, run and interpret, and its implementation ensures reproducibility by including both package management and containerization. Results: We evaluated ARPEGGIO in two ways. First, we tested EAGLE-RC’s performance with publicly available datasets given a ground truth, and we show that EAGLE-RC decreases the error rate by 3 to 4 times compared to standard approaches. Second, using the same initial dataset, we show agreement between ARPEGGIO’s output and published results. Compared to other similar workflows, ARPEGGIO is the only one supporting polyploid data. Conclusions: The goal of ARPEGGIO is to promote, support and improve polyploid research with a reproducible and automated set of analyses in a convenient implementation. ARPEGGIO is available at https://github.com/supermaxiste/ARPEGGIO., MC genomics, 22 (1), ISSN:1471-2164

Summary Polyploidization is pervasive in plants, but little is known about the niche divergence of wild allopolyploids (species that harbor polyploid genomes originating from different diploid species) relative to their diploid progenitor species and the gene expression patterns that may underlie such ecological divergence. We conducted a fine‐scale empirical study on habitat and gene expression of an allopolyploid and its diploid progenitors.We quantified soil properties and light availability of habitats of an allotetraploid Cardamine flexuosa and its diploid progenitors Cardamine amara and Cardamine hirsuta in two seasons. We analyzed expression patterns of genes and homeologs (homeologous gene copies in allopolyploids) using RNA sequencing.We detected niche divergence between the allopolyploid and its diploid progenitors along water availability gradient at a fine scale: the diploids in opposite extremes and the allopolyploid in a broader range between diploids, with limited overlap with diploids at both ends. Most of the genes whose homeolog expression ratio changed among habitats in C. flexuosa varied spatially and temporally.These findings provide empirical evidence for niche divergence between an allopolyploid and its diploid progenitor species at a fine scale and suggest that divergent expression patterns of homeologs in an allopolyploid may underlie its persistence in diverse habitats.

Contemporary speciation provides a unique opportunity to directly observe the traits and environmental responses of a new species. Cardamine insueta is an allotriploid species that appeared within the past 150 years in a Swiss village, Urnerboden. In contrast to its two progenitor species, C. amara and C. rivularis that live in wet and open habitats, respectively, C. insueta is found in-between their habitats with temporal water level fluctuation. This triploid species propagates clonally and serves as a triploid bridge to form higher ploidy species. Although niche separation is observed in field studies, the mechanisms underlying the environmental robustness of C. insueta are not clear. To characterize responses to a fluctuating environment, we performed a time-course analysis of homeolog gene expression in C. insueta in response to submergence treatment. For this purpose, the two parental (C. amara and C. rivularis) genome sequences were assembled with a reference-guided approach, and homeolog-specific gene expression was quantified by using HomeoRoq software. We found that C. insueta and C. rivularis initiated vegetative propagation by forming ectopic meristems on leaves, while C. amara did not. We examined homeolog-specific gene expression of three species at nine time points during the treatment. The genome-wide expression ratio of homeolog pairs was 2:1 over the time-course, consistent with the ploidy number. By searching the genes with high coefficient of variation of expression over time-course transcriptome data, we found many known key transcriptional factors related to meristem development and formation upregulated in both C. rivularis and rivularis-homeolog of C. insueta, but not in C. amara. Moreover, some amara-homeologs of these genes were also upregulated in the triploid, suggesting trans-regulation. In turn, Gene Ontology analysis suggested that the expression pattern of submergence tolerant genes in the triploid was inherited from C. amara. These results suggest that the triploid C. insueta combined advantageous patterns of parental transcriptomes to contribute to its establishment in a new niche along a water-usage gradient.

To study the genetic structure of clonal plant populations, genotyping and genet detection using genetic markers are necessary to assign ramets to corresponding genets. Assignment is difficult as it involves setting a robust threshold of genetic distance for genet distinction as neighbouring genets in a plant population are often genetically related. Here, we used restriction site-associated DNA sequencing (RAD-seq) for a rhizomatous clonal herb, Cardamine leucantha [Brassicaceae] to accurately determine genet structure in a natural population. We determined a draft genome sequence of this species for the first time, which resulted in 66 617 scaffolds with N50 = 6086 bp and an estimated genome size of approximately 253 Mbp. Using genetic distances based on the RAD-seq analysis, we successfully distinguished ramets that belonged to distinct genets even from a half-sib family. We applied these methods to 372 samples of C. leucantha collected at 1-m interval grids within a 20 × 20 m plot in a natural population in Hokkaido, Japan. From these samples, we identified 61 genets with high inequality in terms of genet size and patchy distribution. Spatial autocorrelation analyses indicated significant aggregation within 7 and 4 m at ramet and genet levels, respectively. An analysis of parallel DNA microsatellite loci (simple sequence repeats) suggested that RAD-seq can provide data that allows robust genet assignment. It remains unclear whether the large genets identified here became dominant stochastically or deterministically. Precise identification of genets will assist further study and characterization of dominant genets., Clonal plants produce genetically identical offspring through vegetative growth. A ‘ramet’ refers to a single physiological individual produced by clonal propagation; a ‘genet’ refers to a group of ramets that originate from a single seed. Here we investigated how ramets that belong to different genets distribute in a natural population of Carrdamine leucantha using a genome-wide SNP markers. We applied the method to 372 samples collected at 1-m interval grids within a 20 × 20 m plot, and identified 61 genets with high inequality in terms of genet size and patchy distribution.

Population genetics studies of allopolyploid species lag behind those of diploid species because of practical difficulties in analysis of homeologs-duplicated gene copies originating from hybridized parental species. Pool-Seq, i.e. massive parallel sequencing of pooled individuals, has high potential for detecting nucleotide polymorphisms within and among multiple populations; however, its use has been limited to diploid species. We applied Pool-Seq to an allopolyploid species by developing a bioinformatic pipeline that assigns reads to each homeolog as well as to each polymorphic allele within each homeolog. We simultaneously sequenced eight genes from twenty individuals from each of 24 populations, and found over 100 polymorphic sites in each homeolog. For two sites, we estimated allele frequencies using the number of reads and then validated these estimations by making individual-based estimations. Pool-Seq using our bioinformatic pipeline allows efficient evaluation of nucleotide polymorphisms in a large number of individuals, even in allopolyploid species.

Background:Genetic variation in plants alters insect abundance and community structure in the field; however, little is known about the importance of a single gene among diverse plant genotypes. In this context,Arabidopsistrichomes provide an excellent system to discern the roles of natural variation and a key gene,GLABRA1, in shaping insect communities. In this study, we transplanted two independent glabrous mutants (gl1-1andgl1-2) and 17 natural accessions ofArabidopsis thalianato two localities in Switzerland and Japan.Results:Fifteen insect species inhabited plant accessions, with 10–30% broad-sense heritability of community indices being detected, such as species richness and diversity. The total abundance of leaf-chewing herbivores was negatively correlated with trichome density at both the field sites, while glucosinolates had variable effects on leaf chewers between the two sites. Interestingly, there was a parallel tendency for the abundance of leaf chewers to be higher ongl1-1andgl1-2than for their different parental accessions, Ler-1 and Col-0, respectively. Furthermore, the loss of function in theGLABRA1gene significantly decreased the resistance of plants to the two predominant chewers, flea beetles and turnip sawflies.Conclusions:Overall, our results indicate that insect community composition onA. thalianais heritable across two distant field sites, withGLABRA1playing a key role in altering the abundance of leaf-chewing herbivores. Given that such a trichome variation is widely observed in Brassicaceae plants, the present study exemplifies the community-wide impact of a single plant gene on crucifer-feeding insects in the field.

Mouse APOBEC3 (mA3) inhibits murine leukemia virus (MuLV) replication by a deamination-independent mechanism in which the reverse transcription is considered the main target process. However, other steps in virus replication that can be targeted by mA3 have not been examined. We have investigated the possible effect of mA3 on MuLV protease-mediated processes and found that mA3 binds both mature viral protease and Pr180gag-pol precursor polyprotein. Using replication-competent MuLVs, we also show that mA3 inhibits the processing of Pr65 Gag precursor. Furthermore, we demonstrate that the autoprocessing of Pr180gag-pol is impeded by mA3, resulting in reduced production of mature viral protease. This reduction appears to link with the above inefficient Pr65gag processing in the presence of mA3. Two major isoforms of mA3, exon 5-containing and -lacking ones, equally exhibit this antiviral activity. Importantly, physiologically expressed levels of mA3 impedes both Pr180gag-pol autocatalysis and Pr65gag processing. This blockade is independent of the deaminase activity and requires the C-terminal region of mA3. These results suggest that the above impairment of Pr180gag-pol autoprocessing may significantly contribute to the deaminase-independent antiretroviral activity exerted by mA3., Author summary Soon after the identification of the polynucleotide cytidine deaminase APOBEC3 as a host restriction factor against vif-deficient HIV, it was noticed that deamination-independent mechanisms are involved in the inhibition of viral replication in addition to the deaminase-dependent mechanism. We previously showed that mouse APOBEC3 (mA3) physiologically restricted mouse retrovirus replication in their natural hosts without causing significant G-to-A hypermutations. Inhibition of reverse transcription is reported to be the most plausible mechanism for the deamination-independent antiretroviral function. However, it remains unknown whether the inhibition of reverse transcription is the only way to explain the whole picture of deamination-independent antiviral activity exerted by APOBEC3. Here we show that mA3 targets the autoprocessing of Pr180gag-pol polyprotein. This activity does not require the deaminase catalytic center and mainly exerted by the C-terminal half of mA3. mA3 physically interacts with murine retroviral protease and its precursor Pr180gag-pol. mA3-induced disruption of the autocatalytic Pr180gag-pol cleavage leads to a significant reduction of mature viral protease, resulting in the inhibition of Pr65gag processing to mature Gag proteins. As the Pr180gag-pol autoprocessing is necessary for the maturation of other viral enzymes including the reverse transcriptase, its inhibition by host APOBEC3 may precede the previously described impairment of reverse transcription. Our discovery may lead to the development of novel antiretroviral drugs through the future identification of detailed molecular interfaces between retroviral Gag-Pol polyprotein and APOBEC3.

Recently, increasing attempts have been made to understand how plant genes function in natura. In this context, transcriptional profiles represent plant physiological status in response to environmental stimuli. Herein, we combined high-throughput RNA-Seq with insect survey data on 19 accessions of Arabidopsis thaliana grown at a field site in Switzerland. We found that genes with the gene ontology (GO) annotations of “glucosinolate biosynthetic process” and “response to insects” were most significantly enriched, and the expression of these genes was highly variable among plant accessions. Nearly half of the total expression variation in the glucosinolate biosynthetic genes (AOPs, ESM1, ESP, and TGG1) was explained by among-accession variation. Of these genes, the expression level of AOP3 differed among Col-0 accession individuals depending on the abundance of the mustard aphid (Lipaphis erysimi). We also found that the expression of the major cis-jasmone activated gene CYP81D11 was positively correlated with the number of flea beetles (Phyllotreta striolata and Phyllotreta atra). Combined with the field RNA-Seq data, bioassays confirmed that AOP3 was up-regulated in response to attack by mustard aphids. The combined results from RNA-Seq and our ecological survey illustrate the feasibility of using field transcriptomics to detect an inducible defense, providing a first step towards an in natura understanding of biotic interactions involving phenotypic plasticity.

The self-incompatible species Arabidopsis halleri is a close relative of the self-compatible model plant Arabidopsis thaliana. The broad European and Asian distribution and heavy metal hyperaccumulation ability make A. halleri a useful model for ecological genomics studies. We used long-insert mate-pair libraries to improve the genome assembly of the A. halleri ssp. gemmifera Tada mine genotype (W302) collected from a site with high contamination by heavy metals in Japan. After five rounds of forced selfing, heterozygosity was reduced to 0.04%, which facilitated subsequent genome assembly. Our assembly now covers 196 Mb or 78% of the estimated genome size and achieved scaffold N50 length of 712 kb. To validate assembly and annotation, we used synteny of A. halleri Tada mine with a previously published high-quality reference assembly of a closely related species, Arabidopsis lyrata. Further validation of the assembly quality comes from synteny and phylogenetic analysis of the HEAVY METAL ATPASE4 (HMA4) and METAL TOLERANCE PROTEIN1 (MTP1) regions using published sequences from European A. halleri for comparison. Three tandemly duplicated copies of HMA4, key gene involved in cadmium and zinc hyperaccumulation, were assembled on a single scaffold. The assembly will enhance the genomewide studies of A. halleri as well as the allopolyploid Arabidopsis kamchatica derived from A. lyrata and A. halleri.

The habitats of polyploid species are generally distinct from their parental species. Stebbins described polyploids as 'general purpose genotypes', which can tolerate a wide range of environmental conditions. However, little is known about its molecular basis because of the complexity of polyploid genomes. We hypothesized that allopolyploid species might utilize the expression patterns of both parents depending on environments (polyploid plasticity hypothesis). We focused on hydrological niche segregation along fine-scale soil moisture and waterlogging gradients. Two diploid species, Cardamine amara and Cardamine hirsuta, grew best in submerged and unsubmerged conditions, respectively, consistent with their natural habitats. Interestingly, the allotetraploid Cardamine flexuosa derived from them grew similarly in fluctuating as well as submerged and unsubmerged conditions, consistent with its wide environmental tolerance. A similar pattern was found in another species trio: allotetraploid Cardamine scutata and its parents. Using the close relatedness of Cardamine and Arabidopsis, we quantified genomewide expression patterns following dry and wet treatments using an Arabidopsis microarray. Hierarchical clustering analysis revealed that the expression pattern of C. flexuosa clustered with C. hirsuta in the dry condition and with C. amara in the wet condition, supporting our hypothesis. Furthermore, the induction levels of most genes in the allopolyploid were lower than in a specialist diploid species. This reflects a disadvantage of being allopolyploid arising from fixed heterozygosity. We propose that recurrent allopolyploid speciation along soil moisture and waterlogging gradients confers niche differentiation and reproductive isolation simultaneously and serves as a model for studying the molecular basis of ecological speciation and adaptive radiation.

Whole genome sequences, which can be provided even for non-model organisms owing to high-throughput sequencers, are valuable in enhancing the understanding of adaptive evolution. Metrosideros polymorpha, a tree species endemic to the Hawaiian Islands, occupies a wide range of ecological habitats and shows remarkable polymorphism in phenotypes among/within populations. The biological functions of genetic variations observed within this species could provide significant insights into the adaptive radiation found in a single species. Here de novo assembled genome sequences of M. polymorpha are presented to reveal basic genomic parameters about this species and to develop our knowledge of ecological divergences. The assembly yielded 304-Mbp genome sequences, half of which were covered by 19 scaffolds with >5 Mbp, and contained 30 K protein-coding genes. Demographic history inferred from the genome-wide heterozygosity indicated that this species experienced a dramatic rise and fall in the effective population size, possibly owing to past geographic or climatic changes in the Hawaiian Islands. This M. polymorpha genome assembly represents a high-quality genome resource useful for future functional analyses of both intra- and interspecies genetic variations or comparative genomics.

Allopolyploids possess complete sets of genomes derived from different parental species and exhibit a range of variation in various traits. Reproductive traits may play a key role in the reproductive isolation between allopolyploids and their parental species, thus affecting the thriving of allopolyploids. However, empirical data, especially in natural habitats, comparing reproductive trait variation between allopolyploids and their parental species remain rare. Here, we documented the flowering phenology and floral morphology of the allopolyploid wild plantCardamine flexuosaand its diploid parentsC. amaraandC. hirsutain their native range in Switzerland. The flowering ofC. flexuosastarted at an intermediate time compared with those of the parents and the flowering period ofC. flexuosaoverlapped with those of the parents.Cardamine flexuosaresembledC. hirsutain the size of flowers and petals and the length/width ratio of petals, while it resembledC. amarain the length/width ratio of flowers. These results provide empirical evidence of the trait-dependent variation of allopolyploid phenotypes in natural habitats at the local scale. They also suggest that the variation in some reproductive traits inC. flexuosais associated with self-fertilization. Therefore, it is helpful to consider the mating system in furthering the understanding of the processes that may have shaped trait variation in polyploids in nature.

Genome duplication is widespread in wild and crop plants. However, little is known about genome-wide selection in polyploids due to the complexity of duplicated genomes. In polyploids, the patterns of purifying selection and adaptive substitutions may be affected by masking owing to duplicated genes or homeologs as well as effective population size. Here, we resequence 25 accessions of the allotetraploid Arabidopsis kamchatica, which is derived from the diploid species A. halleri and A. lyrata. We observe a reduction in purifying selection compared with the parental species. Interestingly, proportions of adaptive non-synonymous substitutions are significantly positive in contrast to most plant species. A recurrent pattern observed in both frequency and divergence–diversity neutrality tests is that the genome-wide distributions of both subgenomes are similar, but the correlation between homeologous pairs is low. This may increase the opportunity of different evolutionary trajectories such as in the HMA4 gene involved in heavy metal hyperaccumulation., Nature Communications, 9, ISSN:2041-1723

Although genome duplication is widespread in wild and crop plants, little is known about genome-wide selection due to the complexity of polyploid genomes. In allopolyploid species, the patterns of purifying selection and adaptive substitutions would be affected by masking owing to duplicated genes or ‘homeologs’ as well as by effective population size. We resequenced 25 distribution-wide accessions of the allotetraploid Arabidopsis kamchatica, which has a relatively small genome size (450 Mb) derived from the diploid species A. halleri and A. lyrata. The level of nucleotide polymorphism and linkage disequilibrium decay were comparable to A. thaliana, indicating the feasibility of association studies. A reduction in purifying selection compared with parental species was observed. Interestingly, the proportion of adaptive substitutions (α) was significantly positive in contrast to the majority of plant species. A recurrent pattern observed in both frequency and divergence-based neutrality tests is that the genome-wide distributions of both subgenomes were similar, but the correlation between homeologous pairs was low. This may increase the opportunity of different evolutionary trajectories such as in the HMA4 gene involved in heavy metal hyperaccumulation.

Finger millet (Eleusine coracana (L.) Gaertn) is an important crop for food security because of its tolerance to drought, which is expected to be exacerbated by global climate changes. Nevertheless, it is often classified as an orphan/underutilized crop because of the paucity of scientific attention. Among several small millets, finger millet is considered as an excellent source of essential nutrient elements, such as iron and zinc; hence, it has potential as an alternate coarse cereal. However, high-quality genome sequence data of finger millet are currently not available. One of the major problems encountered in the genome assembly of this species was its polyploidy, which hampers genome assembly compared with a diploid genome. To overcome this problem, we sequenced its genome using diverse technologies with sufficient coverage and assembled it via a novel multiple hybrid assembly workflow that combines next-generation with single-molecule sequencing, followed by whole-genome optical mapping using the Bionano Irys® system. The total number of scaffolds was 1,897 with an N50 length >2.6 Mb and detection of 96% of the universal single-copy orthologs. The majority of the homeologs were assembled separately. This indicates that the proposed workflow is applicable to the assembly of other allotetraploid genomes., DNA Research, 25 (1), ISSN:1340-2838, ISSN:1756-1663

Pitchers are modified leaves used by carnivorous plants for trapping prey. Their fluids contain digestive enzymes from the plant and they harbor abundant microbes. In this study, the diversity of bacterial communities was assessed in Nepenthes pitcher fluids and the composition of the bacterial community was compared to that in other environments, including the phyllosphere of Arabidopsis, animal guts and another pitcher plant, Sarracenia. Diversity was measured by 454 pyrosequencing of 16S rRNA gene amplicons. A total of 232,823 sequences were obtained after chimera and singleton removal that clustered into 3260 distinct operational taxonomic units (OTUs) (3% dissimilarity), which were taxonomically distributed over 17 phyla, 25 classes, 45 orders, 100 families, and 195 genera. Pyrosequencing and fluorescence in situ hybridization yielded similar estimates of community composition. Most pitchers contained high proportions of unique OTUs, and only 22 OTUs (

Background and Aims The coexistence of hermaphrodites and female-sterile individuals, or androdioecy, has been documented in only a handful of plants and animals. This study reports its existence in the plant species Cardamine amara (Brassicaceae), in which female-sterile individuals have shorter pistils than seed-producing hermaphrodites. Methods Morphological analysis, in situ manual pollination, microsatellite genotyping and differential gene expression analysis using Arabidopsis microarrays were used to delimit variation between female-sterile individuals and hermaphrodites. Key Results Female sterility in C. amara appears to be caused by disrupted ovule development. It was associated with a 2.4- to 2.9-fold increase in clonal propagation. This made the pollen number of female-sterile genets more than double that of hermaphrodite genets, which fulfils a condition of co-existence predicted by simple androdioecy theories. When female-sterile individuals were observed in wild androdioecious populations, their ramet frequencies ranged from 5 to 54 %; however, their genet frequencies ranged from 11 to 29 %, which is consistent with the theoretically predicted upper limit of 50 %. Conclusions The results suggest that a combination of sexual reproduction and increased asexual proliferation by female-sterile individuals probably explains the invasion and maintenance of female sterility in otherwise hermaphroditic populations. To our knowledge, this is the first report of the coexistence of female sterility and hermaphrodites in the Brassicaceae.

Beneficial utilization of yuzu (Citrus junos) peel was investigated to improve the commercial value of cultured yellowtail fish. The antioxidant effects of grated yuzu peel from pomace (yuzu paste) on prevention of dark muscle discoloration during storage in sliced yellowtail fish were tested. Four diets were prepared, containing 0, 10, 50, and 100 g yuzu paste in a 1,450-g diet (wet basis). Fish were fed one of the four diets for 10 weeks. Growth performance was not significantly different among the dietary groups, but fish fed a diet containing 100-g yuzu paste showed the least average body weight and a significant increase in fecal amount and fat levels in the feces, which suggest low digestibility and low absorption. Dark muscle discoloration was significantly reduced in all yuzu dietary groups compared to the control group during storage at 4°C. The optimal supplementation of yuzu paste in the fish diet is very effective at maintaining fish flesh freshness as a finishing diet without causing addition...

This article describes the use of cytogenomic and molecular approaches to explore the origin and evolution of Cardamine schulzii, a textbook example of a recent allopolyploid, in its ;110-year history of human-induced hybridization and allopolyploidy in the Swiss Alps. Triploids are typically viewed as bridges between diploids and tetraploids but rarely as parental genomes of high-level hybrids and polyploids. The genome of the triploid semifertile hybrid Cardamine 3 insueta (2n = 24, RRA) was shown to combine the parental genomes of two diploid (2n =2 x = 16) species, Cardamine amara (AA) and Cardamine rivularis (RR). These parental genomes have remained structurally stable within the triploid genome over the >100 years since its origin. Furthermore, we provide compelling evidence that the alleged recent polyploid C. schulzii is not an autohexaploid derivative of C. 3 insueta. Instead, at least two hybridization events involving C. 3 insueta and the hypotetraploid Cardamine pratensis (PPPP, 2n =4 x22 = 30) have resulted in the origin of the trigenomic hypopentaploid (2n =5 x22 = 38, PPRRA) and hypohexaploid (2n =6 x22 = 46, PPPPRA). These data show that the semifertile triploid hybrid can promote a merger of three different genomes and demonstrate how important it is to reexamine the routinely repeated textbook examples using modern techniques.

Allopolyploidization combines parental genomes and often confers broader species distribution. However, little is known about parentally transmitted gene expression underlying quantitative traits following allopolyploidization because of the complexity of polyploid genomes. The allopolyploid species Arabidopsis kamchatica is a natural hybrid of the zinc hyperaccumulator Arabidopsis halleri and of the nonaccumulator Arabidopsis lyrata. We found that A. kamchatica retained the ability to hyperaccumulate zinc from A. halleri and grows in soils with both low and high metal content. Hyperaccumulation of zinc by A. kamchatica was reduced to about half of A. halleri, but is 10-fold greater than A. lyrata. Homeologs derived from A. halleri had significantly higher levels of expression of genes such as HEAVY METAL ATPASE4 (HMA4), METAL TRANSPORTER PROTEIN1 and other metal ion transporters than those derived from A. lyrata, which suggests cis-regulatory differences. A. kamchatica has on average about half the expression of these genes compared with A. halleri due to fixed heterozygosity inherent in allopolyploids. Zinc treatment significantly changed the ratios of expression of 1% of homeologous pairs, including genes putatively involved in metal homeostasis. Resequencing data showed a significant reduction in genetic diversity over a large genomic region (290 kb) surrounding the HMA4 locus derived from the A. halleri parent compared with the syntenic A. lyrata-derived region, which suggests different evolutionary histories. We also estimated that three A. halleri-derived HMA4 copies are present in A. kamchatica. Our findings support a transcriptomic model in which environment-related transcriptional patterns of both parents are conserved but attenuated in the allopolyploids., Molecular Biology and Evolution, 33 (11), ISSN:0737-4038, ISSN:1537-1719

We determined the environmental requirements for Lepus brachyurus lyoni on Sado Island, Japan, during the winter when food resources are low. Hare track survey sites were classified by vegetation type. The movement distances of hares, which indicated habitat preference, were evaluated as the relative length of tracks in a 1-ha area. In addition, we examined browse marks on plants and the protein content of plants to determine the preferred plant species. Movement distances of hares decreased significantly with increasing tree stand height, mean tree diameter, basal area of the tree stand, and forest canopy closure, but increased with branch-and-stem density. The percentage of browse marks on the branches of tree/shrub species decreased significantly with increasing branch-and-stem diameter. The rate of feeding increased significantly with increased protein content in the branches. These findings suggest that hares prefer shrub stands with high branch-and-stem density to mature forests, which hav...

Polyploidization, or genome duplication, has played a critical role in the diversification of animals, fungi and plants. Little is known about the population structure and multiple origins of polyploid species because of the difficulty in identifying multiple homeologous nuclear genes. The allotetraploid species Arabidopsis kamchatica is closely related to the model species Arabidopsis thaliana and is distributed in a broader climatic niche than its parental species. Here, we performed direct sequencing of homeologous pairs of the low-copy nuclear genes WER and CHS by designing homeolog-specific primers, and obtained also chloroplast and ribosomal internal transcribed spacer sequences. Phylogenetic analysis showed that 50 individuals covering the distribution range including North America are allopolyploids derived from Arabidopsis lyrata and Arabidopsis halleri. Three major clusters within A. kamchatica were detected using Bayesian clustering. One cluster has widespread distribution. The other two are restricted to the southern part of the distribution range including Japan, where the parent A. lyrata is not currently distributed. This suggests that the mountains in Central Honshu and surrounding areas in Japan served as refugia during glacial-interglacial cycles and retained this diversity. We also found that multiple haplotypes of nuclear and chloroplast sequences of A. kamchatica are identical to those of their parental species. This indicates that multiple diploid individuals contributed to the origin of A. kamchatica. The haplotypes of low-copy nuclear genes in Japan suggest independent polyploidization events rather than introgression. Our findings suggest that self-compatibility and gene silencing occurred independently in different origins.

We investigated the substrate specificity of farnesyl diphosphate (FPP) synthase derived from Bacillus stearothermophilus and porcine liver by examining the reactivities of two cyclic substrate homologs, cyclohexylideneethyl diphosphate and cyclohexenylethyl diphosphate. Reaction of geranyl diphosphate with 2-cyclohexenylethyl diphosphate using bacterial or porcine liver FPP synthase produced (S)-geranylcyclohexylideneethyl diphosphate, with relative yields of 13.6% for the bacterial enzyme and 42.2% for the porcine liver enzyme. Reaction of cyclohexylideneethyl diphosphate with isopentenyl diphosphate produced 10-cyclohexyliden-3,7-dimethyldeca-2,6-dien-1-ol as a double condensation product, with relative yields of 23.1% (bacterial enzyme) and 3.0% (porcine liver enzyme). Reaction of cyclohexylideneethyl diphosphate with 2-cyclohexenylethyl diphosphate using bacterial enzyme produced (cyclohexylideneethyl)-cyclohexylideneethyl diphosphate (0.8% yield).

The evolution from outcrossing based on self-incompatibility (SI) to a selfing system is one of the most prevalent transitions in flowering plants. It has been suggested that the loss of SI in Arabidopsis thaliana is associated with pseudogene formation at the SCR male component of the S locus. Recent work, however, suggests that alternative alleles with large deletions at the S locus are also present and may be responsible for the evolution of self-compatibility in this species. We demonstrate that most of these deletion alleles are evolutionarily derived from an S haplotype (haplogroups A) that already possessed the SCR pseudogene. This haplotype and its deletion variants are nearly fixed in Europe. Together with previous transgenic data, these results suggest that the pseudogenization of PsiSCR1 gene changed the SI phenotype in the majority of A. thaliana accessions, and was a critical step in the evolution of selfing in this species. Two other haplogroups (B and C) were also identified, the former of which contains a novel and possibly functional SCR allele. In contrast to haplogroups A, these two haplogroups are found primarily in Africa and Asia. These results suggest that self-compatibility, which appears to be fixed in this species, arose multiple times with different genetic bases, and indicates that a species-specific trait is associated with parallel evolution at the molecular level.

Genomewide markers enable us to study genetic differentiation within a species and the factors underlying it at a much higher resolution than before, which advances our understanding of adaptation in organisms. We investigated genomic divergence in Metrosideros polymorpha, a woody species that occupies a wide range of ecological habitats across the Hawaiian Islands and shows remarkable phenotypic variation. Using 1659 single nucleotide polymorphism (SNP) markers annotated with the genome assembly, we examined the population genetic structure and demographic history of nine populations across five elevations and two ages of substrates on Mauna Loa, the island of Hawaii. The nine populations were differentiated into two genetic clusters distributed on the lower and higher elevations and were largely admixed on the middle elevation. Demographic modelling revealed that the two genetic clusters have been maintained in the face of gene flow, and the effective population size of the high-altitude cluster was much smaller. A FST -based outlier search among the 1659 SNPs revealed that 34 SNPs (2.05%) were likely to be under divergent selection and the allele frequencies of 21 of them were associated with environmental changes along elevations, such as temperature and precipitation. This study shows a genomic mosaic of M. polymorpha, in which contrasting divergence patterns were found. While most genomic polymorphisms were shared among populations, a small fraction of the genome was significantly differentiated between populations in diverse environments and could be responsible for the dramatic adaptation to a wide range of environments.

Ralstonia eutropha is a facultative chemolithoautotrophic bacterium that uses the Calvin–Benson–Bassham (CBB) cycle for CO2 fixation. This study showed that R. eutropha strain H16G incorporated 13CO2, emitted by the oxidative decarboxylation of [1-13C1]-glucose, into key metabolites of the CBB cycle and finally into poly(3-hydroxybutyrate) [P(3HB)] with up to 5.6% 13C abundance. The carbon yield of P(3HB) produced from glucose by the strain H16G was 1.2 times higher than that by the CBB cycle-inactivated mutants, in agreement with the possible fixation of CO2 estimated from the balance of energy and reducing equivalents through sugar degradation integrated with the CBB cycle. The results proved that the ‘gratuitously’ functional CBB cycle in R. eutropha under aerobic heterotrophic conditions participated in the reutilization of CO2 emitted during sugar degradation, leading to an advantage expressed as increased carbon yield of the storage compound. This is a new insight into the role of the CBB cycle and may be applicable for more efficient utilization of biomass resources.

Self-incompatibility (SI) systems in flowering plants distinguish self- and non-self pollen to prevent inbreeding. While other SI systems rely on the self-recognition between specific male- and female-determinants, the Solanaceae family has a non-self recognition system resulting in the detoxification of female-determinants of S-ribonucleases (S-RNases), expressed in pistils, by multiple male-determinants of S-locus F-box proteins (SLFs), expressed in pollen. It is not known how many SLF components of this non-self recognition system there are in Solanaceae species, or how they evolved. We identified 16?20 SLFs in each S-haplotype in SI Petunia, from a total of 168 SLF sequences using large-scale next-generation sequencing and genomic polymerase chain reaction (PCR) techniques. We predicted the target S-RNases of SLFs by assuming that a particular S-allele must not have a conserved SLF that recognizes its own S-RNase, and validated these predictions by transformation experiments. A simple mathematical model confirmed that 16?20 SLF sequences would be adequate to recognize the vast majority of target S-RNases. We found evidence of gene conversion events, which we suggest are essential to the constitution of a non-self recognition system and also contribute to self-compatible mutations., [プレスリリース]植物が自分以外の花粉を選ぶための遺伝子セットを特定！ ～生物の非自己認識システムの進化の謎解明へ～, [Author's version title]Gene duplication and genetic exchange drive the evolution of S-RNase based self-incompatibility

To evaluate the effect of low-molecular-weight glutenin subunit (LMW-GS) genes on dough strength, locus-specific primers of LMW-GS genes and gliadin bands tightly linked to LMW-GS genes were analyzed in common wheat. Segregation analysis of the F2 progeny from a cross between ‘Haruhikari’, a good bread-making quality cultivar, and ‘Asakaze-komugi’, a poor bread-making quality cultivar, showed that dough strength significantly correlated with one amplified LMW-GS gene located at the Glu-B3 locus from ‘Haruhikari’. There was no specific reference to the gliadin bands identified as promising markers in the cross under study. The LMW-GS gene of ‘Haruhikari’ had a seven amino-acid deletion in the repetitive domain relative to ‘Asakaze-komugi’, as well as six amino-acid substitutions, three of which would be expected to cause changes in hydrophilicity. The presence of the LMW-GS gene and other LMW-GS genes tightly linked to it may affect the dough strength of wheat.

Finding causal relationships between genotypic and phenotypic variation is a key focus of evolutionary biology, human genetics and plant breeding. To identify genome-wide patterns underlying trait diversity, we assembled a high-quality reference genome of Cardamine hirsuta, a close relative of the model plant Arabidopsis thaliana. We combined comparative genome and transcriptome analyses with the experimental tools available in C. hirsuta to investigate gene function and phenotypic diversification. Our findings highlight the prevalent role of transcription factors and tandem gene duplications in morphological evolution. We identified a specific role for the transcriptional regulators PLETHORA5/7 in shaping leaf diversity and link tandem gene duplication with differential gene expression in the explosive seed pod of C. hirsuta. Our work highlights the value of comparative approaches in genetically tractable species to understand the genetic basis for evolutionary change.

The genus Nepenthes, a carnivorous plant, has a pitcher to trap insects and digest them in the contained fluid to gain nutrient. A distinctive character of the pitcher fluid is the digestive enzyme activity that may be derived from plants and dwelling microbes. However, little is known about in situ digestive enzymes in the fluid. Here we examined the pitcher fluid from four species of Nepenthes. High bacterial density was observed within the fluids, ranging from 7×10(6) to 2.2×10(8) cells ml(-1). We measured the activity of three common enzymes in the fluid: acid phosphatases, β-D-glucosidases, and β-D-glucosaminidases. All the tested enzymes detected in the liquid of all the pitcher species showed activity that considerably exceeded that observed in aquatic environments such as freshwater, seawater, and sediment. Our results indicate that high enzyme activity within a pitcher could assist in the rapid decomposition of prey to maximize efficient nutrient use. In addition, we filtered the fluid to distinguish between dissolved enzyme activity and particle-bound activity. As a result, filtration treatment significantly decreased the activity in all enzymes, while pH value and Nepenthes species did not affect the enzyme activity. It suggested that enzymes bound to bacteria and other organic particles also would significantly contribute to the total enzyme activity of the fluid. Since organic particles are themselves usually colonized by attached and highly active bacteria, it is possible that microbe-derived enzymes also play an important role in nutrient recycling within the fluid and affect the metabolism of the Nepenthes pitcher plant.

Ever since Darwin's pioneering research, the evolution of self-fertilisation (selfing) has been regarded as one of the most prevalent evolutionary transitions in flowering plants. A major mechanism to prevent selfing is the self-incompatibility (SI) recognition system, which consists of male and female specificity genes at the S-locus and SI modifier genes. Under conditions that favour selfing, mutations disabling the male recognition component are predicted to enjoy a relative advantage over those disabling the female component, because male mutations would increase through both pollen and seeds whereas female mutations would increase only through seeds. Despite many studies on the genetic basis of loss of SI in the predominantly selfing plant Arabidopsis thaliana, it remains unknown whether selfing arose through mutations in the female specificity gene (S-receptor kinase, SRK), male specificity gene (S-locus cysteine-rich protein, SCR; also known as S-locus protein 11, SP11) or modifier genes, and whether any of them rose to high frequency across large geographic regions. Here we report that a disruptive 213-base-pair (bp) inversion in the SCR gene (or its derivative haplotypes with deletions encompassing the entire SCR-A and a large portion of SRK-A) is found in 95% of European accessions, which contrasts with the genome-wide pattern of polymorphism in European A. thaliana. Importantly, interspecific crossings using Arabidopsis halleri as a pollen donor reveal that some A. thaliana accessions, including Wei-1, retain the female SI reaction, suggesting that all female components including SRK are still functional. Moreover, when the 213-bp inversion in SCR was inverted and expressed in transgenic Wei-1 plants, the functional SCR restored the SI reaction. The inversion within SCR is the first mutation disrupting SI shown to be nearly fixed in geographically wide samples, and its prevalence is consistent with theoretical predictions regarding the evolutionary advantage of mutations in male components.

Genome duplication with hybridization, or allopolyploidization, occurs commonly in plants, and is considered to be a strong force for generating new species. However, genome-wide quantification of homeolog expression ratios was technically hindered because of the high homology between homeologous gene pairs. To quantify the homeolog expression ratio using RNA-seq obtained from polyploids, a new method named HomeoRoq was developed, in which the genomic origin of sequencing reads was estimated using mismatches between the read and each parental genome. To verify this method, we first assembled the two diploid parental genomes of Arabidopsis halleri subsp. gemmifera and Arabidopsis lyrata subsp. petraea (Arabidopsis petraea subsp. umbrosa), then generated a synthetic allotetraploid, mimicking the natural allopolyploid Arabidopsis kamchatica. The quantified ratios corresponded well to those obtained by Pyrosequencing. We found that the ratios of homeologs before and after cold stress treatment were highly correlated (r = 0.870). This highlights the presence of nonstochastic polyploid gene regulation despite previous research identifying stochastic variation in expression. Moreover, our new statistical test incorporating overdispersion identified 226 homeologs (1.11% of 20 369 expressed homeologs) with significant ratio changes, many of which were related to stress responses. HomeoRoq would contribute to the study of the genes responsible for polyploid-specific environmental responses.

Background Ralstonia eutropha H16 is well known to produce polyhydroxyalkanoates (PHAs), which are potential bio-based biodegradable plastics, in an efficient manner as an energy storage material under unbalanced growth conditions. To obtain further knowledge of PHA biosynthesis, this study performed a quantitative transcriptome analysis based on deep sequencing of the complementary DNA generated from the RNA (RNA-seq) of R. eutropha H16. Results Total RNAs were extracted from R. eutropha cells in growth, PHA production, and stationary phases on fructose. rRNAs in the preparation were removed by repeated treatments with magnetic beads specific to bacterial rRNAs, and then the 36 bp sequences were determined using an Illumina high-throughput sequencer. The RNA-seq results indicated the induction of gene expression for transcription, translation, cell division, peptidoglycan biosynthesis, pilus and flagella assembly, energy conservation, and fatty acid biosynthesis in the growth phase; and the repression trends of genes involved in central metabolisms in the PHA production phase. Interestingly, the transcription of genes for Calvin-Benson-Bassham (CBB) cycle and several genes for β-oxidation were significantly induced in the PHA production phase even when the cells were grown on fructose. Moreover, incorporation of 13C was observed in poly(3-hydroxybutyrate) synthesized by R. eutropha H16 from fructose in the presence of NaH13CO3, and further gene deletion analyses revealed that both of the two ribulose 1,5-bisphosphate carboxylase (Rubiscos) in CBB cycle were actually functional in CO2 fixation under the heterotrophic condition. Conclusions The results revealed the phase-dependent transcriptomic changes and a CO2 fixation capability under heterotrophic conditions by PHA-producing R. eutropha.

1,2,3,4-Tetra-t-butylnaphthalene (2), a highly crowded naphthalene derivative, has been synthesized. Upon irradiation, 2 underwent photovalence isomerization to afford 1,2,5,6-tetra-t-butyl-3,4-benzo-3-tricyclo[3.1.0.02,6]-hexene (3). The synthesis of 1,4,5,6-tetra-t-butyl-2,3-benzo-2,5-bicyclo[2.2.0]hexadiene (1), which was used for the synthesis of 2, is also described.

Two forms of erythropoietin, EPO-bi and EPO-tetra, with different biological activities were isolated from the culture medium of a recombinant Chinese hamster ovary cell line, B8-300, into which the human erythropoietin gene had been introduced. EPO-bi, an unusual form, showed only one-seventh the in vivo activity and 3 times higher in vitro activity of the previously described recombinant human EPO (standard EPO). In contrast, EPO-tetra showed both in vivo and in vitro activities comparable to those of the standard EPO. EPO-bi, EPO-tetra, and the standard EPO had the same amino acid composition and immunoreactivity. However, structural analyses of their N-linked sugar chains revealed that EPO-bi contains the biantennary complex type as the major sugar chain, while EPO-tetra and the standard EPO contain the tetraantennary complex type as the major sugar chain. From examination of various preparations of recombinant human EPO, we found a positive correlation between the in vivo activity of EPO and the ratio of tetraantennary to biantennary oligosaccharides. These results suggest that higher branching of the N-linked sugar chains is essential for effective expression of in vivo biological activity of EPO.