Ancient duplication of cereal genomes

Transcriptomic and QTL Analysis of Seed Germination Vigor under Low Temperature in Weedy Rice WR04-6

Low temperature is one of the major factors affecting rice germination, and low temperature germination (LTG) is an important agronomic trait. Although significant progress has been made in the study of rice LTG, the molecular mechanism of LTG remains poorly understood. To explore more rice LTG gene resources, we first demonstrated that weedy rice WR04-6 (Oryza sativa f. spontanea) had significantly higher LTG ability at 10 °C than the cultivated rice Qishanzhan (QSZ Oryza sativa L. ssp. indica). RNA-seq was used to investigate the gene expression of WR04-6 and QSZ at 10 °C for 10, 12 and 14 days after imbibition (DAI) of seed germination. The results of Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment revealed that the differentially expressed genes (DEGs) between WR04-6 and QSZ were mainly concentrated on the response to starch catabolic processes and the response to abscisic acid (ABA). This is consistent with the results of α-amylase activity, ABA and gibberellins (GA) treatment. A recombinant inbred line (RIL) population derived from a cross between WR04-6 and QSZ and its high-density SNP genetic map were used to detect quantitative trait loci (QTL) for LTG rates. The results showed that two new QTLs were located on chromosome 3 and chromosome 12. Combined with the mapped QTLs and RNA-seq DEGs, sixteen candidate genes potentially associated with LTG were identified. Validation of the expression of the candidates by qRT-PCR were consistent with the RNA-seq data. These results will enable us to understand the genetic basis of LTG in weedy rice and provide new genetic resources for the generation of rice germplasm with improved LTG.

Differences between diploid donors are the main contributing factor for subgenome asymmetry measured in either gene ratio or relative diversity in allopolyploids

Subgenome asymmetry (SA) has routinely been attributed to different responses between the subgenomes of a polyploid to various stimuli during evolution. Here, we compared subgenome differences in gene ratio and relative diversity between artificial and natural genotypes of several allopolyploid species. Surprisingly, consistent differences were not detected between these two types of polyploid genotypes, although they differ in times exposed to evolutionary selection. The estimated ratio of shared genes between a subgenome and its diploid donor was invariably higher for the artificial allopolyploid genotypes than those for the natural genotypes, which is expected as it is now well-known that many genes in a species are not shared among all individuals. As the exact diploid parent for a given subgenome is unknown, the estimated ratios of shared genes for the natural genotypes would also include difference among individual genotypes of the diploid donor species. Further, we detected the presence of SA in genotypes before the completion of the polyploidization events as well as in those which were not formed via polyploidization. These results indicate that SA may, to a large degree, reflect differences between its diploid donors or that changes occurred during polyploid evolution are defined by their donor genomes.

Paleo-polyploidization in Lycophytes

Lycophytes and seed plants constitute the typical vascular plants. Lycophytes have been thought to have no paleo-polyploidization although the event is known to be critical for the fast expansion of seed plants. Here, genomic analyses including the homologous gene dot plot analysis detected multiple paleo-polyploidization events, with one occurring approximately 13-15 million years ago (MYA) and another about 125-142 MYA, during the evolution of the genome of Selaginella moellendorffii, a model lycophyte. In addition, comparative analysis of reconstructed ancestral genomes of lycophytes and angiosperms suggested that lycophytes were affected by more paleo-polyploidization events than seed plants. Results from the present genomic analyses indicate that paleo-polyploidization has contributed to the successful establishment of both lineages-lycophytes and seed plants-of vascular plants.

Sequential Paleotetraploidization shaped the carrot genome

BACKGROUND

Carrot (Daucus carota subsp. carota L.) is an important root crop with an available high-quality genome. The carrot genome is thought to have undergone recursive paleo-polyploidization, but the extent, occurrences, and nature of these events are not clearly defined.

RESULTS

Using a previously published comparative genomics pipeline, we reanalysed the carrot genome and characterized genomic fractionation, as well as gene loss and retention, after each of the two tetraploidization events and inferred a dominant and sensitive subgenome for each event. In particular, we found strong evidence of two sequential tetraploidization events, with one (Dc-α) approximately 46-52 million years ago (Mya) and the other (Dc-β) approximately 77-87 Mya, both likely allotetraploidization in nature. The Dc-β event was likely common to all Apiales plants, occurring around the divergence of Apiales-Bruniales and after the divergence of Apiales-Asterales, likely playing an important role in the derivation and divergence of Apiales species. Furthermore, we found that rounds of polyploidy events contributed to the expansion of gene families responsible for plastidial methylerythritol phosphate (MEP), the precursor of carotenoid accumulation, and shaped underlying regulatory pathways. The alignment of orthologous and paralogous genes related to different events of polyploidization and speciation constitutes a comparative genomics platform for studying Apiales, Asterales, and many other related species.

CONCLUSIONS

Hierarchical inference of homology revealed two tetraploidization events that shaped the carrot genome, which likely contributed to the successful establishment of Apiales plants and the expansion of MEP, upstream of the carotenoid accumulation pathway.

QTLs and analysis of the candidate gene for amylose, protein, and moisture content in rice (Oryza sativa L.)

In this study, we determined using NIRS the heritability percentage of amylose, protein, and moisture content in polished and unpolished rice in a CNDH population derived from a cross between Cheongcheong and Nagdong rice varieties. The results revealed a higher heritability percentage for the amylose content and compromised heritability for protein and moisture contents. We also conducted QTL analysis of rice for these major components and identified their chromosomal locations on a physical map. We found a total of four QTLs affecting the amylose, protein, and moisture contents of grain on chromosome 7. We constructed physical maps of seven DNA markers responsible for amylose content, six responsible for protein content, and three responsible for moisture content. Furthermore, we classified these genes according to their functions and found 17 genes (over 77%) to be involved in secondary metabolite synthesis, two genes (about 9%), each related to cell function and abiotic stress, and one gene (about 5%) involved in redox signaling.

Genome Alignment Spanning Major Poaceae Lineages Reveals Heterogeneous Evolutionary Rates and Alters Inferred Dates for Key Evolutionary Events

Multiple comparisons among genomes can clarify their evolution, speciation, and functional innovations. To date, the genome sequences of eight grasses representing the most economically important Poaceae (grass) clades have been published, and their genomic-level comparison is an essential foundation for evolutionary, functional, and translational research. Using a formal and conservative approach, we aligned these genomes. Direct comparison of paralogous gene pairs all duplicated simultaneously reveal striking variation in evolutionary rates among whole genomes, with nucleotide substitution slowest in rice and up to 48% faster in other grasses, adding a new dimension to the value of rice as a grass model. We reconstructed ancestral genome contents for major evolutionary nodes, potentially contributing to understanding the divergence and speciation of grasses. Recent fossil evidence suggests revisions of the estimated dates of key evolutionary events, implying that the pan-grass polyploidization occurred ∼96 million years ago and could not be related to the Cretaceous-Tertiary mass extinction as previously inferred. Adjusted dating to reflect both updated fossil evidence and lineage-specific evolutionary rates suggested that maize subgenome divergence and maize-sorghum divergence were virtually simultaneous, a coincidence that would be explained if polyploidization directly contributed to speciation. This work lays a solid foundation for Poaceae translational genomics.

PolyCat: a resource for genome categorization of sequencing reads from allopolyploid organisms

Read mapping is a fundamental part of next-generation genomic research but is complicated by genome duplication in many plants. Categorizing DNA sequence reads into their respective genomes enables current methods to analyze polyploid genomes as if they were diploid. We present PolyCat-a pipeline for mapping and categorizing all types of next-generation sequence data produced from allopolyploid organisms. PolyCat uses GSNAP's single-nucleotide polymorphism (SNP)-tolerant mapping to minimize the mapping efficiency bias caused by SNPs between genomes. PolyCat then uses SNPs between genomes to categorize reads according to their respective genomes. Bisulfite-treated reads have a significant reduction in nucleotide complexity because nucleotide conversion events are confounded with transition substitutions. PolyCat includes special provisions to properly handle bisulfite-treated data. We demonstrate the functionality of PolyCat on allotetraploid cotton, Gossypium hirsutum, and create a functional SNP index for efficiently mapping sequence reads to the D-genome sequence of G. raimondii. PolyCat is appropriate for all allopolyploids and all types of next-generation genome analysis, including differential expression (RNA sequencing), differential methylation (bisulfite sequencing), differential DNA-protein binding (chromatin immunoprecipitation sequencing), and population diversity.

Molecular aspects of flower development in grasses

The grass family (Poaceae) of the monocotyledons includes about 10,000 species and represents one of the most important taxa among angiosperms. Their flower morphology is remarkably different from those of other monocotyledons and higher eudicots. The peculiar floral structure of grasses is the floret, which contains carpels and stamens, like eudicots, but lacks petals and sepals. The reproductive organs are surrounded by two lodicules, which correspond to eudicot petals, and by a palea and lemma, whose correspondence to eudicot organs remains controversial. The molecular and genetic analysis of floral morphogenesis and organ specification, primarily performed in eudicot model species, led to the ABCDE model of flower development. Several genes required for floral development in grasses correspond to class A, B, C, D, and E genes of eudicots, but others appear to have unique and diversified functions. In this paper, we outline the present knowledge on the evolution and diversification of grass genes encoding MIKC-type MADS-box transcription factors, based on information derived from studies in rice, maize, and wheat. Moreover, we review recent advances in studying the genes involved in the control of flower development and the extent of structural and functional conservation of these genes between grasses and eudicots.

Gene conversion in angiosperm genomes with an emphasis on genes duplicated by polyploidization

Angiosperm genomes differ from those of mammals by extensive and recursive polyploidizations. The resulting gene duplication provides opportunities both for genetic innovation, and for concerted evolution. Though most genes may escape conversion by their homologs, concerted evolution of duplicated genes can last for millions of years or longer after their origin. Indeed, paralogous genes on two rice chromosomes duplicated an estimated 60–70 million years ago have experienced gene conversion in the past 400,000 years. Gene conversion preserves similarity of paralogous genes, but appears to accelerate their divergence from orthologous genes in other species. The mutagenic nature of recombination coupled with the buffering effect provided by gene redundancy, may facilitate the evolution of novel alleles that confer functional innovations while insulating biological fitness of affected plants. A mixed evolutionary model, characterized by a primary birth-and-death process and occasional homoeologous recombination and gene conversion, may best explain the evolution of multigene families.

The flowering world: a tale of duplications

Flowering plants contain many genes, most of which were created during the past 200 or so million years through small- and large-scale duplications. Paleo-polyploidy events, in particular, have been the subject of much recent research. There is a growing consensus that one or more genome doubling or merging events occurred early during the evolution of the flowering plants, and that many lineages have since undergone additional, independent and more recent duplication events. Here, we review the difficulties in determining the number of genome duplications and discuss how the completion of some additional genome sequences of species occupying key phylogenetic positions has led to a better understanding of the timing of certain duplication events. This is important if we want to demonstrate the significance of genome duplications for the evolution and radiation of (different groups of) flowering plants.

Expression pattern divergence of duplicated genes in rice

Background

Genome-wide duplication is ubiquitous during diversification of the angiosperms, and gene duplication is one of the most important mechanisms for evolutionary novelties. As an indicator of functional evolution, the divergence of expression patterns following duplication events has drawn great attention in recent years. Using large-scale whole-genome microarray data, we systematically analyzed expression divergence patterns of rice genes from block, tandem and dispersed duplications.

Results

We found a significant difference in expression divergence patterns for the three types of duplicated gene pairs. Expression correlation is significantly higher for gene pairs from block and tandem duplications than those from dispersed duplications. Furthermore, a significant correlation was observed between the expression divergence and the synonymous substitution rate which is an approximate proxy of divergence time. Thus, both duplication types and divergence time influence the difference in expression divergence. Using a linear model, we investigated the influence of these two variables and found that the difference in expression divergence between block and dispersed duplicates is attributed largely to their different divergence time. In addition, the difference in expression divergence between tandem and the other two types of duplicates is attributed to both divergence time and duplication type.

Conclusion

Consistent with previous studies on Arabidopsis, our results revealed a significant difference in expression divergence between the types of duplicated genes and a significant correlation between expression divergence and synonymous substitution rate. We found that the attribution of duplication mode to the expression divergence implies a different evolutionary course of duplicated genes.

Phylogeny and evolution of Rab7 and Rab9 proteins

Background

An important role in the evolution of intracellular trafficking machinery in eukaryotes played small GTPases belonging to the Rab family known as pivotal regulators of vesicle docking, fusion and transport. The Rab family is very diversified and divided into several specialized subfamilies. We focused on the VII functional group comprising Rab7 and Rab9, two related subfamilies, and analysed 210 sequences of these proteins. Rab7 regulates traffic from early to late endosomes and from late endosome to vacuole/lysosome, whereas Rab9 participates in transport from late endosomes to the trans-Golgi network.

Results

Although Rab7 and Rab9 proteins are quite small and show heterogeneous rates of substitution in different lineages, we found a phylogenetic signal and inferred evolutionary relationships between them. Rab7 proteins evolved before radiation of main eukaryotic supergroups while Rab9 GTPases diverged from Rab7 before split of choanoflagellates and metazoans. Additional duplication of Rab9 and Rab7 proteins resulting in several isoforms occurred in the early evolution of vertebrates and next in teleost fishes and tetrapods. Three Rab7 lineages emerged before divergence of monocots and eudicots and subsequent duplications of Rab7 genes occurred in particular angiosperm clades. Interestingly, several Rab7 copies were identified in some representatives of excavates, ciliates and amoebozoans. The presence of many Rab copies is correlated with significant differences in their expression level. The diversification of analysed Rab subfamilies is also manifested by non-conserved sequences and structural features, many of which are involved in the interaction with regulators and effectors. Individual sites discriminating different subgroups of Rab7 and Rab9 GTPases have been identified.

Tracking the evolutionary history of polyploidy in Fragaria L. (strawberry): new insights from phylogenetic analyses of low-copy nuclear genes

Phylogenetic utility of two nuclear genes (GBSSI-2 and DHAR) was explored in genus Fragaria in order to clarify phylogenetic relationships among taxa and to elucidate the origin of the polyploid species. Orthology of the amplified products was assessed by several methods. Our results strongly suggest the loss of one GBSSI duplicated copy (GBSSI-1) in the Fragariinae subtribe. Phylogenetic analyses provided new insights into the evolutionary history of Fragaria, such as evidence supporting the presence of three main diploid genomic pools in the genus and demonstrating the occurrence of independent events of polyploidisation. In addition, the data provide evidence supporting an allopolyploid origin of the hexaploid F. moschata, and the octoploids F. chiloensis, F. iturupensis and F. virginiana. Accordingly, a new pattern summarizing our present knowledge on the Fragaria evolutionary history is proposed. Additionally, sequence analyses also revealed relaxed constraints on homoeologous copies at high ploidy level, as demonstrated by deletion events within DHAR coding sequences of some allo-octoploid haplotypes.

How to usefully compare homologous plant genes and chromosomes as DNA sequences

There are four sequenced and publicly available plant genomes to date. With many more slated for completion, one challenge will be to use comparative genomic methods to detect novel evolutionary patterns in plant genomes. This research requires sequence alignment algorithms to detect regions of similarity within and among genomes. However, different alignment algorithms are optimized for identifying different types of homologous sequences. This review focuses on plant genome evolution and provides a tutorial for using several sequence alignment algorithms and visualization tools to detect useful patterns of conservation: conserved non-coding sequences, false positive noise, subfunctionalization, synteny, annotation errors, inversions and local duplications. Our tutorial encourages the reader to experiment online with the reviewed tools as a companion to the text.

Molecular evolution study in China: progress and future promise

China has a large land area with highly diverse topography, climate and vegetation, and animal resources and is ranked eighth in the world and first in the Northern Hemisphere on richness of biodiversity. Even though little work on molecular evolution had been reported a decade ago, studies on both the evolution of macromolecules and the molecular phylogeny have become active in China in recent years. This review highlights some of the interesting and important developments in molecular evolution study in China. Chinese scientists have made significant contribution on the methods inferring phylogeny and biogeography of animals and plants in East Asia using molecular data. Studies on population and conservation genetics of animals and plants, such as Golden monkey and Chinese sturgeon, provided useful information for conserving the endangered species. East and South Asia has been demonstrated to be one of the centres of domestication. Origin and evolution of genes and gene families have been explored, which shed new insight on the genetic mechanism of adaptation. In the genomic era, Chinese researchers also made a transition from single-gene to a genomic investigation approach. Considering the fact that amazing progress has been made in the past few years, and more and more talented young scientists are entering field, the future of molecular evolution study in China holds much promise.

Biochemical Characterization of Plant Ethylene Receptors Following Transgenic Expression in Yeast

The ethylene receptors of plants are related to and originated from bacterial histidine kinases. As such they represent a system by which one can study not only how the ethylene signal is perceived and its signal transduced, but also how bacterial two-component systems have been adapted for signal transduction in a eukaryote. Much of the biochemical characterization of the ethylene receptors, including the demonstration of kinase activity, ethylene binding, and interaction with other signaling components, has relied on the ability of the receptors to be functionally expressed in transgenic yeast. This chapter describes some of the key approaches used for such work, with a special emphasis on techniques employed to analyze ethylene binding. In many cases the approaches used in transgenic yeast may also be used for studies of the receptors in the native plant.

Statistical inference of chromosomal homology based on gene colinearity and applications to Arabidopsis and rice

Background

The identification of chromosomal homology will shed light on such mysteries of genome evolution as DNA duplication, rearrangement and loss. Several approaches have been developed to detect chromosomal homology based on gene synteny or colinearity. However, the previously reported implementations lack statistical inferences which are essential to reveal actual homologies.

Results

In this study, we present a statistical approach to detect homologous chromosomal segments based on gene colinearity. We implement this approach in a software package ColinearScan to detect putative colinear regions using a dynamic programming algorithm. Statistical models are proposed to estimate proper parameter values and evaluate the significance of putative homologous regions. Statistical inference, high computational efficiency and flexibility of input data type are three key features of our approach.

Conclusion

We apply ColinearScan to the Arabidopsis and rice genomes to detect duplicated regions within each species and homologous fragments between these two species. We find many more homologous chromosomal segments in the rice genome than previously reported. We also find many small colinear segments between rice and Arabidopsis genomes.

Gene-balanced duplications, like tetraploidy, provide predictable drive to increase morphological complexity

Controversy surrounds the apparent rising maximums of morphological complexity during eukaryotic evolution, with organisms increasing the number and nestedness of developmental areas as evidenced by morphological elaborations reflecting area boundaries. No "predictable drive" to increase this sort of complexity has been reported. Recent genetic data and theory in the general area of gene dosage effects has engendered a robust "gene balance hypothesis," with a theoretical base that makes specific predictions as to gene content changes following different types of gene duplication. Genomic data from both chordate and angiosperm genomes fit these predictions: Each type of duplication provides a one-way injection of a biased set of genes into the gene pool. Tetraploidies and balanced segments inject bias for those genes whose products are the subunits of the most complex biological machines or cascades, like transcription factors (TFs) and proteasome core proteins. Most duplicate genes are removed after tetraploidy. Genic balance is maintained by not removing those genes that are dose-sensitive, which tends to leave duplicate "functional modules" as the indirect products (spandrels) of purifying selection. Functional modules are the likely precursors of coadapted gene complexes, a unit of natural selection. The result is a predictable drive mechanism where "drive" is used rigorously, as in "meiotic drive." Rising morphological gain is expected given a supply of duplicate functional modules. All flowering plants have survived at least three large-scale duplications/diploidizations over the last 300 million years (Myr). An equivalent period of tetraploidy and body plan evolution may have ended for animals 500 million years ago (Mya). We argue that "balanced gene drive" is a sufficient explanation for the trend that the maximums of morphological complexity have gone up, and not down, in both plant and animal eukaryotic lineages.

Following tetraploidy in an Arabidopsis ancestor, genes were removed preferentially from one homeolog leaving clusters enriched in dose-sensitive genes

Approximately 90% of Arabidopsis' unique gene content is found in syntenic blocks that were formed during the most recent whole-genome duplication. Within these blocks, 28.6% of the genes have a retained pair; the remaining genes have been lost from one of the homeologs. We create a minimized genome by condensing local duplications to one gene, removing transposons, and including only genes within blocks defined by retained pairs. We use a moving average of retained and non-retained genes to find clusters of retention and then identify the types of genes that appear in clusters at frequencies above expectations. Significant clusters of retention exist for almost all chromosomal segments. Detailed alignments show that, for 85% of the genome, one homeolog was preferentially (1.6x) targeted for fractionation. This homeolog fractionation bias suggests an epigenetic mechanism. We find that islands of retention contain "connected genes," those genes predicted-by the gene balance hypothesis-to be resistant to removal because the products they encode interact with other products in a dose-sensitive manner, creating a web of dependency. Gene families that are overrepresented in clusters include those encoding components of the proteasome/protein modification complexes, signal transduction machinery, ribosomes, and transcription factor complexes. Gene pair fractionation following polyploidy or segmental duplication leaves a genome enriched for "connected" genes. These clusters of duplicate genes may help explain the evolutionary origin of coregulated chromosomal regions and new functional modules.

Leafing through the genomes of our major crop plants: strategies for capturing unique information

Crop plants not only have economic significance, but also comprise important botanical models for evolution and development. This is reflected by the recent increase in the percentage of publicly available sequence data that are derived from angiosperms. Further genome sequencing of the major crop plants will offer new learning opportunities, but their large, repetitive, and often polyploid genomes present challenges. Reduced-representation approaches - such as EST sequencing, methyl filtration and Cot-based cloning and sequencing - provide increased efficiency in extracting key information from crop genomes without full-genome sequencing. Combining these methods with phylogenetically stratified sampling to allow comparative genomic approaches has the potential to further accelerate progress in angiosperm genomics.

Grains of knowledge: genomics of model cereals

The economic and scientific importance of the cereals has motivated a rich history of research into their genetics, development, and evolution. The nearly completed sequence of the rice genome is emblematic of a transition to high-throughput genomics and computational biology that has also pervaded study of many other cereals. The relatively close (ca. <50 million years old) relationships among morphologically diverse cereals native to environments that sample much of global geographic diversity make the cereals particularly attractive for comparative studies of plant genome evolution. Extensive germplasm resources, largely a byproduct of their economic importance, together with growing collections of defined mutants, provide foundations for a host of post-genomic studies to shed more light on the relationship between sequence and function in this important group. Using the rapidly growing capabilities of several informatics resources, genomic data from model cereals are likely to be leveraged tremendously in the study and improvement of a wide range of crop plants that sustain much of the world's population, including many which still lack primary genomic resources.