Importance of gene duplication in the evolution of genomic imprinting revealed by molecular evolutionary analysis of the type I MADS-box gene family in Arabidopsis species

Structural evidence for MADS-box type I family expansion seen in new assemblies of Arabidopsis arenosa and A. lyrata

Arabidopsis thaliana diverged from A. arenosa and A. lyrata at least 6 million years ago. The three species differ by genome-wide polymorphisms and morphological traits. The species are to a high degree reproductively isolated, but hybridization barriers are incomplete. A special type of hybridization barrier is based on the triploid endosperm of the seed, where embryo lethality is caused by endosperm failure to support the developing embryo. The MADS-box type I family of transcription factors is specifically expressed in the endosperm and has been proposed to play a role in endosperm-based hybridization barriers. The gene family is well known for its high evolutionary duplication rate, as well as being regulated by genomic imprinting. Here we address MADS-box type I gene family evolution and the role of type I genes in the context of hybridization. Using two de-novo assembled and annotated chromosome-level genomes of A. arenosa and A. lyrata ssp. petraea we analyzed the MADS-box type I gene family in Arabidopsis to predict orthologs, copy number, and structural genomic variation related to the type I loci. Our findings were compared to gene expression profiles sampled before and after the transition to endosperm cellularization in order to investigate the involvement of MADS-box type I loci in endosperm-based hybridization barriers. We observed substantial differences in type-I expression in the endosperm of A. arenosa and A. lyrata ssp. petraea, suggesting a genetic cause for the endosperm-based hybridization barrier between A. arenosa and A. lyrata ssp. petraea.

Genome-Wide Analysis and the Expression Pattern of the MADS-Box Gene Family in Bletilla striata

Bletilla striata (Thunb. ex A. Murray) Rchb. f., a species of the perennial herb Orchidaceae, has potent anti-inflammatory and antiviral biological activities. MADS-box transcription factors play critical roles in the various developmental processes of plants. Although this gene family has been extensively investigated in many species, it has not been analyzed for B. striata. In total, 45 MADS-box genes were identified from B. striata in this study, which were classified into five subfamilies (Mδ, MIKC, Mα, Mβ, and Mγ). Meanwhile, the highly correlated protein domains, motif compositions, and exon-intron structures of BsMADSs were investigated according to local B. striata databases. Chromosome distribution and synteny analyses revealed that segmental duplication and homologous exchange were the main BsMADSs expansion mechanisms. Further, RT-qPCR analysis revealed that BsMADSs had different expression patterns in response to various stress treatments. Our results provide a potential theoretical basis for further investigation of the functions of MADS genes during the growth of B. striata.

Polyploidization events shaped the transcription factor repertoires in legumes (Fabaceae)

Transcription factors (TFs) are essential for plant growth and development. Several legumes (e.g. soybean) are rich sources of protein and oil and have great economic importance. Here we report a phylogenomic analysis of TF families in legumes and their potential association with important traits (e.g. nitrogen fixation). We used TF DNA-binding domains to systematically screen the genomes of 15 leguminous and five non-leguminous species. Transcription factor orthologous groups (OGs) were used to estimate OG sizes in ancestral nodes using a gene birth-death model, which allowed the identification of lineage-specific expansions. The OG analysis and rate of synonymous substitutions show that major TF expansions are strongly associated with whole-genome duplication (WGD) events in the legume (approximately 58 million years ago) and Glycine (approximately 13 million years ago) lineages, which account for a large fraction of the Phaseolus vulgaris and Glycine max TF repertoires. Of the 3407 G. max TFs, 1808 and 676 have homeologs within single syntenic regions in Phaseolus vulgaris and Vitis vinifera, respectively. We found a trend for TFs expanded in legumes to be preferentially transcribed in roots and nodules, supporting their recruitment early in the evolution of nodulation in the legume clade. Some families also showed count differences between G. max and the wild soybean Glycine soja, including genes located within important quantitative trait loci. Our findings strongly support the roles of two WGDs in shaping the TF repertoires in the legume and Glycine lineages, and these are probably related to important aspects of legume and soybean biology.

Identification and expression analysis of the small auxin-up RNA (SAUR) gene family in apple by inducing of auxin

The small auxin-up RNA (SAUR) family plays a vital role in the regulation of plant growth and development. We identified 80 MdSAUR genes in this study. Phylogenetic analysis indicated that the SAUR proteins from Arabidopsis, rice, and apple were divided into six groups. Of the 80 MdSAURs, 71 were randomly distributed along the 17 chromosomes, while the remaining genes were located along unassigned scafoolds. Among them, a comprehensive overview of SAUR gene family is presented, including gene structures, chromosome locations, duplication and selection pressure analyses, synteny and promoter analyses, and protein interaction. The expression profiles based on microarray data found that 80 genes showed increased expression levels in at least one tissue including seed, seedling, root, stem, leaf, flower, fruit 100daa, and harvested fruit. MdSAUR7 possibly regulate the development of flower organs, and MdSAUR15, MdSAUR24, and MdSAUR80 promote the growth of fruits by regulating cell division. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis indicated the expression levels of 79 MdSAUR genes in leaves under exogenous IAA treatment. MdSAUR4, MdSAUR22, MdSAUR37, MdSAUR38, MdSAUR49, and MdSAUR54 were up-regulated after IAA treatment compared with the control, indicating that they may play specific roles in the IAA signaling transduction pathway. This work provided a foundation for further investigations for the functional analyses of SAURs in apple.

Genetic variation and temperature affects hybrid barriers during interspecific hybridization

Genomic imprinting regulates parent-specific transcript dosage during seed development and is mainly confined to the endosperm. Elucidation of the function of many imprinted genes has been hampered by the lack of corresponding mutant phenotypes, and the role of imprinting is mainly associated with genome dosage regulation or allocation of resources. Disruption of imprinted genes has also been suggested to mediate endosperm-based post-zygotic hybrid barriers depending on genetic variation and gene dosage. Here, we have analyzed the conservation of a clade from the MADS-box type I class transcription factors in the closely related species Arabidopsis arenosa, A. lyrata, and A. thaliana, and show that AGL36-like genes are imprinted and maternally expressed in seeds of Arabidopsis species and in hybrid seeds between outbreeding species. In hybridizations between outbreeding and inbreeding species the paternally silenced allele of the AGL36-like gene is reactivated in the hybrid, demonstrating that also maternally expressed imprinted genes are perturbed during hybridization and that such effects on imprinted genes are specific to the species combination. Furthermore, we also demonstrate a quantitative effect of genetic diversity and temperature on the strength of the post-zygotic hybridization barrier. Markedly, a small decrease in temperature during seed development increases the survival of hybrid F1 seeds, suggesting that abiotic and genetic parameters play important roles in post-zygotic species barriers, pointing at evolutionary scenarios favoring such effects. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA562212. All sequences generated in this study have been deposited in the National Center for Biotechnology Information Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra/) with project number PRJNA562212.

Genome-Wide Analysis of Parent-of-Origin Allelic Expression in Endosperms of Brassicaceae Species, Brassica rapa

Uniparental gene expression, observed in both animals and plants, is termed genomic imprinting. Genomic imprinting is a well-known epigenetic phenomenon regulated through epigenetic modifications such as DNA methylation and histone modifications. Recent genome-wide studies of endosperm transcription have revealed the rapid change of imprinted genes between species, suggesting the flexibility of this phenomenon. Although the functional significance and evolutionary trends of imprinted genes are still obscure, it can be clarified by inter-species comparisons. In this study, we analyzed the pattern of genomic imprinting in Brassica rapa, a species related to Arabidopsis thaliana. Compared with the ancient karyotype of A. thaliana and B. rapa, B. rapa has a triplicated genome. Many imprinted genes, beyond the estimated number previously reported in other species, were observed. Several imprinted genes have been conserved among species in Brassicaceae. We also observed rapid molecular evolution of imprinted genes compared to non-imprinted genes in B. rapa. Especially, imprinted gene overlapping between species showed more rapid molecular evolution and preferential expression in endosperms. It may imply that a small number of imprinted genes have retained functional roles among diverged species and have been the target of natural selection.

Genome-wide analysis of poplar SAUR gene family and expression profiles under cold, polyethylene glycol and indole-3-acetic acid treatments

Small auxin-up RNA (SAUR) proteins play an important role in the regulation of plant growth and development. Here, we identified 105 SAUR genes and comprehensively analyzed them in Populus trichocarpa. Based on the phylogenetic relationships, the PtSAURs were classified into ten subfamilies. Of the 105 PtSAURs, 100 were randomly distributed along the nineteen chromosomes, while the remaining genes were located along unassigned scafoolds. These genes mainly evolved through segmental duplications. In total, 94 PtSAURs contained no introns, and each group had a similar conserved motif structure. A promoter analysis revealed various cis-elements related to growth, development and stress responses, and a synteny analysis established orthologous relationships among SAURs in Arabidopsis, rice, grape and poplar. The qRT-PCR and tissue expression analyses indicated that PtSAURs show different expression levels in various tissues in response to different treatments. PtSAUR53 was located on the nuclear and plasma membrane by conducting subcellular localization analysis. This study provides a comprehensive overview of poplar SAUR proteins and a foundation for further investigations for functional analysis of SAURs in poplar growth and development. At the same time, it will be valuable to further study the poplar SAUR genes to reveal their biological effects.

Gene duplication, silencing and expression alteration govern the molecular evolution of PRC2 genes in plants

PRC2 genes were analyzed for their number of gene duplications, d_N/d_S ratios and expression patterns among Brassicaceae and Gramineae species. Although both amino acid sequences and copy number of the PRC2 genes were generally well conserved in both Brassicaceae and Gramineae species, we observed that some rapidly evolving genes experienced duplications and expression pattern changes. After multiple duplication events, all but one or two of the duplicated copies tend to be silenced. Silenced copies were reactivated in the endosperm and showed ectopic expression in developing seeds. The results indicated that rapid evolution of some PRC2 genes is initially caused by a relaxation of selective constraint following the gene duplication events. Several loci could become maternally expressed imprinted genes and acquired functional roles in the endosperm.

Frequent changes in expression profile and accelerated sequence evolution of duplicated imprinted genes in arabidopsis

Eukaryotic genomes have large numbers of duplicated genes that can evolve new functions or expression patterns by changes in coding and regulatory sequences, referred to as neofunctionalization. In flowering plants, some duplicated genes are imprinted in the endosperm, where only one allele is expressed depending on its parental origin. We found that 125 imprinted genes in Arabidopsis arose from gene duplication events during the evolution of the Brassicales. Analyses of 46 gene pairs duplicated by an ancient whole-genome duplication (alpha WGD) indicated that many imprinted genes show an accelerated rate of amino acid changes compared with their paralogs. Analyses of microarray expression data from 63 organ types and developmental stages indicated that many imprinted genes have expression patterns restricted to flowers and/or seeds in contrast to their broadly expressed paralogs. Assays of expression in orthologs from outgroup species revealed that some imprinted genes have acquired an organ-specific expression pattern restricted to flowers and/or seeds. The changes in expression pattern and the accelerated sequence evolution in the imprinted genes suggest that some of them may have undergone neofunctionalization. The imprinted genes MPC, HOMEODOMAIN GLABROUS6 (HDG6), and HDG3 are particularly interesting cases that have different functions from their paralogs. This study indicates that a large number of imprinted genes in Arabidopsis are evolutionarily recent duplicates and that many of them show changes in expression profiles and accelerated sequence evolution. Acquisition of imprinting is a mode of duplicate gene divergence in plants that is more common than previously thought.

Possible roles for polycomb repressive complex 2 in cereal endosperm

The polycomb repressive complex 2 (PRC2) is an evolutionarily conserved multimeric protein complex in both plants and animals. In contrast to animals, plants have evolved a range of different components of PRC2 and form diverse complexes that act in the control of key regulatory genes at many stages of development during the life cycle. A number of studies, particularly in the model species Arabidopsis thaliana, have highlighted the role of PRC2 and of epigenetic controls via parent-of-origin specific gene expression for endosperm development. However, recent research in cereal plants has revealed that although some components of PRC2 show evolutionary conservation with respect to parent-of-origin specific gene expression patterns, the identity of the imprinted genes encoding PRC2 components is not conserved. This disparity may reflect the facts that cereal plant genomes have undergone different patterns of duplication during evolution compared to A. thaliana and that the endosperm development program is not identical in monocots and eudicots. In this context, we focus this review on the expression of imprinted PRC2 genes and their roles in endosperm development in cereals.