Search

Elucidating the cold tolerance mechanism of Paeonia lactiflora, which is one of the most valuable ornamental and medicinal plants in Asia, fundamentally impacts its breeding and production. The glycerol-3-phosphate acyltransferase (GPAT) gene plays a pivotal role in cold resistance in a variety of plant species. Here, we cloned the P. lactiflora GPAT gene, determined its expression pattern, and tested its role in cold resistance. We obtained the full-length P. lactiflora GPAT gene using tissue-cultured seedlings and real-time polymerase chain reaction and rapid amplification of cDNA ends analyses. We named this gene PlGPAT in P. lactiflora. Phylogenetic analysis indicates that the PlGPAT gene is closely related with the GPAT genes in core eudicots. The phylogenetic tree containing 31 angiosperm species based on GPAT protein sequences is largely consistent with the known phylogeny in flowering plants. We conducted a time-course PlGPAT expression analysis and demonstrated that PlGPAT expression is correlated with low-temperature stress. Our results suggest that the PlGPAT gene plays an important role in regulating cold resistance in P. lactiflora.

The genus Oryza is a model system for the study of molecular evolution over time scales ranging from a few thousand to 15 million years. Using 13 reference genomes spanning the Oryza species tree, we show that despite few large-scale chromosomal rearrangements rapid species diversification is mirrored by lineage-specific emergence and turnover of many novel elements, including transposons, and potential new coding and noncoding genes. Our study resolves controversial areas of the Oryza phylogeny, showing a complex history of introgression among different chromosomes in the young 'AA' subclade containing the two domesticated species. This study highlights the prevalence of functionally coupled disease resistance genes and identifies many new haplotypes of potential use for future crop protection. Finally, this study marks a milestone in modern rice research with the release of a complete long-read assembly of IR 8 'Miracle Rice', which relieved famine and drove the Green Revolution in Asia 50 years ago.

In an effort to identify newly evolved genes in rice, we searched the genomes of Asian-cultivated rice Oryza sativa ssp. japonica and its wild progenitors, looking for lineage-specific genes. Using genome pairwise comparison of approximately 20-Mb DNA sequences from the chromosome 3 short arm (Chr3s) in six rice species, O. sativa, O. nivara, O. rufipogon, O. glaberrima, O. barthii, and O. punctata, combined with synonymous substitution rate tests and other evidence, we were able to identify potential recently duplicated genes, which evolved within the last 1 Myr. We identified 28 functional O. sativa genes, which likely originated after O. sativa diverged from O. glaberrima. These genes account for around 1% (28/3,176) of all annotated genes on O. sativa's Chr3s. Among the 28 new genes, two recently duplicated segments contained eight genes. Fourteen of the 28 new genes consist of chimeric gene structure derived from one or multiple parental genes and flanking targeting sequences. Although the majority of these 28 new genes were formed by single or segmental DNA-based gene duplication and recombination, we found two genes that were likely originated partially through exon shuffling. Sequence divergence tests between new genes and their putative progenitors indicated that new genes were most likely evolving under natural selection. We showed all 28 new genes appeared to be functional, as suggested by Ka/Ks analysis and the presence of RNA-seq, cDNA, expressed sequence tag, massively parallel signature sequencing, and/or small RNA data. The high rate of new gene origination and of chimeric gene formation in rice may demonstrate rice's broad diversification, domestication, its environmental adaptation, and the role of new genes in rice speciation.

Oryza (23 species; 10 genome types) contains the world's most important food crop — rice. Although the rice genome serves as an essential tool for biological research, little is known about the evolution of the other Oryza genome types. They contain a historical record of genomic changes that led to diversification of this genus around the world as well as an untapped reservoir of agriculturally important traits. To investigate the evolution of the collective Oryza genome, we sequenced and compared nine orthologous genomic regions encompassing the Adh1-Adh2 genes (from six diploid genome types) with the rice reference sequence. Our analysis revealed the architectural complexities and dynamic evolution of this region that have occurred over the past ∼15 million years. Of the 46 intact genes and four pseudogenes in the japonica genome, 38 (76%) fell into eight multigene families. Analysis of the evolutionary history of each family revealed independent and lineage-specific gain and loss of gene family members as frequent causes of synteny disruption. Transposable elements were shown to mediate massive replacement of intergenic space (>95%), gene disruption, and gene/gene fragment movement. Three cases of long-range structural variation (inversions/deletions) spanning several hundred kilobases were identified that contributed significantly to genome diversification.

Tandem gene duplication is one of the major gene duplication mechanisms in eukaryotes, as illustrated by the prevalence of gene family clusters. Tandem duplicated paralogs usually share the same regulatory element, and as a consequence, they are likely to perform similar biological functions. Here, we provide an example of a newly evolved tandem duplicate acquiring novel functions, which were driven by positive selection. CG32708, CG32706, and CG6999 are 3 clustered genes residing in the X chromosome of Drosophila melanogaster. CG6999 and CG32708 have been examined for their molecular population genetic properties (Thornton and Long 2005). We further investigated the evolutionary forces acting on these genes with greater sample sizes and a broader approach that incorporate between-species divergence, using more variety of statistical methods. We explored the possible functional implications by characterizing the tissue-specific and developmental expression patterns of these genes. Sequence comparison of species within D. melanogaster subgroup reveals that this 3-gene cluster was created by 2 rounds of tandem gene duplication in the last 5 Myr. Based on phylogenetic analysis, CG32708 is clearly the parental copy that is shared by all species. CG32706 appears to have originated in the ancestor of Drosophila simulans and D. melanogaster about 5 Mya, and CG6999 is the newest duplicate that is unique to D. melanogaster. All 3 genes have different expression profiles, and CG6999 has in addition acquired a novel transcript. Biased polymorphism frequency spectrum, linkage disequilibrium, nucleotide substitution, and McDonald–Kreitman analyses suggested that the evolution of CG6999 and CG32706 were driven by positive Darwinian selection.

A genomic pattern of new gene origination is often dependent on a genomic method that can efficiently identify a statistically adequate number of recently originated genes. The heterochromatic regions have often been viewed as genomic deserts with low coding potential and thus a low flux of new genes. However, increasing reports revealed unexpected roles of heterochromatic regions in the evolution of genes and genomes. We identified recently retroposed genes that originated in heterochromatic regions in Drosophila, by developing microarray-based comparative genomic hybridization (CGH) with multiple species. This new gene family, named Ifc-2h, originated in the common ancestor of the clade of D. simulans, D. mauritiana, and D. sechellia. The sequence features and phylogenetic distribution indicated that Ifc-2h resulted from the retroposition from its parental gene, Infertile crescent (Ifc), and integrated into heterochromatic region of common ancestor of the three sibling species 2 million years ago. Expression analysis revealed that Ifc-2h had developed a new expression pattern by recruiting a putative regulatory element from its target sequence. The distribution of indel variation in Ifc-2h of D. simulans and D. mauritiana revealed a significant sequence constraint, suggesting that the Ifc-2h gene may be functional. These analyses cast fresh insight into the evolution of heterochromatin and the origin of its coding regions.

Anthocyanin is a major pigment in vegetative and floral organs of most plants and plays an important role in plant evolution. The anthocyanin regulatory genes are responsible for regulating transcription of genes in the anthocyanin synthetic pathway. To assess evolutionary significance of sequence variation and evaluate the phylogenetic utility of an anthocyanin regulatory gene, we compared nucleotide sequences of the myc-like anthocyanin regulatory gene in the genus of dogwoods (Cornus: Cornaceae). Phylogenetic analyses demonstrate that the myc-like anthocyanin regulatory gene has potential as an informative phylogenetic marker at different taxonomic levels, depending on the data set considered (DNA or protein sequences) and regions applied (exons or introns). Pairwise nonsynonymous and synonymous substitution rate tests and codon-based substitution models were applied to characterize variation and to identify sites under diversifying selection. Mosaic evolution and heterogeneous rates among different domains and sites were detected.

Despite knowledge that polyploidy is widespread and a major evolutionary force in flowering plant diversification, detailed comparative molecular studies on polyploidy have been confined to only a few species and families. The genus Oryza is composed of 23 species that are classified into ten distinct 'genome types' (six diploid and four polyploid), and is emerging as a powerful new model system to study polyploidy. Here we report the identification, sequence and comprehensive comparative annotation of eight homoeologous genomes from a single orthologous region (Adh1-Adh2) from four allopolyploid species representing each of the known Oryza genome types (BC, CD, HJ and KL). Detailed comparative phylogenomic analyses of these regions within and across species and ploidy levels provided several insights into the spatio-temporal dynamics of genome organization and evolution of this region in 'natural' polyploids of Oryza. The major findings of this study are that: (i) homoeologous genomic regions within the same nucleus experience both independent and parallel evolution, (ii) differential lineage-specific selection pressures do not occur between polyploids and their diploid progenitors, (iii) there have been no dramatic structural changes relative to the diploid ancestors, (iv) a variation in the molecular evolutionary rate exists between the two genomes in the BC complex species even though the BC and CD polyploid species appear to have arisen2 million years ago, and (v) there are no clear distinctions in the patterns of genome evolution in the diploid versus polyploid species.

Despite general observations of non-random genomic distribution of new genes, it is unclear whether or not new genes preferentially occur in certain genomic regions driven by related molecular mechanisms. Using 1.5 Mb of genomic sequences from short arms of chromosome 3 of Oryza glaberrima and O. punctata, we conducted a comparative genomic analysis with the reference O. sativa ssp. japonica genome. We identified a 60-kb segment located in the middle of the subtelomeric region of chromosome 3, which is unique to the species O. sativa. The region contained gene duplicates that occurred in Asian cultivated rice species that diverged from the ancestor of Asian and African cultivated rice one million years ago (MYA). For the 12 genes and one complete retrotransposon identified in this segment in O. sativa ssp. japonica, we searched for their parental genes. The high similarity between duplicated paralogs further supports the recent origination of these genes. We found that this segment was recently generated through multiple independent gene recombination and transposon insertion events. Among the 12 genes, we found that five had chimeric gene structures derived from multiple parental genes. Nine out of the 12 new genes seem to be functional, as suggested by Ka/Ks analysis and the presence of cDNA and/or MPSS data. Furthermore, for the eight transcribed genes, at least two genes could be classified as defense or stress response-related genes. Given these findings, and the fact that subtelomeres are associated with high rates of recombination and transcription, it is likely that subtelomeres may facilitate gene recombination and transposon insertions and serve as hot spots for new gene origination in rice genomes.

Most eukaryotic centromeres contain large quantities of repetitive DNA, such as satellite repeats and retrotransposons. Unlike most transposons in plant genomes, the centromeric retrotransposon (CR) family is conserved over long evolutionary periods among a majority of the grass species. CR elements are highly concentrated in centromeres, and are likely to play a role in centromere function. In order to study centromere evolution in the Oryza (rice) genus, we sequenced the orthologous region to centromere 8 of Oryza sativa from a related species, Oryza brachyantha. We found that O. brachyantha does not have the canonical CRR (CR of rice) found in the centromeres of all other Oryza species. Instead, a new Ty3-gypsy (Metaviridae) retroelement (FRetro3) was found to colonize the centromeres of this species. This retroelement is found in high copy numbers in the O. brachyantha genome, but not in other Oryza genomes, and based on the dating of long terminal repeats (LTRs) of FRetro3 it was amplified in the genome in the last few million years. Interestingly, there is a high level of removal of FRetro3 based on solo-LTRs to full-length elements, and this rapid turnover may have played a role in the replacement of the canonical CRR with the new element by active deletion. Comparison with previously described ChIP cloning data revealed that FRetro3 is found in CENH3-associated chromatin sequences. Thus, within a single lineage of the Oryza genus, the canonical component of grass centromeres has been replaced with a new retrotransposon that has all the hallmarks of a centromeric retroelement.

PEX11 gene family has been shown to be involved in peroxisome biogenesis but very little is known about this gene family in rice. Here we show that five putative PEX11 genes (OsPEX11-1–5) present in rice genome and each contain three conserved motifs. The PEX11 sequences from rice and other species can be classified into three major groups. Among the five rice PEX11 genes, OsPEX11-2 and -3 are most likely duplicated. Expression profile and RT-PCR analysis suggested that the members of PEX11 family in rice had differential expression patterns: OsPEX11-1 and OsPEX11-4 had higher expression levels in leaf tissues than in the other tissues, OsPEX11-2 was detected only in germinated seeds, OsPEX11-3 was expressed predominantly in endosperm and germinated seeds, and OsPEX11-5 was expressed in all the tissues investigated. We also observed that the rice PEX11 genes had differential expression patterns under different abiotic stresses. OsPEX11-1 and OsPEX11-4 were induced by abscisic acid (ABA), hydrogen peroxide (H2O2), salt and low nitrogen stress conditions. OsPEX11-3 was responsive to ABA and H2O2 treatments, and OsPEX11-5 was responsive to ABA, H2O2, and salt treatments. However, OsPEX11-2 had no response to any of the stresses. Our results suggest that the rice PEX11 genes have diversification not only in sequences but also in expression patterns under normal and various stress conditions.

The evolutionary pattern of the myc-like anthocyanin regulatory gene antR-Cor was examined in the dwarf dogwood species complex (Cornus Subgenus Arctocrania) that contains two diploid species (C. canadensis and C. suecica), their putative hybrids with intermediate phenotypes, and a tetraploid derivative (C. unalaschkensis). Full-length sequences of this gene ( approximately 4 kb) were sequenced and characterized for 47 dwarf dogwood samples representing all taxa categories from 43 sites in the Pacific Northwest. Analysis of nucleotide diversity indicated departures from neutral evolution, due most likely to local population structure. Neighbor-joining and haplotype network analyses show that sequences from the tetraploid and diploid intermediates are much more strongly diverged from C. suecica than from C. canadensis, and that the intermediate phenotypes may represent an ancestral group to C. canadensis rather than interspecific hybrids. Seven amino acid mutations that are potentially linked to myc-like anthocyanin regulatory gene function correlate with petal colors differences that characterize the divergence between two diploid species and the tetraploid species in this complex. The evidence provides a working hypothesis for testing the role of the gene in speciation and its link to the petal coloration. Sequencing and analysis of additional nuclear genes will be necessary to resolve questions about the evolution of the dwarf dogwood complex.

Protein disulfide isomerases (PDIs) are molecular chaperones that contain thioredoxin (TRX) domains and aid in the formation of proper disulfide bonds during protein folding. To identify plant PDI-like (PDIL) proteins, a genome-wide search of Arabidopsis (Arabidopsis thaliana) was carried out to produce a comprehensive list of 104 genes encoding proteins with TRX domains. Phylogenetic analysis was conducted for these sequences using Bayesian and maximum-likelihood methods. The resulting phylogenetic tree showed that evolutionary relationships of TRX domains alone were correlated with conserved enzymatic activities. From this tree, we identified a set of 22 PDIL proteins that constitute a well-supported clade containing orthologs of known PDIs. Using the Arabidopsis PDIL sequences in iterative BLAST searches of public and proprietary sequence databases, we further identified orthologous sets of 19 PDIL sequences in rice (Oryza sativa) and 22 PDIL sequences in maize (Zea mays), and resolved the PDIL phylogeny into 10 groups. Five groups (I–V) had two TRX domains and showed structural similarities to the PDIL proteins in other higher eukaryotes. The remaining five groups had a single TRX domain. Two of these (quiescin-sulfhydryl oxidase-like and adenosine 5′-phosphosulfate reductase-like) had putative nonisomerase enzymatic activities encoded by an additional domain. Two others (VI and VIII) resembled small single-domain PDIs from Giardia lamblia, a basal eukaryote, and from yeast. Mining of maize expressed sequence tag and RNA-profiling databases indicated that members of all of the single-domain PDIL groups were expressed throughout the plant. The group VI maize PDIL ZmPDIL5-1 accumulated during endoplasmic reticulum stress but was not found within the intracellular membrane fractions and may represent a new member of the molecular chaperone complement in the cell.

Chitinase proteins have evolved and diversified almost in all organisms ranging from prokaryotes to eukaryotes. During evolution, internal repeats may appear in amino acid sequences of proteins which alter the structural and functional features. Here we deciphered the internal repeats from Chitinase and characterized the structural similarities between them. Out of 24 diverse Chitinase sequences selected, six sequences (2CJL, 2DSK, 2XVP, 2Z37, 3EBV and 3HBE) did not contain any internal repeats of amino acid sequences. Ten sequences contained repeats of length

Gene duplication is an important mechanism for the origination of functional novelties in organisms. We performed a comparative genome analysis to systematically estimate recent lineage specific gene duplication events in Arabidopsis thaliana and further investigate whether and how these new duplicate genes (NDGs) play a functional role in the evolution and adaption of A. thaliana. We accomplished this using syntenic relationship among four closely related species, A. thaliana, A. lyrata, Capsella rubella and Brassica rapa. We identified 100 NDGs, showing clear origination patterns, whose parental genes are located in syntenic regions and/or have clear orthologs in at least one of three outgroup species. All 100 NDGs were transcribed and under functional constraints, while 24% of the NDGs have differential expression patterns compared to their parental genes. We explored the underlying evolutionary forces of these paralogous pairs through conducting neutrality tests with sequence divergence and polymorphism data. Evolution of about 15% of NDGs appeared to be driven by natural selection. Moreover, we found that 3 NDGs not only altered their expression patterns when compared with parental genes, but also evolved under positive selection. We investigated the underlying mechanisms driving the differential expression of NDGs and their parents, and found a number of NDGs had different cis-elements and methylation patterns from their parental genes. Overall, we demonstrated that NDGs acquired divergent cis-elements and methylation patterns and may experience sub-functionalization or neo-functionalization influencing the evolution and adaption of A. thaliana.

Data from four DNA regions (rbcL, matK, 26S rDNA, and ITS) as well as extant and fossil morphology were used to reconstruct the phylogeny and biogeographic history of an intercontinentally disjunct plant group, the cornelian cherries of Cornus (dogwoods). The study tests previous hypotheses on the relative roles of two Tertiary land bridges, the North Atlantic land bridge (NALB) and the Bering land bridge (BLB), in plant migration across continents. Three approaches, the Bayesian, nonparametric rate smoothing (NPRS), and penalized likelihood (PL) methods, were employed to estimate the times of geographic isolations of species. Dispersal and vicariance analysis (DIVA) was performed to infer the sequence and directionality of biogeographic pathways. Results of phylogenetic analyses suggest that among the six living species, C. sessilis from western North America represents the oldest lineage, followed by C. volkensii from Africa. The four Eurasian species form a clade consisting of two sister pairs, C. mas-C. officinalis and C. chinensis-C. eydeana. Results of DIVA and data from fossils and molecular dating indicate that the cornelian cherry subgroup arose in Europe as early as the Paleocene. Fossils confirm that the group was present in North America by the late Paleocene, consistent with the DIVA predictions that, by the end of the Eocene, it had diversified into several species and expanded its distribution to North America via the NALB and to Africa via the last direct connection between Eurasia and Africa prior to the Miocene, or via long-distance dispersal. The cornelian cherries in eastern Asia appear to be derived from two independent dispersal events from Europe. These events are inferred to have occurred during the Oligocene and Miocene. This study supports the hypothesis that the NALB served as an important land bridge connecting the North American and European floras, as well as connecting American and African floras via Europe during the early Tertiary.