Simultaneous targeting of duplicated genes in Petunia protoplasts for flower color modification via CRISPR-Cas9 ribonucleoproteins.

KEY MESSAGE: We obtained a complete mutant line of Petunia having mutations in both F3H genes via Cas9-ribonucleoproteins delivery, which exhibited a pale purplish pink flower color. The CRISPR-Cas system is now revolutionizing agriculture by allowing researchers to generate various desired mutations in plants at will. In particular, DNA-free genome editing via Cas9-ribonucleoproteins (RNPs) delivery has many advantages in plants; it does not require codon optimization or specific promoters for expression in plant cells; furthermore, it can bypass GMO regulations in some countries. Here, we have performed site-specific mutagenesis in Petunia to engineer flower color modifications. We determined that the commercial Petunia cultivar 'Madness Midnight' has two F3H coding genes and designed one guide RNA that targets both F3H genes at once. Among 67 T0 plants regenerated from Cas9-RNP transfected protoplasts, we obtained seven mutant lines that contain mutations in either F3HA or F3HB gene and one complete mutant line having mutations in both F3H genes without any selectable markers. It is noteworthy that only the f3ha f3hb exhibited a clearly modified, pale purplish pink flower color (RHS 69D), whereas the others, including the single copy gene knock-out plants, displayed purple violet (RHS 93A) flowers similar to the wild-type Petunia. To the best of our knowledge, we demonstrated a precedent of ornamental crop engineering by DNA-free CRISPR method for the first time, which will greatly accelerate a transition from a laboratory to a farmer's field.

MiR1885 Regulates Disease Tolerance Genes in Brassica rapa during Early Infection with Plasmodiophora brassicae.

Clubroot caused by Plasmodiophora brassicae is a severe disease of cruciferous crops that decreases crop quality and productivity. Several clubroot resistance-related quantitative trait loci and candidate genes have been identified. However, the underlying regulatory mechanism, the interrelationships among genes, and how genes are regulated remain unexplored. MicroRNAs (miRNAs) are attracting attention as regulators of gene expression, including during biotic stress responses. The main objective of this study was to understand how miRNAs regulate clubroot resistance-related genes in P. brassicae-infected Brassica rapa. Two Brassica miRNAs, Bra-miR1885a and Bra-miR1885b, were revealed to target TIR-NBS genes. In non-infected plants, both miRNAs were expressed at low levels to maintain the balance between plant development and basal immunity. However, their expression levels increased in P. brassicae-infected plants. Both miRNAs down-regulated the expression of the TIR-NBS genes Bra019412 and Bra019410, which are located at a clubroot resistance-related quantitative trait locus. The Bra-miR1885-mediated down-regulation of both genes was detected for up to 15 days post-inoculation in the clubroot-resistant line CR Shinki and in the clubroot-susceptible line 94SK. A qRT-PCR analysis revealed Bra019412 expression was negatively regulated by miR1885. Both Bra019412 and Bra019410 were more highly expressed in CR Shinki than in 94SK; the same expression pattern was detected in multiple clubroot-resistant and clubroot-susceptible inbred lines. A 5' rapid amplification of cDNA ends analysis confirmed the cleavage of Bra019412 by Bra-miR1885b. Thus, miR1885s potentially regulate TIR-NBS gene expression during P. brassicae infections of B. rapa.

Site-directed mutagenesis in Petunia × hybrida protoplast system using direct delivery of purified recombinant Cas9 ribonucleoproteins.

KEY MESSAGE: Site-directed mutagenesis of nitrate reductase genes using direct delivery of purified Cas9 protein preassembled with guide RNA produces mutations efficiently in Petunia × hybrida protoplast system. The clustered, regularly interspaced, short palindromic repeat (CRISPR)-CRISPR associated endonuclease 9 (CRISPR/Cas9) system has been recently announced as a powerful molecular breeding tool for site-directed mutagenesis in higher plants. Here, we report a site-directed mutagenesis method targeting Petunia nitrate reductase (NR) gene locus. This method could create mutations efficiently using direct delivery of purified Cas9 protein and single guide RNA (sgRNA) into protoplast cells. After transient introduction of RNA-guided endonuclease (RGEN) ribonucleoproteins (RNPs) with different sgRNAs targeting NR genes, mutagenesis at the targeted loci was detected by T7E1 assay and confirmed by targeted deep sequencing. T7E1 assay showed that RGEN RNPs induced site-specific mutations at frequencies ranging from 2.4 to 21 % at four different sites (NR1, 2, 4 and 6) in the PhNR gene locus with average mutation efficiency of 14.9 ± 2.2 %. Targeted deep DNA sequencing revealed mutation rates of 5.3-17.8 % with average mutation rate of 11.5 ± 2 % at the same NR gene target sites in DNA fragments of analyzed protoplast transfectants. Further analysis from targeted deep sequencing showed that the average ratio of deletion to insertion produced collectively by the four NR-RGEN target sites (NR1, 2, 4, and 6) was about 63:37. Our results demonstrated that direct delivery of RGEN RNPs into protoplast cells of Petunia can be exploited as an efficient tool for site-directed mutagenesis of genes or genome editing in plant systems.

Proteome and phosphoproteome characterization reveals new response and defense mechanisms of Brachypodium distachyon leaves under salt stress.

Salinity is a major abiotic stress affecting plant growth and development. Understanding the molecular mechanisms of salt response and defense in plants will help in efforts to improve the salt tolerance of crops. Brachypodium distachyon is a new model plant for wheat, barley, and several potential biofuel grasses. In the current study, proteome and phosphoproteome changes induced by salt stress were the focus. The Bd21 leaves were initially treated with salt in concentrations ranging from 80 to 320 mm and then underwent a recovery process prior to proteome analysis. A total of 80 differentially expressed protein spots corresponding to 60 unique proteins were identified. The sample treated with a median salt level of 240 mm and the control were selected for phosphopeptide purification using TiO2 microcolumns and LC-MS/MS for phosphoproteome analysis to identify the phosphorylation sites and phosphoproteins. A total of 1509 phosphoproteins and 2839 phosphorylation sites were identified. Among them, 468 phosphoproteins containing 496 phosphorylation sites demonstrated significant changes at the phosphorylation level. Nine phosphorylation motifs were extracted from the 496 phosphorylation sites. Of the 60 unique differentially expressed proteins, 14 were also identified as phosphoproteins. Many proteins and phosphoproteins, as well as potential signal pathways associated with salt response and defense, were found, including three 14-3-3s (GF14A, GF14B, and 14-3-3A) for signal transduction and several ABA signal-associated proteins such as ABF2, TRAB1, and SAPK8. Finally, a schematic salt response and defense mechanism in B. distachyon was proposed.

Dynamic development of starch granules and the regulation of starch biosynthesis in Brachypodium distachyon: comparison with common wheat and Aegilops peregrina.

BACKGROUND: Thorough understanding of seed starch biosynthesis and accumulation mechanisms is of great importance for agriculture and crop improvement strategies. We conducted the first comprehensive study of the dynamic development of starch granules and the regulation of starch biosynthesis in Brachypodium distachyon and compared the findings with those reported for common wheat (Chinese Spring, CS) and Aegilops peregrina. RESULTS: Only B-granules were identified in Brachypodium Bd21, and the shape variation and development of starch granules were similar in the B-granules of CS and Bd21. Phylogenetic analysis showed that most of the Bd21 starch synthesis-related genes were more similar to those in wheat than in rice. Early expression of key genes in Bd21 starch biosynthesis mediate starch synthesis in the pericarp; intermediate-stage expression increases the number and size of starch granules. In contrast, these enzymes in CS and Ae. peregrina were mostly expressed at intermediate stages, driving production of new B-granules and increasing the granule size, respectively. Immunogold labeling showed that granule-bound starch synthase (GBSSI; related to amylose synthesis) was mainly present in starch granules: at lower levels in the B-granules of Bd21 than in CS. Furthermore, GBSSI was phosphorylated at threonine 183 and tyrosine 185 in the starch synthase catalytic domain in CS and Ae. peregrina, but neither site was phosphorylated in Bd21, suggesting GBSSI phosphorylation could improve amylose biosynthesis. CONCLUSIONS: Bd21 contains only B-granules, and the expression of key genes in the three studied genera is consistent with the dynamic development of starch granules. GBSSI is present in greater amounts in the B-granules of CS than in Bd21; two phosphorylation sites (Thr183 and Tyr185) were found in Triticum and Aegilops; these sites were not phosphorylated in Bd21. GBSSI phosphorylation may reflect its importance in amylose synthesis.

Molecular Evaluation of the Effects of FLC Homologs and Coordinating Regulators on the Flowering Responses to Vernalization in Cabbage (Brassica oleracea var. capitata) Genotypes.

The flowering loci of cabbage must be understood to boost their productivity. In this study, to clarify the flowering mechanisms of cabbage, we examined the three flowering repressors BoFLC1, 2 and 3, and the flowering regulators BoGI, BoCOOLAIR, and BoVIN3 of early (CAB1), middle (CAB3), and late (CAB5) flowering cabbage genotypes. Analysis of allele-specifically amplified genomic DNA and various sequence alignments demonstrated that maximal insertions and deletions influenced cabbage flowering behavior, notably in CAB3 and CAB5. Phylogenetic studies showed that BoFLC1, 2, and 3 in the CAB1, 3, and 5 genotypes had the highest homologies to other Brassica species, with CAB3 and 5 the most similar. Although CAB3 and CAB5 have comparable genetic patterns, flowering repressors and flowering regulators were investigated individually with and without vernalization to determine their minor flowering differences. The expression investigation revealed that vernalized CAB5 downregulated all BoFLC genes compared to CAB3 and, in contrast, CAB3 exhibited upregulated BoCOOLAIR. We hypothesized that the CAB3 BoFLC locus' additional insertions may have led to BoCOOLAIR overexpression and BoFLC downregulation. This study sheds light on cabbage genotypes-particularly those of CAB1 and CAB5-and suggests that structural variations in BoFLC2 and 3 bind flowering regulators, such as COOLAIR, which may affect cabbage flowering time.

iTRAQ-based quantitative proteomic analysis reveals new metabolic pathways of wheat seedling growth under hydrogen peroxide stress.

As an abundant ROS, hydrogen peroxide (H2 O2 ) plays pivotal roles in plant growth and development. In this work, we conducted for the first time an iTRAQ-based quantitative proteomic analysis of wheat seedling growth under different exogenous H2 O2 treatments. The growth of seedlings and roots was significantly restrained by increased H2 O2 concentration stress. Malondialdehyde, soluble sugar, and proline contents as well as peroxidase activity increased with increasing H2 O2 levels. A total of 3,425 proteins were identified by iTRAQ, of which 157 showed differential expression and 44 were newly identified H2 O2 -responsive proteins. H2 O2 -responsive proteins were mainly involved in stress/defense/detoxification, signal transduction, and carbohydrate metabolism. It is clear that up-regulated expression of signal transduction and stress/defence/detoxification-related proteins under H2 O2 stress, such as plasma membrane intrinsic protein 1, fasciclin-like arabinogalactan protein, and superoxide dismutase, could contribute to H2 O2 tolerance of wheat seedlings. Increased gluconeogenesis (phosphoenol-pyruvate carboxykinase) and decreased pyruvate kinase proteins are potentially related to the higher H2 O2 tolerance of wheat seedlings. A metabolic pathway of wheat seedling growth under H2 O2 stress is presented.

Establishment of targeted mutagenesis in soybean protoplasts using CRISPR/Cas9 RNP delivery via electro-transfection.

The soybean (Glycine max L.) is an important crop with high agronomic value. The improvement of agronomic traits through gene editing techniques has broad application prospects in soybean. The polyethylene glycol (PEG)-mediated cell transfection has been successfully used to deliver the CRISPR/Cas9-based ribonucleoprotein (RNP) into soybean protoplasts. However, several downstream analyses or further cell regeneration protocols might be hampered by PEG contamination within the samples. Here in this study, we attempted to transfect CRISPR/Cas9 RNPs into trifoliate leaf-derived soybean protoplasts using Neon electroporation to overcome the need for PEG transfection for the first time. We investigated different electroporation parameters including pulsing voltage (V), strength and duration of pulses regarding protoplast morphology, viability, and delivery of CRISPR/Cas9. Electroporation at various pulsing voltages with 3 pulses and 10 ms per pulse was found optimal for protoplast electro-transfection. Following electro-transfection at various pulsing voltages (500 V, 700 V, 1,000 V, and 1,300 V), intact protoplasts were observed at all treatments. However, the relative frequency of cell viability and initial cell divisions decreased with increasing voltages. Confocal laser scanning microscopy (CLSM) confirmed that the green fluorescent protein (GFP)-tagged Cas9 was successfully internalized into the protoplasts. Targeted deep sequencing results revealed that on-target insertion/deletion (InDel) frequencies were increased with increasing voltages in protoplasts electro-transfected with CRISPR/Cas9 RNPs targeting constitutive pathogen response 5 (CPR5). InDel patterns ranged from +1 bp to -6 bp at three different target sites in CPR5 locus with frequencies ranging from 3.8% to 8.1% following electro-transfection at 1,300 V and 2.1% to 3.8% for 700 V and 1,000 V, respectively. Taken together, our results demonstrate that the CRISPR/Cas9 RNP system can be delivered into soybean protoplasts by the Neon electroporation system for efficient and effective gene editing. The electro-transfection system developed in this study would also further facilitate and serve as an alternative delivery method for DNA-free genome editing of soybean and other related species for genetic screens and potential trait improvement.

Proteome characterization of developing grains in bread wheat cultivars (Triticum aestivum L.).

BACKGROUND: The analyses of protein synthesis, accumulation and regulation during grain development in wheat are more complex because of its larger genome size compared to model plants such as Arabidopsis and rice. In this study, grains from two wheat cultivars Jimai 20 and Zhoumai 16 with different gluten quality properties were harvested at five development stages, and were used to displayed variable expression patterns of grain proteins. RESULTS: Proteome characterization during grain development in Chinese bread wheat cultivars Jimai 20 and Zhoumai 16 with different quality properties was investigated by 2-DE and tandem MALDI-TOF/TOF-MS. Identification of 117 differentially accumulated protein spots representing 82 unique proteins and five main expression patterns enabled a chronological description of wheat grain formation. Significant proteome expression differences between the two cultivars were found; these included 14 protein spots that accumulated in both cultivars but with different patterns and 27 cultivar-different spots. Among the cultivar-different protein spots, 14 accumulated in higher abundance in Jimai 20 than in Zhoumai 16, and included NAD-dependent isocitrate dehydrogenase, triticin precursor, LMW-s glutenin subunit and replication factor C-like protein. These proteins are likely to be associated with superior gluten quality. In addition, some proteins such as class II chitinase and peroxidase 1 with isoforms in developing grains were shown to be phosphorylated by Pro-Q Diamond staining and phosphorprotein site prediction. Phosphorylation could have important roles in wheat grain development. qRT-PCR analysis demonstrated that transcriptional and translational expression patterns of many genes were significantly different. CONCLUSIONS: Wheat grain proteins displayed variable expression patterns at different developmental stages and a considerable number of protein spots showed differential accumulation between two cultivars. Differences in seed storage proteins were considered to be related to different quality performance of the flour from these wheat cultivars. Some proteins with isoforms were phosphorylated, and this may reflect their importance in grain development. Our results provide new insights into proteome characterization during grain development in different wheat genotypes.

Molecular characterization and dynamic expression patterns of two types of γ-gliadin genes from Aegilops and Triticum species.

Gliadins were the major components of wheat storage proteins and determine the extensibility properties of gluten dough. In this work, 19 new full-length γ-gliadin genes were isolated from various Aegilops and Triticum species. Sequence characterization showed that a specific octapeptide and celiac disease (CD)-toxic epitope Gliγ-3 (VQGQGIIQPQQPAQL) were present in the rich glutamine domain and C-terminal non-repetitive domain, respectively. Based on the sequence features of both peptides, a new classification system for γ-gliadin gene family was established, in which γ-gliadins were classified into two types (types I and II) with each consisting of two groups. An uneven distribution of different types and groups of γ-gliadin genes was exhibited among 11 Aegilops and Triticum genomes. Phylogenetic analysis revealed that types I and II genes diverged at about 14 MYA while the divergence of 4 γ-gliadin group genes occurred at around 10 MYA almost simultaneously. The γ-gliadin genes from S(l) and B genomes displayed a different transcriptional expression pattern during grain development, and rapid increasing of gliadin mRNA and proteins occurred at 15-20 DPA. In addition, genome-specific variations of CD-toxic epitopes among Aegilops and Triticum genomes were found. The A genome and its related progenitor genomes A(u) and A(m) had fewer CD epitopes than other genomes, suggesting that these genomes might be valuable gene resources to remove CD toxic peptides for wheat quality improvement.

A DNA-Free Editing Platform for Genetic Screens in Soybean via CRISPR/Cas9 Ribonucleoprotein Delivery.

CRISPR/Cas9-based ribonucleoprotein (RNP)-mediated system has the property of minimizing the effects related to the unwanted introduction of vector DNA and random integration of recombinant DNA. Here, we describe a platform based on the direct delivery of Cas9 RNPs to soybean protoplasts for genetic screens in knockout gene-edited soybean lines without the transfection of DNA vectors. The platform is based on the isolation of soybean protoplasts and delivery of Cas RNP complex. To empirically test our platform, we have chosen a model gene from the soybean genetic toolbox. We have used five different guide RNA (gRNA) sequences that targeted the constitutive pathogen response 5 (CPR5) gene associated with the growth of trichomes in soybean. In addition, efficient protoplast transformation, concentration, and ratio of Cas9 and gRNAs were optimized for soybean for the first time. Targeted mutagenesis insertion and deletion frequency and sequences were analyzed using both Sanger and targeted deep sequencing strategies. We were able to identify different mutation patterns within insertions and deletions (InDels) between + 5 nt and -30 bp and mutation frequency ranging from 4.2 to 18.1% in the GmCPR5 locus. Our results showed that DNA-free delivery of Cas9 complexes to protoplasts is a useful approach to perform early-stage genetic screens and anticipated analysis of Cas9 activity in soybeans.

Molecular Characterization and Identification of Calnexin 1 As a Radiation Biomarker from Tradescantia BNL4430.

Calnexin (CNX) is an integral membrane protein that functions as a chaperone in the endoplasmic reticulum for the correct folding of proteins under stress conditions, rendering organisms tolerant under adverse conditions. Studies have investigated the cytogenetic effects of gamma irradiation (Ɣ-IR) on plants, but information on the molecular response under Ɣ-IR remains limited. Previously, we constructed a cDNA library of an irradiation-sensitive bioindicator plant, Tradescantia BNL4430 (T-4430) under Ɣ-IR, in which the Calnexin-1 gene was highly upregulated at 50 mGy treatment. TrCNX1 encodes a 61.4 kDa protein with conserved signature motifs similar to already reported CNX1s. TrCNX1 expression was evaluated by semiquantitative reverse transcriptase PCR and quantitative real-time PCR and was ubiquitously expressed in various tissues and highly upregulated in flower petals under 50 mGy Ɣ-IR stress. The protective function of TrCNX1 was investigated by overexpression of TrCNX1 in an Escherichia coli BL21(DE3) heterologous system. Using plate assay, we showed that TrCNX1 increased the viability of E. coli transformants under both UV-B and Ɣ-IR compared with the control, demonstrating that TrCNX1 functions under irradiation stress. TrCNX1 may enhance irradiation stress tolerance in crops and act as a radio marker gene to monitor the effects of radiation.

A Genome-Wide Analysis of the Pentatricopeptide Repeat (PPR) Gene Family and PPR-Derived Markers for Flesh Color in Watermelon (Citrullus lanatus).

Watermelon (Citrullus lanatus) is an economically important fruit crop grown for consumption of its large edible fruit flesh. Pentatricopeptide-repeat (PPR) encoding genes, one of the large gene families in plants, are important RNA-binding proteins involved in the regulation of plant growth and development by influencing the expression of organellar mRNA transcripts. However, systematic information regarding the PPR gene family in watermelon remains largely unknown. In this comprehensive study, we identified and characterized a total of 422 C. lanatus PPR (ClaPPR) genes in the watermelon genome. Most ClaPPRs were intronless and were mapped across 12 chromosomes. Phylogenetic analysis showed that ClaPPR proteins could be divided into P and PLS subfamilies. Gene duplication analysis suggested that 11 pairs of segmentally duplicated genes existed. In-silico expression pattern analysis demonstrated that ClaPPRs may participate in the regulation of fruit development and ripening processes. Genotyping of 70 lines using 4 single nucleotide polymorphisms (SNPs) from 4 ClaPPRs resulted in match rates of over 0.87 for each validated SNPs in correlation with the unique phenotypes of flesh color, and could be used in differentiating red, yellow, or orange watermelons in breeding programs. Our results provide significant insights for a comprehensive understanding of PPR genes and recommend further studies on their roles in watermelon fruit growth and ripening, which could be utilized for cultivar development of watermelon.

Whole genome resequencing of watermelons to identify single nucleotide polymorphisms related to flesh color and lycopene content.

Cultivated watermelon (Citrullus lanatus) is one of the most important food crops in the Cucurbitaceae family. Diversification after domestication has led cultivated watermelons to exhibit diverse fruit flesh colors, including red, yellow, and orange. Recently, there has been increased interest in red-fleshed watermelons because they contain the antioxidant cis-isomeric lycopene. We performed whole genome resequencing (WGRS) of 24 watermelons with different flesh colors to identify single-nucleotide polymorphisms (SNPs) related to high lycopene content. The resequencing data revealed 203,894-279,412 SNPs from read mapping between inbred lines and the 97103 reference genome. In total, 295,065 filtered SNPs were identified, which had an average polymorphism information content of 0.297. Most of these SNPs were intergenic (90.1%) and possessed a transversion (Tv) rate of 31.64%. Overall, 2,369 SNPs were chosen at 0.5 Mb physical intervals to analyze genetic diversity across the 24 inbred lines. A neighbor-joining dendrogram and principal coordinate analysis (PCA) based on the 2,369 SNPs revealed that the 24 inbred lines could be grouped into high and low lycopene-type watermelons. In addition, we analyzed SNPs that could discriminate high lycopene content, red-fleshed watermelon from low lycopene, yellow or orange watermelon inbred lines. For validation, 19 SNPs (designated as WMHL1-19) were chosen randomly, and cleavage amplified polymorphic sequence (CAPS) markers were designed. Genotyping of the above 24 lines and 12 additional commercial cultivars using WMHL1-19 CAPS markers resulted in match rates of over 0.92 for most validated markers in correlation with the flesh color phenotypes. Our results provide valuable genomic information regarding the high lycopene content phenotype of red-fleshed cultivated watermelons, and the identified SNPs will be useful for the development of molecular markers in the marker-assisted breeding of watermelons with high lycopene content.

Molecular Characterization and Expression Profiling of Brachypodium distachyon L. Cystatin Genes Reveal High Evolutionary Conservation and Functional Divergence in Response to Abiotic Stress.

Cystatin is a class of proteins mainly involved in cysteine protease inhibition and plant growth and development, as well as tolerance under various abiotic stresses. In this study, we performed the first comprehensive analysis of the molecular characterization and expression profiling in response to various abiotic stresses of the cystatin gene family in Brachypodium distachyon, a novel model plant for Triticum species with huge genomes. Comprehensive searches of the Brachypodium genome database identified 25 B. distachyon cystatin (BdC) genes that are distributed unevenly on chromosomes; of these, nine and two were involved in tandem and segmental duplication events, respectively. All BdC genes had similar exon/intron structural organization, with three conserved motifs similar to those from other plant species, indicating their high evolutionary conservation. Expression profiling of 10 typical BdC genes revealed ubiquitous expression in different organs at varying expression levels. BdC gene expression in seedling leaves was particularly highly induced by various abiotic stresses, including the plant hormone abscisic acid and various environmental cues (cold, H2O2, CdCl2, salt, and drought). Interestingly, most BdC genes were significantly upregulated under multiple abiotic stresses, including BdC15 under all stresses, BdC7-2 and BdC10 under five stresses, and BdC7-1, BdC2-1, BdC14, and BdC12 under four stresses. The putative metabolic pathways of cytastin genes in response to various abiotic stresses mainly involve the aberrant protein degradation pathway and reactive oxygen species (ROS)-triggered programmed cell death signaling pathways. These observations provide a better understanding of the structural and functional characteristics of the plant cystatin gene family.

Phylogenetic Analysis, Lineage-Specific Expansion and Functional Divergence of seed dormancy 4-Like Genes in Plants.

The rice gene seed dormancy 4 (OsSdr4) functions in seed dormancy and is a major factor associated with pre-harvest sprouting (PHS). Although previous studies of this protein family were reported for rice and other species, knowledge of the evolution of genes homologous to OsSdr4 in plants remains inadequate. Fifty four Sdr4-like (hereafter designated Sdr4L) genes were identified in nine plant lineages including 36 species. Phylogenetic analysis placed these genes in eight subfamilies (I-VIII). Genes from the same lineage clustered together, supported by analysis of conserved motifs and exon-intron patterns. Segmental duplications were present in both dicot and monocot clusters, while tandemly duplicated genes occurred only in monocot clusters indicating that both tandem and segmental duplications contributed to expansion of the grass I and II subfamilies. Estimation of the approximate ages of the duplication events indicated that ancestral Sdr4 genes evolved from a common angiosperm ancestor, about 160 million years ago (MYA). Moreover, diversification of Sdr4L genes in mono and dicot plants was mainly associated with genome-wide duplication and speciation events. Functional divergence was observed in all subfamily pairs, except IV/VIIIa. Further analysis indicated that functional constraints between subfamily pairs I/II, I/VIIIb, II/VI, II/VIIIb, II/IV, and VI/VIIIb were statistically significant. Site and branch-site model analyses of positive selection suggested that these genes were under strong adaptive selection pressure. Critical amino acids detected for both functional divergence and positive selection were mostly located in the loops, pointing to functional importance of these regions in this protein family. In addition, differential expression studies by transcriptome atlas of 11 Sdr4L genes showed that the duplicated genes may have undergone divergence in expression between plant species. Our findings showed that Sdr4L genes are functionally divergent and positively selected. These may contribute to further functional analysis and molecular evolution of Sdr4L gene families in land plants.

Molecular characterization and evolutionary origins of farinin genes in Brachypodium distachyon L.

Farinins are one of the oldest members of the gluten family in wheat and Aegilops species, and they influence dough properties. Here, we performed the first detailed molecular genetic study on farinin genes in Brachypodium distachyon L., the model species for Triticum aestivum. A total of 51 b-type farinin genes were cloned and characterized, including 27 functional and 24 non-functional pseudogenes from 14 different B. distachyon accessions. All genes were highly similar to those previously reported from wheat and Aegilops species. The identification of deduced amino acid sequences showed that b-type farinins across Triticeae genomes could be classified as b1-, b2-, b3-, and b4-type farinins; however, B. distachyon had only b3- and b4-type farinins. Real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR) revealed that farinin genes are transcribed into mRNA in B. distachyon at much lower levels than in Triticeae, despite the presence of cis-acting elements in promoter regions. Phylogenetic analysis suggested that Brachypodium farinins may have closer relationships with common wheat and further confirmed four different types of b-type farinins in Triticeae and Brachypodium genomes, corresponding to b1, b2, b3 (group 1), and b4 (group 2). A putative evolutionary origin model of farinin genes in Brachypodium, Triticum, and the related species suggests that all b-type farinins diverged from their common ancestor ~3.2 million years ago (MYA). The b3 and b4 types could be considered older in the farinin family. The results explain the loss of b1- and b2-type farinin alleles in Brachypodium.

Molecular characterisation and evolution of HMW glutenin subunit genes in Brachypodium distachyon L.

Brachypodium distachyon, a small wild grass within the Pooideae family, is a new model organism for exploring the functional genomics of cereal crops. It was shown to have close relationships to wheat, barley and rice. Here, we describe the molecular characterisation and evolutionary relationships of high molecular weight glutenin subunits (HMW-GS) genes from B. distachyon. Sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE), high performance capillary electrophoresis (HPCE) and liquid chromatography-tandem mass spectrometry (LC-MS/MS) analyses demonstrated that there was no HMW-GS expression in the Brachypodium grains due to the silencing of their encoding genes. Through allele-specific polymerase chain reaction (AS-PCR) amplification and cloning, a total of 13 HMW-GS encoding genes from diploid, tetraploid and hexaploid Brachypodium species were obtained, and all of them had typical structural features of y-type HMW-GS genes from common wheat and related species, particularly more similar to the 1Dy12 gene. However, the presence of an in-frame premature stop codon (TAG) at position 1521 in the coding region resulted in the conversion of all the genes to pseudogenes. Further, quantitative real-time PCR (qRT-PCR) analysis revealed that HMW-GS genes in B. distachyon displayed a similar trend, but with a low transcriptional expression profile during grain development due to the occurrence of the stop codon. Phylogenetic analysis showed that the highly conserved Glu-1-2 loci were presented in B. distachyon, which displayed close phylogenetic evolutionary relationships with Triticum and related species.