Seed dormancy 4 like1 of Arabidopsis is a key regulator of phase transition from embryo to vegetative development.

Seed dormancy is an adaptive trait that enables plants to survive adverse conditions and restart growth in a season and location suitable for vegetative and reproductive growth. Control of seed dormancy is also important for crop production and food quality because it can help induce uniform germination and prevent preharvest sprouting. Rice preharvest sprouting quantitative trait locus analysis has identified Seed dormancy 4 (Sdr4) as a positive regulator of dormancy development. Here, we analyzed the loss-of-function mutant of the Arabidopsis ortholog, Sdr4 Like1 (SFL1), and found that the sfl1-1 seeds showed precocious germination at the mid- to late-maturation stage similar to rice sdr4 mutant, but converted to become more dormant than the wild type during maturation drying. Coordinated with the dormancy levels, expression levels of the seed maturation and dormancy master regulator genes, ABI3, FUS3, and DOG1 in sfl1-1 seeds were lower than in wild type at early- and mid-maturation stages, but higher at the late-maturation stage. In addition to the seed dormancy phenotype, sfl1-1 seedlings showed a growth arrest phenotype and heterochronic expression of LAFL (LEC1, ABI3, FUS3, LEC2) and DOG1 in the seedlings. These data suggest that SFL1 is a positive regulator of initiation and termination of the seed dormancy program. We also found genetic interaction between SFL1 and the SFL2, SFL3, and SFL4 paralogs of SFL1, which impacts on the timing of the phase transition from embryo maturation to seedling growth.

Root angle modifications by the DRO1 homolog improve rice yields in saline paddy fields.

The root system architecture (RSA) of crops can affect their production, particularly in abiotic stress conditions, such as with drought, waterlogging, and salinity. Salinity is a growing problem worldwide that negatively impacts on crop productivity, and it is believed that yields could be improved if RSAs that enabled plants to avoid saline conditions were identified. Here, we have demonstrated, through the cloning and characterization of qSOR1 (quantitative trait locus for SOIL SURFACE ROOTING 1), that a shallower root growth angle (RGA) could enhance rice yields in saline paddies. qSOR1 is negatively regulated by auxin, predominantly expressed in root columella cells, and involved in the gravitropic responses of roots. qSOR1 was found to be a homolog of DRO1 (DEEPER ROOTING 1), which is known to control RGA. CRISPR-Cas9 assays revealed that other DRO1 homologs were also involved in RGA. Introgression lines with combinations of gain-of-function and loss-of-function alleles in qSOR1 and DRO1 demonstrated four different RSAs (ultra-shallow, shallow, intermediate, and deep rooting), suggesting that natural alleles of the DRO1 homologs could be utilized to control RSA variations in rice. In saline paddies, near-isogenic lines carrying the qSOR1 loss-of-function allele had soil-surface roots (SOR) that enabled rice to avoid the reducing stresses of saline soils, resulting in increased yields compared to the parental cultivars without SOR. Our findings suggest that DRO1 homologs are valuable targets for RSA breeding and could lead to improved rice production in environments characterized by abiotic stress.

Mutations within the miR172 target site of wheat AP2 homoeologs regulate lodicule size and rachis internode length.

Closed fertilization in flowers, or cleistogamy, reduces the risk of fungal infection in Triticeae crops. In barley (Hordeum vulgare), cleistogamy is determined by a single recessive gene, cly1, which results from a single nucleotide polymorphism within the microRNA172 target site of the Apetala2 (AP2) transcription factor gene. The recessive cly1 allele negatively regulates the development of lodicules, keeping florets closed at anthesis. However, cleistogamy is not evident in hexaploid wheat (Triticum aestivum) cultivars. This study aimed at identifying mutations in wheat AP2 orthologs by ethyl methane sulfonate-induced mutagenesis and high-resolution melt analysis. Although flowers of AP2 mutants induced in the A and D genomes opened at anthesis, their lodicule size was significantly smaller, especially in the direction of depth, than that of wild-type plants. One of the mutants that carried a nucleotide replacement in AP2 from the D genome produced a compact spike caused by a substantial decrease in rachis internode length, analogous to the barley dense spike. Cleistogamous hexaploid wheat might be generated by combining effective mutant alleles of AP2-homoeologous genes.

Development of 12 sets of chromosome segment substitution lines that enhance allele mining in Asian cultivated rice.

Many agronomic traits that are important in rice breeding are controlled by multiple genes. The extensive time and effort devoted so far to identifying and selecting such genes are still not enough to target multiple agronomic traits in practical breeding in Japan because of a lack of suitable plant materials in which to efficiently detect and validate beneficial alleles from diverse genetic resources. To facilitate the comprehensive analysis of genetic variation in agronomic traits among Asian cultivated rice, we developed 12 sets of chromosome segment substitution lines (CSSLs) with the japonica background, 11 of them in the same genetic background, using donors representing the genetic diversity of Asian cultivated rice. Using these materials, we overviewed the chromosomal locations of 1079 putative QTLs for seven agronomic traits and their allelic distribution in Asian cultivated rice through multiple linear regression analysis. The CSSLs will allow the effects of putative QTLs in the highly homogeneous japonica background to be validated.

Deficiency in alcohol dehydrogenase 2 reduces arsenic in rice grains by suppressing silicate transporters.

Paddy fields are anaerobic and facilitate arsenite (As(III)) elution from the soil. Paddy-field rice accumulates arsenic (As) in its grains because silicate transporters actively assimilate As(III) during the reproductive stage. Reducing the As level in rice grains is an important challenge for agriculture. Using a forward genetic approach, we isolated a rice (Oryza sativa) mutant, low arsenic line 3 (las3), whose As levels were decreased in aerial tissues, including grains. The low-As phenotype was not observed in young plants before heading (emergence of the panicle). Genetic analyses revealed that a deficiency in alcohol dehydrogenase (ADH) 2 by mutation is responsible for the phenotype. Among the three rice ADH paralogues, ADH2 was the most efficiently produced in root tissue under anaerobic conditions. In wild-type (WT), silicon and As concentrations in aerial tissues increased with growth. However, the increase was suppressed in las3 during the reproductive stage. Accordingly, the gene expression of two silicate transporters, Lsi1 and Lsi2, was increased in WT around the time of heading, whereas the increase was suppressed in las3. These results indicate that the low-As phenotype in las3 is due to silicate transporter suppression. Measurement of intracellular pH by 31P-nuclear magnetic resonance revealed intracellular acidification of las3 roots under hypoxia, suggesting that silicate transporter suppression in las3 might arise from an intracellular pH decrease, which is known to be facilitated by a deficiency in ADH activity under anaerobic conditions. This study provides valuable insight into reducing As levels in rice grains.

A weak allele of OsNRAMP5 confers moderate cadmium uptake while avoiding manganese deficiency in rice.

Decreasing cadmium (Cd) concentrations in rice grains can effectively reduce potential risks to human health because rice is the major contributor to Cd intake in many diets. Among several genes involved in rice Cd accumulation, the loss of function of OsNRAMP5 is known to be effective in reducing grain concentration by inhibiting root uptake. However, disruption of this gene simultaneously decreases manganese (Mn) uptake because OsNRAMP5 is a major Mn transporter. With the aim of improving Mn uptake in OsNRAMP5 mutants while still restricting the grain Cd concentration below the upper limit of international standards, we identified a novel OsNRAMP5 allele encoding a protein in which glutamine (Q) at position 337 was replaced by lysine (K). The mutant carrying the OsNRAMP5-Q337K allele showed intermediate Cd and Mn accumulation between that of the wild-type and OsNRAMP5-knockout lines, and exhibited more resistance to Mn deficiency than the knockout lines. Different amino acid substitutions at position Q337 significantly affected the Cd and Mn transport activity in yeast cells, indicating that it is one of the crucial sites for OsNRAMP5 function. Our results suggest that the OsNRAMP5-Q337K allele might be useful for reducing grain Cd concentrations without causing severe Mn deficiency in rice cultivars through DNA marker-assisted breeding.

Unleashing floret fertility in wheat through the mutation of a homeobox gene.

Floret fertility is a key determinant of the number of grains per inflorescence in cereals. During the evolution of wheat (Triticum sp.), floret fertility has increased, such that current bread wheat (Triticum aestivum) cultivars set three to five grains per spikelet. However, little is known regarding the genetic basis of floret fertility. The locus Grain Number Increase 1 (GNI1) is shown here to be an important contributor to floret fertility. GNI1 evolved in the Triticeae through gene duplication. The gene, which encodes a homeodomain leucine zipper class I (HD-Zip I) transcription factor, was expressed most abundantly in the most apical floret primordia and in parts of the rachilla, suggesting that it acts to inhibit rachilla growth and development. The level of GNI1 expression has decreased over the course of wheat evolution under domestication, leading to the production of spikes bearing more fertile florets and setting more grains per spikelet. Genetic analysis has revealed that the reduced-function allele GNI-A1 contributes to the increased number of fertile florets per spikelet. The RNAi-based knockdown of GNI1 led to an increase in the number of both fertile florets and grains in hexaploid wheat. Mutants carrying an impaired GNI-A1 allele out-yielded WT allele carriers under field conditions. The data show that gene duplication generated evolutionary novelty affecting floret fertility while mutations favoring increased grain production have been under selection during wheat evolution under domestication.

Investigation of the Genetic Diversity of a Rice Core Collection of Japanese Landraces using Whole-Genome Sequencing.

The Rice Core Collection of Japanese Landraces (JRC) consisting of 50 accessions was developed by the genebank at the National Agriculture and Food Research Organization (NARO) in 2008. As a Japanese landrace core collection, the JRC has been used for many research projects, including screening for different phenotypes and allele mining for target genes. To understand the genetic diversity of Japanese Landraces, we performed whole-genome resequencing of these 50 accessions and obtained a total of 2,145,095 single nucleotide polymorphism (SNPs) and 317,832 insertion-deletions (indels) by mapping against the Oryza sativa ssp. japonica Nipponbare genome. A JRC phylogenetic tree based on 1,394 representative SNPs showed that JRC accessions were divided into two major groups and one small group. We used the multiple genome browser, TASUKE+, to examine the haplotypes of flowering genes and detected new mutations in these genes. Finally, we performed genome-wide association studies (GWAS) for agronomical traits using the JRC and another core collection, the World Rice Core Collection (WRC), comprising 69 accessions also provided by the NARO genebank. In leaf blade width, a strong peak close to NAL1, a key gene for the regulation of leaf width, and, in heading date, a peak near HESO1 involved in flowering regulation were observed in GWAS using the JRC. They were also detected in GWAS using the combined JRC + WRC. Thus, JRC and JRC + WRC are suitable populations for GWAS of particular traits.

NDH-PSI Supercomplex Assembly Precedes Full Assembly of the NDH Complex in Chloroplast.

The chloroplast NADH dehydrogenase-like (NDH) complex is structurally similar to respiratory complex I and mediates PSI cyclic electron flow. In Arabidopsis (Arabidopsis thaliana), chloroplast NDH is composed of at least 29 subunits and associates with two copies of PSI to form the NDH-PSI supercomplex. Here, we found that CHLORORESPIRATORY REDUCTION3 (CRR3) is an assembly factor required for the accumulation of subcomplex B (SubB) of chloroplast NDH. In Suc density gradient centrifugation, CRR3 was detected in three protein complexes. Accumulation of the largest peak III complex was impaired in mutants defective in the SubB subunits PnsB2-PnsB5. The oligomeric form of CRR3 likely functions to assemble the core of SubB to form the peak III complex as an assembly intermediate. A defect in the PnsL3 subunit increased the level of the peak III complex, suggesting that CRR3 was released from the assembly intermediate after PnsL3 binding. Unlike PnsB2-PnsB5 and PnsL3, PnsB1 was not absolutely necessary for stabilizing SubB. PnsB1 is likely incorporated into the intermediate at the final step during SubB assembly. Lhca6 is a linker protein mediating NDH-PSI supercomplex formation, and its site of contact with NDH was suggested to be SubB. In the lhca6 mutant, accumulation of the peak III complex was impaired, suggesting that SubB interacted with Lhca6 during the step of SubB assembly. The process of supercomplex formation was triggered before the completion of the NDH assembly. Consistent with its predicted function, CRR3 accumulated in young leaves, where the NDH complex was assembled.

N-Glycoproteomic Characterization of Mannosidase and Xylosyltransferase Mutant Strains of Chlamydomonasreinhardtii.

At present, only little is known about the enzymatic machinery required for N-glycosylation in Chlamydomonas reinhardtii, leading to the formation of N-glycans harboring Xyl and methylated Man. This machinery possesses new enzymatic features, as C. reinhardtii N-glycans are independent of ß1,2-N-acetylglucosaminyltransferase I. Here we have performed comparative N-glycoproteomic analyses of insertional mutants of mannosidase 1A (IM Man1A ) and xylosyltransferase 1A (IM XylT1A ). The disruption of man1A affected methylation of Man and the addition of terminal Xyl. The absence of XylT1A led to shorter N-glycans compared to the wild type. The use of a IM Man1A xIM XylT1A double mutant revealed that the absence of Man1A suppressed the IM XylT1A phenotype, indicating that the increased N-glycan trimming is regulated by core ß1,2-Xyl and is dependent on Man1A activity. These data point toward an enzymatic cascade in the N-glycosylation pathway of C. reinhardtii with interlinked roles of Man1A and XylT1A. The results described herein represent the first step toward a functional characterization of the enzymatic N-glycosylation machinery in C. reinhardtii.

Two patients with MIRAGE syndrome lacking haematological features: role of somatic second-site reversion SAMD9 mutations.

BACKGROUND: Myelodysplasia, infection, restriction of growth, adrenal hypoplasia, genital phenotypes and enteropathy (MIRAGE) syndrome is a recently described congenital disorder caused by heterozygous SAMD9 mutations. The phenotypic spectrum of the syndrome remains to be elucidated. METHODS AND RESULTS: We describe two unrelated patients who showed manifestations compatible with MIRAGE syndrome, with the exception of haematological features. Leucocyte genomic DNA samples were analysed with next-generation sequencing and Sanger sequencing, revealing the patients to have two de novoSAMD9 mutations on the same allele (patient 1 p.[Gln695*; Ala722Glu] and patient 2 p.[Gln39*; Asp769Gly]). In patient 1, p.Gln695* was absent in genomic DNA extracted from hair follicles, implying that the non-sense mutation was acquired somatically. In patient 2, with the 46,XX karyotype, skewed X chromosome inactivation pattern was found in leucocyte DNA, suggesting monoclonality of cells in the haematopoietic system. In vitro expression experiments confirmed the growth-restricting capacity of the two missense mutant SAMD9 proteins that is a characteristic of MIRAGE-associated SAMD9 mutations. CONCLUSIONS: Acquisition of a somatic nonsense SAMD9 mutation in the cells of the haematopoietic system might revert the cellular growth repression caused by the germline SAMD9 mutations (ie, second-site reversion mutations). Unexpected lack of haematological features in the two patients would be explained by the reversion mutations.

Genetic dissection of pre-harvest sprouting resistance in an upland rice cultivar.

Seed dormancy is important in rice breeding because it confers resistance to pre-harvest sprouting (PHS). To detect quantitative trait loci (QTLs) for pre-harvest sprouting resistance, we used chromosome segment substitution lines (CSSLs) derived from a cross between the Japanese upland rice cultivar 'Owarihatamochi' and the lowland rice cultivar 'Koshihikari'. In the CSSLs, several chromosomal regions were associated with PHS resistance. Among these, the chromosome 9 segment from 'Owarihatamochi' had the greatest association with increased PHS resistance. Further QTL analysis using an advanced backcross population (BC4F2) derived from a 'Koshihikari' × 'Owarihatamochi' cross revealed two putative QTLs, here designated qSDR9.1 (Seed dormancy 9.1) and qSDR9.2, on chromosome 9. The 'Owarihatamochi' alleles of the two QTLs reduced germination. Further fine mapping revealed that qSDR9.1 and qSDR9.2 were located within 4.1-Mb and 2.3-Mb intervals (based on the 'Nipponbare' reference genome sequence) defined by the simple sequence repeat marker loci RM24039 and RM24260 and Indel_2 and RM24540, respectively. We thus identified two QTLs for PHS resistance in 'Owarihatamochi', even though resistance levels are relatively low in this cultivar. This unexpected finding suggests the advantages of using CSSLs for QTL detection.

A Defect in DNA Ligase4 Enhances the Frequency of TALEN-Mediated Targeted Mutagenesis in Rice.

We have established methods for site-directed mutagenesis via transcription activator-like effector nucleases (TALENs) in the endogenous rice (Oryza sativa) waxy gene and demonstrated stable inheritance of TALEN-induced somatic mutations to the progeny. To analyze the role of classical nonhomologous end joining (cNHEJ) and alternative nonhomologous end joining (altNHEJ) pathways in TALEN-induced mutagenesis in plant cells, we investigated whether a lack of DNA Ligase4 (Lig4) affects the kinetics of TALEN-induced double-strand break repair in rice cells. Deep-sequencing analysis revealed that the frequency of all types of mutations, namely deletion, insertion, combination of insertion with deletion, and substitution, in lig4 null mutant calli was higher than that in a lig4 heterozygous mutant or the wild type. In addition, the ratio of large deletions (greater than 10 bp) and deletions repaired by microhomology-mediated end joining (MMEJ) to total deletion mutations in lig4 null mutant calli was higher than that in the lig4 heterozygous mutant or wild type. Furthermore, almost all insertions (2 bp or greater) were shown to be processed via copy and paste of one or more regions around the TALENs cleavage site and rejoined via MMEJ regardless of genetic background. Taken together, our findings indicate that the dysfunction of cNHEJ leads to a shift in the repair pathway from cNHEJ to altNHEJ or synthesis-dependent strand annealing.

CHLORORESPIRATORY REDUCTION 9 is a Novel Factor Required for Formation of Subcomplex A of the Chloroplast NADH Dehydrogenase-Like Complex.

In vascular plants, the chloroplast NADH dehydrogenase-like (NDH) complex, a homolog of respiratory NADH:quinone oxidoreductase (Complex I), mediates plastoquinone reduction using ferredoxin as an electron donor in cyclic electron transport around PSI in the thylakoid membrane. In angiosperms, chloroplast NDH is composed of five subcomplexes and forms a supercomplex with PSI. The modular assembly of stroma-protruded subcomplex A, which corresponds to the Q module of Complex I, was recently reported. However, the factors involved in the specific assembly steps have not been completely identified. Here, we isolated an Arabidopsis mutant, chlororespiratory reduction 9 (crr9), defective in NDH activity. The CRR9 gene encodes a novel stromal protein without any known functional domains or motifs. CRR9 is highly conserved in cyanobacteria and land plants but not in green algae, which do not have chloroplast NDH. Blue native-PAGE and immunoblot analyses of thylakoid proteins indicated that formation of subcomplex A was impaired in crr9 CRR9 was specifically required for the accumulation of NdhK, a subcomplex A subunit, in NDH assembly intermediates in the stroma. Furthermore, two-dimensional clear native/SDS-PAGE analysis of the stroma fraction indicated that incorporation of NdhM into NDH assembly intermediate complex 400 was impaired in crr9 These results suggest that CRR9 is a novel factor required for the formation of NDH subcomplex A.

RICE SALT SENSITIVE3 forms a ternary complex with JAZ and class-C bHLH factors and regulates jasmonate-induced gene expression and root cell elongation.

Plasticity of root growth in response to environmental cues and stresses is a fundamental characteristic of land plants. However, the molecular basis underlying the regulation of root growth under stressful conditions is poorly understood. Here, we report that a rice nuclear factor, RICE SALT SENSITIVE3 (RSS3), regulates root cell elongation during adaptation to salinity. Loss of function of RSS3 only moderately inhibits cell elongation under normal conditions, but it provokes spontaneous root cell swelling, accompanied by severe root growth inhibition, under saline conditions. RSS3 is preferentially expressed in the root tip and forms a ternary complex with class-C basic helix-loop-helix (bHLH) transcription factors and JASMONATE ZIM-DOMAIN proteins, the latter of which are the key regulators of jasmonate (JA) signaling. The mutated protein arising from the rss3 allele fails to interact with bHLH factors, and the expression of a significant portion of JA-responsive genes is upregulated in rss3. These results, together with the known roles of JAs in root growth regulation, suggest that RSS3 modulates the expression of JA-responsive genes and plays a crucial role in a mechanism that sustains root cell elongation at appropriate rates under stressful conditions.

Genetic architecture of variation in heading date among Asian rice accessions.

BACKGROUND: Heading date, a crucial factor determining regional and seasonal adaptation in rice (Oryza sativa L.), has been a major selection target in breeding programs. Although considerable progress has been made in our understanding of the molecular regulation of heading date in rice during last two decades, the previously isolated genes and identified quantitative trait loci (QTLs) cannot fully explain the natural variation for heading date in diverse rice accessions. RESULTS: To genetically dissect naturally occurring variation in rice heading date, we collected QTLs in advanced-backcross populations derived from multiple crosses of the japonica rice accession Koshihikari (as a common parental line) with 11 diverse rice accessions (5 indica, 3 aus, and 3 japonica) that originate from various regions of Asia. QTL analyses of over 14,000 backcrossed individuals revealed 255 QTLs distributed widely across the rice genome. Among the detected QTLs, 128 QTLs corresponded to genomic positions of heading date genes identified by previous studies, such as Hd1, Hd6, Hd3a, Ghd7, DTH8, and RFT1. The other 127 QTLs were detected in different chromosomal regions than heading date genes. CONCLUSIONS: Our results indicate that advanced-backcross progeny allowed us to detect and confirm QTLs with relatively small additive effects, and the natural variation in rice heading date could result from combinations of large- and small-effect QTLs. We also found differences in the genetic architecture of heading date (flowering time) among maize, Arabidopsis, and rice.

Chromosomal locations of a gene underlying heat-accelerated brown spot formation and its suppressor genes in rice.

Brown spots on mature leaves from the heading to ripening stages in rice are considered to be lesions induced by heat stress. However, there are few studies of lesions that are induced by heat stress rather than by pathogen infections. To understand the genetic background underlying such lesions, we used the chromosome segment substitution line (CSSL) SL518, derived from a distant cross between rice cultivars Koshihikari (japonica) and Nona Bokra (indica). We observed brown spots on mature leaf blades of the CSSL, although the parents barely showed any spots. Spot formation in SL518 was accelerated by high temperature, suggesting that the candidate gene for spot formation is related to heat stress response. Using progeny derived from a cross between SL518 and Koshihikari, we mapped the causative gene, BROWN-SPOTTED LEAF 1 (BSPL1), on chromosome 5. We speculated that one or more Nona Bokra genes suppress spot formation caused by BSPL1 and identified candidate genomic regions on chromosomes 2 and 9 using a cross between a near-isogenic line for BSPL1 and other CSSLs possessing Nona Bokra segments in the Koshihikari genetic background. In conclusion, our data support the concept that multiple genes are complementarily involved in brown spot formation induced by heat stress and will be useful for cloning of the novel gene(s) related to the spot formation.

RiceXPro version 3.0: expanding the informatics resource for rice transcriptome.

A wide range of resources on gene expression profiling enhance various strategies in plant molecular biology particularly in characterization of gene function. We have updated our gene expression profile database, RiceXPro (http://ricexpro.dna.affrc.go.jp/), to provide more comprehensive information on the transcriptome of rice encompassing the entire growth cycle and various experimental conditions. The gene expression profiles are currently grouped into three categories, namely, 'field/development' with 572 data corresponding to 12 data sets, 'plant hormone' with 143 data corresponding to 13 data sets and 'cell- and tissue-type' comprising of 38 microarray data. In addition to the interface for retrieving expression information of a gene/genes in each data set, we have incorporated an interface for a global approach in searching an overall view of the gene expression profiles from multiple data sets within each category. Furthermore, we have also added a BLAST search function that enables users to explore expression profile of a gene/genes with similarity to a query sequence. Therefore, the updated version of RiceXPro can be used more efficiently to survey the gene expression signature of rice in sufficient depth and may also provide clues on gene function of other cereal crops.

RiceFREND: a platform for retrieving coexpressed gene networks in rice.

Similarity of gene expression across a wide range of biological conditions can be efficiently used in characterization of gene function. We have constructed a rice gene coexpression database, RiceFREND (http://ricefrend.dna.affrc.go.jp/), to identify gene modules with similar expression profiles and provide a platform for more accurate prediction of gene functions. Coexpression analysis of 27 201 genes was performed against 815 microarray data derived from expression profiling of various organs and tissues at different developmental stages, mature organs throughout the growth from transplanting until harvesting in the field and plant hormone treatment conditions, using a single microarray platform. The database is provided with two search options, namely, 'single guide gene search' and 'multiple guide gene search' to efficiently retrieve information on coexpressed genes. A user-friendly web interface facilitates visualization and interpretation of gene coexpression networks in HyperTree, Cytoscape Web and Graphviz formats. In addition, analysis tools for identification of enriched Gene Ontology terms and cis-elements provide clue for better prediction of biological functions associated with the coexpressed genes. These features allow users to clarify gene functions and gene regulatory networks that could lead to a more thorough understanding of many complex agronomic traits.

Composition and physiological function of the chloroplast NADH dehydrogenase-like complex in Marchantia polymorpha.

The chloroplast NADH dehydrogenase-like (NDH) complex mediates cyclic electron transport and chloro-respiration and consists of five sub-omplexes, which in angiosperms further associate with photosystem I (PSI) to form a super-complex. In Marchantia polymorpha, 11 plastid-encoded subunits and all the nuclear-encoded subunits of the A, B, membrane and ferredoxin-binding sub-complexes are conserved. However, it is unlikely that the genome of this liverwort encodes Lhca5 and Lhca6, both of which mediate NDH-PSI super-complex formation. It is also unlikely that the subunits of the lumen sub-complex, PnsL1-L4, are encoded by the genome. Consistent with this in silico prediction, the results of blue-native gel electrophoresis showed that NDH subunits were detected in a protein complex with lower molecular mass in Marchantia than the NDH-PSI super-complex in Arabidopsis. Using the plastid transformation technique, we knocked out the ndhB gene in Marchantia. Although the wild-type genome copies were completely segregated out, the ΔndhB lines grew like the wild-type photoautotrophically. A post-illumination transient increase in chlorophyll fluorescence, which reflects NDH activity in vivo in angiosperms, was absent in the thalli of the ΔndhB lines. In ruptured chloroplasts, antimycin A-insensitive, and ferredoxin-dependent plastoquinone reduction was impaired, suggesting that chloroplast NDH mediates similar electron transport in Marchantia and Arabidopsis, despite its possible difference in structure. As in angiosperms, linear electron transport was not strongly affected in the ΔndhB lines. However, the plastoquinone pool was slightly more reduced at low light intensity, suggesting that chloroplast NDH functions in redox balancing of the inter system, especially under low light conditions.