The hidden duplication past of Arabidopsis thaliana

Structural and Phylogenetic In Silico Characterization of Vitis vinifera PRR Protein as Potential Target for Plasmopara viticola Infection

Fungi infection, especially derived from Plasmopara viticola, causes severe grapevine economic losses worldwide. Despite the availability of chemical treatments, looking for eco-friendly ways to control Vitis vinifera infection is gaining much more attention. When a plant is infected, multiple disease-control molecular mechanisms are activated. PRRs (Pattern Recognition Receptors) and particularly RLKs (receptor-like kinases) take part in the first barrier of the immune system, and, as a consequence, the kinase signaling cascade is activated, resulting in an immune response. In this context, discovering new lectin-RLK (LecRLK) membrane-bounded proteins has emerged as a promising strategy. The genome-wide localization of potential LecRLKs involved in disease defense was reported in two grapevine varieties of great economic impact: Chardonnay and Pinot Noir. A total of 23 potential amino acid sequences were identified, exhibiting high-sequence homology and evolution related to tandem events. Based on the domain architecture, a carbohydrate specificity ligand assay was conducted with docking, revealing two sequences as candidates for specific Vitis vinifera-Plasmopara viticola host-pathogen interaction. This study confers a starting point for designing new effective antifungal treatments directed at LecRLK targets in Vitis vinifera.

Genomic insights into the chromosomal elongation in a family of Collembola

Collembola is a highly diverse and abundant group of soil arthropods with chromosome numbers ranging from 5 to 11. Previous karyotype studies indicated that the Tomoceridae family possesses an exceptionally long chromosome. To better understand chromosome size evolution in Collembola, we obtained a chromosome-level genome of Yoshiicerus persimilis with a size of 334.44 Mb and BUSCO completeness of 97.0% (n = 1013). Both genomes of Y. persimilis and Tomocerus qinae (recently published) have an exceptionally large chromosome (ElChr greater than 100 Mb), accounting for nearly one-third of the genome. Comparative genomic analyses suggest that chromosomal elongation occurred independently in the two species approximately 10 million years ago, rather than in the ancestor of the Tomoceridae family. The ElChr elongation was caused by large tandem and segmental duplications, as well as transposon proliferation, with genes in these regions experiencing weaker purifying selection (higher dN/dS) than conserved regions. Moreover, inter-genomic synteny analyses indicated that chromosomal fission/fusion events played a crucial role in the evolution of chromosome numbers (ranging from 5 to 7) within Entomobryomorpha. This study provides a valuable resource for investigating the chromosome evolution of Collembola.

Structural and functional characterization of genes PYL-PP2C-SnRK2s in the ABA signalling pathway of Cucurbita pepo

BACKGROUND

The core regulation of the abscisic acid (ABA) signalling pathway comprises the multigenic families PYL, PP2C, and SnRK2. In this work, we conducted a genome-wide study of the components of these families in Cucurbita pepo.

RESULTS

The bioinformatic analysis of the C. pepo genome resulted in the identification of 19 CpPYL, 102 CpPP2C and 10 CpSnRK2 genes. The investigation of gene structure and protein motifs allowed to define 4 PYL, 13 PP2C and 3 SnRK2 subfamilies. RNA-seq analysis was used to determine the expression of these gene families in different plant organs, as well as to detect their differential gene expression during germination, and in response to ABA and cold stress in leaves. The specific tissue expression of some gene members indicated the relevant role of some ABA signalling genes in plant development. Moreover, their differential expression under ABA treatment or cold stress revealed those ABA signalling genes that responded to ABA, and those that were up- or down-regulated in response to cold stress. A reduced number of genes responded to both treatments. Specific PYL-PP2C-SnRK2 genes that had potential roles in germination were also detected, including those regulated early during the imbibition phase, those regulated later during the embryo extension and radicle emergence phase, and those induced or repressed during the whole germination process.

CONCLUSIONS

The outcomes of this research open new research lines for agriculture and for assessing gene function in future studies.

Comprehensive Evolutionary Analysis of the SMXL Gene Family in Rosaceae: Further Insights into Its Origin, Expansion, Diversification, and Role in Regulating Pear Branching

SMXL genes constitute a conserved gene family that is ubiquitous in angiosperms and involved in regulating various plant processes, including branching, leaf elongation, and anthocyanin biosynthesis, but little is known about their molecular functions in pear branching. Here, we performed genome-wide identification and investigation of the SMXL genes in 16 angiosperms and analyzed their phylogenetics, structural features, conserved motifs, and expression patterns. In total, 121 SMXLs genes were identified and were classified into four groups. The number of non-redundant SMXL genes in each species varied from 3 (Amborella trichopoda Baill.) to 18 (Glycine max Merr.) and revealed clear gene expansion events over evolutionary history. All the SMXL genes showed conserved structures, containing no more than two introns. Three-dimensional protein structure prediction revealed distinct structures between but similar structures within groups. A quantitative real-time PCR analysis revealed different expressions of 10 SMXL genes from pear branching induced by fruit-thinning treatment. Overall, our study provides a comprehensive investigation of SMXL genes in the Rosaceae family, especially pear. The results offer a reference for understanding the evolutionary history of SMXL genes and provide excellent candidates for studying fruit tree branching regulation, and in facilitating pear pruning and planting strategies.

The genetic and molecular basis of haploinsufficiency in flowering plants

In diploid organisms, haploinsufficiency can be defined as the requirement for more than one fully functional copy of a gene. In contrast to most genes, whose loss-of-function alleles are recessive, loss-of-function alleles of haploinsufficient genes are dominant. However, forward and reverse genetic screens are biased toward obtaining recessive, loss-of-function mutations, and therefore, dominant mutations of all types are underrepresented in mutant collections. Despite this underrepresentation, haploinsufficient loci have intriguing implications for studies of genome evolution, gene dosage, stability of protein complexes, genetic redundancy, and gene expression. Here we review examples of haploinsufficiency in flowering plants and describe the underlying molecular mechanisms and evolutionary forces driving haploinsufficiency. Finally, we discuss the masking of haploinsufficiency by genetic redundancy, a widespread phenomenon among angiosperms.

Dynamics of cis-regulatory sequences and transcriptional divergence of duplicated genes in soybean

Transcriptional divergence of duplicated genes after whole genome duplication (WGD) has been described in many plant lineages and is often associated with subgenome dominance, a genome-wide mechanism. However, it is unknown what underlies the transcriptional divergence of duplicated genes in polyploid species that lack subgenome dominance. Soybean is a paleotetraploid with a WGD that occurred 5 to 13 Mya. Approximately 50% of the duplicated genes retained from this WGD exhibit transcriptional divergence. We developed accessible chromatin region (ACR) datasets from leaf, flower, and seed tissues using MNase-hypersensitivity sequencing. We validated enhancer function of several ACRs associated with known genes using CRISPR/Cas9-mediated genome editing. The ACR datasets were used to examine and correlate the transcriptional patterns of 17,111 pairs of duplicated genes in different tissues. We demonstrate that ACR dynamics are correlated with divergence of both expression level and tissue specificity of individual gene pairs. Gain or loss of flanking ACRs and mutation of cis-regulatory elements (CREs) within the ACRs can change the balance of the expression level and/or tissue specificity of the duplicated genes. Analysis of DNA sequences associated with ACRs revealed that the extensive sequence rearrangement after the WGD reshaped the CRE landscape, which appears to play a key role in the transcriptional divergence of duplicated genes in soybean. This may represent a general mechanism for transcriptional divergence of duplicated genes in polyploids that lack subgenome dominance.

Genomic imbalance modulates transposable element expression in maize

Genomic imbalance refers to the more severe phenotypic consequences of changing part of a chromosome compared with the whole genome set. Previous genome imbalance studies in maize have identified prevalent inverse modulation of genes on the unvaried chromosomes (trans) with both the addition or subtraction of chromosome arms. Transposable elements (TEs) comprise a substantial fraction of the genome, and their reaction to genomic imbalance is therefore of interest. Here, we analyzed TE expression using RNA-seq data of aneuploidy and ploidy series and found that most aneuploidies showed an inverse modulation of TEs, but reductions in monosomy and increases in disomy and trisomy were also common. By contrast, the ploidy series showed little TE modulation. The modulation of TEs and genes in the same experimental group were compared, and TEs showed greater modulation than genes, especially in disomy. Class I and II TEs were differentially modulated in most aneuploidies, and some superfamilies in each TE class also showed differential modulation. Finally, the significantly upregulated TEs in three disomies (TB-7Lb, TB9Lc, and TB-10L19) did not increase the proportion of adjacent gene expression when compared with non-differentially expressed TEs, indicating that modulations of TEs do not compound the effect on genes. These results suggest that the prevalent inverse TE modulation in aneuploidy results from stoichiometric upset of the regulatory machinery used by TEs, similar to the response of core genes to genomic imbalance.

Contrasting modes of macro and microsynteny evolution in a eukaryotic subphylum

Examination of the changes in order and arrangement of homologous genes is key for understanding the mechanisms of genome evolution in eukaryotes. Previous comparisons between eukaryotic genomes have revealed considerable conservation across species that diverged hundreds of millions of years ago (e.g., vertebrates,1,2,3 bilaterian animals,4,5 and filamentous fungi6). However, understanding how genome organization evolves within and between eukaryotic major lineages remains underexplored. We analyzed high-quality genomes of 120 representative budding yeast species (subphylum Saccharomycotina) spanning ∼400 million years of eukaryotic evolution to examine how their genome organization evolved and to compare it with the evolution of animal and plant genome organization.7 We found that the decay of both macrosynteny (the conservation of homologous chromosomes) and microsynteny (the conservation of local gene content and order) was strongly associated with evolutionary divergence across budding yeast major clades. However, although macrosynteny decayed very fast, within ∼100 million years, the microsynteny of many genes-especially genes in metabolic clusters (e.g., in the GAL gene cluster8)-was much more deeply conserved both within major clades and across the subphylum. We further found that when genomes with similar evolutionary divergence times were compared, budding yeasts had lower macrosynteny conservation than animals and filamentous fungi but higher conservation than angiosperms. In contrast, budding yeasts had levels of microsynteny conservation on par with mammals, whereas angiosperms exhibited very low conservation. Our results provide new insight into the tempo and mode of the evolution of gene and genome organization across an entire eukaryotic subphylum.

Comparative genomics of orobanchaceous species with different parasitic lifestyles reveals the origin and stepwise evolution of plant parasitism

Orobanchaceae is the largest family of parasitic plants, containing autotrophic and parasitic plants with all degrees of parasitism. This makes it by far the best family for studying the origin and evolution of plant parasitism. Here we provide three high-quality genomes of orobanchaceous plants, the autotrophic Lindenbergia luchunensis and the holoparasitic plants Phelipanche aegyptiaca and Orobanche cumana. Phylogenomic analysis of these three genomes together with those previously published and the transcriptomes of other orobanchaceous species created a robust phylogenetic framework for Orobanchaceae. We found that an ancient whole-genome duplication (WGD; about 73.48 million years ago), which occurred earlier than the origin of Orobanchaceae, might have contributed to the emergence of parasitism. However, no WGD events occurred in any lineage of orobanchaceous parasites except for Striga after divergence from their autotrophic common ancestor, suggesting that, in contrast with previous speculations, WGD is not associated with the emergence of holoparasitism. We detected evident convergent gene loss in all parasites within Orobanchaceae and between Orobanchaceae and dodder Cuscuta australis. The gene families in the orobanchaceous parasites showed a clear pattern of recent gains and expansions. The expanded gene families are enriched in functions related to the development of the haustorium, suggesting that recent gene family expansions may have facilitated the adaptation of orobanchaceous parasites to different hosts. This study illustrates a stepwise pattern in the evolution of parasitism in the orobanchaceous parasites and will facilitate future studies on parasitism and the control of parasitic plants in agriculture.

The multiple fates of gene duplications: Deletion, hypofunctionalization, subfunctionalization, neofunctionalization, dosage balance constraints, and neutral variation

Gene duplications have long been recognized as a contributor to the evolution of genes with new functions. Multiple copies of genes can result from tandem duplication, from transposition to new chromosomes, or from whole-genome duplication (polyploidy). The most common fate is that one member of the pair is deleted to return the gene to the singleton state. Other paths involve the reduced expression of both copies (hypofunctionalization) that are held in duplicate to maintain sufficient quantity of function. The two copies can split functions (subfunctionalization) or can diverge to generate a new function (neofunctionalization). Retention of duplicates resulting from doubling of the whole genome occurs for genes involved with multicomponent interactions such as transcription factors and signal transduction components. In contrast, these classes of genes are underrepresented in small segmental duplications. This complementary pattern suggests that the balance of interactors affects the fate of the duplicate pair. We discuss the different mechanisms that maintain duplicated genes, which may change over time and intersect.

Exploring the SiCCT Gene Family and Its Role in Heading Date in Foxtail Millet

CCT transcription factors are involved in the regulation of photoperiod and abiotic stress in Arabidopsis and rice. It is not clear that how CCT gene family expand and regulate heading date in foxtail millet. In this study, we conducted a systematic analysis of the CCT gene family in foxtail millet. Thirty-nine CCT genes were identified and divided into four subfamilies based on functional motifs. Analysis showed that dispersed duplication played a predominant role in the expansion of CCT genes during evolution. Nucleotide diversity analysis suggested that genes in CONSTANS (COL)-like, CCT MOTIF FAMILY (CMF)-like, and pseudoresponse response regulator (PRR)-like subfamilies were subjected to selection. Fifteen CCT genes were colocalized with previous heading date quantitative trait loci (QTL) and genome-wide association analysis (GWAS) signals. Transgenic plants were then employed to confirm that overexpression of the CCT gene SiPRR37 delayed the heading date and increased plant height. Our study first investigated the characterization and expansion of the CCT family in foxtail millet and demonstrated the role of SiPRR37. These results lay a significant foundation for further research on the function of CCT genes and provide a cue for the regulation of heading date.

A Population-Genetic Lens into the Process of Gene Loss Following Whole-Genome Duplication

Whole-genome duplications (WGDs) have occurred in many eukaryotic lineages. However, the underlying evolutionary forces and molecular mechanisms responsible for the long-term retention of gene duplicates created by WGDs are not well understood. We employ a population-genomic approach to understand the selective forces acting on paralogs and investigate ongoing duplicate-gene loss in multiple species of Paramecium that share an ancient WGD. We show that mutations that abolish protein function are more likely to be segregating in retained WGD paralogs than in single-copy genes, most likely because of ongoing nonfunctionalization post-WGD. This relaxation of purifying selection occurs in only one WGD paralog, accompanied by the gradual fixation of nonsynonymous mutations and reduction in levels of expression, and occurs over a long period of evolutionary time, "marking" one locus for future loss. Concordantly, the fitness effects of new nonsynonymous mutations and frameshift-causing indels are significantly more deleterious in the highly expressed copy compared with their paralogs with lower expression. Our results provide a novel mechanistic model of gene duplicate loss following WGDs, wherein selection acts on the sum of functional activity of both duplicate genes, allowing the two to wander in expression and functional space, until one duplicate locus eventually degenerates enough in functional efficiency or expression that its contribution to total activity is too insignificant to be retained by purifying selection. Retention of duplicates by such mechanisms predicts long times to duplicate-gene loss, which should not be falsely attributed to retention due to gain/change in function.

Dosage-sensitive miRNAs trigger modulation of gene expression during genomic imbalance in maize

The genomic imbalance caused by varying the dosage of individual chromosomes or chromosomal segments (aneuploidy) has more detrimental effects than altering the dosage of complete chromosome sets (ploidy). Previous analysis of maize (Zea mays) aneuploids revealed global modulation of gene expression both on the varied chromosome (cis) and the remainder of the genome (trans). However, little is known regarding the role of microRNAs (miRNAs) under genomic imbalance. Here, we report the impact of aneuploidy and polyploidy on the expression of miRNAs. In general, cis miRNAs in aneuploids present a predominant gene-dosage effect, whereas trans miRNAs trend toward the inverse level, although other types of responses including dosage compensation, increased effect, and decreased effect also occur. By contrast, polyploids show less differential miRNA expression than aneuploids. Significant correlations between expression levels of miRNAs and their targets are identified in aneuploids, indicating the regulatory role of miRNAs on gene expression triggered by genomic imbalance.

Genome-wide identification and expression analysis of CPP-like gene family in Triticum aestivum L. under different hormone and stress conditions

The CPP-like plant‐specific transcription factor has a prominent role in plant development and growth through cell division and differential activities. However, little information is available about the CPP gene family in Triticum aestivum L. Herein, we identified 37 and 11 CPP genes in the wheat and rice genome databases, respectively. The phylogeny of the CPP protein-like family members was further divided into five subfamilies based on structural similarities and phenotypic functional diversities. The in silico expression analysis showed that CPP genes are highly expressed in some tissues, such as shoot apex, shoot, leaf, leaf sheath, and microspore. Furthermore, the qRT-PCR found higher expression for TaCPP gene family members in leaf, leaf blade, young spike, mature spike, and differential expression patterns under abiotic stresses, including heat, drought, salt, and hormonal treatment, such as indole acetic acid and 1-aminocyclopropane-1 carboxylic acid. We found that CPP gene family members are mostly located in the nucleus after infiltrating the CPP5-1B-GFP and TaCPP11-3B-GFP into tobacco leaves. The overexpression of the TaCPP5-1D gene revealed that the CPP gene positively regulates the germanium, shoot, and root activities in Arabidopsis. The TaCPP5-1D-overexpressed plants showed less anti-oxidative sensitivity under drought stress conditions. These results demonstrated that TaCPP5-1D protein has a crucial contribution by interacting with TaCPP11-3B protein in maintaining stress homeostasis under the natural and unfavorable environmental conditions for growth, development, and stress resistance activities. Therefore, this study could be used as pioneer knowledge to further investigate the function of CPP genes in plant growth and development.

Illegitimate Recombination between Duplicated Genes Generated from Recursive Polyploidizations Accelerated the Divergence of the Genus Arachis

The peanut (Arachis hypogaea L.) is the leading oil and food crop among the legume family. Extensive duplicate gene pairs generated from recursive polyploidizations with high sequence similarity could result from gene conversion, caused by illegitimate DNA recombination. Here, through synteny-based comparisons of two diploid and three tetraploid peanut genomes, we identified the duplicated genes generated from legume common tetraploidy (LCT) and peanut recent allo-tetraploidy (PRT) within genomes. In each peanut genome (or subgenomes), we inferred that 6.8–13.1% of LCT-related and 11.3–16.5% of PRT-related duplicates were affected by gene conversion, in which the LCT-related duplicates were the most affected by partial gene conversion, whereas the PRT-related duplicates were the most affected by whole gene conversion. Notably, we observed the conversion between duplicates as the long-lasting contribution of polyploidizations accelerated the divergence of different Arachis genomes. Moreover, we found that the converted duplicates are unevenly distributed across the chromosomes and are more often near the ends of the chromosomes in each genome. We also confirmed that well-preserved homoeologous chromosome regions may facilitate duplicates’ conversion. In addition, we found that these biological functions contain a higher number of preferentially converted genes, such as catalytic activity-related genes. We identified specific domains that are involved in converted genes, implying that conversions are associated with important traits of peanut growth and development.

One Hundred Years of Gene Balance: How Stoichiometric Issues Affect Gene Expression, Genome Evolution, and Quantitative Traits

A century ago experiments with the flowering plant Datura stramonium and the fruit fly Drosophila melanogaster revealed that adding an extra chromosome to a karyotype was much more detrimental than adding a whole set of chromosomes. This phenomenon was referred to as gene balance and has been recapitulated across eukaryotic species. Here, we retrace some developments in this field. Molecular studies suggest that the basis of balance involves stoichiometric relationships of multi-component interactions. This concept has implication for the mechanisms controlling gene expression, genome evolution, sex chromosome evolution/dosage compensation, speciation mechanisms, and the underlying genetics of quantitative traits.

An anchored chromosome-scale genome assembly of spinach improves annotation and reveals extensive gene rearrangements in euasterids

Spinach (Spinacia oleracea L.) is a member of the Caryophyllales family, a basal eudicot asterid that consists of sugar beet (Beta vulgaris L. subsp. vulgaris), quinoa (Chenopodium quinoa Willd.), and amaranth (Amaranthus hypochondriacus L.). With the introduction of baby leaf types, spinach has become a staple food in many homes. Production issues focus on yield, nitrogen-use efficiency and resistance to downy mildew (Peronospora effusa). Although genomes are available for the above species, a chromosome-level assembly exists only for quinoa, allowing for proper annotation and structural analyses to enhance crop improvement. We independently assembled and annotated genomes of the cultivar Viroflay using short-read strategy (Illumina) and long-read strategies (Pacific Biosciences) to develop a chromosome-level, genetically anchored assembly for spinach. Scaffold N50 for the Illumina assembly was 389 kb, whereas that for Pacific BioSciences was 4.43 Mb, representing 911 Mb (93% of the genome) in 221 scaffolds, 80% of which are anchored and oriented on a sequence-based genetic map, also described within this work. The two assemblies were 99.5% collinear. Independent annotation of the two assemblies with the same comprehensive transcriptome dataset show that the quality of the assembly directly affects the annotation with significantly more genes predicted (26,862 vs. 34,877) in the long-read assembly. Analysis of resistance genes confirms a bias in resistant gene motifs more typical of monocots. Evolutionary analysis indicates that Spinacia is a paleohexaploid with a whole-genome triplication followed by extensive gene rearrangements identified in this work. Diversity analysis of 75 lines indicate that variation in genes is ample for hypothesis-driven, genomic-assisted breeding enabled by this work.

Genomic imbalance determines positive and negative modulation of gene expression in diploid maize

Genomic imbalance caused by changing the dosage of individual chromosomes (aneuploidy) has a more detrimental effect than varying the dosage of complete sets of chromosomes (ploidy). We examined the impact of both increased and decreased dosage of 15 distal and 1 interstitial chromosomal regions via RNA-seq of maize (Zea mays) mature leaf tissue to reveal new aspects of genomic imbalance. The results indicate that significant changes in gene expression in aneuploids occur both on the varied chromosome (cis) and the remainder of the genome (trans), with a wider spread of modulation compared with the whole-ploidy series of haploid to tetraploid. In general, cis genes in aneuploids range from a gene-dosage effect to dosage compensation, whereas for trans genes the most common effect is an inverse correlation in that expression is modulated toward the opposite direction of the varied chromosomal dosage, although positive modulations also occur. Furthermore, this analysis revealed the existence of increased and decreased effects in which the expression of many genes under genome imbalance are modulated toward the same direction regardless of increased or decreased chromosomal dosage, which is predicted from kinetic considerations of multicomponent molecular interactions. The findings provide novel insights into understanding mechanistic aspects of gene regulation.

Predominantly inverse modulation of gene expression in genomically unbalanced disomic haploid maize

The phenotypic consequences of the addition or subtraction of part of a chromosome is more severe than changing the dosage of the whole genome. By crossing diploid trisomies to a haploid inducer, we identified 17 distal segmental haploid disomies that cover ∼80% of the maize genome. Disomic haploids provide a level of genomic imbalance that is not ordinarily achievable in multicellular eukaryotes, allowing the impact to be stronger and more easily studied. Transcriptome size estimates revealed that a few disomies inversely modulate most of the transcriptome. Based on RNA sequencing, the expression levels of genes located on the varied chromosome arms (cis) in disomies ranged from being proportional to chromosomal dosage (dosage effect) to showing dosage compensation with no expression change with dosage. For genes not located on the varied chromosome arm (trans), an obvious trans-acting effect can be observed, with the majority showing a decreased modulation (inverse effect). The extent of dosage compensation of varied cis genes correlates with the extent of trans inverse effects across the 17 genomic regions studied. The results also have implications for the role of stoichiometry in gene expression, the control of quantitative traits, and the evolution of dosage-sensitive genes.

The Landscape of the Genomic Distribution and the Expression of the F-Box Genes Unveil Genome Plasticity in Hexaploid Wheat during Grain Development and in Response to Heat and Drought Stress

FBX proteins are subunits of the SCF complex (Skp1-cullin-FBX) belonging to the E3 ligase family, which is involved in the ubiquitin-proteasome 26S (UPS) pathway responsible for the post-translational protein turnover. By targeting, in a selective manner, key regulatory proteins for ubiquitination and 26S proteasome degradation, FBX proteins play a major role in plant responses to diverse developmental and stress conditions. Although studies on the genomic organization of the FBX gene family in various species have been reported, knowledge related to bread wheat (Triticum aestivum) is scarce and needs to be broadened. Using the latest assembly of the wheat genome, we identified 3670 TaFBX genes distributed non-homogeneously within the three subgenomes (A, B and D) and between the 21 chromosomes, establishing it as one of the richest gene families among plant species. Based on the presence of the five different chromosomal regions previously identified, the present study focused on the genomic distribution of the TaFBX family and the identification of differentially expressed genes during the embryogenesis stages and in response to heat and drought stress. Most of the time, when comparing the expected number of genes (taking into account the formal gene distribution on the entire wheat genome), the TaFBX family harbors a different pattern at the various stratum of observation (subgenome, chromosome, chromosomal regions). We report here that the local gene expansion of the TaFBX family must be the consequence of multiple and complex events, including tandem and small-scale duplications. Regarding the differentially expressed TaFBX genes, while the majority of the genes are localized in the distal chromosomal regions (R1 and R3), differentially expressed genes are more present in the interstitial regions (R2a and R2b) than expected, which could be an indication of the preservation of major genes in those specific chromosomal regions.

The Gene Balance Hypothesis: Epigenetics and Dosage Effects in Plants

Dosage effects in plants are caused by changes in the copy number of chromosomes, segments of chromosomes, or multiples of individual genes. Genes often exhibit a dosage effect in which the amount of product is closely correlated with the number of copies present. However, when larger segments of chromosomes are varied, there are trans-acting effects across the genome that are unleashed that modulate gene expression in cascading effects. These appear to be mediated by the stoichiometric relationship of gene regulatory machineries. There are both positive and negative modulations of target gene expression, but the latter is the plurality effect. When this inverse effect is combined with a dosage effect, compensation for a gene can occur in which its expression is similar to the normal diploid regardless of the change in chromosomal dosage. In contrast, changing the whole genome in a polyploidy series has fewer relative effects as the stoichiometric relationship is not disrupted. Together, these observations suggest that the stoichiometry of gene regulation is important as a reflection of the mode of assembly of the individual subunits involved in the effective regulatory macromolecular complexes. This principle has implications for gene expression mechanisms, quantitative trait genetics, and the evolution of genes depending on the mode of duplication, either segmentally or via whole-genome duplication.

A chromosome-scale reference genome of Lobularia maritima, an ornamental plant with high stress tolerance

Lobularia maritima (L.) Desv. is an ornamental plant cultivated across the world. It belongs to the family Brassicaceae and can tolerate dry, poor and contaminated habitats. Here, we present a chromosome-scale, high-quality genome assembly of L. maritima based on integrated approaches combining Illumina short reads and Hi-C chromosome conformation data. The genome was assembled into 12 pseudochromosomes with a 197.70 Mb length, and it includes 25,813 protein-coding genes. Approximately 41.94% of the genome consists of repetitive sequences, with abundant long terminal repeat transposable elements. Comparative genomic analysis confirmed that L. maritima underwent a species-specific whole-genome duplication (WGD) event ~22.99 million years ago. We identified ~1900 species-specific genes, 25 expanded gene families, and 50 positively selected genes in L. maritima. Functional annotations of these genes indicated that they are mainly related to stress tolerance. These results provide new insights into the stress tolerance of L. maritima, and this genomic resource will be valuable for further genetic improvement of this important ornamental plant.

DNA Content Variation and SNP Diversity Within a Single Population of Asexual Snails

A growing body of research suggests that many clonal populations maintain genetic diversity even without occasional sexual reproduction. The purpose of our study was to document variation in single-nucleotide polymorphism (SNP) diversity, DNA content, and pathogen susceptibility in clonal lineages of the New Zealand freshwater snail, Potamopyrgus antipodarum. We studied snails that were collected from multiple field sites around a single lake (Lake Alexandrina), as well as isofemale clonal lineages that had been isolated and maintained in the laboratory. We used the kompetitive allele specific PCR (KASP) method to genotype our samples at 46 nuclear SNP sites, and we used flow cytometry to estimate DNA content. We found high levels of SNP diversity, both in our field samples and in our clonal laboratory lines. We also found evidence of high variation in DNA content among clones, even among clones with identical genotypes across all SNP sites. Controlled pathogen exposures of the laboratory populations revealed variation in susceptibility among distinct clonal genotypes, which was independent of DNA content. Taken together, these results show high levels of diversity among asexual snails, especially for DNA content, and they suggest rapid genome evolution in asexuals.

Evolutionary strategies drive a balance of the interacting gene products for the CBL and CIPK gene families

Genes encoding interacting proteins tend to be co-retained after whole-genome duplication (WGD). The preferential retention after WGD has been explained by the gene balance hypothesis (GBH). However, small-scale duplications could independently occur in the connected gene families. Certain evolutionary strategies might keep the dosage balanced. Here, we examined the gene duplication, interaction and expression patterns of calcineurin B-like (CBL) and CBL-interacting protein kinase (CIPK) gene families to understand the underlying principles. The ratio of the CBL and CIPK gene numbers evolved from 5 : 7 in Physcomitrella to 10 : 26 in Arabidopsis, and retrotransposition, tandem duplication, and WGDs contributed to the expansion. Two pairs of CBLs and six pairs of CIPKs were retained after the α WGD in Arabidopsis, in which specific interaction patterns were identified. In some cases, two retained CBLs (CIPKs) might compete to interact with a sole CIPK (CBL). Results of gene expression analyses indicated that the relatively over-retained duplicates tend to show asymmetric expression, thus avoiding competition. In conclusion, our results suggested that the highly specific interaction, together with the differential gene expression pattern, jointly maintained the balanced dosage for the interacting CBL and CIPK proteins.

Multi-faceted analysis provides little evidence for recurrent whole-genome duplications during hexapod evolution

BACKGROUND

Gene duplication events play an important role in the evolution and adaptation of organisms. Duplicated genes can arise through different mechanisms, including whole-genome duplications (WGDs). Recently, WGD was suggested to be an important driver of evolution, also in hexapod animals.

RESULTS

Here, we analyzed 20 high-quality hexapod genomes using whole-paranome distributions of estimated synonymous distances (KS), patterns of within-genome co-linearity, and phylogenomic gene tree-species tree reconciliation methods. We observe an abundance of gene duplicates in the majority of these hexapod genomes, yet we find little evidence for WGD. The majority of gene duplicates seem to have originated through small-scale gene duplication processes. We did detect segmental duplications in six genomes, but these lacked the within-genome co-linearity signature typically associated with WGD, and the age of these duplications did not coincide with particular peaks in KS distributions. Furthermore, statistical gene tree-species tree reconciliation failed to support all but one of the previously hypothesized WGDs.

CONCLUSIONS

Our analyses therefore provide very limited evidence for WGD having played a significant role in the evolution of hexapods and suggest that alternative mechanisms drive gene duplication events in this group of animals. For instance, we propose that, along with small-scale gene duplication events, episodes of increased transposable element activity could have been an important source for gene duplicates in hexapods.

Peculiar genomic traits in the stress-adapted cryptoendolithic Antarctic fungus Friedmanniomyces endolithicus

Friedmanniomyces endolithicus is a highly melanized fungus endemic to the Antarctic, occurring exclusively in endolithic communities of the ice-free areas of the Victoria Land, including the McMurdo Dry Valleys, the coldest and most hyper-arid desert on Earth and accounted as the Martian analog on our planet. F. endolithicus is highly successful in these inhospitable environments, the most widespread and commonly isolated species from these peculiar niches, indicating a high degree of adaptation. The nature of its extremo tolerance has not been previously studied. To investigate this, we sequenced genome of F. endolithicus CCFEE 5311 to explore gene content and genomic patterns that could be attributed to its specialization. The predicted functional potential of the genes was assigned by similarity to InterPro and CAZy domains. The genome was compared to phylogenetically close relatives which are also melanized fungi occurring in extreme environments including Friedmanniomyces simplex, Baudoinia panamericana, Acidomyces acidophilus, Hortaea thailandica and Hortaea werneckii. We tested if shared genomic traits existed among these species and hyper-extremotolerant fungus F. endolithicus. We found that some characters for stress tolerance such as meristematic growth and cold tolerance are enriched in F. endolithicus that may be triggered by the exposure to Antarctic prohibitive conditions.

Expansion and Evolutionary Patterns of Glycosyltransferase Family 8 in Gramineae Crop Genomes and Their Expression under Salt and Cold Stresses in Oryza sativa ssp. japonica

Plant cell walls play a fundamental role in several ways, providing structural support for cells, resistance against pathogens and facilitating the communication between cells. The glycosyltransferase family 8 (GT8) is involved in the formation of the plant cell wall. However, the evolutionary relationship and the functional differentiation of this important gene family remain obscure in Gramineae crop genomes. In the present investigation, we identified 269 GT8 genes in the seven Gramineae representative crop genomes, namely, 33 in Hordeum vulgare, 37 in Brachypodium distachyon, 40 in Oryza sativa ssp. japonica, 41 in Oryza rufipogon, 36 in Setaria italica, 37 in Sorghum bicolor, and 45 in Zea mays. Phylogenetic analysis suggested that all identified GT8 proteins belonged to seven subfamilies: galacturonosyltransferase (GAUT), galacturonosyltransferase-like (GATL), GATL-related (GATR), galactinol synthase (GolS), and plant glycogenin-like starch initiation proteins A (PGSIP-A), PGSIP-B, and PGSIP-C. We estimated that the GAUT subfamily might be further divided into four subgroups (I-IV) due to differentiation of gene structures and expression patterns. Our orthogroup analysis identified 22 orthogroups with different sizes. Of these orthogroups, several orthogroups were lost in some species, such as S. italica and Z. mays. Moreover, lots of duplicate pairs and collinear pairs were discovered among these species. These results indicated that multiple duplication modes led to the expansion of this important gene family and unequal loss of orthogroups and subfamilies might have happened during the evolutionary process. RNA-seq, microarray analysis, and qRT-PCR analyses indicated that GT8 genes are critical for plant growth and development, and for stresses responses. We found that OsGolS1 was significantly up-regulated under salt stress, while OsGAUT21, OsGATL2, and OsGATL5 had remarkable up-regulation under cold stress. The current study highlighted the expansion and evolutionary patterns of the GT8 gene family in these seven Gramineae crop genomes and provided potential candidate genes for future salt- and cold- resistant molecular breeding studies in O. sativa.

Genome-wide investigation of the MADS gene family and dehulling genes in tartary buckwheat (Fagopyrum tataricum)

Genome-wide identification, expression analysis and potential functional characterization of previously uncharacterized MADS family of tartary buckwheat, emphasized the importance of this gene family in plant growth and development. The MADS transcription factor is a key regulatory factor in the development of most plants. The MADS gene in plants controls all aspects of tissue and organ growth and reproduction and can be used to regulate plant seed cracking. However, there has been little research on the MADS genes of tartary buckwheat (Fagopyrum tataricum), which is an important edible and medicinal crop. The recently published whole genome sequence of tartary buckwheat allows us to study the tissue and expression profiles of the MADS gene in tartary buckwheat at a genome-wide level. In this study, 65 MADS genes of tartary buckwheat were identified and renamed according to the chromosomal distribution of the FtMADS genes. Here, we provide a complete overview of the gene structure, gene expression, genomic mapping, protein motif organization, and phylogenetic relationships of each member of the gene family. According to the phylogenetic relationship of MADS genes, the transcription factor family was divided into two subfamilies, the M subfamily (28 genes) and the MIKC subfamily (37 genes). The results showed that the FtMADS genes belonged to related sister pairs and the chromosomal map showed that the replication of FtMADSs was related to the replication of chromosome blocks. In different tissues and at different fruit development stages, the FtMADS genes obtained by real-time quantitative PCR (RT-qPCR) showed obvious expression patterns. A comprehensive analysis of the MADS genes in tartary buckwheat was conducted. Through systematic analysis, the potential genes that may regulate the growth and development of tartary buckwheat and the genes that may regulate the easy dehulling of tartary buckwheat fruit were screened, which laid a solid foundation for improving the quality of tartary buckwheat.

Genome-wide analysis of the pentatricopeptide repeat gene family in different maize genomes and its important role in kernel development

BACKGROUND

The pentatricopeptide repeat (PPR) gene family is one of the largest gene families in land plants (450 PPR genes in Arabidopsis, 477 PPR genes in rice and 486 PPR genes in foxtail millet) and is important for plant development and growth. Most PPR genes are encoded by plastid and mitochondrial genomes, and the gene products regulate the expression of the related genes in higher plants. However, the functions remain largely unknown, and systematic analysis and comparison of the PPR gene family in different maize genomes have not been performed.

RESULTS

In this study, systematic identification and comparison of PPR genes from two elite maize inbred lines, B73 and PH207, were performed. A total of 491 and 456 PPR genes were identified in the B73 and PH207 genomes, respectively. Basic bioinformatics analyses, including of the classification, gene structure, chromosomal location and conserved motifs, were conducted. Examination of PPR gene duplication showed that 12 and 15 segmental duplication gene pairs exist in the B73 and PH207 genomes, respectively, with eight duplication events being shared between the two genomes. Expression analysis suggested that 53 PPR genes exhibit qualitative variations in the different genetic backgrounds. Based on analysis of the correlation between PPR gene expression in kernels and kernel-related traits, four PPR genes are significantly negatively correlated with hundred kernel weight, 12 are significantly negatively correlated with kernel width, and eight are significantly correlated with kernel number. Eight of the 24 PPR genes are also located in metaQTL regions associated with yield and kernel-related traits in maize. Two important PPR genes (GRMZM2G353195 and GRMZM2G141202) might be regarded as important candidate genes associated with maize kernel-related traits.

CONCLUSIONS

Our results provide a more comprehensive understanding of PPR genes in different maize inbred lines and identify important candidate genes related to kernel development for subsequent functional validation in maize.

Global impacts of chromosomal imbalance on gene expression in Arabidopsis and other taxa

Changes in dosage of part of the genome (aneuploidy) have long been known to produce much more severe phenotypic consequences than changes in the number of whole genomes (ploidy). To examine the basis of these differences, global gene expression in mature leaf tissue for all five trisomies and in diploids, triploids, and tetraploids of Arabidopsis thaliana was studied. The trisomies displayed a greater spread of expression modulation than the ploidy series. In general, expression of genes on the varied chromosome ranged from compensation to dosage effect, whereas genes from the remainder of the genome ranged from no effect to reduced expression approaching the inverse level of chromosomal imbalance (2/3). Genome-wide DNA methylation was examined in each genotype and found to shift most prominently with trisomy 4 but otherwise exhibited little change, indicating that genetic imbalance is generally mechanistically unrelated to DNA methylation. Independent analysis of gene functional classes demonstrated that ribosomal, proteasomal, and gene body methylated genes were less modulated compared with all classes of genes, whereas transcription factors, signal transduction components, and organelle-targeted protein genes were more tightly inversely affected. Comparing transcription factors and their targets in the trisomies and in expression networks revealed considerable discordance, illustrating that altered regulatory stoichiometry is a major contributor to genetic imbalance. Reanalysis of published data on gene expression in disomic yeast and trisomic mouse cells detected similar stoichiometric effects across broad phylogenetic taxa, and indicated that these effects reflect normal gene regulatory processes.

Redundancy, Feedback, and Robustness in the Arabidopsis thaliana BZR/BEH Gene Family

Organismal development is remarkably robust, tolerating stochastic errors to produce consistent, so-called canalized adult phenotypes. The mechanistic underpinnings of developmental robustness are poorly understood, but recent studies implicate certain features of genetic networks such as functional redundancy, connectivity, and feedback. Here, we examine the BZR/BEH gene family, whose function contributes to embryonic stem development in the plant Arabidopsis thaliana, to test current assumptions on functional redundancy and trait robustness. Our analyses of BZR/BEH gene mutants and mutant combinations revealed that functional redundancy among these gene family members is not necessary for trait robustness. Connectivity is another commonly cited determinant of robustness; however, we found no correlation between connectivity among gene family members or their connectivity with other transcription factors and effects on developmental robustness. Instead, our data suggest that BEH4, the earliest diverged family member, modulates developmental robustness. We present evidence indicating that regulatory cross-talk among gene family members is integrated by BEH4 to promote wild-type levels of developmental robustness. Further, the chaperone HSP90, a known determinant of developmental robustness, appears to act via BEH4 in maintaining robustness of embryonic stem length. In summary, we demonstrate that even among closely related transcription factors, trait robustness can arise through the activity of a single gene family member, challenging common assumptions about the molecular underpinnings of robustness.

Evolution by duplication: paleopolyploidy events in plants reconstructed by deciphering the evolutionary history of VOZ transcription factors

BACKGROUND

Facilitated by the rapid progress of sequencing technology, comparative genomic studies in plants have unveiled recurrent whole genome duplication (i.e. polyploidization) events throughout plant evolution. The evolutionary past of plant genes should be analyzed in a background of recurrent polyploidy events in distinctive plant lineages. The Vascular Plant One Zinc-finger (VOZ) gene family encode transcription factors associated with a number of important traits including control of flowering time and photoperiodic pathways, but the evolutionary trajectory of this gene family remains uncharacterized.

RESULTS

In this study, we deciphered the evolutionary history of the VOZ gene family by analyses of 107 VOZ genes in 46 plant genomes using integrated methods: phylogenic reconstruction, Ks-based age estimation and genomic synteny comparisons. By scrutinizing the VOZ gene family phylogeny the core eudicot γ event was well circumscribed, and relics of the precommelinid τ duplication event were detected by incorporating genes from oil palm and banana. The more recent T and ρ polyploidy events, closely coincident with the species diversification in Solanaceae and Poaceae, respectively, were also identified. Other important polyploidy events captured included the "salicoid" event in poplar and willow, the "early legume" and "soybean specific" events in soybean, as well as the recent polyploidy event in Physcomitrella patens. Although a small transcription factor gene family, the evolutionary history of VOZ genes provided an outstanding record of polyploidy events in plants. The evolutionary past of VOZ gene family demonstrated a close correlation with critical plant polyploidy events which generated species diversification and provided answer to Darwin's "abominable mystery".

CONCLUSIONS

We deciphered the evolutionary history of VOZ transcription factor family in plants and ancestral polyploidy events in plants were recapitulated simultaneously. This analysis allowed for the generation of an idealized plant gene tree demonstrating distinctive retention and fractionation patterns following polyploidy events.

Spatio-temporal expression dynamics differ between homologues of flowering time genes in the allopolyploid Brassica napus

Polyploidy is a recurrent feature of eukaryotic evolution and has been linked to increases in complexity, adaptive radiation and speciation. Within angiosperms such events have occurred repeatedly in many plant lineages. Here we investigate the retention and spatio-temporal expression dynamics of duplicated genes predicted to regulate the floral transition in Brassica napus (oilseed rape, OSR). We show that flowering time genes are preferentially retained relative to other genes in the OSR genome. Using a transcriptome time series in two tissues (leaf and shoot apex) across development we show that 67% of these retained flowering time genes are expressed. Furthermore, between 64% (leaf) and 74% (shoot apex) of the retained gene homologues show diverged expression patterns relative to each other across development, suggesting neo- or subfunctionalization. A case study of homologues of the shoot meristem identity gene TFL1 reveals differences in cis-regulatory elements that could explain this divergence. Such differences in the expression dynamics of duplicated genes highlight the challenges involved in translating gene regulatory networks from diploid model systems to more complex polyploid crop species.

Microsynteny analysis to understand evolution and impact of polyploidization on MIR319 family within Brassicaceae

The availability of a large number of whole-genome sequences allows comparative genomic analysis to reveal and understand evolution of regulatory regions and elements. The role played by events such as whole-genome and segmental duplications followed by genome fractionation in shaping genomic landscape and in expansion of gene families is crucial toward developing insights into evolutionary trends and consequences such as sequence and functional diversification. Members of Brassicaceae are known to have experienced several rounds of whole-genome duplication (WGD) that have been termed as paleopolyploidy, mesopolyploidy, and neopolyploidy. Such repeated events led to the creation and expansion of a large number of gene families. MIR319 is reported to be one of the most ancient and conserved plant MIRNA families and plays a role in growth and development including leaf development, seedling development, and embryo patterning. We have previously reported functional diversification of members of miR319 in Brassica oleracea affecting leaf architecture; however, the evolutionary history of the MIR319 gene family across Brassicaceae remains unknown and requires investigation. We therefore identified homologous and homeologous segments of ca. 100 kb, with or without MIR319, performed comparative synteny analysis and genome fractionation studies. We detected variable rates of gene retention across members of Brassicaceae when genomic blocks of MIR319a, MIR319b, and MIR319c were compared either between themselves or against Arabidopsis thaliana genome which was taken as the base genome. The highest levels of shared genes were found between A. thaliana and Capsella rubella in both MIR319b- and MIR319c-containing genomic segments, and with the closest species of A. thaliana, A. lyrata, only in MIR319a-containing segment. Synteny analysis across 12 genomes (with 30 sub-genomes) revealed MIR319c to be the most conserved MIRNA loci (present in 27 genomes/sub-genomes) followed by MIR319a (present in 23 genomes/sub-genomes); MIR319b was found to be frequently lost (present in 20 genomes/sub-genomes) and thus is under least selection pressure for retention. Genome fractionation revealed extensive and differential loss of MIRNA homeologous loci and flanking genes from various sub-genomes of Brassica species that is in accordance with their older history of polyploidy when compared to Camelina sativa, a recent neopolyploid, where the effect of genome fractionation was least. Finally, estimation of phylogenetic relationship using precursor sequences of MIR319 reveals MIR319a and MIR319b form sister clades, with MIR319c forming a separate clade. An intra-species synteny analysis between MIR319a-, MIR319b-, and MIR319c-containing genomic segments suggests segmental duplications at the base of Brassicaceae to be responsible for the origin of MIR319a and MIR319b.

Genome-Wide Identification, Classification, and Expression Divergence of Glutathione-Transferase Family in Brassica rapa under Multiple Hormone Treatments

The GSTs is one of the most important multifunctional protein families which has been playing a crucial role in the different aspects of plant growth. This extensive study about GSTs may establish a solid foundation for the brief functional analysis of BraGSTs in future. In this study, a total of 75 genes were identified in B. rapa. Phylogenetic analysis characterized them into eight different subclasses, while Tau and Phi subclasses were the most numerous. The exon-intron structure and the motif composition of BraGSTs were exhibited accordingly to their subclasses. Notably, we also investigated 15 tandem paralogous pairs of genes, which highlighted that all the pairs were purifying in nature as their synonymous values were lower than 1.00. Duplication analysis indicated that about 45.33% of genes mainly occurred through tandem duplication in B. rapa. Predominately, the tandem cluster of genes in subclass Tau was greater than the other subclasses. Furthermore, among eight multiple hormonal treatments (ABA, GA, BR, ETH, IAA, IBA, NPA, and JA), most number of BraGSTs was activated by NPA, BR, and ABA treatments. This analysis has provided comprehensive information about GSTs family which may assist in elucidating their exact functions in B. rapa.

The fate of ribosomal RNA genes in spontaneous polyploid dogrose hybrids [Rosa L. sect. Caninae (DC.) Ser.] exhibiting non-symmetrical meiosis

Dogroses represent an exceptional system for studying the effects of genome doubling and hybridization: their asymmetrical meiosis enables recombination in bi-parentally inherited chromosomes but prevents it in maternally inherited ones. We employed fluorescent in situ hybridization, genome skimming, amplicon sequencing of genomic and cDNA as well as conventional cloning of nuclear ribosomal DNA in two phylogenetically distinct pentaploid (2n = 5x = 35) species, Rosa canina and Rosa inodora, and their naturally occurring reciprocal hybrids, Rosa dumalis (5x) and Rosa agrestis (5x, 6x). Both progenitor species differed in composition, meiotic behaviour and expression of rDNA loci: R. canina (five 18S and 5-8 5S loci) was dominated by the Canina ribotypes, but R. inodora (four 18S loci and 7-8 5S loci) by the Rubiginosa ribotype. The co-localized 5S/18S loci occurred on either bivalent-forming (R. canina) or univalent-forming (R. inodora) chromosomes. Ribosomal DNA loci were additively inherited; however, the Canina ribotypes were dominantly expressed, even in genotypes with relatively low copy number of these genes. Moreover, we observed rDNA homogenization towards the paternally transmitted Canina ribotype in 6x R. agrestis. The here-observed variation in arrangement and composition of rDNA types between R. canina and R. inodora suggests the involvement of different genomes in bivalent formation. This results supports the hypothesis that the asymmetrical meiosis arose at least twice by independent ancient hybridization events.

Widespread Whole Genome Duplications Contribute to Genome Complexity and Species Diversity in Angiosperms

Gene duplications provide evolutionary potentials for generating novel functions, while polyploidization or whole genome duplication (WGD) doubles the chromosomes initially and results in hundreds to thousands of retained duplicates. WGDs are strongly supported by evidence commonly found in many species-rich lineages of eukaryotes, and thus are considered as a major driving force in species diversification. We performed comparative genomic and phylogenomic analyses of 59 public genomes/transcriptomes and 46 newly sequenced transcriptomes covering major lineages of angiosperms to detect large-scale gene duplication events by surveying tens of thousands of gene family trees. These analyses confirmed most of the previously reported WGDs and provided strong evidence for novel ones in many lineages. The detected WGDs supported a model of exponential gene loss during evolution with an estimated half-life of approximately 21.6 million years, and were correlated with both the emergence of lineages with high degrees of diversification and periods of global climate changes. The new datasets and analyses detected many novel WGDs widely spread during angiosperm evolution, uncovered preferential retention of gene functions in essential cellular metabolisms, and provided clues for the roles of WGD in promoting angiosperm radiation and enhancing their adaptation to environmental changes.

Evolutionary history of the cobalamin-independent methionine synthase gene family across the land plants

Plants are successful paleopolyploids. The wide diversity of land plants is driven strongly by their gene duplicates undergoing distinct evolutionary fates after duplication. We used genomic resources from 35 model plant species to unravel the evolutionary fate of gene copies (paralogs) of the cobalamin-independent methionine synthase (metE) gene family across the land plants. To explore genealogical relationships and characterize positive selection as a driving force in the evolution of metE paralogs within a single species, we carried out complementary analyses on genomic data of 32 genotypes of soybean. The size of the metE gene family remained small across the land plants; most of the studied species possessed 1-6 paralogs. Gene products were either cytosolic or chloroplastic; this dual subcellular distribution arose early during the divergence of the land plants and reached all extant lineages. Biased gene loss and gene retention events took place multiple times; recurrent evolution remodeled redundant metE paralogs to recover and maintain the dual subcellular distribution of MetE. Shared whole-genome duplication events gave rise to the metE paralogs of both soybean and Medicago truncatula. In soybean, the ancestral paralog pair GlymaPP2A encoded a cytosolic isoform of MetE, was under strong purifying selection, and retained high levels of expression across eight RNA-seq expression libraries. The daughters GlymaPP1 and GlymaPP2B showed accelerated rates of evolution, accumulated many sites predicted to be under positive selection, and possessed low levels of expression. Our results suggest that the metE paralogs of soybean follow Ohno's neofunctionalization model of gene duplicate evolution.

Genetic redundancy of senescence-associated transcription factors in Arabidopsis

Leaf senescence is a genetically programmed process that constitutes the last stage of leaf development, and involves massive changes in gene expression. As a result of the intensive efforts that have been made to elucidate the molecular genetic mechanisms underlying leaf senescence, 184 genes that alter leaf senescence phenotypes when mutated or overexpressed have been identified in Arabidopsis thaliana over the past two decades. Concurrently, experimental evidence on functional redundancy within senescence-associated genes (SAGs) has increased. In this review, we focus on transcription factors that play regulatory roles in Arabidopsis leaf senescence, and describe the relationships among gene duplication, gene expression level, and senescence phenotypes. Previous findings and our re-analysis demonstrate the widespread existence of duplicate SAG pairs and a correlation between gene expression levels in duplicate genes and senescence-related phenotypic severity of the corresponding mutants. We also highlight effective and powerful tools that are available for functional analyses of redundant SAGs. We propose that the study of duplicate SAG pairs offers a unique opportunity to understand the regulation of leaf senescence and can guide the investigation of the functions of redundant SAGs via reverse genetic approaches.

Evolutionary conservation and functional divergence of the LFK gene family play important roles in the photoperiodic flowering pathway of land plants

ZEITLUPE (ZTL), LOV KELCH PROTEIN 2 (LKP2), and FLAVIN-BINDING KELCH REPEAT F-BOX 1 (FKF1)-blue-light photoreceptors-play important roles in regulating the circadian clock and photoperiodic flowering pathway in plants. In this study, phylogenetic analysis revealed that the LOV (Light, Oxygen, or Voltage) and Kelch repeat-containing F-box (LFK) gene family can be classified into two clades, ZTL/LKP2 and FKF1, with clear differentiation between monocots and dicots within each clade. The LFK family genes underwent strong purifying selection; however, signatures of positive selection to adapt to local conditions still existed in 18 specific codons. In 87 diverse maize inbred lines, significant differences were identified (P ≤ 0.01) for days to female flowering between the haplotypes consisting of eight positive selection sites at ZmFKF1b corresponding to tropical and temperate maize groups of the phylogenetic tree, indicating a key role of ZmFKF1b in maize adaptive evolution. In addition, positive coevolution was detected in the domains of the LFK family for long-term cooperation to targets. The Type-I and Type-II functional divergence analysis revealed subfunctionalization or neofunctionalization of the LFKs, and the ZTL subfamily is most likely to maintain the ancestral function of LFKs. Over 50% of critical amino acid sites involved in the functional divergence were identified in the Kelch repeat domain, resulting in the distinction of substrates for ubiquitination and degradation. These results suggest that evolutionary conservation contributes to the maintenance of critical physiological functions, whereas functional divergence after duplication helps to generate diverse molecular regulation mechanisms.

Intrinsic protein disorder reduces small-scale gene duplicability

Whereas the rate of gene duplication is relatively high, only certain duplications survive the filter of natural selection and can contribute to genome evolution. However, the reasons why certain genes can be retained after duplication whereas others cannot remain largely unknown. Many proteins contain intrinsically disordered regions (IDRs), whose structures fluctuate between alternative conformational states. Due to their high flexibility, IDRs often enable protein–protein interactions and are the target of post-translational modifications. Intrinsically disordered proteins (IDPs) have characteristics that might either stimulate or hamper the retention of their encoding genes after duplication. On the one hand, IDRs may enable functional diversification, thus promoting duplicate retention. On the other hand, increased IDP availability is expected to result in deleterious unspecific interactions. Here, we interrogate the proteomes of human, Drosophila melanogaster, Caenorhabditis elegans, Saccharomyces cerevisiae, Arabidopsis thaliana and Escherichia coli, in order to ascertain the impact of protein intrinsic disorder on gene duplicability. We show that, in general, proteins encoded by duplicated genes tend to be less disordered than those encoded by singletons. The only exception is proteins encoded by ohnologs, which tend to be more disordered than those encoded by singletons or genes resulting from small-scale duplications. Our results indicate that duplication of genes encoding IDPs outside the context of whole-genome duplication (WGD) is often deleterious, but that IDRs facilitate retention of duplicates in the context of WGD. We discuss the potential evolutionary implications of our results.

Evolutionary history and functional divergence of the cytochrome P450 gene superfamily between Arabidopsis thaliana and Brassica species uncover effects of whole genome and tandem duplications

Background

The cytochrome P450 monooxygenase (P450) superfamily is involved in the biosynthesis of various primary and secondary metabolites. However, little is known about the effects of whole genome duplication (WGD) and tandem duplication (TD) events on the evolutionary history and functional divergence of P450s in Brassica after splitting from a common ancestor with Arabidopsis thaliana.

Results

Using Hidden Markov Model search and manual curation, we detected that Brassica species have nearly 1.4-fold as many P450 members as A. thaliana. Most P450s in A. thaliana and Brassica species were located on pseudo-chromosomes. The inferred phylogeny indicated that all P450s were clustered into two different subgroups. Analysis of WGD event revealed that different P450 gene families had appeared after evolutionary events of species. For the TD event analyses, the P450s from TD events in Brassica species can be divided into ancient and recent parts. Our comparison of influence of WGD and TD events on the P450 gene superfamily between A. thaliana and Brassica species indicated that the family-specific evolution in the Brassica lineage can be attributed to both WGD and TD, whereas WGD was recognized as the major mechanism for the recent evolution of the P450 super gene family. Expression analysis of P450s from A. thaliana and Brassica species indicated that WGD-type P450s showed the same expression pattern but completely different expression with TD-type P450s across different tissues in Brassica species. Selection force analysis suggested that P450 orthologous gene pairs between A. thaliana and Brassica species underwent negative selection, but no significant differences were found between P450 orthologous gene pairs in A. thaliana–B. rapa and A. thaliana–B. oleracea lineages, as well as in different subgenomes in B. rapa or B. oleracea compared with A. thaliana.

Conclusions

This study is the first to investigate the effects of WGD and TD on the evolutionary history and functional divergence of P450 gene families in A. thaliana and Brassica species. This study provides a biology model to study the mechanism of gene family formation, particularly in the context of the evolutionary history of angiosperms, and offers novel insights for the study of angiosperm genomes.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4094-7) contains supplementary material, which is available to authorized users.

Genome-Wide Analysis of Soybean LATERAL ORGAN BOUNDARIES Domain-Containing Genes: A Functional Investigation of GmLBD12

Plant-specific () genes play critical roles in various plant growth and development processes. However, the number and characteristics of genes in soybean [ (L.) Merr.] remain unknown. Here, we identified 90 homologous genes in the soybean genome that phylogenetically clustered into two classes (I and II). The majority of the genes were evenly distributed across all 20 soybean chromosomes, and 77 (81.11%) of them were detected in segmental duplicated regions. Furthermore, the exon-intron organization and motif composition for each were analyzed. A close phylogenetic relationship was identified between the soybean genes and 41 previously reported genes of different plants in the same group, providing insights into their putative functions. Expression analysis indicated that more than half of the genes were expressed, with the two gene classes showing differential tissue expression characteristics; in addition, they were differentially induced by biotic and abiotic stresses. To further explore the functions of genes in soybean, was selected for functional characterization. GmLBD12 was mainly localized to the nucleus and showed high expression in root and seed tissues. Overexpressing in (L.) Heynh resulted in increases in lateral root (LR) number and plant height. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis demonstrated that was induced by drought, salt, cold, indole acetic acid (IAA), abscisic acid (ABA), and salicylic acid SA treatments. This study provides the first comprehensive analysis of the soybean gene family and a valuable foundation for future functional studies of genes.

Both mechanism and age of duplications contribute to biased gene retention patterns in plants

BACKGROUND

All extant seed plants are successful paleopolyploids, whose genomes carry duplicate genes that have survived repeated episodes of diploidization. However, the survival of gene duplicates is biased with respect to gene function and mechanism of duplication. Transcription factors, in particular, are reported to be preferentially retained following whole-genome duplications (WGDs), but disproportionately lost when duplicated by tandem events. An explanation for this pattern is provided by the Gene Balance Hypothesis (GBH), which posits that duplicates of highly connected genes are retained following WGDs to maintain optimal stoichiometry among gene products; but such connected gene duplicates are disfavored following tandem duplications.

RESULTS

We used genomic data from 25 taxonomically diverse plant species to investigate the roles of duplication mechanism, gene function, and age of duplication in the retention of duplicate genes. Enrichment analyses were conducted to identify Gene Ontology (GO) functional categories that were overrepresented in either WGD or tandem duplications, or across ranges of divergence times. Tandem paralogs were much younger, on average, than WGD paralogs and the most frequently overrepresented GO categories were not shared between tandem and WGD paralogs. Transcription factors were overrepresented among ancient paralogs regardless of mechanism of origin or presence of a WGD. Also, in many cases, there was no bias toward transcription factor retention following recent WGDs.

CONCLUSIONS

Both the fixation and the retention of duplicated genes in plant genomes are context-dependent events. The strong bias toward ancient transcription factor duplicates can be reconciled with the GBH if selection for optimal stoichiometry among gene products is strongest following the earliest polyploidization events and becomes increasingly relaxed as gene families expand.

Multiple Polyploidization Events across Asteraceae with Two Nested Events in the Early History Revealed by Nuclear Phylogenomics

Biodiversity results from multiple evolutionary mechanisms, including genetic variation and natural selection. Whole-genome duplications (WGDs), or polyploidizations, provide opportunities for large-scale genetic modifications. Many evolutionarily successful lineages, including angiosperms and vertebrates, are ancient polyploids, suggesting that WGDs are a driving force in evolution. However, this hypothesis is challenged by the observed lower speciation and higher extinction rates of recently formed polyploids than diploids. Asteraceae includes about 10% of angiosperm species, is thus undoubtedly one of the most successful lineages and paleopolyploidization was suggested early in this family using a small number of datasets. Here, we used genes from 64 new transcriptome datasets and others to reconstruct a robust Asteraceae phylogeny, covering 73 species from 18 tribes in six subfamilies. We estimated their divergence times and further identified multiple potential ancient WGDs within several tribes and shared by the Heliantheae alliance, core Asteraceae (Asteroideae-Mutisioideae), and also with the sister family Calyceraceae. For two of the WGD events, there were subsequent great increases in biodiversity; the older one proceeded the divergence of at least 10 subfamilies within 10 My, with great variation in morphology and physiology, whereas the other was followed by extremely high species richness in the Heliantheae alliance clade. Our results provide different evidence for several WGDs in Asteraceae and reveal distinct association among WGD events, dramatic changes in environment and species radiations, providing a possible scenario for polyploids to overcome the disadvantages of WGDs and to evolve into lineages with high biodiversity.

The functional divergence between SPA1 and SPA2 in Arabidopsis photomorphogenesis maps primarily to the respective N-terminal kinase-like domain

BACKGROUND

Plants have evolved complex mechanisms to adapt growth and development to the light environment. The COP1/SPA complex is a key repressor of photomorphogenesis in dark-grown Arabidopsis plants and acts as an E3 ubiquitin ligase to ubiquitinate transcription factors involved in the light response. In the light, COP1/SPA activity is inhibited by photoreceptors, thereby allowing accumulation of these transcription factors and a subsequent light response. Previous results have shown that the four members of the SPA family exhibit partially divergent functions. In particular, SPA1 and SPA2 strongly differ in their responsiveness to light, while they have indistinguishable activities in darkness. The much higher light-responsiveness of SPA2 is partially explained by the much stronger light-induced degradation of SPA2 when compared to SPA1. Here, we have conducted SPA1/SPA2 domain swap experiments to identify the protein domain(s) responsible for the functional divergence between SPA1 and SPA2.

RESULTS

We have individually swapped the three domains between SPA1 and SPA2 - the N-terminal kinase-like domain, the coiled-coil domain and the WD-repeat domain - and expressed them in spa mutant Arabidopsis plants. The phenotypes of transgenic seedlings show that the respective N-terminal kinase-like domain is primarily responsible for the respective light-responsiveness of SPA1 and SPA2. Furthermore, the most divergent part of the N-terminal domain was sufficient to confer a SPA1- or SPA2-like activity to the respective SPA protein. The stronger light-induced degradation of SPA2 when compared to SPA1 was also primarily conferred by the SPA2 N-terminal domain. At last, the different affinities of SPA1 and SPA2 for cryptochrome 2 are defined by the N-terminal domain of the respective SPA protein. In contrast, both SPA1 and SPA2 similarly interacted with COP1 in light-grown seedlings.

CONCLUSIONS

Our results show that the distinct activities and protein stabilities of SPA1 and SPA2 in light-grown seedlings are primarily encoded by their N-terminal kinase-like domains. Similarly, the different affinities of SPA1 and SPA2 for cry2 are explained by their respective N-terminal domain. Hence, after a duplication event during evolution, the N-terminal domains of SPA1 and SPA2 underwent subfunctionalization, possibly to allow optimal adaptation of growth and development to a changing light environment.

Genome of Plant Maca (Lepidium meyenii) Illuminates Genomic Basis for High-Altitude Adaptation in the Central Andes

Maca (Lepidium meyenii Walp, 2n = 8x = 64), belonging to the Brassicaceae family, is an economic plant cultivated in the central Andes sierra in Peru (4000-4500 m). Considering that the rapid uplift of the central Andes occurred 5-10 million years ago (Ma), an evolutionary question arises regarding how plants such as maca acquire high-altitude adaptation within a short geological period. Here, we report the high-quality genome assembly of maca, in which two closely spaced maca-specific whole-genome duplications (WGDs; ∼6.7 Ma) were identified. Comparative genomic analysis between maca and closely related Brassicaceae species revealed expansions of maca genes and gene families involved in abiotic stress response, hormone signaling pathway, and secondary metabolite biosynthesis via WGDs. The retention and subsequent functional divergence of many duplicated genes may account for the morphological and physiological changes (i.e., small leaf shape and self-fertility) in maca in a high-altitude environment. In addition, some duplicated maca genes were identified with functions in morphological adaptation (i.e., LEAF CURLING RESPONSIVENESS) and abiotic stress response (i.e., GLYCINE-RICH RNA-BINDING PROTEINS and DNA-DAMAGE-REPAIR/TOLERATION 2) under positive selection. Collectively, the maca genome provides useful information to understand the important roles of WGDs in the high-altitude adaptation of plants in the Andes.

Major Chromosomal Rearrangements Distinguish Willow and Poplar After the Ancestral "Salicoid" Genome Duplication

Populus (poplar) and Salix (willow) are sister genera in the Salicaceae family. In both lineages extant species are predominantly diploid. Genome analysis previously revealed that the two lineages originated from a common tetraploid ancestor. In this study, we conducted a syntenic comparison of the corresponding 19 chromosome members of the poplar and willow genomes. Our observations revealed that almost every chromosomal segment had a parallel paralogous segment elsewhere in the genomes, and the two lineages shared a similar syntenic pinwheel pattern for most of the chromosomes, which indicated that the two lineages diverged after the genome reorganization in the common progenitor. The pinwheel patterns showed distinct differences for two chromosome pairs in each lineage. Further analysis detected two major interchromosomal rearrangements that distinguished the karyotypes of willow and poplar. Chromosome I of willow was a conjunction of poplar chromosome XVI and the lower portion of poplar chromosome I, whereas willow chromosome XVI corresponded to the upper portion of poplar chromosome I. Scientists have suggested that Populus is evolutionarily more primitive than Salix. Therefore, we propose that, after the “salicoid” duplication event, fission and fusion of the ancestral chromosomes first give rise to the diploid progenitor of extant Populus species. During the evolutionary process, fission and fusion of poplar chromosomes I and XVI subsequently give rise to the progenitor of extant Salix species. This study contributes to an improved understanding of genome divergence after ancient genome duplication in closely related lineages of higher plants.