Posttranscriptional tuning of gene expression over a large dynamic range in synthetic tobacco chloroplast operons.

Achieving optimally balanced gene expression within synthetic operons requires regulatory elements capable of providing a spectrum of expression levels. In this study, we investigate the expression of gfp reporter gene in tobacco chloroplasts, guided by variants of the plastid atpH 5' UTR, which harbors a binding site for PPR10, a protein that activates atpH at the posttranscriptional level. Our findings reveal that endogenous tobacco PPR10 confers distinct levels of reporter activation when coupled with the tobacco and maize atpH 5' UTRs in different design contexts. Notably, high GFP expression was not coupled to the stabilization of monocistronic gfp transcripts in dicistronic reporter lines, adding to the evidence that PPR10 activates translation via a mechanism that is independent of its stabilization of monocistronic transcripts. Furthermore, the incorporation of a tRNA upstream of the UTR nearly abolishes gfp mRNA (and GFP protein), presumably by promoting such rapid RNA cleavage and 5' exonucleolytic degradation that PPR10 had insufficient time to bind and protect gfp RNA, resulting in a substantial reduction in GFP accumulation. When combined with a mutant atpH 5' UTR, the tRNA leads to an exceptionally low level of transgene expression. Collectively, this approach allows for tuning of reporter gene expression across a wide range, spanning from a mere 0.02-25% of the total soluble cellular protein. These findings highlight the potential of employing cis-elements from heterologous species and expand the toolbox available for plastid synthetic biology applications requiring multigene expression at varying levels.

Turnover of mammal sex chromosomes in the Sry-deficient Amami spiny rat is due to male-specific upregulation of Sox9.

Mammalian sex chromosomes are highly conserved, and sex is determined by SRY on the Y chromosome. Two exceptional rodent groups in which some species lack a Y chromosome and Sry offer insights into how novel sex genes can arise and replace Sry, leading to sex chromosome turnover. However, intensive study over three decades has failed to reveal the identity of novel sex genes in either of these lineages. We here report our discovery of a male-specific duplication of an enhancer of Sox9 in the Amami spiny rat Tokudaia osimensis, in which males and females have only a single X chromosome (XO/XO) and the Y chromosome and Sry are completely lost. We performed a comprehensive survey to detect sex-specific genomic regions in the spiny rat. Sex-related genomic differences were limited to a male-specific duplication of a 17-kb unit located 430 kb upstream of Sox9 on an autosome. Hi-C analysis using male spiny rat cells showed the duplicated region has potential chromatin interaction with Sox9. The duplicated unit harbored a 1,262-bp element homologous to mouse enhancer 14 (Enh14), a candidate Sox9 enhancer that is functionally redundant in mice. Transgenic reporter mice showed that the spiny rat Enh14 can function as an embryonic testis enhancer in mice. Embryonic gonads of XX mice in which Enh14 was replaced by the duplicated spiny rat Enh14 showed increased Sox9 expression and decreased Foxl2 expression. We propose that male-specific duplication of this Sox9 enhancer substituted for Sry function, defining a novel Y chromosome in the spiny rat.

Genome-wide identification of ZmMYC2 binding sites and target genes in maize.

BACKGROUND: Jasmonate (JA) is the important phytohormone to regulate plant growth and adaption to stress signals. MYC2, an bHLH transcription factor, is the master regulator of JA signaling. Although MYC2 in maize has been identified, its function remains to be clarified. RESULTS: To understand the function and regulatory mechanism of MYC2 in maize, the joint analysis of DAP-seq and RNA-seq is conducted to identify the binding sites and target genes of ZmMYC2. A total of 3183 genes are detected both in DAP-seq and RNA-seq data, potentially as the directly regulating genes of ZmMYC2. These genes are involved in various biological processes including plant growth and stress response. Besides the classic cis-elements like the G-box and E-box that are bound by MYC2, some new motifs are also revealed to be recognized by ZmMYC2, such as nGCATGCAnn, AAAAAAAA, CACGTGCGTGCG. The binding sites of many ZmMYC2 regulating genes are identified by IGV-sRNA. CONCLUSIONS: All together, abundant target genes of ZmMYC2 are characterized with their binding sites, providing the basis to construct the regulatory network of ZmMYC2 and better understanding for JA signaling in maize.

Genome-wide identification and phylogenetic and expression pattern analyses of EPF/EPFL family genes in the Rye (Secale cereale L.).

Rye (Secale cereale L.) is one of the major cereal crop species in the Triticeae family and is known to be most tolerant to diverse abiotic stresses, such as cold, heat, osmotic, and salt stress. The EPIDERMAL PATTERNING FACTOR (EPF) and EPF-LIKE (EPFL) families of small secreted peptides act to regulate many aspects of plant growth and development; however, their functions are not widely characterized in rye. In this study, we identified 12 ScEPF/EPFL genes, which can be divided into six groups and are evenly distributed on six rye chromosomes. Further examination of the gene structure and protein conservation motifs of EPF/EPFL family members demonstrated the high conservation of the ScEPF/EPFL sequence. Interactions between ScEPF/EPFL proteins and promoters containing hormone- and stress-responsive cis-acting elements suggest that the regulation of ScEPF/EPFL expression is complex. Expression profiling analyses revealed that ScEPF/EPFL genes exhibited tissue-specific expression patterns. Notably, ScEPFL1,ScEPFL7, ScEPFL9, and ScEPFL10 displayed significantly higher expression levels in spikelets compared to other tissues. Moreover, fluorescence quantification experiments demonstrated that these genes exhibited distinct expression patterns in response to various stress conditions, suggesting that each gene plays a unique role in stress signaling pathways. Our research findings provide a solid basis for further investigation into the functions of ScEPF/EPFLs. Furthermore, these genes can serve as potential candidates for breeding stress-resistant rye varieties and improving production yields.

Unraveling the role of PlARF2 in regulating deed formancy in Paeonia lactiflora.

MAIN CONCLUSION: PlARF2 can positively regulate the seed dormancy in Paeonia lactiflora Pall. and bind the RY cis-element. Auxin, a significant phytohormone influencing seed dormancy, has been demonstrated to be regulated by auxin response factors (ARFs), key transcriptional modulators in the auxin signaling pathway. However, the role of this class of transcription factors (TFs) in perennials with complex seed dormancy mechanisms remains largely unexplored. Here, we cloned and characterized an ARF gene from Paeonia lactiflora, named PlARF2, which exhibited differential expression levels in the seeds during the process of seed dormancy release. The deduced amino acid sequence of PlARF2 had high homology with those of other plants and contained typical conserved Auxin_resp domain of the ARF family. Phylogenetic analysis revealed that PlARF2 was closely related to VvARF3 in Vitis vinifera. The subcellular localization and transcriptional activation assay showed that PlARF2 is a nuclear protein possessing transcriptional activation activity. The expression levels of dormancy-related genes in transgenic callus indicated that PlARF2 was positively correlated with the contents of PlABI3 and PlDOG1. The germination assay showed that PlARF2 promoted seed dormancy. Moreover, TF Centered Yeast one-hybrid assay (TF-Centered Y1H), electrophoretic mobility shift assay (EMSA) and dual-luciferase reporter assay analysis (Dual-Luciferase) provided evidence that PlARF2 can bind to the 'CATGCATG' motif. Collectively, our findings suggest that PlARF2, as TF, could be involved in the regulation of seed dormancy and may act as a repressor of germination.

A rapid and scalable approach to build synthetic repetitive hormone-responsive promoters.

Advancement of DNA-synthesis technologies has greatly facilitated the development of synthetic biology tools. However, high-complexity DNA sequences containing tandems of short repeats are still notoriously difficult to produce synthetically, with commercial DNA synthesis companies usually rejecting orders that exceed specific sequence complexity thresholds. To overcome this limitation, we developed a simple, single-tube reaction method that enables the generation of DNA sequences containing multiple repetitive elements. Our strategy involves commercial synthesis and PCR amplification of padded sequences that contain the repeats of interest, along with random intervening sequence stuffers that include type IIS restriction enzyme sites. GoldenBraid molecular cloning technology is then employed to remove the stuffers, rejoin the repeats together in a predefined order, and subclone the tandem(s) in a vector using a single-tube digestion-ligation reaction. In our hands, this new approach is much simpler, more versatile and efficient than previously developed solutions to this problem. As a proof of concept, two different phytohormone-responsive, synthetic, repetitive proximal promoters were generated and tested in planta in the context of transcriptional reporters. Analysis of transgenic lines carrying the synthetic ethylene-responsive promoter 10x2EBS-S10 fused to the GUS reporter gene uncovered several developmentally regulated ethylene response maxima, indicating the utility of this reporter for monitoring the involvement of ethylene in a variety of physiologically relevant processes. These encouraging results suggest that this reporter system can be leveraged to investigate the ethylene response to biotic and abiotic factors with high spatial and temporal resolution.

T-hairpin structure found in the RNA element involved in piRNA biogenesis.

PIWI-interacting RNAs (piRNAs) repress transposons to protect the germline genome from DNA damage caused by transposon transposition. In Drosophila, the Traffic jam (Tj) mRNA is consumed to produce piRNA in its 3'-UTR. A cis element located within the 3'-UTR, Tj-cis, is necessary for piRNA biogenesis. In this study, we analyzed the structure of the Tj-cis RNA, a 100-nt RNA corresponding to the Tj-cis element, by the SHAPE and NMR analyses and found that a stable hairpin structure formed in the 5' half of the Tj-cis RNA. The tertiary structure of the 16-nt stable hairpin was analyzed by NMR, and a novel stem-loop structure, the T-hairpin, was found. In the T-hairpin, four uridine residues are exposed to the solvent, suggesting that this stem-loop is the target of Yb protein, a Tudor domain-containing piRNA biogenesis factor. The piRNA biogenesis assay showed that both the T-hairpin and the 3' half are required for the function of the Tj-cis element, suggesting that both the T-hairpin and the 3' half are recognized by Yb protein.

A novel natural variation in the promoter of GmCHX1 regulates conditional gene expression to improve salt tolerance in soybean.

Identification and characterization of soybean germplasm and gene(s)/allele(s) for salt tolerance is an effective way to develop improved varieties for saline soils. Previous studies identified GmCHX1 (Glyma03g32900) as a major salt tolerance gene in soybean, and two main functional variations were found in the promoter region (148/150 bp insertion) and the third exon with a retrotransposon insertion (3.78 kb). In the current study, we identified four salt-tolerant soybean lines, including PI 483460B (Glycine soja), carrying the previously identified salt-sensitive variations at GmCHX1, suggesting new gene(s) or new functional allele(s) of GmCHX1 in these soybean lines. Subsequently, we conducted quantitative trait locus (QTL) mapping in a recombinant-inbred line population (Williams 82 (salt-sensitive) × PI 483460B) to identify the new salt tolerance loci/alleles. A new locus, qSalt_Gm18, was mapped on chromosome 18 associated with leaf scorch score. Another major QTL, qSalt_Gm03, was identified to be associated with chlorophyll content ratio and leaf scorch score in the same chromosomal region of GmCHX1 on chromosome 3. Novel variations in a STRE (stress response element) cis-element in the promoter region of GmCHX1 were found to regulate the salt-inducible expression of the gene in these four newly identified salt-tolerant lines including PI 483460B. This new allele of GmCHX1 with salt-inducible expression pattern provides an energy cost efficient (conditional gene expression) strategy to protect soybean yield in saline soils without yield penalty under non-stress conditions. Our results suggest that there might be no other major salt tolerance locus similar to GmCHX1 in soybean germplasm, and further improvement of salt tolerance in soybean may rely on gene-editing techniques instead of looking for natural variations.

The transcription factor ThDOF8 binds to a novel cis-element and mediates molecular responses to salt stress in Tamarix hispida.

Salt stress is a common abiotic factor that restricts plant growth and development. As a halophyte, Tamarix hispida is a good model plant for exploring salt-tolerance genes and regulatory mechanisms. DNA-binding with one finger (DOF) is an important transcription factor (TF) that influences and controls various signaling substances involved in diverse biological processes related to plant growth and development, but the regulatory mechanisms of DOF TFs in response to salt stress are largely unknown in T. hispida. In the present study, a newly identified Dof gene, ThDOF8, was cloned from T. hispida, and its expression was found to be induced by salt stress. Transient overexpression of ThDOF8 enhanced T. hispida salt tolerance by enhancing proline levels, and increasing the activities of the antioxidant enzymes superoxide dismutase (SOD) and peroxidase (POD). These results were also verified in stably transformed Arabidopsis. Results from TF-centered yeast one-hybrid (Y1H) assays and EMSAs showed that ThDOF8 binds to a newly identified cis-element (TGCG). Expression profiling by gene chip analysis identified four potential direct targets of ThDOF8, namely the cysteine-rich receptor-like kinases genes, CRK10 and CRK26, and two glutamate decarboxylase genes, GAD41, and GAD42, and these were further verified by ChIP-quantitative-PCR, EMSAs, Y1H assays, and ß-glucuronidase enzyme activity assays. ThDOF8 can bind to the TGCG element in the promoter regions of its target genes, and transient overexpression of ThCRK10 also enhanced T. hispida salt tolerance. On the basis of our results, we propose a new regulatory mechanism model, in which ThDOF8 binds to the TGCG cis-element in the promoter of the target gene CRK10 to regulate its expression and improve salt tolerance in T. hispida. This study provides a basis for furthering our understanding the role of DOF TFs and identifying other downstream candidate genes that have the potential for improving plant salt tolerance via molecular breeding.

Two independent cis-acting elements are required for the guard cell-specific expression of SCAP1, which is essential for late stomatal development.

Regulating the stomatal aperture to adapt to environmental changes is critical for plants as stomatal guard cells are responsible for gas exchange between plants and the atmosphere. We previously showed that a plant-specific DNA-binding with one finger (Dof)-type transcription factor, SCAP1, functions as a key regulator in the final stages of guard cell differentiation. In the present study, we performed deletion and gain-of-function analyses with the 5' flanking region of SCAP1 to identify the regulatory region controlling the guard cell-specific expression of SCAP1. The results revealed that two cis-acting elements, 5'-CACGAGA-3' and 5'-CACATGTTTCCC-3', are crucial for the guard cell-specific expression of SCAP1. Consistently, when an 80-bp promoter region including these two cis-elements was fused to a gene promoter that is not active in guard cells, it functioned as a promoter that directed gene expression in guard cells. Furthermore, the promoter region of HT1 encoding the central regulator of stomatal CO2 signaling was also found to contain a 5'-CACGAGA-3' sequence, which was confirmed to function as a cis-element necessary for guard cell-specific expression of HT1. These findings suggest the existence of a novel transcriptional regulatory mechanism that synchronously promotes the expression of multiple genes required for the stomatal maturation and function.

Genome-wide identification and characterization of the lettuce GASA family in response to abiotic stresses.

BACKGROUND: Lettuce is one of the most extensively farmed vegetables in the world, and it prefers cool growing conditions. High temperatures promote premature bolt formation, reducing quality and yield. The gibberellic acid-stimulated Arabidopsis (GASA) family genes play critical roles in plant growth, development, and stress responses. However, the biological functions of GASA proteins in lettuce have yet to be thoroughly investigated. RESULTS: Using genome-wide analysis, 20 GASAs were identified in lettuce including, three groups of LsGASA proteins based on the phylogenetic analysis. Except for one, all GASA proteins included a conserved GASA domain with 12 cysteine residues. Cis-element analysis showed that LsGASAs were closely associated with light, phytohormones, and stress resistance. Five segmental and three tandem duplication events were observed in the LsGASA family based on duplication analysis. GASA synteny analysis among lettuce, Arabidopsis, tobacco, and rice revealed that LsGASA5 is highly collinear with all species. Six of the 20 LsGASA showed increased expression patterns at specific time points in the shoot apical meristem when subjected to heat stress. According to gene expression analysis, the majority of GASA were highly expressed in flowers compared to other organs, and six GASA exhibited highly increased expression levels in response to NaCl, abscisic acid, and gibberellin treatment. Furthermore, LsGASA proteins are predominantly found in the plasma membrane and/or the cytosol. CONCLUSIONS: This study provides a comprehensive characterization of LsGASA genes for their diversity and biological functions. Moreover, our results will be useful for further studies on the function of lettuce GASA in abiotic stress- and heat-induced bolting signaling.

Genome-Wide Identification, Characterization and Expression Analysis of the TaDUF724 Gene Family in Wheat (Triticum aestivum).

Unknown functional domain (DUF) proteins constitute a large number of functionally uncharacterized protein families in eukaryotes. DUF724s play crucial roles in plants. However, the insight understanding of wheat TaDUF724s is currently lacking. To explore the possible function of TaDUF724s in wheat growth and development and stress response, the family members were systematically identified and characterized. In total, 14 TaDUF724s were detected from a wheat reference genome; they are unevenly distributed across the 11 chromosomes, and, according to chromosome location, they were named TaDUF724-1 to TaDUF724-14. Evolution analysis revealed that TaDUF724s were under negative selection, and fragment replication was the main reason for family expansion. All TaDUF724s are unstable proteins; most TaDUF724s are acidic and hydrophilic. They were predicted to be located in the nucleus and chloroplast. The promoter regions of TaDUF724s were enriched with the cis-elements functionally associated with growth and development, as well as being hormone-responsive. Expression profiling showed that TaDUF724-9 was highly expressed in seedings, roots, leaves, stems, spikes and grains, and strongly expressed throughout the whole growth period. The 12 TaDUF724 were post-transcription regulated by 12 wheat MicroRNA (miRNA) through cleavage and translation. RT-qPCR showed that six TaDUF724s were regulated by biological and abiotic stresses. Conclusively, TaDUF724s were systematically analyzed using bioinformatics methods, which laid a theoretical foundation for clarifying the function of TaDUF724s in wheat.

Mechanisms and consequences of subcellular RNA localization across diverse cell types.

Essentially all cells contain a variety of spatially restricted regions that are important for carrying out specialized functions. Often, these regions contain specialized transcriptomes that facilitate these functions by providing transcripts for localized translation. These transcripts play a functional role in maintaining cell physiology by enabling a quick response to changes in the cellular environment. Here, we review how RNA molecules are trafficked within cells, with a focus on the subcellular locations to which they are trafficked, mechanisms that regulate their transport and clinical disorders associated with misregulation of the process.

The R2R3-MYB gene family in Cicer arietinum: genome-wide identification and expression analysis leads to functional characterization of proanthocyanidin biosynthesis regulators in the seed coat.

MAIN CONCLUSION: We identified 119 typical CaMYB encoding genes and reveal the major components of the proanthocyanidin regulatory network. CaPARs emerged as promising targets for genetic engineering toward improved agronomic traits in C. arietinum. Chickpea (Cicer arietinum) is among the eight oldest crops and has two main types, i.e., desi and kabuli, whose most obvious difference is the color of their seeds. We show that this color difference is due to differences in proanthocyanidin content of seed coats. Using a targeted approach, we performed in silico analysis, metabolite profiling, molecular, genetic, and biochemical studies to decipher the transcriptional regulatory network involved in proanthocyanidin biosynthesis in the seed coat of C. arietinum. Based on the annotated C. arietinum reference genome sequence, we identified 119 typical CaMYB encoding genes, grouped in 32 distinct clades. Two CaR2R3-MYB transcription factors, named CaPAR1 and CaPAR2, clustering with known proanthocyanidin regulators (PARs) were identified and further analyzed. The expression of CaPAR genes correlated well with the expression of the key structural proanthocyanidin biosynthesis genes CaANR and CaLAR and with proanthocyanidin levels. Protein-protein interaction studies suggest the in vivo interaction of CaPAR1 and CaPAR2 with the bHLH-type transcription factor CaTT8. Co-transfection analyses using Arabidopsis thaliana protoplasts showed that the CaPAR proteins form a MBW complex with CaTT8 and CaTTG1, able to activate the promoters of CaANR and CaLAR in planta. Finally, transgenic expression of CaPARs in the proanthocyanidin-deficient A. thaliana mutant tt2-1 leads to complementation of the transparent testa phenotype. Taken together, our results reveal main components of the proanthocyanidin regulatory network in C. arietinum and suggest that CaPARs are relevant targets of genetic engineering toward improved agronomic traits.

Interplay between R2R3 MYB-type activators and repressors regulates proanthocyanidin biosynthesis in banana (Musa acuminata).

Proanthocyanidins are oligomeric flavonoids that promote plant disease resistance and benefit human health. Banana is one of the world's most extensively farmed crops and its fruit pulp contain proanthocyanidins. However, the transcriptional regulatory network that fine tunes proanthocyanidin biosynthesis in banana remains poorly understood. We characterised two proanthocyanidin-specific R2R3 MYB activators (MaMYBPA1-MaMYBPA2) and four repressors (MaMYBPR1-MaMYBPR4) to elucidate the mechanisms underlying the transcriptional regulation of proanthocyanidin biosynthesis in banana. Heterologous expression of MaMYBPA1 and MaMYBPA2 partially complemented the Arabidopsis thaliana proanthocyanidin-deficient transparent testa2 mutant. MaMYBPA1 and MaMYBPA2 interacted physically with MaMYCs to transactivate anthocyanin synthase, leucoanthocyanidin reductase, and anthocyanidin reductase genes in vitro and form functional MYB-bHLH-WD Repeat (MBW) complexes with MaTTG1 to transactivate these promoters in vivo. Overexpression of MaMYBPAs alone or with MaMYC in banana fruits induced proanthocyanidin accumulation and transcription of proanthocyanidin biosynthesis-related genes. MaMYBPR repressors are also shown to interact with MaMYCs forming repressing MBW complexes, and diminished proanthocyanidin accumulation. Interestingly overexpression of MaMYBPA induces the expression of MaMYBPR, indicating an agile regulation of proanthocyanidin biosynthesis through the formation of competitive MBW complexes. Our results reveal regulatory modules of R2R3 MYB- that fine tune proanthocyanidin biosynthesis and offer possible targets for genetic manipulation for nutritional improvement of banana.

Genome-Wide Identification and Characterization of the RCI2 Gene Family in Allotetraploid Brassica napus Compared with Its Diploid Progenitors.

Brassica napus and its diploid progenitors (B. rapa and B. oleracea) are suitable for studying the problems associated with polyploidization. As an important anti-stress protein, RCI2 proteins widely exist in various tissues of plants, and are crucial to plant growth, development, and stress response. In this study, the RCI2 gene family was comprehensively identified and analyzed, and 9, 9, and 24 RCI2 genes were identified in B. rapa, B. oleracea, and B. napus, respectively. Phylogenetic analysis showed that all of the identified RCI2 genes were divided into two groups, and further divided into three subgroups. Ka/Ks analysis showed that most of the identified RCI2 genes underwent a purifying selection after the duplication events. Moreover, gene structure analysis showed that the structure of RCI2 genes is largely conserved during polyploidization. The promoters of the RCI2 genes in B. napus contained more cis-acting elements, which were mainly involved in plant development and growth, plant hormone response, and stress responses. Thus, B. napus might have potential advantages in some biological aspects. In addition, the changes of RCI2 genes during polyploidization were also discussed from the aspects of gene number, gene structure, gene relative location, and gene expression, which can provide reference for future polyploidization analysis.

Cis-element amplified polymorphism (CEAP), a novel promoter- and gene-targeted molecular marker of plants.

In this study, we selected eight cis-elements: AAAG, ACGTG, CCGA, ACTCAT, GGTCA, TATCC, TGAC and GATAA, which are closely related to plant growth and development, signal transduction and stress response. The CEAP primers were 18 nucleotides long and consisted of a central cis-element nucleotide core flanked by a filler sequence at the 5' end and di- or tri-nucleotides at the 3' end. A total of two hundred and twenty-four primers were developed, and the PCR procedure consisted of 5 cycles of low-temperature annealing and 35 subsequent cycles of annealing at 50°C. The PCR products are electrophoretically separated by 1.8-2.3% agarose. The polymorphism of the CEAP marker was amplified in eight mango (Mangifera indica L.) species. The results showed that the CEAP primers could amplify clear, repeatable bands in mango and combine at least four cis-elements from which a large number of bands were amplified and six highly polymorphic primers for each cis-element can reach an accurate clustering result. The results of CEAP marker assays compared with ISSR, CBDP and iPBS marker assays showed that CEAP marker was better than the other three markers in the number of fragment bands, H and I indexes. In addition, we also tested the CEAP markers in rice, tomato, potato, wax gourd, citrus and longan and the results showed that the CEAP marker assay could amplify clear polymorphic bands in different species. Our results indicate that the CEAP markers could be universally used in different species for genetic diversity analysis, relationship analysis, and marker-assisted selection for breeding. Supplementary Information: The online version contains supplementary material available at 10.1007/s12298-022-01212-5.

Dual Eigen-modules of Cis-Element Regulation Profiles and Selection of Cognition-Language Eigen-direction along Evolution in Hominidae.

To understand the genomic basis accounting for the phenotypic differences between human and apes, we compare the matrices consisting of the cis-element frequencies in the proximal regulatory regions of their genomes. One such frequency matrix is represented by a robust singular value decomposition. For each singular value, the negative and positive ends of the sorted motif eigenvector correspond to the dual ends of the sorted gene eigenvector, respectively, comprising a dual eigen-module defined by cis-regulatory element frequencies (CREF). The CREF eigen-modules at levels 1, 2, 3, and 6 are highly conserved across humans, chimpanzees, and orangutans. The key biological processes embedded in the top three CREF eigen-modules are reproduction versus embryogenesis, fetal maturation versus immune system, and stress responses versus mitosis. Although the divergence at the nucleotide level between the chimpanzee and human genome was small, their cis-element frequency matrices crossed a singularity point, at which the fourth and fifth singular values were identical. The CREF eigen-modules corresponding to the fourth and fifth singular values were reorganized along the evolution from apes to human. Interestingly, the fourth sorted gene eigenvector encodes the phenotypes unique to human such as long-term memory, language development, and social behavior. The number of motifs present on Alu elements increases substantially at the fourth level. The motif analysis together with the cases of human-specific Alu insertions suggests that mutations related to Alu elements play a critical role in the evolution of the human-phenotypic gene eigenvector.

Characteristic and evolution of HAT and HDAC genes in Gramineae genomes and their expression analysis under diverse stress in Oryza sativa.

MAIN CONCLUSION: Comprehensive characterization of Gramineae HATs and HDACs reveals their conservation and variation. The recent WGD/SD gene pairs in the CBP and RPD/HDA1 gene family may confer specific adaptive evolutionary changes. Expression of OsHAT and OsHDAC genes provides a new vision in different aspects of development and response to diverse stress. The histone acetylase (HAT) and histone deacetylase (HDAC) have been proven to be tightly linked to play a crucial role in plant growth, development and response to abiotic stress by regulating histone acetylation levels. However, the evolutionary dynamics and functional differentiation of HATs and HDACs in Gramineae remain largely unclear. In the present study, we identified 37 HAT genes and 110 HDAC genes in seven Gramineae genomes by a detailed analysis. Phylogenetic trees of these HAT and HDAC proteins were constructed to illustrate evolutionary relationship in Gramineae. Gene structure, protein property and protein motif composition illustrated the conservation and variation of HATs and HDACs in Gramineae. Gene duplication analysis suggested that recent whole genome duplication (WGD)/segmental duplication (SD) events contributed to the diversification of the CBP and RPD3/HDA1 gene family in Gramineae. Furthermore, promoter cis-element prediction indicated that OsHATs and OsHDACs were likely functional proteins and involved in various signaling pathways. Expression analysis by RNA-seq data showed that all OsHAT and OsHDAC genes were expressed in different tissues or development stages, revealing that they were ubiquitously expressed. In addition, we found that their expression patterns were altered in response to cold, drought, salt, light, abscisic acid (ABA), and indole-3-acetic acid (IAA) treatments. These findings provide the basis for further identification of candidate OsHAT and OsHDAC genes that may be utilized in regulating growth and development and improving crop tolerance to abiotic stress.

Genome-wide characterization and analysis of the CCT motif family genes in soybean (Glycine max).

MAIN CONCLUSION: Soybean possesses 19 CMF genes which mainly arose from duplication events. Their features and motifs are highly conserved but transcriptional data indicated functional diversity in metabolism and stress responses. CCT [for CONSTANS, CONSTANS-like (CO-like), and timing of CAB expression1 (TOC1)] domain-containing genes play important roles in regulating flowering, plant growth, and grain yield and are also involved in stress responses. The CMF (CCT motif family) genes, included in the CCT family, contain a single CCT domain as the only identifiable domain in their predicted protein sequence and are interesting targets for breeding programs. In this study, we identified 19 putative GmCMF genes, based on the latest soybean (Glycine max) genome annotation. The predicted GmCMF proteins were characterized based on conserved structural features, and a phylogenetic tree was constructed including all CMF proteins from rice and Arabidopsis as representative examples of the monocotyledonous (monocot) and dicotyledonous (dicot) plants, respectively. High similarities in the conserved motifs of the protein sequences and the gene structures were found. In addition, by analyzing the CMF gene family in soybean, we identified seven pairs of genes that originated from segmental chromosomal duplication events attributable to the most recent whole-genome duplication (WGD) event in the Glycine lineage. Expression analysis of GmCMF genes in various tissues and after specific treatments demonstrated tissue and stress-response specific differential expression. Gene expression analysis was complemented by the identification of putative cis-elements present in the promoter regions of the genes through a bioinformatics approach, using the existing soybean reference genome sequence and gene models. Co-functional networks inferred from distinct types of genomics data-including microarrays and RNA-seq samples from soybean-revealed that GmCMF genes might play crucial roles in metabolism and transport processes. The results of this study, the first systematic analysis of the soybean CCT gene family, can serve as a strong foundation for further elucidation of their physiological functions and biological roles.