Duplication of a well-conserved homeodomain-leucine zipper transcription factor gene in barley generates a copy with more specific functions

HOMEOBOX2, the paralog of SIX-ROWED SPIKE1/HOMEOBOX1, is dispensable for barley spikelet development

The HD-ZIP class I transcription factor Homeobox 1 (HvHOX1), also known as Vulgare Row-type Spike 1 (VRS1) or Six-rowed Spike 1, regulates lateral spikelet fertility in barley (Hordeum vulgare L.). It was shown that HvHOX1 has a high expression only in lateral spikelets, while its paralog HvHOX2 was found to be expressed in different plant organs. Yet, the mechanistic functions of HvHOX1 and HvHOX2 during spikelet development are still fragmentary. Here, we show that compared with HvHOX1, HvHOX2 is more highly conserved across different barley genotypes and Hordeum species, hinting at a possibly vital but still unclarified biological role. Using bimolecular fluorescence complementation, DNA-binding, and transactivation assays, we validate that HvHOX1 and HvHOX2 are bona fide transcriptional activators that may potentially heterodimerize. Accordingly, both genes exhibit similar spatiotemporal expression patterns during spike development and growth, albeit their mRNA levels differ quantitatively. We show that HvHOX1 delays the lateral spikelet meristem differentiation and affects fertility by aborting the reproductive organs. Interestingly, the ancestral relationship of the two genes inferred from their co-expressed gene networks suggested that HvHOX1 and HvHOX2 might play a similar role during barley spikelet development. However, CRISPR-derived mutants of HvHOX1 and HvHOX2 demonstrated the suppressive role of HvHOX1 on lateral spikelets, while the loss of HvHOX2 does not influence spikelet development. Collectively, our study shows that through the suppression of reproductive organs, lateral spikelet fertility is regulated by HvHOX1, whereas HvHOX2 is dispensable for spikelet development in barley.

Boosting Triticeae crop grain yield by manipulating molecular modules to regulate inflorescence architecture: insights and knowledge from other cereal crops

One of the challenges for global food security is to reliably and sustainably improve the grain yield of cereal crops. One solution is to modify the architecture of the grain-bearing inflorescence to optimize for grain number and size. Cereal inflorescences are complex structures, with determinacy, branching patterns, and spikelet/floret growth patterns that vary by species. Recent decades have witnessed rapid advancements in our understanding of the genetic regulation of inflorescence architecture in rice, maize, wheat, and barley. Here, we summarize current knowledge on key genetic factors underlying the different inflorescence morphologies of these crops and model plants (Arabidopsis and tomato), focusing particularly on the regulation of inflorescence meristem determinacy and spikelet meristem identity and determinacy. We also discuss strategies to identify and utilize these superior alleles to optimize inflorescence architecture and, ultimately, improve crop grain yield.

GRASSY TILLERS1 (GT1) and SIX-ROWED SPIKE1 (VRS1) homologs share conserved roles in growth repression

Crop engineering and de novo domestication using gene editing are new frontiers in agriculture. However, outside of well-studied crops and model systems, prioritizing engineering targets remains challenging. Evolution can guide us, revealing genes with deeply conserved roles that have repeatedly been selected in the evolution of plant form. Homologs of the transcription factor genes GRASSY TILLERS1 (GT1) and SIX-ROWED SPIKE1 (VRS1) have repeatedly been targets of selection in domestication and evolution, where they repress growth in many developmental contexts. This suggests a conserved role for these genes in regulating growth repression. To test this, we determined the roles of GT1 and VRS1 homologs in maize (Zea mays) and the distantly related grass brachypodium (Brachypodium distachyon) using gene editing and mutant analysis. In maize, gt1; vrs1-like1 (vrl1) mutants have derepressed growth of floral organs. In addition, gt1; vrl1 mutants bore more ears and more branches, indicating broad roles in growth repression. In brachypodium, Bdgt1; Bdvrl1 mutants have more branches, spikelets, and flowers than wild-type plants, indicating conserved roles for GT1 and VRS1 homologs in growth suppression over ca. 59 My of grass evolution. Importantly, many of these traits influence crop productivity. Notably, maize GT1 can suppress growth in arabidopsis (Arabidopsis thaliana) floral organs, despite ca. 160 My of evolution separating the grasses and arabidopsis. Thus, GT1 and VRS1 maintain their potency as growth regulators across vast timescales and in distinct developmental contexts. This work highlights the power of evolution to inform gene editing in crop improvement.

Evolutionary Consequences of Functional and Regulatory Divergence of HD-Zip I Transcription Factors as a Source of Diversity in Protein Interaction Networks in Plants

The HD superfamily has been studied in detail for several decades. The plant-specific HD-Zip I subfamily attracts the most attention because of its involvement in plant development and stress responses. In this review, we provide a comprehensive insight into the evolutionary events responsible for the functional redundancy and diversification of the HD-Zip I genes in regulating various biological processes. We summarized the evolutionary history of the HD-Zip family, highlighting the important role of WGDs in its expansion and divergence of retained duplicates in the genome. To determine the relationship between the evolutionary origin and functional conservation of HD-Zip I in different species, we performed a phylogenetic analysis, compared their expression profiles in different tissues and under stress and traced the role of orthologs and paralogs in regulating developmental processes. We found that HD-Zip I from different species have similar gene structures with a highly conserved HD and Zip, bind to the same DNA sequences and are involved in similar biological processes. However, they exhibit a functional diversity, which is manifested in altered expression patterns. Some of them are involved in the regulation of species-specific leaf morphology and phenotypes. Here, we discuss the role of changes in functional domains involved in DNA binding and protein interaction of HD-Zip I and in cis-regulated regions of its target genes in promoting adaptive innovations through the formation of de novo regulatory systems. Understanding the role of the HD-Zip I subfamily in organism-environment interactions remains a challenge for evolutionary developmental biology (evo-devo).

Grain Disarticulation in Wild Wheat and Barley

Our industrial-scale crop monocultures, which are necessary to provide grain for large-scale food and feed production, are highly vulnerable to biotic and abiotic stresses. Crop wild relatives have adapted to harsh environmental conditions over millennia; thus, they are an important source of genetic variation and crop diversification. Despite several examples where significant yield increases have been achieved through the introgression of genomic regions from wild relatives, more detailed understanding of the differences between wild and cultivated species for favorable and unfavorable traits is still required to harness these valuable resources. Recently, as an alternative to the introgression of beneficial alleles from the wild into domesticated species, a radical suggestion is to domesticate wild relatives to generate new crops. A first and critical step for the domestication of cereal wild relatives would be to prevent grain disarticulation from the inflorescence at maturity. Discovering the molecular mechanisms and understanding the network of interactions behind grain retention/disarticulation would enable the implementation of approaches to select for this character in targeted species. Brittle rachis 1 and Brittle rachis 2 are major genes responsible for grain disarticulation in the wild progenitors of wheat and barley that were the target of mutations during domestication. These two genes are only found in the Triticeae tribe and are hypothesized to have evolved by a duplication followed by neo-functionalization. Current knowledge gaps include the molecular mechanisms controlling grain retention in cereals and the genomic consequences of strong selection for this essential character.

Global Role of Crop Genomics in the Face of Climate Change

The development of climate change resilient crops is necessary if we are to meet the challenge of feeding the growing world's population. We must be able to increase food production despite the projected decrease in arable land and unpredictable environmental conditions. This review summarizes the technological and conceptual advances that have the potential to transform plant breeding, help overcome the challenges of climate change, and initiate the next plant breeding revolution. Recent developments in genomics in combination with high-throughput and precision phenotyping facilitate the identification of genes controlling critical agronomic traits. The discovery of these genes can now be paired with genome editing techniques to rapidly develop climate change resilient crops, including plants with better biotic and abiotic stress tolerance and enhanced nutritional value. Utilizing the genetic potential of crop wild relatives (CWRs) enables the domestication of new species and the generation of synthetic polyploids. The high-quality crop plant genome assemblies and annotations provide new, exciting research targets, including long non-coding RNAs (lncRNAs) and cis-regulatory regions. Metagenomic studies give insights into plant-microbiome interactions and guide selection of optimal soils for plant cultivation. Together, all these advances will allow breeders to produce improved, resilient crops in relatively short timeframes meeting the demands of the growing population and changing climate.

Interaction between row-type genes in barley controls meristem determinacy and reveals novel routes to improved grain

Hordeum species develop a central spikelet flanked by two lateral spikelets at each inflorescence node. In 'two-rowed' spikes, the central spikelet alone is fertile and sets grain, while in 'six-rowed' spikes, lateral spikelets can also produce grain. Induced loss-of-function alleles of any of five Six-rowed spike (VRS) genes (VRS1-5) cause complete to intermediate gains of lateral spikelet fertility. Current six-row cultivars contain natural defective vrs1 and vrs5 alleles. Little information is known about VRS mechanism(s). We used comparative developmental, expression and genetic analyses on single and double vrs mutants to learn more about how VRS genes control development and assess their agronomic potential. We show that all VRS genes repress fertility at carpel and awn emergence in developing lateral spikelets. VRS4, VRS3 and VRS5 work through VRS1 to suppress fertility, probably by inducing VRS1 expression. Pairing vrs3, vrs4 or vrs5 alleles increased lateral spikelet fertility, despite the presence of a functional VRS1 allele. The vrs3 allele caused loss of spikelet identity and determinacy, improved grain homogeneity and increased tillering in a vrs4 background, while with vrs5, decreased tiller number and increased grain weight. Interactions amongst VRS genes control spikelet infertility, determinacy and outgrowth, and novel routes to improving six-row grain.

Genetic pathways controlling inflorescence architecture and development in wheat and barley

Modifications of inflorescence architecture have been crucial for the successful domestication of wheat and barley, which are central members of the Triticeae tribe that provide essential grains for the human diet. Investigation of the genes and alleles that underpin domestication-related traits has provided valuable insights into the molecular regulation of inflorescence development of the Triticeae, and further investigation of modified forms of architecture are proving to be equally fruitful. The identified genes are involved in diverse biological processes, including transcriptional regulation, hormone biosynthesis and metabolism, post-transcriptional and post-translational regulation, which alter inflorescence architecture by modifying the development and fertility of lateral organs, called spikelets and florets. Recent advances in sequencing capabilities and the generation of mutant populations are accelerating the identification of genes that influence inflorescence development, which is important given that genetic variation for this trait promises to be a valuable resource for optimizing grain production. This review assesses recent advances in our understanding of the genes controlling inflorescence development in wheat and barley, with the aim of highlighting the importance of improvements in developmental biology for optimizing the agronomic performance of staple crop plants.

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.

Functional analysis of the HD-Zip transcription factor genes Oshox12 and Oshox14 in rice

The homeodomain-leucine zipper (HD-Zip) transcription factor family plays vital roles in plant development and morphogenesis as well as responses to biotic and abiotic stresses. In barley, a recessive mutation in Vrs1 (HvHox1) changes two-rowed barley to six-rowed barley, which improves yield considerably. The Vrs1 gene encodes an HD-Zip subfamily I transcription factor. Phylogenetic analysis has shown that the rice HD-Zip I genes Oshox12 and Oshox14 are the closest homologues of Vrs1. Here, we show that Oshox12 and Oshox14 are ubiquitously expressed with higher levels in developing panicles. Trans-activation assays in yeast and rice protoplasts demonstrated that Oshox12 and Oshox14 can bind to a specific DNA sequence, AH1 (CAAT(A/T)ATTG), and activate reporter gene expression. Overexpression of Oshox12 and Oshox14 in rice resulted in reduced panicle length and a dwarf phenotype. In addition, Oshox14 overexpression lines showed a deficiency in panicle exsertion. Our findings suggest that Oshox12 and Oshox14 may be involved in the regulation of panicle development. This study provides a significant advancement in understanding the functions of HD-Zip transcription factors in rice.

Extreme Suppression of Lateral Floret Development by a Single Amino Acid Change in the VRS1 Transcription Factor

Increasing grain yield is an endless challenge for cereal crop breeding. In barley (Hordeum vulgare), grain number is controlled mainly by Six-rowed spike 1 (Vrs1), which encodes a homeodomain leucine zipper class I transcription factor. However, little is known about the genetic basis of grain size. Here, we show that extreme suppression of lateral florets contributes to enlarged grains in deficiens barley. Through a combination of fine-mapping and resequencing of deficiens mutants, we have identified that a single amino acid substitution at a putative phosphorylation site in VRS1 is responsible for the deficiens phenotype. deficiens mutant alleles confer an increase in grain size, a reduction in plant height, and a significant increase in thousand grain weight in contemporary cultivated germplasm. Haplotype analysis revealed that barley carrying the deficiens allele (Vrs1.t1) originated from two-rowed types carrying the Vrs1.b2 allele, predominantly found in germplasm from northern Africa. In situ hybridization of histone H4, a marker for cell cycle or proliferation, showed weaker expression in the lateral spikelets compared with central spikelets in deficiens Transcriptome analysis revealed that a number of histone superfamily genes were up-regulated in the deficiens mutant, suggesting that enhanced cell proliferation in the central spikelet may contribute to larger grains. Our data suggest that grain yield can be improved by suppressing the development of specific organs that are not positively involved in sink/source relationships.

Leaf primordium size specifies leaf width and vein number among row-type classes in barley

Exploring genes with impact on yield-related phenotypes is the preceding step to accomplishing crop improvements while facing a growing world population. A genome-wide association scan on leaf blade area (LA) in a worldwide spring barley collection (Hordeum vulgare L.), including 125 two- and 93 six-rowed accessions, identified a gene encoding the homeobox transcription factor, Six-rowed spike 1 (VRS1). VRS1 was previously described as a key domestication gene affecting spike development. Its mutation converts two-rowed (wild-type VRS1, only central fertile spikelets) into six-rowed spikes (mutant vrs1, fully developed fertile central and lateral spikelets). Phenotypic analyses of mutant and wild-type leaves revealed that mutants had an increased leaf width with more longitudinal veins. The observed significant increase of LA and leaf nitrogen (%) during pre-anthesis development in vrs1 mutants also implies a link between wider leaf and grain number, which was validated from the association of vrs1 locus with wider leaf and grain number. Histological and gene expression analyses indicated that VRS1 might influence the size of leaf primordia by affecting cell proliferation of leaf primordial cells. This finding was supported by the transcriptome analysis of mutant and wild-type leaf primordia where in the mutant transcriptional activation of genes related to cell proliferation was detectable. Here we show that VRS1 has an independent role on barley leaf development which might influence the grain number.

Major genes determining yield-related traits in wheat and barley

Current development of advanced biotechnology tools allows us to characterize the role of key genes in plant productivity. The implementation of this knowledge in breeding strategies might accelerate the progress in obtaining high-yielding cultivars. The achievements of the Green Revolution were based on a specific plant ideotype, determined by a single gene involved in gibberellin signaling or metabolism. Compared with the 1950s, an enormous increase in our knowledge about the biological basis of plant productivity has opened new avenues for novel breeding strategies. The large and complex genomes of diploid barley and hexaploid wheat represent a great challenge, but they also offer a large reservoir of genes that can be targeted for breeding. We summarize examples of productivity-related genes/mutants in wheat and barley, identified or characterized by means of modern biology. The genes are classified functionally into several groups, including the following: (1) transcription factors, regulating spike development, which mainly affect grain number; (2) genes involved in metabolism or signaling of growth regulators-cytokinins, gibberellins, and brassinosteroids-which control plant architecture and in consequence stem hardiness and grain yield; (3) genes determining cell division and proliferation mainly impacting grain size; (4) floral regulators influencing inflorescence architecture and in consequence seed number; and (5) genes involved in carbohydrate metabolism having an impact on plant architecture and grain yield. The implementation of selected genes in breeding programs is discussed, considering specific genotypes, agronomic and climate conditions, and taking into account that many of the genes are members of multigene families.

Mutations in Barley Row Type Genes Have Pleiotropic Effects on Shoot Branching

Cereal crop yield is determined by different yield components such as seed weight, seed number per spike and the tiller number and spikes. Negative correlations between these traits are often attributed to resource limitation. However, recent evidence suggests that the same genes or regulatory modules can regulate both inflorescence branching and tillering. It is therefore important to explore the role of genetic correlations between different yield components in small grain cereals. In this work, we studied pleiotropic effects of row type genes on seed size, seed number per spike, thousand grain weight, and tillering in barley to better understand the genetic correlations between individual yield components. Allelic mutants of nine different row type loci (36 mutants), in the original spring barley varieties Barke, Bonus and Foma and introgressed in the spring barley cultivar Bowman, were phenotyped under greenhouse and outdoor conditions. We identified two main mutant groups characterized by their relationships between seed and tillering parameters. The first group comprises all mutants with an increased number of seeds and significant change in tiller number at early development (group 1a) or reduced tillering only at full maturity (group 1b). Mutants in the second group are characterized by a reduction in seeds per spike and tiller number, thus exhibiting positive correlations between seed and tiller number. Reduced tillering at full maturity (group 1b) is likely due to resource limitations. In contrast, altered tillering at early development (groups 1a and 2) suggests that the same genes or regulatory modules affect inflorescence and shoot branching. Understanding the genetic bases of the trade-offs between these traits is important for the genetic manipulation of individual yield components.

Transcriptomic and phytochemical analysis of the biosynthesis of characteristic constituents in tea (Camellia sinensis) compared with oil tea (Camellia oleifera)

BACKGROUND

Tea plants (Camellia sinensis) are used to produce one of the most important beverages worldwide. The nutritional value and healthful properties of tea are closely related to the large amounts of three major characteristic constituents including polyphenols (mainly catechins), theanine and caffeine. Although oil tea (Camellia oleifera) belongs to the genus Camellia, this plant lacks these three characteristic constituents. Comparative analysis of tea and oil tea via RNA-Seq would help uncover the genetic components underlying the biosynthesis of characteristic metabolites in tea.

RESULTS

We found that 3,787 and 3,359 bud genes, as well as 4,042 and 3,302 leaf genes, were up-regulated in tea and oil tea, respectively. High-performance liquid chromatography (HPLC) analysis revealed high levels of all types of catechins, theanine and caffeine in tea compared to those in oil tea. Activation of the genes involved in the biosynthesis of these characteristic compounds was detected by RNA-Seq analysis. In particular, genes encoding enzymes involved in flavonoid, theanine and caffeine pathways exhibited considerably different expression levels in tea compared to oil tea, which were also confirmed by quantitative RT-PCR (qRT-PCR).

CONCLUSION

We assembled 81,826 and 78,863 unigenes for tea and oil tea, respectively, based on their differences at the transcriptomic level. A potential connection was observed between gene expression and content variation for catechins, theanine and caffeine in tea and oil tea. The results demonstrated that the metabolism was activated during the accumulation of characteristic metabolites in tea, which were present at low levels in oil tea. From the molecular biological perspective, our comparison of the transcriptomes and related metabolites revealed differential regulatory mechanisms underlying secondary metabolic pathways in tea versus oil tea.

Iron deficiency in barley plants: phytosiderophore release, iron translocation, and DNA methylation

All living organisms require iron (Fe) to carry out many crucial metabolic pathways. Despite its high concentrations in the geosphere, Fe bio-availability to plant roots can be very scarce. To cope with Fe shortage, plants can activate different strategies. For these reasons, we investigated Fe deficient Hordeum vulgare L. plants by monitoring growth, phytosiderophores (PS) release, iron content, and translocation, and DNA methylation, with respect to Fe sufficient ones. Reductions of plant growth, roots to shoots Fe translocation, and increases in PS release were found. Experiments on DNA methylation highlighted significant differences between fully and hemy-methylated sequences in Fe deficient plants, with respect to Fe sufficient plants. Eleven DNA bands differently methylated were found in starved plants. Of these, five sequences showed significant alignment to barley genes encoding for a glucosyltransferase, a putative acyl carrier protein, a peroxidase, a β-glucosidase and a transcription factor containing a Homeodomin. A resupply experiment was carried out on starved barley re-fed at 13 days after sowing (DAS), and it showed that plants did not recover after Fe addition. In fact, Fe absorption and root to shoot translocation capacities were impaired. In addition, resupplied barley showed DNA methylation/demethylation patterns very similar to that of barley grown in Fe deprivation. This last finding is very encouraging because it indicates as these variations/modifications could be transmitted to progenies.

Plant genetics. A Y-chromosome-encoded small RNA acts as a sex determinant in persimmons

In plants, multiple lineages have evolved sex chromosomes independently, providing a powerful comparative framework, but few specific determinants controlling the expression of a specific sex have been identified. We investigated sex determinants in the Caucasian persimmon, Diospyros lotus, a dioecious plant with heterogametic males (XY). Male-specific short nucleotide sequences were used to define a male-determining region. A combination of transcriptomics and evolutionary approaches detected a Y-specific sex-determinant candidate, OGI, that displays male-specific conservation among Diospyros species. OGI encodes a small RNA targeting the autosomal MeGI gene, a homeodomain transcription factor regulating anther fertility in a dosage-dependent fashion. This identification of a feminizing gene suppressed by a Y-chromosome-encoded small RNA contributes to our understanding of the evolution of sex chromosome systems in higher plants.

Plant homeodomain-leucine zipper I transcription factors exhibit different functional AHA motifs that selectively interact with TBP or/and TFIIB

Different members of the HD-Zip I family of transcription factors exhibit differential AHA-like activation motifs, able to interact with proteins of the basal transcriptional machinery. Homeodomain-leucine zipper proteins are transcription factors unique to plants, classified in four subfamilies. Subfamily I members have been mainly associated to abiotic stress responses. Several ones have been characterized using knockout or overexpressors plants, indicating that they take part in different signal transduction pathways even when their expression patterns are similar and they bind the same DNA sequence. A bioinformatic analysis has revealed the existence of conserved motifs outside the HD-Zip domain, including transactivation AHA motifs. Here, we demonstrate that these putative activation motifs are functional. Four members of the Arabidopsis family were chosen: AtHB1, AtHB7, AtHB12 and AtHB13. All of them exhibited activation activity in yeast and in plants but with different degrees. The protein segment necessary for such activation was different for these four transcription factors as well as the role of the tryptophans they present. When interaction with components of the basal transcription machinery was tested, AtHB1 was able to interact with TBP, AtHB12 interacted with TFIIB, AtHB7 interacted with both, TBP and TFIIB while AtHB13 showed weak interactions with any of them, in yeast two-hybrid as well as in pull-down assays. Transient transformation of Arabidopsis seedlings confirmed the activation capacity and specificity of these transcription factors and showed some differences with the results obtained in yeast. In conclusion, the differential activation functionality of these transcription factors adds an important level of functional divergence of these proteins, and together with their expression patterns, these differences could explain, at least in part, their functional divergence.

Arabidopsis AtHB7 and AtHB12 evolved divergently to fine tune processes associated with growth and responses to water stress

BACKGROUND

Arabidopsis AtHB7 and AtHB12 transcription factors (TFs) belong to the homeodomain-leucine zipper subfamily I (HD-Zip I) and present 62% amino acid identity. These TFs have been associated with the control of plant development and abiotic stress responses; however, at present it is not completely understood how AtHB7 and AtHB12 regulate these processes.

RESULTS

By using different expression analysis approaches, we found that AtHB12 is expressed at higher levels during early Arabidopsis thaliana development whereas AtHB7 during later developmental stages. Moreover, by analysing gene expression in single and double Arabidopsis mutants and in transgenic plants ectopically expressing these TFs, we discovered a complex mechanism dependent on the plant developmental stage and in which AtHB7 and AtHB12 affect the expression of each other. Phenotypic analysis of transgenic plants revealed that AtHB12 induces root elongation and leaf development in young plants under standard growth conditions, and seed production in water-stressed plants. In contrast, AtHB7 promotes leaf development, chlorophyll levels and photosynthesis and reduces stomatal conductance in mature plants. Moreover AtHB7 delays senescence processes in standard growth conditions.

CONCLUSIONS

We demonstrate that AtHB7 and AtHB12 have overlapping yet specific roles in several processes related to development and water stress responses. The analysis of mutant and transgenic plants indicated that the expression of AtHB7 and AtHB12 is regulated in a coordinated manner, depending on the plant developmental stage and the environmental conditions. The results suggested that AtHB7 and AtHB12 evolved divergently to fine tune processes associated with development and responses to mild water stress.

Genetic structure and linkage disequilibrium in a diverse, representative collection of the C4 model plant, Sorghum bicolor

To facilitate the mapping of genes in sorghum [Sorghum bicolor (L.) Moench] underlying economically important traits, we analyzed the genetic structure and linkage disequilibrium in a sorghum mini core collection of 242 landraces with 13,390 single-nucleotide polymorphims. The single-nucleotide polymorphisms were produced using a highly multiplexed genotyping-by-sequencing methodology. Genetic structure was established using principal component, Neighbor-Joining phylogenetic, and Bayesian cluster analyses. These analyses indicated that the mini-core collection was structured along both geographic origin and sorghum race classification. Examples of the former were accessions from Southern Africa, East Asia, and Yemen. Examples of the latter were caudatums with widespread geographical distribution, durras from India, and guineas from West Africa. Race bicolor, the most primitive and the least clearly defined sorghum race, clustered among other races and formed only one clear bicolor-centric cluster. Genome-wide linkage disequilibrium analyses showed linkage disequilibrium decayed, on average, within 10-30 kb, whereas the short arm of SBI-06 contained a linkage disequilibrium block of 20.33 Mb, confirming a previous report of low recombination on this chromosome arm. Four smaller but equally significant linkage disequilibrium blocks of 3.5-35.5 kb were detected on chromosomes 1, 2, 9, and 10. We examined the genes encoded within each block to provide a first look at candidates such as homologs of GS3 and FT that may indicate a selective sweep during sorghum domestication.

Site-directed nucleases: a paradigm shift in predictable, knowledge-based plant breeding

Conventional plant breeding exploits existing genetic variability and introduces new variability by mutagenesis. This has proven highly successful in securing food supplies for an ever-growing human population. The use of genetically modified plants is a complementary approach but all plant breeding techniques have limitations. Here, we discuss how the recent evolution of targeted mutagenesis and DNA insertion techniques based on tailor-made site-directed nucleases (SDNs) provides opportunities to overcome such limitations. Plant breeding companies are exploiting SDNs to develop a new generation of crops with new and improved traits. Nevertheless, some technical limitations as well as significant uncertainties on the regulatory status of SDNs may challenge their use for commercial plant breeding.

Divergence of expression pattern contributed to neofunctionalization of duplicated HD-Zip I transcription factor in barley

Barley (Hordeum vulgare) spikes are developmentally switched from two-rowed to six-rowed by a single recessive gene, six-rowed spike 1 (vrs1), which encodes a homeodomain-leucine zipper I class transcription factor. Vrs1 is a paralog of HvHox2 and both were generated by duplication of an ancestral gene. HvHox2 is conserved among cereals, whereas Vrs1 acquired its current function during the evolution of barley. It was unclear whether divergence of expression pattern or protein function accounted for the functionalization of Vrs1. Here, we conducted a comparative analysis of protein functions and gene expression between HvHox2 and Vrs1 to clarify the functionalization mechanism. We revealed that the transcriptional activation activity of HvHOX2 and VRS1 was conserved. In situ hybridization analysis showed that HvHox2 is localized in vascular bundles in developing spikes, whereas Vrs1 is expressed exclusively in the pistil, lemma, palea and lodicule of lateral spikelets. The transcript abundance of Vrs1 was > 10-fold greater than that of HvHox2 during the pistil developmental stage, suggesting that the essential function of Vrs1 is to inhibit gynoecial development. We demonstrated the quantitative function of Vrs1 using RNAi transgenic plants and Vrs1 expression variants. Expression analysis of six-rowed spike mutants that are nonallelic to vrs1 showed that Vrs1 expression was up-regulated by Vrs4, whereas HvHox2 expression was not. These data demonstrate that the divergence of gene expression pattern contributed to the neofunctionalization of Vrs1.

Genetic variation and metabolic pathway intricacy govern the active compound content and quality of the Chinese medicinal plant Lonicera japonica thunb

Background

Traditional Chinese medicine uses various herbs for the treatment of various diseases for thousands of years and it is now time to assess the characteristics and effectiveness of these medicinal plants based on modern genetic and molecular tools. The herb Flos Lonicerae Japonicae (FLJ or Lonicera japonica Thunb.) is used as an anti-inflammatory agent but the chemical quality of FLJ and its medicinal efficacy has not been consistent. Here, we analyzed the transcriptomes and metabolic pathways to evaluate the active medicinal compounds in FLJ and hope that this approach can be used for a variety of medicinal herbs in the future.

Results

We assess transcriptomic differences between FLJ and L. japonica Thunb. var. chinensis (Watts) (rFLJ), which may explain the variable medicinal effects. We acquired transcriptomic data (over 100 million reads) from the two herbs, using RNA-seq method and the Illumina GAII platform. The transcriptomic profiles contain over 6,000 expressed sequence tags (ESTs) for each of the three flower development stages from FLJ, as well as comparable amount of ESTs from the rFLJ flower bud. To elucidate enzymatic divergence on biosynthetic pathways between the two varieties, we correlated genes and their expression profiles to known metabolic activities involving the relevant active compounds, including phenolic acids, flavonoids, terpenoids, and fatty acids. We also analyzed the diversification of genes that process the active compounds to distinguish orthologs and paralogs together with the pathways concerning biosynthesis of phenolic acid and its connections with other related pathways.

Conclusions

Our study provides both an initial description of gene expression profiles in flowers of FLJ and its counterfeit rFLJ and the enzyme pool that can be used to evaluate FLJ quality. Detailed molecular-level analyses allow us to decipher the relationship between metabolic pathways involved in processing active medicinal compounds and gene expressions of their processing enzymes. Our evolutionary analysis revealed specific functional divergence of orthologs and paralogs, which lead to variation in gene functions that govern the profile of active compounds.

The domestication syndrome genes responsible for the major changes in plant form in the Triticeae crops

The process of crop domestication began 10,000 years ago in the transition of early humans from hunter/gatherers to pastoralists/farmers. Recent research has revealed the identity of some of the main genes responsible for domestication. Two of the major domestication events in barley were (i) the failure of the spike to disarticulate and (ii) the six-rowed spike. The former mutation increased grain yield by preventing grain loss after maturity, while the latter resulted in an up to 3-fold increase in yield potential. Here we provide an overview of the disarticulation systems and inflorescence characteristics, along with the genes underlying these traits, occurring in the Triticeae tribe.

Uncharacterized conserved motifs outside the HD-Zip domain in HD-Zip subfamily I transcription factors; a potential source of functional diversity

BACKGROUND

Plant HD-Zip transcription factors are modular proteins in which a homeodomain is associated to a leucine zipper. Of the four subfamilies in which they are divided, the tested members from subfamily I bind in vitro the same pseudopalindromic sequence CAAT(A/T)ATTG and among them, several exhibit similar expression patterns. However, most experiments in which HD-Zip I proteins were over or ectopically expressed under the control of the constitutive promoter 35S CaMV resulted in transgenic plants with clearly different phenotypes. Aiming to elucidate the structural mechanisms underlying such observation and taking advantage of the increasing information in databases of sequences from diverse plant species, an in silico analysis was performed. In addition, some of the results were also experimentally supported.

RESULTS

A phylogenetic tree of 178 HD-Zip I proteins together with the sequence conservation presented outside the HD-Zip domains allowed the distinction of six groups of proteins. A motif-discovery approach enabled the recognition of an activation domain in the carboxy-terminal regions (CTRs) and some putative regulatory mechanisms acting in the amino-terminal regions (NTRs) and CTRs involving sumoylation and phosphorylation. A yeast one-hybrid experiment demonstrated that the activation activity of ATHB1, a member of one of the groups, is located in its CTR. Chimerical constructs were performed combining the HD-Zip domain of one member with the CTR of another and transgenic plants were obtained with these constructs. The phenotype of the chimerical transgenic plants was similar to the observed in transgenic plants bearing the CTR of the donor protein, revealing the importance of this module inside the whole protein.

CONCLUSIONS

The bioinformatical results and the experiments conducted in yeast and transgenic plants strongly suggest that the previously poorly analyzed NTRs and CTRs of HD-Zip I proteins play an important role in their function, hence potentially constituting a major source of functional diversity among members of this subfamily.

Molecular analysis of phosphomannomutase (PMM) genes reveals a unique PMM duplication event in diverse Triticeae species and the main PMM isozymes in bread wheat tissues

BACKGROUND

Phosphomannomutase (PMM) is an essential enzyme in eukaryotes. However, little is known about PMM gene and function in crop plants. Here, we report molecular evolutionary and biochemical analysis of PMM genes in bread wheat and related Triticeae species.

RESULTS

Two sets of homologous PMM genes (TaPMM-1 and 2) were found in bread wheat, and two corresponding PMM genes were identified in the diploid progenitors of bread wheat and many other diploid Triticeae species. The duplication event yielding PMM-1 and 2 occurred before the radiation of diploid Triticeae genomes. The PMM gene family in wheat and relatives may evolve largely under purifying selection. Among the six TaPMM genes, the transcript levels of PMM-1 members were comparatively high and their recombinant proteins were all enzymatically active. However, PMM-2 homologs exhibited lower transcript levels, two of which were also inactive. TaPMM-A1, B1 and D1 were probably the main active isozymes in bread wheat tissues. The three isozymes differed from their counterparts in barley and Brachypodium distachyon in being more tolerant to elevated test temperatures.

CONCLUSION

Our work identified the genes encoding PMM isozymes in bread wheat and relatives, uncovered a unique PMM duplication event in diverse Triticeae species, and revealed the main active PMM isozymes in bread wheat tissues. The knowledge obtained here improves the understanding of PMM evolution in eukaryotic organisms, and may facilitate further investigations of PMM function in the temperature adaptability of bread wheat.

Somatic polyploidy promotes cell function under stress and energy depletion: evidence from tissue-specific mammal transcriptome

Polyploid cells show great among-species and among-tissues diversity and relation to developmental mode, suggesting their importance in adaptive evolution and developmental programming. At the same time, excessive polyploidization is a hallmark of functional impairment, aging, growth disorders, and numerous pathologies including cancer and cardiac diseases. To shed light on this paradox and to find out how polyploidy contributes to organ functions, we review here the ploidy-associated shifts in activity of narrowly expressed (tissue specific) genes in human and mouse heart and liver, which have the reciprocal pattern of polyploidization. For this purpose, we use the modular biology approach and genome-scale cross-species comparison. It is evident from this review that heart and liver show similar traits in response to polyploidization. In both organs, polyploidy protects vitality (mainly due to the activation of sirtuin-mediated pathways), triggers the reserve adenosine-5'-triphosphate (ATP) production, and sustains tissue-specific functions by switching them to energy saving mode. In heart, the strongest effects consisted in the concerted up-regulation of contractile proteins and substitution of energy intensive proteins with energy economic ones. As a striking example, the energy intensive alpha myosin heavy chain (providing fast contraction) decreased its expression by a factor of 10, allowing a 270-fold increase of expression of beta myosin heavy chain (providing slow contraction), which has approximately threefold lower ATP-hydrolyzing activity. The liver showed the enhancement of immunity, reactive oxygen species and xenobiotic detoxication, and numerous metabolic adaptations to long-term energy depletion. Thus, somatic polyploidy may be an ingenious evolutionary instrument for fast adaptation to stress and new environments allowing trade-offs between high functional demand, stress, and energy depletion.

Molecular evolution of regulatory genes in spruces from different species and continents: heterogeneous patterns of linkage disequilibrium and selection but correlated recent demographic changes

Genes involved in transcription regulation may represent valuable targets in association genetics studies because of their key roles in plant development and potential selection at the molecular level. Selection and demographic signatures at the sequence level were investigated for five regulatory genes belonging to the knox-I family (KN1, KN2, KN3, KN4) and the HD-Zip III family (HB-3) in three Picea species affected by post-glacial recolonization in North America and Europe. To disentangle neutral and selective forces and estimate linkage disequilibrium (LD) on a gene basis, complete or nearly complete gene sequences were analysed. Nucleotide variation within species, haplotype structure, LD, and neutrality tests, in addition to coalescent simulations based on Tajima’s D and Fay and Wu’s H, were estimated. Nucleotide diversity was generally low in all species (average π = 0.002–0.003) and much heterogeneity was seen in LD and selection signatures among genes and species. Most of the genes harboured an excess of both rare and frequent alleles in the three species. Simulations showed that this excess was significantly higher than that expected under neutrality and a bottleneck during the Last Glacial Maximum followed by population expansion at the Pleistocene/Holocene boundary or shortly after best explains the correlated sequence patterns. These results indicate that despite recent large demographic changes in the three boreal species from two continents, species-specific selection signatures could still be detected from the analysis of nearly complete regulatory gene sequences. Such different signatures indicate differential subfunctionalization of gene family members in the three congeneric species.