The causes and molecular consequences of polyploidy in flowering plants

The link between ancient whole-genome duplications and cold adaptations in the Caryophyllaceae

PREMISE

The Caryophyllaceae (the carnation family) have undergone multiple transitions into colder climates and convergence on cushion plant adaptation, indicating that they may provide a natural system for cold adaptation research. Previous research has suggested that putative ancient whole-genome duplications (WGDs) are correlated with niche shifts into colder climates across the Caryophyllales. Here, we explored the genomic changes potentially involved in one of these discovered shifts in the Caryophyllaceae.

METHODS

We constructed a data set combining 26 newly generated transcriptomes with 45 published transcriptomes, including 11 cushion plant species across seven genera. With this data set, we inferred a dated phylogeny for the Caryophyllaceae and mapped ancient WGDs and gene duplications onto the phylogeny. We also examined functional groups enriched for gene duplications related to the climatic shift.

RESULTS

The ASTRAL topology was mostly congruent with the current consensus of relationships within the family. We inferred 15 putative ancient WGDs in the family, including eight that have not been previously published. The oldest ancient WGD (ca. 64.4-56.7 million years ago), WGD1, was found to be associated with a shift into colder climates by previous research. Gene regions associated with ubiquitination were overrepresented in gene duplications retained after WGD1 and those convergently retained by cushion plants in Colobanthus and Eremogone, along with other functional annotations.

CONCLUSIONS

Gene family expansions induced by ancient WGDs may have contributed to the shifts to cold climatic niches in the Caryophyllaceae. Transcriptomic data are crucial resources that help unravel heterogeneity in deep-time evolutionary patterns in plants.

Genome assembly provides insights into the genome evolution of Baccaurea ramiflora Lour

Baccaurea ramiflora Lour., an evergreen tree of the Baccaurea genus of the Phyllanthaceae family, is primarily distributed in South Asia, Southeast Asia, and southern China, including southern Yunnan Province. It is a wild or semi-cultivated tree species with ornamental, edible, and medicinal value, exhibiting significant development potential. In this study, we present the whole-genome sequencing of B. ramiflora, employing a combination of PacBio SMRT and Illumina HiSeq 2500 sequencing techniques. The assembled genome size was 975.8 Mb, with a contig N50 of 509.33 kb and the longest contig measuring 7.74 Mb. The genome comprises approximately 73.47% highly repetitive sequences, of which 52.1% are long terminal repeat-retrotransposon sequences. A total of 29,172 protein-coding genes were predicted, of which 25,980 (89.06%) have been annotated, Additionally, 3452 non-coding RNAs were identified. Comparative genomic analysis revealed a close relationship between B. ramiflora and the Euphorbiaceae family, with both being sister groups that diverged approximately 59.9 million years ago. During the evolutionary process, B. ramiflora exhibited positive selection in 278 candidate genes. Synonymous substitution rate and collinearity analysis demonstrated that B. ramiflora underwent a single ancient genome-wide triploidization event, without recent genome-wide duplication events. This high-quality B. ramiflora genome provides a valuable resource for basic research and tree improvement programs focusing on the Phyllanthaceae family.

Gamete sex and elevation affect genetically based variation for unreduced gamete production in a mixed-ploidy plant

PREMISE

Unreduced gametes are the primary mechanism of neopolyploid formation. Their production in diploid populations is arguably maladaptive, but the magnitude and patterns of genetically based variation maintained in natural populations are poorly understood.

METHODS

We examined variation in male and female unreduced gamete production among plants from different elevations in fireweed, Chamerion angustifolium, grown in a common environment. Using seeds from three high-elevation and three low-elevation diploid populations in one study, and a single diploid population in another, we estimated realized rates of unreduced male (sperm) and female (egg) gamete production by reciprocally pollinating diploid and tetraploid plants and estimating the incidence of tetraploid seeds using flow cytometry.

RESULTS

Unreduced gamete frequencies per plant were similar in the two studies (0.12% vs. 0.08%). High-elevation populations had a greater percentage of fruit with seeds from unreduced gametes, but a lower percentage of seeds per fruit than low-elevation populations. Female unreduced gamete frequencies differed among elevations, but male frequencies did not, and the gamete sexes were not correlated at the plant level.

CONCLUSIONS

We conclude that genetically based variation for unreduced gametes is maintained within and among natural populations, despite their fitness disadvantages, suggesting that local selection may be ineffective at purging them under some conditions.

Genome resources for the elite bread wheat cultivar Aikang 58 and mining of elite homeologous haplotypes for accelerating wheat improvement

Despite recent progress in crop genomics studies, the genomic changes brought about by modern breeding selection are still poorly understood, thus hampering genomics-assisted breeding, especially in polyploid crops with compound genomes such as common wheat (Triticum aestivum). In this work, we constructed genome resources for the modern elite common wheat variety Aikang 58 (AK58). Comparative genomics between AK58 and the landrace cultivar Chinese Spring (CS) shed light on genomic changes that occurred through recent varietal improvement. We also explored subgenome diploidization and divergence in common wheat and developed a homoeologous locus-based genome-wide association study (HGWAS) approach, which was more effective than single homoeolog-based GWAS in unraveling agronomic trait-associated loci. A total of 123 major HGWAS loci were detected using a genetic population derived from AK58 and CS. Elite homoeologous haplotypes (HHs), formed by combinations of subgenomic homoeologs of the associated loci, were found in both parents and progeny, and many could substantially improve wheat yield and related traits. We built a website where users can download genome assembly sequence and annotation data for AK58, perform blast analysis, and run JBrowse. Our work enriches genome resources for wheat, provides new insights into genomic changes during modern wheat improvement, and suggests that efficient mining of elite HHs can make a substantial contribution to genomics-assisted breeding in common wheat and other polyploid crops.

Genome-Wide Identification and Expression Analysis of Cysteine-Rich Polycomb-like Protein (CPP) Gene Family in Tomato

The cysteine-rich polycomb-like protein (CPP) gene family is a class of transcription factors containing conserved cysteine-rich CRC structural domains that is involved in the regulation of plant growth and stress tolerance to adversity. Relative to other gene families, the CPP gene family has not received sufficient attention. In this study, six SlCPPs were identified for the first time using the most recent genome-wide identification data of tomato. Subsequently, a phylogenetic analysis classified SlCPPs into four subfamilies. The analysis of cis-acting elements in the promoter indicates that SlCPPs are involved in plant growth and development and also stress response. We present for the first time the prediction of the tertiary structure of these SlCPPs proteins using the AlphaFold2 artificial intelligence system developed by the DeepMind team. Transcriptome data analysis showed that SlCPPs were differentially expressed in different tissues. Gene expression profiling showed that all SlCPPs except SlCPP5 were up-regulated under drought stress; SlCPP2, SlCPP3 and SlCPP4 were up-regulated under cold stress; SlCPP2 and SlCPP5 were up-regulated under salt stress; all SlCPPs were up-regulated under inoculation with Cladosporium fulvum; and SlCPP1, SlCPP3, and SlCPP4 were up-regulated under inoculation with Stemphylium lycopersici. We performed a virus-induced gene silencing experiment on SlCPP3, and the results indicated that SlCPP3 was involved in the response to drought stress. Finally, we predicted the interaction network of the key gene SlCPP3, and there was an interaction relationship between SlCPP3 and 10 genes, such as RBR1 and MSI1. The positive outcome showed that SlCPPs responded to environmental stress. This study provides a theoretical and empirical basis for the response mechanisms of tomato in abiotic stresses.

Karyotype evolution of the Asterids insights from the first genome sequences of the family Cornaceae

Cornaceae is a core representative family in Cornales, the earliest branching lineage in the Asterids on the life tree of angiosperms. This family includes the only genus Cornus, a group of ~55 species. These species occur widely in Northern Hemisphere and have been used as resources for horticultural ornaments, medicinal and industrial manufacturing. However, no any genome sequences are available for this family. Here, we reported a chromosome-level genome for Cornus controversa. This was generated using high-fidelity plus Hi-C sequencing, and totally ~771.80 Mb assembled sequences and 39,886 protein-coding genes were obtained. We provided evidence for a whole-genome duplication event (WGD) unique to C. controversa. The evolutionary features of this genome indicated that the expanded and unique genes might have contributed to response to stress, stimulus and defense. By using chromosome-level syntenic blocks shared between eight living genomes, we found high degrees of genomic diversification from the ancestral core-eudicot genome to the present-day genomes, suggesting an important role of WGD in genomic plasticity that leads to speciation and diversification. These results provide foundational insights on the evolutionary history of Cornaceae, as well as on the Asterids diversification.

Pervasive genome duplications across the plant tree of life and their links to major evolutionary innovations and transitions

Whole-genome duplication (WGD) has occurred repeatedly during plant evolution and diversification, providing genetic layers for evolving new functions and phenotypes. Advances in long-read sequencing technologies have enabled sequencing and assembly of over 1000 plant genomes spanning nearly 800 species, in which a large set of ancient WGDs has been uncovered. Here, we review the recently reported WGDs that occurred in major plant lineages and key evolutionary positions, and highlight their contributions to morphological innovation and adaptive evolution. Current gaps and challenges in integrating enormous volumes of sequenced plant genomes, accurately inferring WGDs, and developing web-based analysis tools are emphasized. Looking to the future, ambitious genome sequencing projects and global efforts may substantially recapitulate the plant tree of life based on broader sampling of phylogenetic diversity, reveal much of the timetable of ancient WGDs, and address the biological significance of WGDs in plant adaptation and radiation.

Species Tree Estimation and the Impact of Gene Loss Following Whole-Genome Duplication

Whole-genome duplication (WGD) occurs broadly and repeatedly across the history of eukaryotes and is recognized as a prominent evolutionary force, especially in plants. Immediately following WGD, most genes are present in two copies as paralogs. Due to this redundancy, one copy of a paralog pair commonly undergoes pseudogenization and is eventually lost. When speciation occurs shortly after WGD; however, differential loss of paralogs may lead to spurious phylogenetic inference resulting from the inclusion of pseudoorthologs–paralogous genes mistakenly identified as orthologs because they are present in single copies within each sampled species. The influence and impact of including pseudoorthologs versus true orthologs as a result of gene extinction (or incomplete laboratory sampling) are only recently gaining empirical attention in the phylogenomics community. Moreover, few studies have yet to investigate this phenomenon in an explicit coalescent framework. Here, using mathematical models, numerous simulated data sets, and two newly assembled empirical data sets, we assess the effect of pseudoorthologs on species tree estimation under varying degrees of incomplete lineage sorting (ILS) and differential gene loss scenarios following WGD. When gene loss occurs along the terminal branches of the species tree, alignment-based (BPP) and gene-tree-based (ASTRAL, MP-EST, and STAR) coalescent methods are adversely affected as the degree of ILS increases. This can be greatly improved by sampling a sufficiently large number of genes. Under the same circumstances, however, concatenation methods consistently estimate incorrect species trees as the number of genes increases. Additionally, pseudoorthologs can greatly mislead species tree inference when gene loss occurs along the internal branches of the species tree. Here, both coalescent and concatenation methods yield inconsistent results. These results underscore the importance of understanding the influence of pseudoorthologs in the phylogenomics era. [Coalescent method; concatenation method; incomplete lineage sorting; pseudoorthologs; single-copy gene; whole-genome duplication.]

Exploiting Unreduced Gametes for Improving Ornamental Plants

The formation of gametes with somatic chromosome number or unreduced gametes (2n gametes) is an important process involved in the origin of polyploid plants in nature. Unreduced gametes are the result of meiotic mutations occurring during micro- and mega-sporogenesis. 2n gametes have been identified or artificially induced in a large number of plant species. Breeding of plants through 2n gametes can be advantageous because it combines genetic effects of polyploidy with meiotic recombination and sexual hybridization to produce tremendous genetic variation and heterosis. 2n gametes also occur in ornamental plants, but the potential of using 2n gametes in ornamental plant breeding has not been extensively exploited. Ornamental plants are primarily produced for their esthetic appearance and novelty, not for food and yield, and they can be readily propagated through vegetative means. Triploids, tetraploids, and plants with even higher ploidy levels produced through 2n gametes can be propagated through tissue culture to fix their phenotypes, thus leading to the development of new cultivars. In this review article, we intend to discuss the mechanisms underlying the formation of 2n gametes, techniques for 2n gamete identification, methods for enhancing 2n gamete formation, and the current status in the use of 2n gametes for development of novel ornamental plants. We believe that polyploidy breeding through 2n gametes represents a viable way of developing new cultivars, new species, and even new genera of ornamental plants.

In vitro octaploid induction of Populus hopeiensis with colchicine

BACKGROUND

Autopolyploids, especially artificial lines, provide model systems for understanding the mechanisms of gene dosage effects on trait variation owing to their relatively uniform genetic background. Here, a protocol for in vitro octaploid induction of Populus hopeiensis from leaf blades with colchicine treatment was established through investigation of the effects of different pre-culture durations, colchicine concentrations, and exposure times.

RESULTS

We found that pre-culture duration, colchicine concentration, and exposure time had significant effects on the survival rate, shoot regeneration rate, and octaploid induction rate of P. hopeiensis leaf blades. The highest octaploid induction rate (8.61%) was observed when leaf blades pre-cultured for 9 days were treated for 4 days with 100 μM colchicine. The ploidy level of all regenerated plantlets was analyzed by flow cytometry and further confirmed by chromosome counting. A total of 14 octaploids were obtained. The stomatal length, width, and density of leaf blades significantly differed between tetraploid and octaploid plants. Compared with diploid and tetraploid plants, octaploids had a slower growth rate, smaller leaf blade size, and shorter internodes.

CONCLUSIONS

We established an effective protocol for inducing octaploids in vitro from autotetraploid P. hopeiensis leaf blades by colchicine treatment.

Genetic Divergence of Lineage-Specific Tandemly Duplicated Gene Clusters in Four Diploid Potato Genotypes

Potato (Solanum tuberosum L.) is the most important non-grain food crop. Tandem duplication significantly contributes to genome evolution. The objectives of this study were to (i) identify tandemly duplicated genes and compare their genomic distributions across potato genotypes, (ii) investigate the bias in functional specificities, (iii) explore the relationships among coding sequence, promoter and expression divergences associated with tandemly duplicated genes, (iv) examine the role of tandem duplication in generating and expanding lineage-specific gene families, (v) investigate the evolutionary forces affecting tandemly duplicated genes, and (vi) assess the similarities and differences with respect to above mentioned aspects between cultivated genotypes and their wild-relative. In this study, we used well-annotated and chromosome-scale de novo genome assemblies of multiple potato genotypes. Our results showed that tandemly duplicated genes are abundant and dispersed through the genome. We found that several functional specificities, such as disease resistance, stress-tolerance, and biosynthetic pathways of tandemly duplicated genes were differentially enriched across multiple potato genomes. Our results indicated the existence of a significant correlation among expression, promoter, and protein divergences in tandemly duplicated genes. We found about one fourth of tandemly duplicated gene clusters as lineage-specific among multiple potato genomes, and these tended to localize toward centromeres and revealed distinct selection signatures and expression patterns. Furthermore, our results showed that a majority of duplicated genes were retained through sub-functionalization followed by genetic redundancy, while only a small fraction of duplicated genes was retained though neo-functionalization. The lineage-specific expansion of gene families by tandem duplication coupled with functional bias might have significantly contributed to potato's genotypic diversity, and, thus, to adaption to environmental stimuli.

Intraspecific variation mediates density dependence in a genetically diverse plant species

Interactions between neighboring plants are critical for biodiversity maintenance in plant populations and communities. Intraspecific trait variation and genome duplication are common in plant species and can drive eco-evolutionary dynamics through genotype-mediated plant-plant interactions. However, few studies have examined how species-wide intraspecific variation may alter interactions between neighboring plants. We investigate how subspecies and ploidy variation in a genetically diverse species, big sagebrush (Artemisia tridentata), can alter the demographic outcomes of plant interactions. Using a replicated, long-term common garden experiment that represents range-wide diversity of A. tridentata, we ask how intraspecific variation, environment, and stand age mediate neighbor effects on plant growth and survival. Spatially explicit models revealed that ploidy variation and subspecies identity can mediate plant-plant interactions but that the effect size varied in time and across experimental sites. We found that demographic impacts of neighbor effects were strongest during early stages of stand development and in sites with greater growth rates. Within subspecies, tetraploid populations showed greater tolerance to neighbor crowding compared to their diploid variants. Our findings provide evidence that intraspecific variation related to genome size and subspecies identity impacts spatial demography in a genetically diverse plant species. Accounting for intraspecific variation in studies of conspecific density dependence will improve our understanding of how local populations will respond to novel genotypes and biotic interaction regimes. As introduction of novel genotypes into local populations becomes more common, quantifying demographic processes in genetically diverse populations will help predict long-term consequences of plant-plant interactions.

Solanum steroidal glycoalkaloids: structural diversity, biological activities, and biosynthesis

Chemical structures of typical Solanum steroidal glycoalkaloids from eggplant, tomato, and potato.

Genome-wide survey of soybean papain-like cysteine proteases and their expression analysis in root nodule symbiosis

BACKGROUND

Plant papain-like cysteine proteases (PLCPs) are a large class of proteolytic enzymes and play important roles in root nodule symbiosis (RNS), while the whole-genome studies of PLCP family genes in legume are quite limited, and the roles of Glycine max PLCPs (GmPLCPs) in nodulation, nodule development and senescence are not fully understood.

RESULTS

In the present study, we identified 97 GmPLCPs and performed a genome-wide survey to explore the expansion of soybean PLCP family genes and their relationships to RNS. Nineteen paralogous pairs of genomic segments, consisting of 77 GmPLCPs, formed by whole-genome duplication (WGD) events were identified, showing a high degree of complexity in duplication. Phylogenetic analysis among different species showed that the lineage differentiation of GmPLCPs occurred after family expansion, and large tandem repeat segment were specifically in soybean. The expression patterns of GmPLCPs in symbiosis-related tissues and nodules identified RNS-related GmPLCPs and provided insights into their putative symbiotic functions in soybean. The symbiotic function analyses showed that a RNS-related GmPLCP gene (Glyma.04G190700) really participate in nodulation and nodule development.

CONCLUSIONS

Our findings improved our understanding of the functional diversity of legume PLCP family genes, and provided insights into the putative roles of the legume PLCPs in nodulation, nodule development and senescence.

Genome-wide identification of ATP binding cassette (ABC) transporter and heavy metal associated (HMA) gene families in flax (Linum usitatissimum L.)

Background

The recent release of the reference genome sequence assembly of flax, a self-pollinated crop with 15 chromosome pairs, into chromosome-scale pseudomolecules enables the characterization of gene families. The ABC transporter and HMA gene families are important in the control of cadmium (Cd) accumulation in crops. To date, the genome-wide analysis of these two gene families has been successfully conducted in some plant species, but no systematic evolutionary analysis is available for the flax genome.

Results

Here we describe the ABC transporter and HMA gene families in flax to provide a comprehensive overview of its evolution and some support towards the functional annotation of its members. The 198 ABC transporter and 12 HMA genes identified in the flax genome were classified into eight ABC transporter and four HMA subfamilies based on their phylogenetic analysis and domains’ composition. Nine of these genes, i.e., LuABCC9, LuABCC10, LuABCG58, LuABCG59, LuABCG71, LuABCG72, LuABCG73, LuHMA3, and LuHMA4, were orthologous with the Cd associated genes in Arabidopsis, rice and maize. Ten motifs were identified from all ABC transporter and HMA genes. Also, several motifs were conserved among genes of similar length, but each subfamily each had their own motif structures. Both the ABC transporter and HMA gene families were highly conserved among subfamilies of flax and with those of Arabidopsis. While four types of gene duplication were observed at different frequencies, whole-genome or segmental duplications were the most frequent with 162 genes, followed by 29 dispersed, 14 tandem and 4 proximal duplications, suggesting that segmental duplications contributed the most to the expansion of both gene families in flax. The rates of non-synonymous to synonymous (Ka/Ks) mutations of paired duplicated genes were for the most part lower than one, indicative of a predominant purifying selection. Only five pairs of genes clearly exhibited positive selection with a Ka/Ks ratio greater than one. Gene ontology analyses suggested that most flax ABC transporter and HMA genes had a role in ATP binding, transport, catalytic activity, ATPase activity, and metal ion binding. The RNA-Seq analysis of eight different organs demonstrated diversified expression profiling patterns of the genes and revealed their functional or sub-functional conservation and neo-functionalization.

Conclusion

Characterization of the ABC transporter and HMA gene families will help in the functional analysis of candidate genes in flax and other crop species.

Genome-wide identification of the class III POD gene family and their expression profiling in grapevine (Vitis vinifera L)

BACKGROUND

The class III peroxidases (PODs) are involved in a broad range of physiological activities, such as the formation of lignin, cell wall components, defense against pathogenicity or herbivore, and abiotic stress tolerance. The POD family members have been well-studied and characterized by bioinformatics analysis in several plant species, but no previous genome-wide analysis has been carried out of this gene family in grapevine to date.

RESULTS

We comprehensively identified 47 PODs in the grapevine genome and are further classified into 7 subgroups based on their phylogenetic analysis. Results of motif composition and gene structure organization analysis revealed that PODs in the same subgroup shared similar conjunction while the protein sequences were highly conserved. Intriguingly, the integrated analysis of chromosomal mapping and gene collinearity analysis proposed that both dispersed and tandem duplication events contributed to the expansion of PODs in grapevine. Also, the gene duplication analysis suggested that most of the genes (20) were dispersed followed by (15) tandem, (9) segmental or whole-genome duplication, and (3) proximal, respectively. The evolutionary analysis of PODs, such as Ka/Ks ratio of the 15 duplicated gene pairs were less than 1.00, indicated that most of the gene pairs exhibiting purifying selection and 7 pairs underwent positive selection with value greater than 1.00. The Gene Ontology Enrichment (GO), Kyoto Encyclopedia of Genes Genomics (KEGG) analysis, and cis-elements prediction also revealed the positive functions of PODs in plant growth and developmental activities, and response to stress stimuli. Further, based on the publically available RNA-sequence data, the expression patterns of PODs in tissue-specific response during several developmental stages revealed diverged expression patterns. Subsequently, 30 genes were selected for RT-PCR validation in response to (NaCl, drought, and ABA), which showed their critical role in grapevine.

CONCLUSIONS

In conclusion, we predict that these results will lead to novel insights regarding genetic improvement of grapevine.

Identification, Evolution, and Expression Profiling of Histone Lysine Methylation Moderators in Brassica rapa

Histone modifications, such as methylation and demethylation, are vital for regulating chromatin structure, thus affecting its expression patterns. The objective of this study is to understand the phylogenetic relationships, genomic organization, diversification of motif modules, gene duplications, co-regulatory network analysis, and expression dynamics of histone lysine methyltransferases and histone demethylase in Brassica rapa. We identified 60 SET (HKMTases), 53 JmjC, and 4 LSD (HDMases) genes in B. rapa. The domain composition analysis subcategorized them into seven and nine subgroups, respectively. Duplication analysis for paralogous pairs of SET and JmjC (eight and nine pairs, respectively) exhibited variation. Interestingly, three pairs of SET exhibited Ka/Ks > 1.00 values, signifying positive selection, whereas the remaining underwent purifying selection with values less than 1.00. Furthermore, RT-PCR validation analysis and RNA-sequence data acquired on six different tissues (i.e., leaf, stem, callus, silique, flower, and root) revealed dynamic expression patterns. This comprehensive study on the abundance, classification, co-regulatory network analysis, gene duplication, and responses to heat and cold stress of SET and JmjC provides insights into the structure and diversification of these family members in B. rapa. This study will be helpful to reveal functions of these putative SET and JmjC genes in B. rapa.

Combining genome-wide and transcriptome-wide analyses reveal the evolutionary conservation and functional diversity of aquaporins in cotton

BACKGROUND

Aquaporins (AQPs) are integral membrane proteins from a larger family of major intrinsic proteins (MIPs) and function in a huge variety of processes such as water transport, plant growth and stress response. The availability of the whole-genome data of different cotton species allows us to study systematic evolution and function of cotton AQPs on a genome-wide level.

RESULTS

Here, a total of 53, 58, 113 and 111 AQP genes were identified in G. arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. A comprehensive analysis of cotton AQPs, involved in exon/intron structure, functional domains, phylogenetic relationships and gene duplications, divided these AQPs into five subfamilies (PIP, NIP, SIP, TIP and XIP). Comparative genome analysis among 30 species from algae to angiosperm as well as common tandem duplication events in 24 well-studied plants further revealed the evolutionary conservation of AQP family in the organism kingdom. Combining transcriptome analysis and Quantitative Real-time PCR (qRT-PCR) verification, most AQPs exhibited tissue-specific expression patterns both in G. raimondii and G. hirsutum. Meanwhile, a bias of time to peak expression of several AQPs was also detected after treating G. davidsonii and G. hirsutum with 200 mM NaCl. It is interesting that both PIP1;4 h/i/j and PIP2;2a/e showed the highly conserved tandem structure, but differentially contributed to tissue development and stress response in different cotton species.

CONCLUSIONS

These results demonstrated that cotton AQPs were structural conservation while experienced the functional differentiation during the process of evolution and domestication. This study will further broaden our insights into the evolution and functional elucidation of AQP gene family in cotton.

Genome-Wide Identification and Expression Profiling of the Polygalacturonase (PG) and Pectin Methylesterase (PME) Genes in Grapevine (Vitis vinifera L.)

In pectin regulation, polygalacturonases (PGs) and pectin methylesterases (PMEs) are critical components in the transformation, disassembly network, and remodeling of plant primary cell walls. In the current study, we identified 36 PG and 47 PME genes using the available genomic resources of grapevine. Herein, we provide a comprehensive overview of PGs and PMEs, including phylogenetic and collinearity relationships, motif and gene structure compositions, gene duplications, principal component analysis, and expression profiling during developmental stages. Phylogenetic analysis of PGs and PMEs revealed similar domain composition patterns with Arabidopsis. The collinearity analysis showed high conservation and gene duplications with purifying selection. The type of duplications also varied in terms of gene numbers in PGs (10 dispersed, 1 proximal, 12 tandem, and 13 segmental, respectively) and PMEs (23 dispersed, 1 proximal, 16 tandem, and 7 segmental, respectively). The tissue-specific response of PG and PME genes based on the reported transcriptomic data exhibited diverged expression patterns in various organs during different developmental stages. Among PGs, VvPG8, VvPG10, VvPG13, VvPG17, VvPG18, VvPG19, VvPG20, VvPG22, and VvPG23 showed tissue- or organ-specific expression in majority of the tissues during development. Similarly, in PMEs, VvPME3, VvPME4, VvPME5, VvPME6, VvPME19, VvPME21, VvPME23, VvPME29, VvPME31, and VvPME32 suggested high tissue-specific response. The gene ontology (GO), Kyoto Encyclopedia of Genes and Genomics (KEGG) enrichment, and cis-elements prediction analysis also suggested the putative functions of PGs and PMEs in plant development, such as pectin and carbohydrate metabolism, and stress activities. Moreover, qRT-PCR validation of 32 PG and PME genes revealed their role in various organs of grapevines (i.e., root, stem, tendril, inflorescence, flesh, skins, and leaves). Therefore, these findings will lead to novel insights and encourage cutting-edge research on functional characterization of PGs and PMEs in fruit crop species.

Evolution of Glucosinolate Diversity via Whole-Genome Duplications, Gene Rearrangements, and Substrate Promiscuity

Over several decades, glucosinolates have become a model system for the study of specialized metabolic diversity in plants. The near-complete identification of biosynthetic enzymes, regulators, and transporters has provided support for the role of gene duplication and subsequent changes in gene expression, protein function, and substrate specificity as the evolutionary bases of glucosinolate diversity. Here, we provide examples of how whole-genome duplications, gene rearrangements, and substrate promiscuity potentiated the evolution of glucosinolate biosynthetic enzymes, regulators, and transporters by natural selection. This in turn may have led to the repeated evolution of glucosinolate metabolism and diversity in higher plants.

polyRAD: Genotype Calling with Uncertainty from Sequencing Data in Polyploids and Diploids

Low or uneven read depth is a common limitation of genotyping-by-sequencing (GBS) and restriction site-associated DNA sequencing (RAD-seq), resulting in high missing data rates, heterozygotes miscalled as homozygotes, and uncertainty of allele copy number in heterozygous polyploids. Bayesian genotype calling can mitigate these issues, but previously has only been implemented in software that requires a reference genome or uses priors that may be inappropriate for the population. Here we present several novel Bayesian algorithms that estimate genotype posterior probabilities, all of which are implemented in a new R package, polyRAD. Appropriate priors can be specified for mapping populations, populations in Hardy-Weinberg equilibrium, or structured populations, and in each case can be informed by genotypes at linked markers. The polyRAD software imports read depth from several existing pipelines, and outputs continuous or discrete numerical genotypes suitable for analyses such as genome-wide association and genomic prediction.

Extensive allopolyploidy in the neotropical genus Lachemilla (Rosaceae) revealed by PCR-based target enrichment of the nuclear ribosomal DNA cistron and plastid phylogenomics

PREMISE OF THE STUDY

Polyploidy has been long recognized as an important force in plant evolution. Previous studies had suggested widespread occurrence of polyploidy and the allopolyploid origin of several species in the diverse neotropical genus Lachemilla (Rosaceae). Nonetheless, this evidence has relied mostly on patterns of cytonuclear discordance, and direct evidence from nuclear allelic markers is still needed.

METHODS

Here we used PCR target enrichment in combination with high throughput sequencing to obtain multiple copies of the nuclear ribosomal (nr) DNA cistron and 45 regions of the plastid genome (cpDNA) from 219 accessions representing 48 species of Lachemilla and to explore the allopolyploid origin of species in this group.

KEY RESULTS

We were able to identify multiple nrDNA ribotypes and establish clear evidence of allopolyploidy in 33 species of Lachemilla, showing that this condition is common and widespread in the genus. Additionally, we found evidence for three autopolyploid species. We also established multiple, independent origins of several allopolyploid species. Finally, based solely on the cpDNA phylogeny, we identified that the monotypic genus Farinopsis is the sister group of Lachemilla and allied genera within subtribe Fragariinae.

CONCLUSIONS

Our study demonstrates the utility of the nuclear ribosomal DNA cistron to detect allopolyploidy when concerted evolution of this region is not complete. Additionally, with a robust chloroplast phylogeny in place, the direction of hybridization events can be established, and multiple, independent origins of allopolyploid species can be identified.

Transcriptome-based gene expression profiling of diploid radish (Raphanus sativus L.) and the corresponding autotetraploid

Polyploidy is an important evolutionary factor in most land plant lineages which possess more than two complete sets of chromosomes. Radish (Raphanus sativus L.) is an economically annual/biennial root vegetable crop worldwide. However, the expression patterns of duplicated homologs involved in the autopolyploidization remains unclear. In present study, the autotetraploid radish plants (2n = 4x = 36) were produced with colchicine and exhibited an increase in the size of flowers, leaves, stomata and pollen grains. The differential gene expression (DGE) profiling was performed to investigate the differences in gene expression patterns between diploid and its corresponding autotetraploid by RNA-Sequencing (RNA-Seq). Totally, 483 up-regulated differentially expressed genes (DEGs) and 408 down-regulated DEGs were detected in diploid and autotetraploid radishes, which majorly involved in the pathways of hormones, photosynthesis and stress response. Moreover, the xyloglucan endotransglucosylase/hydrolase (XTH) and pectin methylesterases (PME) family members related to cell enlargement and cell wall construction were found to be enriched in GO enrichment analysis, of which XTH family members enriched in "apoplast" and "cell wall" terms, while PME family members enriched in "cell wall" term. Reverse-transcription quantitative PCR (RT-qPCR) analysis indicated that the expression profile of DEGs were consistent with results from the RNA-Seq analysis. The DEGs involved in cell wall construction and auxin metabolism were predicted to be associated with organs size increase of autotetraploid radishes in the present study. These results could provide valuable information for elucidating the molecular mechanism underlying polyploidization and facilitating further genetic improvements of important traits in radish breeding programs.

Allopolyploidization in Cucumis contributes to delayed leaf maturation with repression of redundant homoeologous genes

The important role of polyploidy in plant evolution is widely recognized. However, many questions remain to be explored to address how polyploidy affects the phenotype of the plant. To shed light on the phenotypic and molecular impacts of allopolyploidy, we investigated the leaf development of a synthesized allotetraploid (Cucumis × hytivus), with an emphasis on chlorophyll development. Delayed leaf maturation was identified in C. × hytivus, based on delayed leaf expansion, initial chlorophyll deficiency in the leaves and disordered sink-source transition. Anatomical observations also revealed disturbed chloroplast development in C. ×hytivus. The determination of chlorophyll biosynthesis intermediates suggested that the chlorophyll biosynthesis pathway of C. × hytivus is blocked at the site at which uroporphyrinogen III is catalysed to coproporphyrinogen III. Three chlorophyll biosynthesis-related genes, HEMA1, HEME2 and POR, were significantly repressed in C. × hytivus. Sequence alignment showed both synonymous and non-synonymous substitutions in the HEMA1, HEME2 and POR genes of the parents. Cloning of the chlorophyll biosynthetic genes suggested the retention of homoeologs. In addition, a chimeric clone of the HEMA1 gene that consisted of homologous genes from the parents was identified in C. × hytivus. Overall, our results showed that allopolyploidization in Cucumis has resulted in disturbed chloroplast development and reduced chlorophyll biosynthesis caused by the repressed expression of duplicated homologous genes, which further led to delayed leaf maturation in the allotetraploid, C. × hytivus. The preferential retention/loss of certain types of genes and non-reciprocal homoeologous recombination were also supported in the present study, which provides new insights into the impact of allopolyploidy.

Different Modes of Gene Duplication Show Divergent Evolutionary Patterns and Contribute Differently to the Expansion of Gene Families Involved in Important Fruit Traits in Pear (Pyrus bretschneideri)

Pear is an important fruit crop of the Rosaceae family and has experienced two rounds of ancient whole-genome duplications (WGDs). However, whether different types of gene duplications evolved differently after duplication remains unclear in the pear genome. In this study, we identified the different modes of gene duplication in pear. Duplicate genes derived from WGD, tandem, proximal, retrotransposed, DNA-based transposed or dispersed duplications differ in genomic distribution, gene features, selection pressure, expression divergence, regulatory divergence and biological roles. Widespread sequence, expression and regulatory divergence have occurred between duplicate genes over the 30-45 million years of evolution after the recent genome duplication in pear. The retrotransposed genes show relatively higher expression and regulatory divergence than other gene duplication modes. In contrast, WGD genes underwent a slower sequence divergence and may be influenced by abundant gene conversion events. Moreover, the different classes of duplicate genes exhibited biased functional roles. We also investigated the evolution and expansion patterns of the gene families involved in sugar and organic acid metabolism pathways, which are closely related to the fruit quality and taste in pear. Single-gene duplications largely account for the extensive expansion of gene families involved in the sorbitol metabolism pathway in pear. Gene family expansion was also detected in the sucrose metabolism pathway and tricarboxylic acid cycle pathways. Thus, this study provides insights into the evolutionary fates of duplicated genes.

Structural, expression and evolutionary analysis of the non-specific phospholipase C gene family in Gossypium hirsutum

BACKGROUND

Nonspecific phospholipase C (NPC), which belongs to a phospholipase C subtype, is a class of phospholipases that hydrolyzes the primary membrane phospholipids, such as phosphatidylcholine, to yield sn-1, 2-diacylglycerol and a phosphorylated head-group. NPC plays multiple physiological roles in lipid metabolism and signaling in plants. To fully understand the putative roles of NPC genes in upland cotton, we cloned NPC genes from Gossypium hirsutum and carried out structural, expression and evolutionary analysis.

RESULTS

Eleven NPC genes were cloned from G. hirsutum, which were found on chromosomes scaffold269.1, D03, A07, D07, A08, D11, and scaffold3511_A13. All GhNPCs had typical phosphoesterase domains and have hydrolase activity that acts on ester bonds. GhNPCs were annotated as phospholipase C, which was involved in glycerophospholipid metabolism, ether lipid metabolism, and biosynthesis of secondary metabolites. These GhNPCs showed differential expression patterns in distinct plant tissues and in response to various types of stress (low-phosphate, salt, drought, and abscisic acid). They also had different types and numbers of cis-element. GhNPCs could be classified into four subfamilies. Four pairs of GhNPCs were generated by whole-genome duplication and they underwent purifying selection.

CONCLUSIONS

Our results suggested that GhNPCs are involved in regulating key abiotic stress responses and ABA signaling transduction, and they may have various functional roles for different members under complex abiotic stress conditions. Functional divergence may be the evolutionary driving force for the retention of four pairs of duplicate NPCs. Our analysis provides a solid foundation for the further functional characterization of the GhNPC gene family, and leads to potential applications in the genetic improvement of cotton cultivars.

Evaluating allopolyploid origins in strawberries (Fragaria) using haplotypes generated from target capture sequencing

BACKGROUND

Hybridization is observed in many eukaryotic lineages and can lead to the formation of polyploid species. The study of hybridization and polyploidization faces challenges both in data generation and in accounting for population-level phenomena such as coalescence processes in phylogenetic analysis. Genus Fragaria is one example of a set of plant taxa in which a range of ploidy levels is observed across species, but phylogenetic origins are unknown.

RESULTS

Here, using 20 diploid and polyploid Fragaria species, we combine approaches from NGS data analysis and phylogenetics to infer evolutionary origins of polyploid strawberries, taking into account coalescence processes. We generate haplotype sequences for 257 low-copy nuclear markers assembled from Illumina target capture sequence data. We then identify putative hybridization events by analyzing gene tree topologies, and further test predicted hybridizations in a coalescence framework. This approach confirms the allopolyploid ancestry of F. chiloensis and F. virginiana, and provides new allopolyploid ancestry hypotheses for F. iturupensis, F. moschata, and F. orientalis. Evidence of gene flow between diploids F. bucharica and F. vesca is also detected, suggesting that it might be appropriate to consider these groups as conspecifics.

CONCLUSIONS

This study is one of the first in which target capture sequencing followed by computational deconvolution of individual haplotypes is used for tracing origins of polyploid taxa. The study also provides new perspectives on the evolutionary history of Fragaria.

Evolutionary Dynamics of Unreduced Gametes

Unreduced gametes, which have the somatic (2n) chromosome number, are an important precursor to polyploid formation and apomixis. The product of irregularities in meiosis, 2n gametes are expected to be rare and deleterious in most natural populations, contrary to their wide taxonomic distribution and the prevalence of polyploidy. To better understand this discrepancy, we review contemporary evidence related to four aspects of 2n gamete dynamics in natural populations: (i) estimates of their frequency; (ii) their environmental and genetic determinants; (iii) adaptive and nonadaptive processes regulating their evolution; and (iv) factors regulating their union and production of polyploids in diploid populations. Aided by high-throughput methods of detection, these foci will advance our understanding of variation in 2n gametes within and among species, and their role in polyploid evolution.

Divergences in hydraulic architecture form an important basis for niche differentiation between diploid and polyploid Betula species in NE China

Habitat differentiation between polyploid and diploid plants are frequently observed, with polyploids usually occupying more stressed environments. In woody plants, polyploidization can greatly affect wood characteristics but knowledge of its influences on xylem hydraulics is scarce. The four Betula species in NE China, representing two diploids and two polyploids with obvious habitat differentiation, provide an exceptional study system for investigating the impact of polyploidization on environmental adaptation of trees from the point view of xylem hydraulics. To test the hypothesis that changes in hydraulic architecture play an important role in determining their niche differentiation, we measured wood structural traits at both the tissue and pit levels and quantified xylem water transport efficiency and safety in these species. The two polyploids had significantly larger hydraulic weighted mean vessel diameters than the two diploids (45.1 and 45.5 vs 25.9 and 24.5 μm) although the polyploids are occupying more stressed environments. As indicated by more negative water potentials corresponding to 50% loss of stem hydraulic conductivities, the two polyploids exhibited significantly higher resistance to drought-induced embolism than the two diploids (-5.23 and -5.05 vs -3.86 and -3.13 MPa) despite their larger vessel diameters. This seeming discrepancy is reconciled by distinct characteristics favoring greater embolism resistance at the pit level in the two polyploid species. Our results showed clearly that the two polyploid species have remarkably different pit-level anatomical traits favoring greater hydraulic safety than their congeneric diploid species, which have likely contributed to the abundance of polyploid birches in more stressed habitats; however, less porous inter-conduit pits together with a reduced leaf to sapwood area may have compromised their competitiveness under more favorable conditions. Contrasts in hydraulic architecture between diploid and polyploid Betula species suggest an important functional basis for their clear habitat differentiation along environmental gradients in Changbai Mountain of NE China.

Heme oxygenase 1 defects lead to reduced chlorophyll in Brassica napus

We previously described a Brassica napus chlorophyll-deficient mutant (ygl) with yellow-green seedling leaves and mapped the related gene, BnaC.YGL, to a 0.35 cM region. However, the molecular mechanisms involved in this chlorophyll defect are still unknown. In this study, the BnaC07.HO1 gene (equivalent to BnaC.YGL) was isolated by the candidate gene approach, and its function was confirmed by genetic complementation. Comparative sequencing analysis suggested that BnaC07.HO1 was lost in the mutant, while a long noncoding-RNA was inserted into the promoter of the homologous gene BnaA07.HO1. This insert was widely present in B. napus cultivars and down-regulated BnaA07.HO1 expression. BnaC07.HO1 was highly expressed in the seedling leaves and encoded heme oxygenase 1, which was localized in the chloroplast. Biochemical analysis showed that BnaC07.HO1 can catalyze heme conversion to form biliverdin IXα. RNA-seq analysis revealed that the loss of BnaC07.HO1 impaired tetrapyrrole metabolism, especially chlorophyll biosynthesis. According, the levels of chlorophyll intermediates were reduced in the ygl mutant. In addition, gene expression in multiple pathways was affected in ygl. These findings provide molecular evidences for the basis of the yellow-green leaf phenotype and further insights into the crucial role of HO1 in B. napus.

Divergent evolutionary patterns of the MAPK cascade genes in Brassica rapa and plant phylogenetics

Mitogen-activated protein kinase (MAPK) cascade signal transduction modules play crucial roles in regulating many biological processes in plants. These cascades are composed of three classes of hierarchically organized protein kinases, MAPKKKs, MAPKKs and MAPKs. Here, we analyzed gene retention, phylogenetic, evolution and expression patterns of MAPK cascade genes in Brassica rapa. We further found that the MAPK branches, classes III and IV, appeared after the split of bryophytes and green algae after analyzing the MAPK cascade genes in 8 species, and their rapid expansion led to the great size of the families of MAPKs. In contrast, the ancestral class I subfamily of MAPKK gene families have been highly conserved from algae to angiosperms. For the MAPKKK family, the MEKK and Raf subfamily share a common evolutionary origin, and Raf plays a major role in the expansion of the MAPKKK gene family. The cis-elements and interaction network analyses showed the important function of MAPK cascade genes in development and stress responses in B. rapa. This study provides a solid foundation for molecular evolution analyses of MAPK cascade genes.

A Simple and Efficient Approach to Elucidate Genomic Contribution of Transcripts to a Target Gene in Polyploids: The Case of Hexaploid Wheat (Triticum aestivum L.)

Common wheat (Triticum aestivum L.) is one of the most economically important crops in the world, however, gene functional studies in this crop have been lagging mainly due to the complexity of its polyploid genome, which is derived through two rounds of intergeneric hybridization events that led to the presence of six copies for each gene. Elucidating the transcript contribution of each genome to the total expression of a target gene in polyploids such as hexaploid wheat has a paramount significance for direct discovery of genes and the associated molecular mechanisms controlling traits of agronomic importance. A polymerase chain reaction approach that involved primers amplifying DNA fragments unique to each homeolog of a target gene and quantitation of the intensity of the resulting fragment bands were able to successfully determine the genomic transcript contributions as a percentage of target gene's total expression in hexaploid wheat. Our results showed that the genomic contributions of transcripts to a target gene vary with genotype and tissue type, suggesting the distinct role of each homeolog in regulating the trait associated with the target gene. The approach described in this study is an effective and economical method to elucidate the genomic transcript contribution to the total expression of individual target genes in hexaploid wheat. It can also be applied to determine the transcript contribution of each genome towards the collective expression of a target gene in other economically important polypoid crop species.

A Genomic Analysis of Factors Driving lincRNA Diversification: Lessons from Plants

Transcriptomic analyses from across eukaryotes indicate that most of the genome is transcribed at some point in the developmental trajectory of an organism. One class of these transcripts is termed long intergenic noncoding RNAs (lincRNAs). Recently, attention has focused on understanding the evolutionary dynamics of lincRNAs, particularly their conservation within genomes. Here, we take a comparative genomic and phylogenetic approach to uncover factors influencing lincRNA emergence and persistence in the plant family Brassicaceae, to which Arabidopsis thaliana belongs. We searched 10 genomes across the family for evidence of > 5000 lincRNA loci from A. thaliana. From loci conserved in the genomes of multiple species, we built alignments and inferred phylogeny. We then used gene tree/species tree reconciliation to examine the duplication history and timing of emergence of these loci. Emergence of lincRNA loci appears to be linked to local duplication events, but, surprisingly, not whole genome duplication events (WGD), or transposable elements. Interestingly, WGD events are associated with the loss of loci for species having undergone relatively recent polyploidy. Lastly, we identify 1180 loci of the 6480 previously annotated A. thaliana lincRNAs (18%) with elevated levels of conservation. These conserved lincRNAs show higher expression, and are enriched for stress-responsiveness and cis-regulatory motifs known as conserved noncoding sequences (CNSs). These data highlight potential functional pathways and suggest that CNSs may regulate neighboring genes at both the genomic and transcriptomic level. In sum, we provide insight into processes that may influence lincRNA diversification by providing an evolutionary context for previously annotated lincRNAs.

Homoeologs: What Are They and How Do We Infer Them?

The evolutionary history of nearly all flowering plants includes a polyploidization event. Homologous genes resulting from allopolyploidy are commonly referred to as 'homoeologs', although this term has not always been used precisely or consistently in the literature. With several allopolyploid genome sequencing projects under way, there is a pressing need for computational methods for homoeology inference. Here we review the definition of homoeology in historical and modern contexts and propose a precise and testable definition highlighting the connection between homoeologs and orthologs. In the second part, we survey experimental and computational methods of homoeolog inference, considering the strengths and limitations of each approach. Establishing a precise and evolutionarily meaningful definition of homoeology is essential for understanding the evolutionary consequences of polyploidization.

Vegetative and generative maintenance of self-incompatibility in six accessions of German chamomile

Self-incompatible (SI) plants are able to form ideal mother lines for hybrid crossing in hermaphroditic plants, assuring fertilization from the desired father line. To find out suitable ways to maintain SI was the aim of this study. Among 220 plants of German chamomile (Matricaria recutita (L.) Rauschert) within six accessions SI-genotypes were selected. SI was determined as staying seedless in three flower heads per plant. Initial SI-plants formed the basic paternal generation (P1) of i) maintaining the same genotypes over six months and repeating seed set analysis (P2) and of ii) conducting crossings in three versions (SI × SI, SI × NSI (not SI evaluated plants) and NSI × SI), thereby producing the F1 population. F1 exhibited 78% SI and P2 62% SI, indicating a higher environmental than genetic influence on SI. But heritability, calculated from the results of SI × SI crossings, showed high values (h(2) = 0.71). Within generative propagation, the influence of generation/crossing version was highly significant (p = 0.001) and the cultivar 'Degumille' explored the highest value of SI (86%) after SI × NSI crossings. Therefore, the intra-cultivar combination of 'Degumille' SI mother plants crossed with NSI father plants can be recommended as the most promising version to maintain SI in chamomile.

Morphological features of different polyploids for adaptation and molecular characterization of CC-NBS-LRR and LEA gene families in Agave L

Polyploidy has been widely described in many Agave L. species, but its influence on environmental response to stress is still unknown. With the objective of knowing the morphological adaptations and regulation responses of genes related to biotic (LEA) and abiotic (NBS-LRR) stress in species of Agave with different levels of ploidy, and how these factors contribute to major response of Agave against environmental stresses, we analyzed 16 morphological trials on five accessions of three species (Agave tequilana Weber, Agave angustifolia Haw. and Agave fourcroydes Lem.) with different ploidy levels (2n=2x=60 2n=3x=90, 2n=5x=150, 2n=6x=180) and evaluated the expression of NBS-LRR and LEA genes regulated by biotic and abiotic stress. It was possible to associate some morphological traits (spines, nuclei, and stomata) to ploidy level. The genetic characterization of stress-related genes NBS-LRR induced by pathogenic infection and LEA by heat or saline stresses indicated that amino acid sequence analysis in these genes showed more substitutions in higher ploidy level accessions of A. fourcroydes Lem. 'Sac Ki' (2n=5x=150) and A. angustifolia Haw. 'Chelem Ki' (2n=6x=180), and a higher LEA and NBS-LRR representativeness when compared to their diploid and triploid counterparts. In all studied Agave accessions expression of LEA and NBS-LRR genes was induced by saline or heat stresses or by infection with Erwinia carotovora, respectively. The transcriptional activation was also higher in A. angustifolia Haw. 'Chelem Ki' (2n=6x=180) and A. fourcroydes 'Sac Ki' (2n=5x=150) than in their diploid and triploid counterparts, which suggests higher adaptation to stress. Finally, the diploid accession A. tequilana Weber 'Azul' showed a differentiated genetic profile relative to other Agave accessions. The differences include similar or higher genetic representativeness and transcript accumulation of LEA and NBS-LRR genes than in polyploid (2n=5x=150 and 2n=6x=180) Agave accessions, thus suggesting a differentiated selection pressure for overcoming the lower ploidy level of the diploid A. tequilana Weber 'Azul'.

Comprehensive analysis of the polygalacturonase and pectin methylesterase genes in Brassica rapa shed light on their different evolutionary patterns

Pectins are fundamental polysaccharides in the plant primary cell wall. Polygalacturonases (PGs) and pectin methylesterases (PMEs), major components of the pectin remodeling and disassembly network, are involved in cell separation processes during many stages of plant development. A comprehensive study of these genes in plants could shed light on the evolution patterns of their structural development. In this study, we conducted whole-genome annotation, molecular evolution and gene expression analyses of PGs and PMEs in Brassica rapa and 8 other plant species. A total of 100 PGs and 110 PMEs were identified in B. rapa; they primarily diverged from 12–18 MYA and PMEs were retained more than PGs. Along with another 305 PGs and 348 PMEs in the 8 species, two different expansion or evolution types were discovered: a new branch of class A PGs appeared after the split of gymnosperms and angiosperms, which led to the rapid expansion of PGs; the pro domain was obtained or lost in the proPMEs through comprehensive analyses among PME genes. In addition, the PGs and PMEs exhibit diverged expression patterns. These findings will lead to novel insight regarding functional divergence and conservation in the gene families and provide more support for molecular evolution analyses.

Evolution of DNA-Binding Sites of a Floral Master Regulatory Transcription Factor

Flower development is controlled by the action of key regulatory transcription factors of the MADS-domain family. The function of these factors appears to be highly conserved among species based on mutant phenotypes. However, the conservation of their downstream processes is much less well understood, mostly because the evolutionary turnover and variation of their DNA-binding sites (BSs) among plant species have not yet been experimentally determined. Here, we performed comparative ChIP (chromatin immunoprecipitation)-seq experiments of the MADS-domain transcription factor SEPALLATA3 (SEP3) in two closely related Arabidopsis species: Arabidopsis thaliana and A. lyrata which have very similar floral organ morphology. We found that BS conservation is associated with DNA sequence conservation, the presence of the CArG-box BS motif and on the relative position of the BS to its potential target gene. Differences in genome size and structure can explain that SEP3 BSs in A. lyrata can be located more distantly to their potential target genes than their counterparts in A. thaliana. In A. lyrata, we identified transposition as a mechanism to generate novel SEP3 binding locations in the genome. Comparative gene expression analysis shows that the loss/gain of BSs is associated with a change in gene expression. In summary, this study investigates the evolutionary dynamics of DNA BSs of a floral key-regulatory transcription factor and explores factors affecting this phenomenon.

Natural variation in ARF18 gene simultaneously affects seed weight and silique length in polyploid rapeseed

Seed weight (SW), which is one of the three major factors influencing grain yield, has been widely accepted as a complex trait that is controlled by polygenes, particularly in polyploid crops. Brassica napus L., which is the second leading crop source for vegetable oil around the world, is a tetraploid (4×) species. In the present study, we identified a major quantitative trait locus (QTL) on chromosome A9 of rapeseed in which the genes for SW and silique length (SL) were colocated. By fine mapping and association analysis, we uncovered a 165-bp deletion in the auxin-response factor 18 (ARF18) gene associated with increased SW and SL. ARF18 encodes an auxin-response factor and shows inhibitory activity on downstream auxin genes. This 55-aa deletion prevents ARF18 from forming homodimers, in turn resulting in the loss of binding activity. Furthermore, reciprocal crossing has shown that this QTL affects SW by maternal effects. Transcription analysis has shown that ARF18 regulates cell growth in the silique wall by acting via an auxin-response pathway. Together, our results suggest that ARF18 regulates silique wall development and determines SW via maternal regulation. In addition, our study reveals the first (to our knowledge) QTL in rapeseed and may provide insights into gene cloning involving polyploid crops.

Transcriptional networks-crops, clocks, and abiotic stress

Several factors affect the yield potential and geographical range of crops including the circadian clock, water availability, and seasonal temperature changes. In order to sustain and increase plant productivity on marginal land in the face of both biotic and abiotic stresses, we need to more efficiently generate stress-resistant crops through marker-assisted breeding, genetic modification, and new genome-editing technologies. To leverage these strategies for producing the next generation of crops, future transcriptomic data acquisition should be pursued with an appropriate temporal design and analyzed with a network-centric approach. The following review focuses on recent developments in abiotic stress transcriptional networks in economically important crops and will highlight the utility of correlation-based network analysis and applications.