The WIP6 transcription factor TOO MANY LATERALS specifies vein type in C4 and C3 grass leaves.

Grass leaves are invariantly strap shaped with an elongated distal blade and a proximal sheath that wraps around the stem. Underpinning this shape is a scaffold of leaf veins, most of which extend in parallel along the proximo-distal leaf axis. Differences between species are apparent both in the vein types that develop and in the distance between veins across the medio-lateral leaf axis. A prominent engineering goal is to increase vein density in leaves of C3 photosynthesizing species to facilitate the introduction of the more efficient C4 pathway. Here, we discover that the WIP6 transcription factor TOO MANY LATERALS (TML) specifies vein rank in both maize (C4) and rice (C3). Loss-of-function tml mutations cause large lateral veins to develop in positions normally occupied by smaller intermediate veins, and TML transcript localization in wild-type leaves is consistent with a role in suppressing lateral vein development in procambial cells that form intermediate veins. Attempts to manipulate TML function in rice were unsuccessful because transgene expression was silenced, suggesting that precise TML expression is essential for shoot viability. This finding may reflect the need to prevent the inappropriate activation of downstream targets or, given that transcriptome analysis revealed altered cytokinin and auxin signaling profiles in maize tml mutants, the need to prevent local or general hormonal imbalances. Importantly, rice tml mutants display an increased occupancy of veins in the leaf, providing a step toward an anatomical chassis for C4 engineering. Collectively, a conserved mechanism of vein rank specification in grass leaves has been revealed.

Genetic Mapping by Sequencing More Precisely Detects Loci Responsible for Anaerobic Germination Tolerance in Rice.

Direct seeded rice (DSR) is a mainstay for planting rice in the Americas, and it is rapidly becoming more popular in Asia. It is essential to develop rice varieties that are suitable for this type of production system. ASD1, a landrace from India, possesses several traits desirable for direct-seeded fields, including tolerance to anaerobic germination (AG). To map the genetic basis of its tolerance, we examined a population of 200 F2:3 families derived from a cross between IR64 and ASD1 using the restriction site-associated DNA sequencing (RAD-seq) technology. This genotyping platform enabled the identification of 1921 single nucleotide polymorphism (SNP) markers to construct a high-resolution genetic linkage map with an average interval of 0.9 cM. Two significant quantitative trait loci (QTLs) were detected on chromosomes 7 and 9, qAG7 and qAG9, with LOD scores of 7.1 and 15.0 and R2 values of 15.1 and 29.4, respectively. Here, we obtained more precise locations of the QTLs than traditional simple sequence repeat and low-density SNP genotyping methods and may help further dissect the genetic factors of these QTLs.

RPW8/HR repeats control NLR activation in Arabidopsis thaliana.

In many plant species, conflicts between divergent elements of the immune system, especially nucleotide-binding oligomerization domain-like receptors (NLR), can lead to hybrid necrosis. Here, we report deleterious allele-specific interactions between an NLR and a non-NLR gene cluster, resulting in not one, but multiple hybrid necrosis cases in Arabidopsis thaliana. The NLR cluster is RESISTANCE TO PERONOSPORA PARASITICA 7 (RPP7), which can confer strain-specific resistance to oomycetes. The non-NLR cluster is RESISTANCE TO POWDERY MILDEW 8 (RPW8) / HOMOLOG OF RPW8 (HR), which can confer broad-spectrum resistance to both fungi and oomycetes. RPW8/HR proteins contain at the N-terminus a potential transmembrane domain, followed by a specific coiled-coil (CC) domain that is similar to a domain found in pore-forming toxins MLKL and HET-S from mammals and fungi. C-terminal to the CC domain is a variable number of 21- or 14-amino acid repeats, reminiscent of regulatory 21-amino acid repeats in fungal HET-S. The number of repeats in different RPW8/HR proteins along with the sequence of a short C-terminal tail predicts their ability to activate immunity in combination with specific RPP7 partners. Whether a larger or smaller number of repeats is more dangerous depends on the specific RPW8/HR autoimmune risk variant.

The complex geography of domestication of the African rice Oryza glaberrima.

While the domestication history of Asian rice has been extensively studied, details of the evolution of African rice remain elusive. The inner Niger delta has been suggested as the center of origin but molecular data to support this hypothesis is lacking. Here, we present a comprehensive analysis of the evolutionary and domestication history of African rice. By analyzing whole genome re-sequencing data from 282 individuals of domesticated African rice Oryza glaberrima and its progenitor O. barthii, we hypothesize a non-centric (i.e. multiregional) domestication origin for African rice. Our analyses showed genetic structure within O. glaberrima that has a geographical association. Furthermore, we have evidence that the previously hypothesized O. barthii progenitor populations in West Africa have evolutionary signatures similar to domesticated rice and carried causal domestication mutations, suggesting those progenitors were either mislabeled or may actually represent feral wild-domesticated hybrids. Phylogeographic analysis of genes involved in the core domestication process suggests that the origins of causal domestication mutations could be traced to wild progenitors in multiple different locations in West and Central Africa. In addition, measurements of panicle threshability, a key early domestication trait for seed shattering, were consistent with the gene phylogeographic results. We suggest seed non-shattering was selected from multiple genotypes, possibly arising from different geographical regions. Based on our evidence, O. glaberrima was not domesticated from a single centric location but was a result of diffuse process where multiple regions contributed key alleles for different domestication traits.

Flax latitudinal adaptation at LuTFL1 altered architecture and promoted fiber production.

After domestication in the Near East around 10,000 years ago several founder crops, flax included, spread to European latitudes. On reaching northerly latitudes the architecture of domesticated flax became more suitable to fiber production over oil, with longer stems, smaller seeds and fewer axillary branches. Latitudinal adaptations in crops typically result in changes in flowering time, often involving the PEBP family of genes that also have the potential to influence plant architecture. Two PEBP family genes in the flax genome, LuTFL1 and LuTFL2, vary in wild and cultivated flax over latitudinal range with cultivated flax receiving LuTFL1 alleles from northerly wild flax populations. Compared to a background of population structure of flaxes over latitude, the LuTFL1 alleles display a level of differentiation that is consistent with selection for an allele III in the north. We demonstrate through heterologous expression in Arabidopsis thaliana that LuTFL1 is a functional homolog of TFL1 in A. thaliana capable of changing both flowering time and plant architecture. We conclude that specialized fiber flax types could have formed as a consequence of a natural adaptation of cultivated flax to higher latitudes.

Comparing Arabidopsis receptor kinase and receptor protein-mediated immune signaling reveals BIK1-dependent differences.

Pattern recognition receptors (PRRs) sense microbial patterns and activate innate immunity against attempted microbial invasions. The leucine-rich repeat receptor kinases (LRR-RK) FLS2 and EFR, and the LRR receptor protein (LRR-RP) receptors RLP23 and RLP42, respectively, represent prototypical members of these two prominent and closely related PRR families. We conducted a survey of Arabidopsis thaliana immune signaling mediated by these receptors to address the question of commonalities and differences between LRR-RK and LRR-RP signaling. Quantitative differences in timing and amplitude were observed for several early immune responses, with RP-mediated responses typically being slower and more prolonged than those mediated by RKs. Activation of RLP23, but not FLS2, induced the production of camalexin. Transcriptomic analysis revealed that RLP23-regulated genes represent only a fraction of those genes differentially expressed upon FLS2 activation. Several positive and negative regulators of FLS2-signaling play similar roles in RLP23 signaling. Intriguingly, the cytoplasmic receptor kinase BIK1, a positive regulator of RK signaling, acts as a negative regulator of RP-type immune receptors in a manner dependent on BIK1 kinase activity. Our study unveiled unexpected differences in two closely related receptor systems and reports a new negative role of BIK1 in plant immunity.

Evolutionary and ecological functional genomics, from lab to the wild.

Plant phenotypes are the result of both genetic and environmental forces that act to modulate trait expression. Over the last few years, numerous approaches in functional genomics and systems biology have led to a greater understanding of plant phenotypic variation and plant responses to the environment. These approaches, and the questions that they can address, have been loosely termed evolutionary and ecological functional genomics (EEFG), and have been providing key insights on how plants adapt and evolve. In particular, by bringing these studies from the laboratory to the field, EEFG studies allow us to gain greater knowledge of how plants function in their natural contexts.

Modulation of ACD6 dependent hyperimmunity by natural alleles of an Arabidopsis thaliana NLR resistance gene.

Plants defend themselves against pathogens by activating an array of immune responses. Unfortunately, immunity programs may also cause unintended collateral damage to the plant itself. The quantitative disease resistance gene ACCELERATED CELL DEATH 6 (ACD6) serves to balance growth and pathogen resistance in natural populations of Arabidopsis thaliana. An autoimmune allele, ACD6-Est, which strongly reduces growth under specific laboratory conditions, is found in over 10% of wild strains. There is, however, extensive variation in the strength of the autoimmune phenotype expressed by strains with an ACD6-Est allele, indicative of genetic modifiers. Quantitative genetic analysis suggests that ACD6 activity can be modulated in diverse ways, with different strains often carrying different large-effect modifiers. One modifier is SUPPRESSOR OF NPR1-1, CONSTITUTIVE 1 (SNC1), located in a highly polymorphic cluster of nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes, which are prototypes for qualitative disease resistance genes. Allelic variation at SNC1 correlates with ACD6-Est activity in multiple accessions, and a common structural variant affecting the NL linker sequence can explain differences in SNC1 activity. Taken together, we find that an NLR gene can mask the activity of an ACD6 autoimmune allele in natural A. thaliana populations, thereby linking different arms of the plant immune system.

Genome-wide analysis of chromatin packing in Arabidopsis thaliana at single-gene resolution.

The three-dimensional packing of the genome plays an important role in regulating gene expression. We have used Hi-C, a genome-wide chromatin conformation capture (3C) method, to analyze Arabidopsis thaliana chromosomes dissected into subkilobase segments, which is required for gene-level resolution in this species with a gene-dense genome. We found that the repressive H3K27me3 histone mark is overrepresented in the promoter regions of genes that are in conformational linkage over long distances. In line with the globally dispersed distribution of RNA polymerase II in A. thaliana nuclear space, actively transcribed genes do not show a strong tendency to associate with each other. In general, there are often contacts between 5' and 3' ends of genes, forming local chromatin loops. Such self-loop structures of genes are more likely to occur in more highly expressed genes, although they can also be found in silent genes. Silent genes with local chromatin loops are highly enriched for the histone variant H3.3 at their 5' and 3' ends but depleted of repressive marks such as heterochromatic histone modifications and DNA methylation in flanking regions. Our results suggest that, different from animals, a major theme of genome folding in A. thaliana is the formation of structural units that correspond to gene bodies.

Species-wide genetic incompatibility analysis identifies immune genes as hot spots of deleterious epistasis.

Intraspecific genetic incompatibilities prevent the assembly of specific alleles into single genotypes and influence genome- and species-wide patterns of sequence variation. A common incompatibility in plants is hybrid necrosis, characterized by autoimmune responses due to epistatic interactions between natural genetic variants. By systematically testing thousands of F1 hybrids of Arabidopsis thaliana strains, we identified a small number of incompatibility hot spots in the genome, often in regions densely populated by nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes. In several cases, these immune receptor loci interact with each other, suggestive of conflict within the immune system. A particularly dangerous locus is a highly variable cluster of NLR genes, DM2, which causes multiple independent incompatibilities with genes that encode a range of biochemical functions, including NLRs. Our findings suggest that deleterious interactions of immune receptors limit the combinations of favorable disease resistance alleles accessible to plant genomes.

Activation of the Arabidopsis thaliana immune system by combinations of common ACD6 alleles.

A fundamental question in biology is how multicellular organisms distinguish self and non-self. The ability to make this distinction allows animals and plants to detect and respond to pathogens without triggering immune reactions directed against their own cells. In plants, inappropriate self-recognition results in the autonomous activation of the immune system, causing affected individuals to grow less well. These plants also suffer from spontaneous cell death, but are at the same time more resistant to pathogens. Known causes for such autonomous activation of the immune system are hyperactive alleles of immune regulators, or epistatic interactions between immune regulators and unlinked genes. We have discovered a third class, in which the Arabidopsis thaliana immune system is activated by interactions between natural alleles at a single locus, ACCELERATED CELL DEATH 6 (ACD6). There are two main types of these interacting alleles, one of which has evolved recently by partial resurrection of a pseudogene, and each type includes multiple functional variants. Most previously studies hybrid necrosis cases involve rare alleles found in geographically unrelated populations. These two types of ACD6 alleles instead occur at low frequency throughout the range of the species, and have risen to high frequency in the Northeast of Spain, suggesting a role in local adaptation. In addition, such hybrids occur in these populations in the wild. The extensive functional variation among ACD6 alleles points to a central role of this locus in fine-tuning pathogen defenses in natural populations.

Large-scale production and evaluation of marker-free indica rice IR64 expressing phytoferritin genes.

Biofortification of rice (Oryza sativa L.) using a transgenic approach to increase the amount of iron in the grain is proposed as a low-cost, reliable, and sustainable solution to help developing countries combat anemia. In this study, we generated and evaluated a large number of rice or soybean ferritin over-accumulators in rice mega-variety IR64, including marker-free events, by introducing soybean or rice ferritin genes into the endosperm for product development. Accumulation of the protein was confirmed by ELISA, in situ immunological detection, and Western blotting. As much as a 37- and 19-fold increase in the expression of ferritin gene in single and co-transformed plants, respectively, and a 3.4-fold increase in Fe content in the grain over the IR64 wild type was achieved using this approach. Agronomic characteristics of a total of 1,860 progenies from 58 IR64 single independent transgenic events and 768 progenies from 27 marker-free transgenic events were evaluated and most trait characteristics did not show a penalty. Grain quality evaluation of high-Fe IR64 transgenic events showed quality similar to that of the wild-type IR64. To understand the effect of transgenes on iron homeostasis, transcript analysis was conducted on a subset of genes involved in iron uptake and loading. Gene expression of the exogenous ferritin gene in grain correlates with protein accumulation and iron concentration. The expression of NAS2 and NAS3 metal transporters increased during the grain milky stage.