Evolutionary Origins of Pseudogenes and Their Association with Regulatory Sequences in Plants

Long-read RNA sequencing can probe organelle genome pervasive transcription

40 years ago, organelle genomes were assumed to be streamlined and, perhaps, unexciting remnants of their prokaryotic past. However, the field of organelle genomics has exposed an unparallel diversity in genome architecture (i.e. genome size, structure, and content). The transcription of these eccentric genomes can be just as elaborate - organelle genomes are pervasively transcribed into a plethora of RNA types. However, while organelle protein-coding genes are known to produce polycistronic transcripts that undergo heavy posttranscriptional processing, the nature of organelle noncoding transcriptomes is still poorly resolved. Here, we review how wet-lab experiments and second-generation sequencing data (i.e. short reads) have been useful to determine certain types of organelle RNAs, particularly noncoding RNAs. We then explain how third-generation (long-read) RNA-Seq data represent the new frontier in organelle transcriptomics. We show that public repositories (e.g. NCBI SRA) already contain enough data for inter-phyla comparative studies and argue that organelle biologists can benefit from such data. We discuss the prospects of using publicly available sequencing data for organelle-focused studies and examine the challenges of such an approach. We highlight that the lack of a comprehensive database dedicated to organelle genomics/transcriptomics is a major impediment to the development of a field with implications in basic and applied science.

Loss to gain: pseudogenes in microorganisms, focusing on eubacteria, and their biological significance

Pseudogenes are defined as "non-functional" copies of corresponding parent genes. The cognition of pseudogenes continues to be refreshed through accumulating and updating research findings. Previous studies have predominantly focused on mammals, but pseudogenes have received relatively less attention in the field of microbiology. Given the increasing recognition on the importance of pseudogenes, in this review, we focus on several aspects of microorganism pseudogenes, including their classification and characteristics, their generation and fate, their identification, their abundance and distribution, their impact on virulence, their ability to recombine with functional genes, the extent to which some pseudogenes are transcribed and translated, and the relationship between pseudogenes and viruses. By summarizing and organizing the latest research progress, this review will provide a comprehensive perspective and improved understanding on pseudogenes in microorganisms. KEY POINTS: • Concept, classification and characteristics, identification and databases, content, and distribution of microbial pseudogenes are presented. • How pseudogenization contribute to pathogen virulence is highlighted. • Pseudogenes with potential functions in microorganisms are discussed.

Balanophora genomes display massively convergent evolution with other extreme holoparasites and provide novel insights into parasite-host interactions

Parasitic plants have evolved to be subtly or severely dependent on host plants to complete their life cycle. To provide new insights into the biology of parasitic plants in general, we assembled genomes for members of the sandalwood order Santalales, including a stem hemiparasite (Scurrula) and two highly modified root holoparasites (Balanophora) that possess chimaeric host-parasite tubers. Comprehensive genome comparisons reveal that hemiparasitic Scurrula has experienced a relatively minor degree of gene loss compared with autotrophic plants, consistent with its moderate degree of parasitism. Nonetheless, patterns of gene loss appear to be substantially divergent across distantly related lineages of hemiparasites. In contrast, Balanophora has experienced substantial gene loss for the same sets of genes as an independently evolved holoparasite lineage, the endoparasitic Sapria (Malpighiales), and the two holoparasite lineages experienced convergent contraction of large gene families through loss of paralogues. This unprecedented convergence supports the idea that despite their extreme and strikingly divergent life histories and morphology, the evolution of these and other holoparasitic lineages can be shaped by highly predictable modes of genome reduction. We observe substantial evidence of relaxed selection in retained genes for both hemi- and holoparasitic species. Transcriptome data also document unusual and novel interactions between Balanophora and host plants at the host-parasite tuber interface tissues, with evidence of mRNA exchange, substantial and active hormone exchange and immune responses in parasite and host.

Unravelling the due importance of pseudogenes and their resurrection in plants

The complexities of a genome are underpinned to the vast expanses of the intergenic region, which constitutes ∼97-98% of the genome. This region is essentially composed of what is colloquially referred to as the "junk DNA" and is composed of various elements like transposons, repeats, pseudogenes, etc. The latter have long been considered as dead elements merely contributing to transcriptional noise in the genome. Many studies now describe the previously unknown regulatory functions of these genes. Recent advances in the Next-generation sequencing (NGS) technologies have allowed unprecedented access to these regions. With the availability of whole genome sequences of more than 788 different plant species in past 20 years, genome annotation has become feasible like never before. Different bioinformatic pipelines are available for the identification of pseudogenes. However, still little is known about their biological functions. The functional validation of these genes remains challenging and research in this area is still in infancy, particularly in plants. CRISPR/Cas-based genome editing could provide solutions to understand the biological roles of these genes by allowing creation of precise edits within these genes. The possibility of pseudogene reactivation or resurrection as has been demonstrated in a few studies might open new avenues of genetic manipulation to yield a desirable phenotype. This review aims at comprehensively summarizing the progress made with regards to the identification of pseudogenes and understanding their biological functions in plants.

High-quality genome of a pioneer mangrove Laguncularia racemosa explains its advantages for intertidal zone reforestation

Ecological restoration of mangrove ecosystems that became susceptible to recent habitat perturbations is crucial for tropical coast conservation. The white mangrove Laguncularia racemosa, a pioneer species inhabiting intertidal environments of the Atlantic East Pacific (AEP) region, has been used for reforestation in China for decades. However, the molecular mechanisms underlying its fast growth and high adaptive potential remain unknown. Using PacBio single-molecule real-time sequencing, we completed a high-quality L. racemosa genome assembly covering 1105 Mb with scaffold N50 of 3.46 Mb. Genomic phylogeny shows that L. racemosa invaded intertidal zones during a period of global warming. Multi-level genomic convergence analyses between L. racemosa and three native dominant mangrove clades show that they experienced convergent changes in genes involved in nutrient absorption and high salinity tolerance. This may explain successful L. racemosa adaptation to stressful intertidal environments after introduction. Without recent whole-genome duplications or activated transposable elements, L. racemosa has retained many tandem gene duplications. Some of them are involved in auxin biosynthesis, intense light stress and cold stress response pathways, associated with L. racemosa's ability to grow fast under high light or cold conditions when used for reforestation. In summary, our study identifies shared mechanisms of intertidal environmental adaptation and unique genetic changes underlying fast growth in mangrove-unfavourable conditions and sheds light on the molecular mechanisms of the white mangrove utility in ecological restoration.

Genome-wide analysis of the serine carboxypeptidase-like protein family reveals Ga09G1039 is involved in fiber elongation in cotton

The Gossypium is a model genus for understanding polyploidy and the evolutionary pattern of inheritance. This study aimed to investigate the characteristics of SCPLs in different cotton species and their role in fiber development. A total of 891 genes from one typical monocot and ten dicot species were naturally divided into three classes based on phylogenetic analysis. The SCPL gene family in cotton has undergone intense purifying selection with some functional variation. Segmental duplication and whole genome duplication were shown to be the two main reasons for the increase in the number of genes during cotton evolution. The identification of Gh_SCPL genes exhibiting differential expression in particular tissues or response to environmental stimuli provides a new measure for the in-depth characterization of selected genes of importance. Ga09G1039 was involved in the developmental process of fibers and ovules, and it is significantly different from proteins from other cotton species in terms of phylogenetic, gene structure, conserved protein motifs and tertiary structure. Overexpression of Ga09G1039 significantly increased the length of stem trichomes. Ga09G1039 may be a serine carboxypeptidase protein with hydrolase activity, according to functional region, prokaryotic expression, and western blotting analysis. The results provide a comprehensive overview of the genetic basis of SCPLs in Gossypium and further our knowledge in understanding the key aspects of SCPLs in cotton with their potential role in fiber development and stress resistance.

Enhanced genome-wide association reveals the role of YABBY11-NGATHA-LIKE1 in leaf serration development of Populus

Leaf margins are complex plant morphological features that contribute to leaf shape diversity, which affects plant structure, yield, and adaptation. Although several leaf margin regulators have been identified to date, the genetic basis of their natural variation has not been fully elucidated. In this study, we profiled two distinct leaf morphology types (serrated and smooth) using the persistent homology mathematical framework (PHMF) in two poplar species (Populus tomentosa and Populus simonii, respectively). A combined genome-wide association study (GWAS) and expression quantitative trait nucleotide (eQTN) mapping were applied to create a leaf morphology control module using data from P. tomentosa and P. simonii populations. Natural variation in leaf margins was associated with YABBY11 (YAB11) transcript abundance in poplar. In P. tomentosa, PtoYAB11 carries a premature stop codon (PtoYAB11PSC), resulting in the loss of its positive regulation of NGATHA-LIKE1 (PtoNGAL-1) and RIBULOSE BISPHOSPHATE CARBOXYLASE LARGE SUBUNIT (PtoRBCL). Overexpression of PtoYAB11PSC promoted serrated leaf margins, enlarged leaves, enhanced photosynthesis, and increased biomass. Overexpression of PsiYAB11 in P. tomentosa promoted smooth leaf margins, higher stomatal density, and greater light damage repair ability. In poplar, YAB11-NGAL1 is sensitive to environmental conditions, acts as a positive regulator of leaf margin serration, and may also link environmental signaling to leaf morphological plasticity.

Differences in pseudogene evolution contributed to the contrasting flavors of turnip and Chiifu, two Brassica rapa subspecies

Pseudogenes are important resources for investigation of genome evolution and genomic diversity because they are nonfunctional but have regulatory effects that influence plant adaptation and diversification. However, few systematic comparative analyses of pseudogenes in closely related species have been conducted. Here, we present a turnip (Brassica rapa ssp. rapa) genome sequence and characterize pseudogenes among diploid Brassica species/subspecies. The results revealed that the number of pseudogenes was greatest in Brassica oleracea (CC genome), followed by B. rapa (AA genome) and then Brassica nigra (BB genome), implying that pseudogene differences emerged after species differentiation. In Brassica AA genomes, pseudogenes were distributed asymmetrically on chromosomes because of numerous chromosomal insertions/rearrangements, which contributed to the diversity among subspecies. Pseudogene differences among subspecies were reflected in the flavor-related glucosinolate (GSL) pathway. Specifically, turnip had the highest content of pungent substances, probably because of expansion of the methylthioalkylmalate synthase-encoding gene family in turnips; these genes were converted into pseudogenes in B. rapa ssp. pekinensis (Chiifu). RNA interference-based silencing of the gene encoding 2-oxoglutarate-dependent dioxygenase 2, which is also associated with flavor and anticancer substances in the GSL pathway, resulted in increased abundance of anticancer compounds and decreased pungency of turnip and Chiifu. These findings revealed that pseudogene differences between turnip and Chiifu influenced the evolution of flavor-associated GSL metabolism-related genes, ultimately resulting in the different flavors of turnip and Chiifu.

Identification of Pseudo-R genes in Vitis vinifera and characterization of their role as immunomodulators in host-pathogen interactions

INTRODUCTION

Duplication events are fundamental to co-evolution in host-pathogen interactions. Pseudogenes (Ψs) are dysfunctional paralogs of functional genes and resistance genes (Rs) in plants are the key to disarming pathogenic invasions. Thus, deciphering the roles of pseudo-R genes in plant defense is momentous.

OBJECTIVES

This study aimed to functionally characterize diverse roles of the resistance Ψs as novel gene footprints and as significant gene regulators in the grapevine genome.

METHODS

PlantPseudo pipeline and HMM-profiling identified whole-genome duplication-derived (WGD) Ψs associated with resistance genes (Ψ-Rs). Further, novel antifungal and antimicrobial peptides were characterized for fungal associations using protein-protein docking with Erysiphe necator proteins. miRNA and tasiRNA target sites and transcription factor (TF) binding sites were predicted in Ψ-Rs. Finally, differential co-expression patterns in Ψ-Rs-lncRNAs-coding genes were identified using the UPGMA method.

RESULTS

2,746 Ψ-Rs were identified from 31,032 WGD Ψs in the genome of grapevine. 69-antimicrobial and 81-antifungal novel peptides were generated from Ψ-Rs. The putative genic potential was predicted for five novel antifungal peptides which were further characterized by docking against E. necator proteins. 395 out of 527 resistance loci-specific Ψ-Rs were acting as parental gene mimics. Further, to explore the diverse roles of Ψ-Rs in plant-defense, we identified 37,026 TF-binding sites, 208 miRNA, and 99 tasiRNA targeting sites on these Ψ-Rs. 194 Ψ-Rs were exhibiting tissue-specific expression patterns. The co-expression network analysis between Ψs-lncRNA-genes revealed six out of 79 pathogen-responsive Ψ-Rs as significant during pathogen invasion.

CONCLUSIONS

Our study provides pathogen responsive Ψ-Rs integral for pathogen invasion, which will offer a useful resource for future experimental validations. In addition, our findings on novel peptide generations from Ψ-Rs offer valuable insights which can serve as a useful resource for predicting novel genes with the futuristic potential of being investigated for their bioactivities in the plant system.

Diversity and Functional Evolution of Terpene Synthases in Rosaceae

Terpenes are organic compounds and play important roles in plant development and stress response. Terpene synthases (TPSs) are the key enzymes for the biosynthesis of terpenes. For Rosaceae species, terpene composition represents a critical quality attribute, but limited information is available regarding the evolution and expansion occurring in the terpene synthases gene family. Here, we selected eight Rosaceae species with sequenced and annotated genomes for the identification of TPSs, including three Prunoideae, three Maloideae, and two Rosoideae species. Our data showed that the TPS gene family in the Rosaceae species displayed a diversity of family numbers and functions among different subfamilies. Lineage and species-specific expansion of the TPSs accompanied by frequent domain loss was widely observed within different TPS clades, which might have contributed to speciation or environmental adaptation in Rosaceae. In contrast to Maloideae and Rosoideae species, Prunoideae species owned less TPSs, with the evolution of Prunoideae species, TPSs were expanded in modern peach. Both tandem and segmental duplication significantly contributed to TPSs expansion. Ka/Ks calculations revealed that TPSs genes mainly evolved under purifying selection except for several pairs, where the divergent time indicated TPS-e clade was diverged relatively anciently. Gene function classification of TPSs further demonstrated the function diversity among clades and species. Moreover, based on already published RNA-Seq data from NCBI, the expression of most TPSs in Malus domestica, Prunus persica, and Fragaria vesca displayed tissue specificity and distinct expression patterns either in tissues or expression abundance between species and TPS clades. Certain putative TPS-like proteins lacking both domains were detected to be highly expressed, indicating the underlying functional or regulatory potentials. The result provided insight into the TPS family evolution and genetic information that would help to improve Rosaceae species quality.

Melatonin synthesis genes N-acetylserotonin methyltransferases evolved into caffeic acid O-methyltransferases and both assisted in plant terrestrialization

Terrestrialization is one of the most momentous events in the history of plant life, which leads to the subsequent evolution of plant diversity. The transition species, in this process, had to acquire a range of adaptive mechanisms to cope with the harsh features of terrestrial environments compared to that of aquatic habitat. As an ancient antioxidant, a leading regulator of ROS signaling or homeostasis, and a presumed plant master regulator, melatonin likely assisted plants transition to land and their adaption to terrestrial ecosystems. N-acetylserotonin methyltransferases (ASMT) and caffeic acid O-methyltransferases (COMT), both in the O-methyltransferase (OMT) family, catalyze the core O-methylation reaction in melatonin biosynthesis. How these two enzymes with close relevance evolved in plant evolutionary history and whether they participated in plant terrestrialization remains unknown. Using combined phylogenetic evidence and protein structure analysis, it is revealed that COMT likely evolved from ASMT by gene duplication and subsequent divergence. Newly emergent COMT gained a significantly higher ASMT activity to produce greater amounts of melatonin for immobile plants to acclimate to the stressful land environments after evolving from the more environmentally-stable aquatic conditions. The COMT genes possess more conserved substrate-binding sites at the amino acid level and more open protein conformation compared to ASMT, and getting a new function to catalyze the lignin biosynthesis. This development directly contributed to the dominance of vascular plants among the Earth's flora and prompted plant colonization of land. Thus, ASMT, together with its descendant COMT, might play key roles in plant transition to land. The current study provides new insights into plant terrestrialization with gene duplication contributing to this process along with well-known horizontal gene transfer.

Patterns and Processes of Diploidization in Land Plants

Most land plants are now known to be ancient polyploids that have rediploidized. Diploidization involves many changes in genome organization that ultimately restore bivalent chromosome pairing and disomic inheritance, and resolve dosage and other issues caused by genome duplication. In this review, we discuss the nature of polyploidy and its impact on chromosome pairing behavior. We also provide an overview of two major and largely independent processes of diploidization: cytological diploidization and genic diploidization/fractionation. Finally, we compare variation in gene fractionation across land plants and highlight the differences in diploidization between plants and animals. Altogether, we demonstrate recent advancements in our understanding of variation in the patterns and processes of diploidization in land plants and provide a road map for future research to unlock the mysteries of diploidization and eukaryotic genome evolution.

Genome-Wide Analysis of Coding and Non-coding RNA Reveals a Conserved miR164-NAC-mRNA Regulatory Pathway for Disease Defense in Populus

MicroRNAs (miRNAs) contribute to plant defense responses by increasing the overall genetic diversity; however, their origins and functional importance in plant defense remain unclear. Here, we employed Illumina sequencing technology to assess how miRNA and messenger RNA (mRNA) populations vary in the Chinese white poplar (Populus tomentosa) during a leaf black spot fungus (Marssonina brunnea) infection. We sampled RNAs from infective leaves at conidia germinated stage [12 h post-inoculation (hpi)], infective vesicles stage (24 hpi), and intercellular infective hyphae stage (48 hpi), three essential stages associated with plant colonization and biotrophic growth in M. brunnea fungi. In total, 8,938 conserved miRNA-target gene pairs and 3,901 Populus-specific miRNA-target gene pairs were detected. The result showed that Populus-specific miRNAs (66%) were more involved in the regulation of the disease resistance genes. By contrast, conserved miRNAs (>80%) target more whole-genome duplication (WGD)-derived transcription factors (TFs). Among the 1,023 WGD-derived TF pairs, 44.9% TF pairs had only one paralog being targeted by a miRNA that could be due to either gain or loss of a miRNA binding site after the WGD. A conserved hierarchical regulatory network combining promoter analyses and hierarchical clustering approach uncovered a miR164–NAM, ATAF, and CUC (NAC) transcription factor–mRNA regulatory module that has potential in Marssonina defense responses. Furthermore, analyses of the locations of miRNA precursor sequences reveal that pseudogenes and transposon contributed a certain proportion (∼30%) of the miRNA origin. Together, these observations provide evolutionary insights into the origin and potential roles of miRNAs in plant defense and functional innovation.

Gauging the trends of pseudogenes in plants

Pseudogenes, the debilitated parts of ancient genes, were previously scrapped off as junk or discarded genes with no functional significance. Pseudogenes have come under scrutiny for their functionality, since recent studies have unveiled their importance in the regulation of their corresponding parent genes and various biological mechanisms. Despite the enormous occurrence of pseudogenes in plants, the lack of experimental validation has contributed toward their unresolved roles in gene regulation. Contrarily, most of the studies associated with gene regulation have been mainly reported for humans, mice, and other mammalian genomes. Consequently, in order to present a cumulative report on plant-based pseudogenes research, an attempt has been made to assemble multiple studies presenting the pseudogene classification, the prediction and the determination of comparative accuracies of various computational pipelines, and recent trends in analyzing their biological functions, and regulatory mechanisms. This review represents the classical, as well as the recent advances on pseudogene identification and their potential roles in transcriptional regulation, which could possibly invigorate the quality of genome annotation, evolutionary analysis, and complexity surrounding the regulatory pathways in plants. Thus, when the ambiguous boundary girdling the pseudogenes eventually recedes on account of their explicit orchestration role, research in flora would no longer saunter compared to that on fauna.

Evolutionary research on the expansin protein family during the plant transition to land provides new insights into the development of Tartary buckwheat fruit

BACKGROUND

Plant transitions to land require robust cell walls for regulatory adaptations and to resist changing environments. Cell walls provide essential plasticity for plant cell division and defense, which are often conferred by the expansin superfamily with cell wall-loosening functions. However, the evolutionary mechanisms of expansin during plant terrestrialization are unclear.

RESULTS

Here, we identified 323 expansin proteins in 12 genomes from algae to angiosperms. Phylogenetic evolutionary, structural, motif gain and loss and Ka/Ks analyses indicated that highly conserved expansin proteins were already present in algae and expanded and purified after plant terrestrialization. We found that the expansion of the FtEXPA subfamily was caused by duplication events and that the functions of certain duplicated genes may have differentiated. More importantly, we generated space-time expression profiles and finally identified five differentially expressed FtEXPs in both large and small fruit Tartary buckwheat that may regulate fruit size by responding to indoleacetic acid.

CONCLUSIONS

A total of 323 expansin proteins from 12 representative plants were identified in our study during terrestrialization, and the expansin family that originated from algae expanded rapidly after the plants landed. The EXPA subfamily has more members and conservative evolution in angiosperms. FtEXPA1, FtEXPA11, FtEXPA12, FtEXPA19 and FtEXPA24 can respond to indole-3-acetic acid (IAA) signals and regulate fruit development. Our study provides a blueprint for improving the agronomic traits of Tartary buckwheat and a reference for defining the evolutionary history of the expansin family during plant transitions to land.

Structural characterization and duplication modes of pseudogenes in plants

We identified and characterized the pseudogene complements of five plant species: four dicots (Arabidopsis thaliana, Vitis vinifera, Populus trichocarpa and Phaseolus vulgaris) and one monocot (Oryza sativa). Retroposition was considered of modest importance for pseudogene formation in all investigated species except V. vinifera, which showed an unusually high number of retro-pseudogenes in non coding genic regions. By using a pipeline for the classification of sequence duplicates in plant genomes, we compared the relative importance of whole genome, tandem, proximal, transposed and dispersed duplication modes in the pseudo and functional gene complements. Pseudogenes showed higher tendencies than functional genes to genomic dispersion. Dispersed pseudogenes were prevalently fragmented and showed high sequence divergence at flanking regions. On the contrary, those deriving from whole genome duplication were proportionally less than expected based on observations on functional loci and showed higher levels of flanking sequence conservation than dispersed pseudogenes. Pseudogenes deriving from tandem and proximal duplications were in excess compared to functional loci, probably reflecting the high evolutionary rate associated with these duplication modes in plant genomes. These data are compatible with high rates of sequence turnover at neutral sites and double strand break repairs mediated duplication mechanisms.

Excision Dominates Pseudogenization During Fractionation After Whole Genome Duplication and in Gene Loss After Speciation in Plants

We take advantage of synteny blocks, the analytical construct enabled at the evolutionary moment of speciation or polyploidization, to follow the independent loss of duplicate genes in two sister species or the loss through fractionation of syntenic paralogs in a doubled genome. By examining how much sequence remains after a contiguous series of genes is deleted, we find that this residue remains at a constant low level independent of how many genes are lost—there are few if any relics of the missing sequence. Pseudogenes are rare or extremely transient in this context. The potential exceptions lie exclusively with a few examples of speciation, where the synteny blocks in some larger genomes tolerate degenerate sequence during genomic divergence of two species, but not after whole genome doubling in the same species where fractionation pressure eliminates virtually all non-coding sequence.

Computational Methods for Pseudogene Annotation Based on Sequence Homology

The number of complete genome sequences explodes more and more with each passing year. Thus, methods for genome annotation need to be honed constantly to handle the deluge of information. Annotation of pseudogenes (i.e., gene copies that appear not to make a functional protein) in genomes is a persistent problem; here, we overview pseudogene annotation methods that are based on the detection of sequence homology in genomic DNA.

Approaching the genetic dissection of indirect adventitious organogenesis process in tomato explants

The screening of 862 T-DNA lines was carried out to approach the genetic dissection of indirect adventitious organogenesis in tomato. Several mutants defective in different phases of adventitious organogenesis, namely callus growth (tdc-1), bud differentiation (tdb-1, -2, -3) and shoot-bud development (tds-1) were identified and characterized. The alteration of the TDC-1 gene blocked callus proliferation depending on the composition of growth regulators in the culture medium. Calli from tds-1 explants differentiated buds but did not develop normal shoots. Histological analysis showed that their abnormal development is due to failure in the organization of normal adventitious shoot meristems. Interestingly, tdc-1 and tds-1 mutant plants were indistinguishable from WT ones, indicating that the respective altered genes play specific roles in cell proliferation from explant cut zones (TDC-1 gene) or in the organization of adventitious shoot meristems (TDS-1 gene). Unlike the previous, plants of the three mutants defective in the differentiation of adventitious shoot-buds (tdb-1, -2, -3) showed multiple changes in vegetative and reproductive traits. Cosegregation analyses revealed the existence of an association between the phenotype of the tdb-3 mutant and a T-DNA insert, which led to the discovery that the SlMAPKKK17 gene is involved in the shoot-bud differentiation process.

Roles of Non-Coding RNAs in Response to Nitrogen Availability in Plants

Nitrogen (N) is an essential nutrient for plant growth and development; therefore, N deficiency is a major limiting factor in crop production. Plants have evolved mechanisms to cope with N deficiency, and the role of protein-coding genes in these mechanisms has been well studied. In the last decades, regulatory non-coding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs), and long ncRNAs (lncRNAs), have emerged as important regulators of gene expression in diverse biological processes. Recent advances in technologies for transcriptome analysis have enabled identification of N-responsive ncRNAs on a genome-wide scale. Characterization of these ncRNAs is expected to improve our understanding of the gene regulatory mechanisms of N response. In this review, we highlight recent progress in identification and characterization of N-responsive ncRNAs in Arabidopsis thaliana and several other plant species including maize, rice, and Populus.

Less Is More, Natural Loss-of-Function Mutation Is a Strategy for Adaptation

Gene gain and loss are crucial factors that shape the evolutionary success of diverse organisms. In the past two decades, more attention has been paid to the significance of gene gain through gene duplication or de novo genes. However, gene loss through natural loss-of-function (LoF) mutations, which is prevalent in the genomes of diverse organisms, has been largely ignored. With the development of sequencing techniques, many genomes have been sequenced across diverse species and can be used to study the evolutionary patterns of gene loss. In this review, we summarize recent advances in research on various aspects of LoF mutations, including their identification, evolutionary dynamics in natural populations, and functional effects. In particular, we discuss how LoF mutations can provide insights into the minimum gene set (or the essential gene set) of an organism. Furthermore, we emphasize their potential impact on adaptation. At the genome level, although most LoF mutations are neutral or deleterious, at least some of them are under positive selection and may contribute to biodiversity and adaptation. Overall, we highlight the importance of natural LoF mutations as a robust framework for understanding biological questions in general.

Expression Partitioning of Duplicate Genes at Single Cell Resolution in Arabidopsis Roots

Gene duplication is a key evolutionary phenomenon, prevalent in all organisms but particularly so in plants, where whole genome duplication (WGD; polyploidy) is a major force in genome evolution. Much effort has been expended in attempting to understand the evolution of duplicate genes, addressing such questions as why some paralog pairs rapidly return to single copy status whereas, in other pairs, both paralogs are retained and may diverge in expression pattern or function. The effect of a gene - its site of expression and thus the initial locus of its function - occurs at the level of a cell comprising a single cell type at a given state of the cell's development. Using Arabidopsis thaliana single cell transcriptomic data we categorized patterns of expression for 11,470 duplicate gene pairs across 36 cell clusters comprising nine cell types and their developmental states. Among these 11,470 pairs, 10,187 (88.8%) had at least one copy expressed in at least one of the 36 cell clusters. Pairs produced by WGD more often had both paralogs expressed in root cells than did pairs produced by small scale duplications. Three quarters of gene pairs expressed in the 36 cell clusters (7,608/10,187) showed extreme expression bias in at least one cluster, including 352 cases of reciprocal bias, a pattern consistent with expression subfunctionalization. More than twice as many pairs showed reciprocal expression bias between cell states than between cell types or between roots and leaves. A group of 33 gene pairs with reciprocal expression bias showed evidence of concerted divergence of gene networks in stele vs. epidermis. Pairs with both paralogs expressed without bias were less likely to have paralogs with divergent mutant phenotypes; such bias-free pairs showed evidence of preservation by maintenance of dosage balance. Overall, we found considerable evidence of shifts in gene expression following duplication, including in >80% of pairs encoding 7,653 genes expressed ubiquitously in all root cell types and states for which we inferred the polarity of change.

Genome-Wide Identification and Characterization of the bHLH Transcription Factor Family in Pepper (Capsicum annuum L.)

Plant basic helix–loop–helix (bHLH) transcription factors are involved in the regulation of various biological processes in plant growth, development, and stress response. However, members of this important transcription factor family have not been systematically identified and analyzed in pepper (Capsicum annuum L.). In this study, we identified 122 CabHLH genes in the pepper genome and renamed them based on their chromosomal locations. CabHLHs were divided into 21 subfamilies according to their phylogenetic relationships, and genes from the same subfamily had similar motif compositions and gene structures. Sixteen pairs of tandem and segmental duplicated genes were detected in the CabHLH family. Cis-elements identification and expression analysis of the CabHLHs revealed that they may be involved in plant development and stress responses. This study is the first comprehensive analysis of the CabHLH genes and will serve as a reference for further characterization of their molecular functions.

Genome-wide identification and comparative analysis of the ADH gene family in Chinese white pear (Pyrus bretschneideri) and other Rosaceae species

Alcohol dehydrogenase (ADH) is essential to the formation of aromatic compounds in fruits. However, the evolutionary history and characteristics of ADH gene expression remain largely unclear in Rosaceae fruit species. In this study, 464 ADH genes were identified in eight Rosaceae fruit species, 68 of the genes were from pear and which were classified into four subgroups. Frequent single gene duplication events were found to have contributed to the formation of ADH gene clusters and the expansion of the ADH gene family in these eight Rosaceae species. Purifying selection was the major force in ADH gene evolution. The younger genes derived from tandem and proximal duplications had evolved faster than those derived from other types of duplication. RNA-Seq and qRT-PCR analysis revealed that the expression levels of three ADH genes were closely correlated with the content of aromatic compounds detected during fruit development.

Genome-wide analysis of MIKC-type MADS-box genes in wheat: pervasive duplications, functional conservation and putative neofunctionalization

Wheat (Triticum aestivum) is one of the most important crops worldwide. Given a growing global population coupled with increasingly challenging cultivation conditions, facilitating wheat breeding by fine-tuning important traits is of great importance. MADS-box genes are prime candidates for this, as they are involved in virtually all aspects of plant development. Here, we present a detailed overview of phylogeny and expression of 201 wheat MIKC-type MADS-box genes. Homoeolog retention is significantly above the average genome-wide retention rate for wheat genes, indicating that many MIKC-type homoeologs are functionally important and not redundant. Gene expression is generally in agreement with the expected subfamily-specific expression pattern, indicating broad conservation of function of MIKC-type genes during wheat evolution. We also found extensive expansion of some MIKC-type subfamilies, especially those potentially involved in adaptation to different environmental conditions like flowering time genes. Duplications are especially prominent in distal telomeric regions. A number of MIKC-type genes show novel expression patterns and respond, for example, to biotic stress, pointing towards neofunctionalization. We speculate that conserved, duplicated and neofunctionalized MIKC-type genes may have played an important role in the adaptation of wheat to a diversity of conditions, hence contributing to the importance of wheat as a global staple food.

Origination and Function of Plant Pseudogenes

Pseudogenes, nonfunctional genomic sequences derived from functional protein-coding genes, form by duplication or retrotransposition, and loss of gene function by disabling mutations. Studies on the evolution and functional aspects of plant pseudogenes are limited, despite their abundance in the plant genome. To date, most researches on pseudogenes focus on mammals. Here, we summarized current knowledge on pseudogenes including the historical and recent progress, analyzes their essential roles in gene regulation in hope of further stimulating researches in plant species for understanding gene regulation and evolution.