Phylogenomic Analysis of micro-RNA Involved in Juvenile to Flowering-Stage Transition in Photophilic Rice and Its Sister Species

Vegetative to reproductive phase transition in phototropic plants is an important developmental process and is sequentially mediated by the expression of micro-RNA MIR172. To obtain insight into the evolution, adaptation, and function of MIR172 in photophilic rice and its wild relatives, we analyzed the genescape of a 100 kb segment harboring MIR172 homologs from 11 genomes. The expression analysis of MIR172 revealed its incremental accumulation from the 2-leaf to 10-leaf stage, with maximum expression coinciding with the flag-leaf stage in rice. Nonetheless, the microsynteny analysis of MIR172s revealed collinearity within the genus Oryza, but a loss of synteny was observed in (i) MIR172A in O. barthii (AA) and O. glaberima (AA); (ii) MIR172B in O. brachyantha (FF); and (iii) MIR172C in O. punctata (BB). Phylogenetic analysis of precursor sequences/region of MIR172 revealed a distinct tri-modal clade of evolution. The genomic information generated in this investigation through comparative analysis of MIRNA, suggests mature MIR172s to have evolved in a disruptive and conservative mode amongst all Oryza species with a common origin of descent. Further, the phylogenomic delineation provided an insight into the adaptation and molecular evolution of MIR172 to changing environmental conditions (biotic and abiotic) of phototropic rice through natural selection and the opportunity to harness untapped genomic regions from rice wild relatives (RWR).

Contrasting modes of macro and microsynteny evolution in a eukaryotic subphylum

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

Fine mapping of Rf5 region for a sorghum fertility restorer gene and microsynteny analysis across grass species

Cytoplasmic male sterility (CMS) is widely used to control pollination in the production of commercial F₁ hybrid seed in sorghum. So far, 6 major fertility restorer genes, Rf1 to Rf6, have been reported in sorghum. Here, we fine-mapped the Rf5 locus on sorghum chromosome 5 using descendant populations of a 'Nakei MS-3A' × 'JN43' cross. The Rf5 locus was narrowed to a 140-kb region in BTx623 genome (161-kb in JN43) with 16 predicted genes, including 6 homologous to the rice fertility restorer Rf1 (PPR.1 to PPR.6). These 6 homologs have tandem pentatricopeptide repeat (PPR) motifs. Many Rf genes encode PPR proteins, which bind RNA transcripts and modulate gene expression at the RNA level. No PPR genes were detected at the Rf5 locus on the corresponding homologous chromosome of rice, foxtail millet, or maize, so this gene cluster may have originated by chromosome translocation and duplication after the divergence of sorghum from these species. Comparison of the sequences of these genes between fertile and CMS lines identified PPR.4 as the most plausible candidate gene for Rf5.

Sawfly Genomes Reveal Evolutionary Acquisitions That Fostered the Mega-Radiation of Parasitoid and Eusocial Hymenoptera

The tremendous diversity of Hymenoptera is commonly attributed to the evolution of parasitoidism in the last common ancestor of parasitoid sawflies (Orussidae) and wasp-waisted Hymenoptera (Apocrita). However, Apocrita and Orussidae differ dramatically in their species richness, indicating that the diversification of Apocrita was promoted by additional traits. These traits have remained elusive due to a paucity of sawfly genome sequences, in particular those of parasitoid sawflies. Here, we present comparative analyses of draft genomes of the primarily phytophagous sawfly Athalia rosae and the parasitoid sawfly Orussus abietinus. Our analyses revealed that the ancestral hymenopteran genome exhibited traits that were previously considered unique to eusocial Apocrita (e.g., low transposable element content and activity) and a wider gene repertoire than previously thought (e.g., genes for CO₂ detection). Moreover, we discovered that Apocrita evolved a significantly larger array of odorant receptors than sawflies, which could be relevant to the remarkable diversification of Apocrita by enabling efficient detection and reliable identification of hosts.

SWEET Transporters and the Potential Functions of These Sequences in Tea (Camellia sinensis)

Tea (Camellia sinensis) is an important economic beverage crop. Its flowers and leaves could be used as healthcare tea for its medicinal value. SWEET proteins were recently identified in plants as sugar transporters, which participate in diverse physiological processes, including pathogen nutrition, seed filling, nectar secretion, and phloem loading. Although SWEET genes have been characterized and identified in model plants, such as Arabidopsis thaliana and Oryza sativa, there is very little knowledge of these genes in C. sinensis. In this study, 28 CsSWEETs were identified in C. sinensis and further phylogenetically divided into four subfamilies with A. thaliana. These identified CsSWEETs contained seven transmembrane helixes (TMHs) which were generated by an ancestral three-TMH unit with an internal duplication experience. Microsynteny analysis revealed that the large-scale duplication events were the main driving forces for members from CsSWEET family expansion in C. sinensis. The expression profiles of the 28 CsSWEETs revealed that some genes were highly expressed in reproductive tissues. Among them, CsSWEET1a might play crucial roles in the efflux of sucrose, and CsSWEET17b could control fructose content as a hexose transporter in C. sinensis. Remarkably, CsSWEET12 and CsSWEET17c were specifically expressed in flowers, indicating that these two genes might be involved in sugar transport during flower development. The expression patterns of all CsSWEETs were differentially regulated under cold and drought treatments. This work provided a systematic understanding of the members from the CsSWEET gene family, which would be helpful for further functional studies of CsSWEETs in C. sinensis.

Bracon brevicornis Genome Showcases the Potential of Linked-Read Sequencing in Identifying a Putative Complementary Sex Determiner Gene

Bracon brevicornis is an ectoparasitoid of a wide range of larval-stage Lepidopterans, including several pests of important crops, such as the corn borer, Ostrinia nubilalis. It is also one of the earliest documented cases of complementary sex determination in Hymenoptera. Here, we present the linked-read-based genome of B. brevicornis, complete with an ab initio-derived annotation and protein comparisons with fellow braconids, Fopius arisanus and Diachasma alloeum. We demonstrate the potential of linked-read assemblies in exploring regions of heterozygosity and search for structural and homology-derived evidence of the complementary sex determiner gene (csd).

Evolutionary History of GLIS Genes Illuminates Their Roles in Cell Reprograming and Ciliogenesis

The GLIS family transcription factors, GLIS1 and GLIS3, potentiate generation of induced pluripotent stem cells (iPSCs). In contrast, another GLIS family member, GLIS2, suppresses cell reprograming. To understand how these disparate roles arose, we examined evolutionary origins and genomic organization of GLIS genes. Comprehensive phylogenetic analysis shows that GLIS1 and GLIS3 originated during vertebrate whole genome duplication, whereas GLIS2 is a sister group to the GLIS1/3 and GLI families. This result is consistent with their opposing functions in cell reprograming. Glis1 evolved faster than Glis3, losing many protein-interacting motifs. This suggests that Glis1 acquired new functions under weakened evolutionary constraints. In fact, GLIS1 induces induced pluripotent stem cells more strongly. Transcriptomic data from various animal embryos demonstrate that glis1 is maternally expressed in some tetrapods, whereas vertebrate glis3 and invertebrate glis1/3 genes are rarely expressed in oocytes, suggesting that vertebrate (or tetrapod) Glis1 acquired a new expression domain and function as a maternal factor. Furthermore, comparative genomic analysis reveals that glis1/3 is part of a bilaterian-specific gene cluster, together with rfx3, ndc1, hspb11, and lrrc42. Because known functions of these genes are related to cilia formation and function, the last common ancestor of bilaterians may have acquired this cluster by shuffling gene order to establish more sophisticated epithelial tissues involving cilia. This evolutionary study highlights the significance of GLIS1/3 for cell reprograming, development, and diseases in ciliated organs such as lung, kidney, and pancreas.

A Chinese White Pear (Pyrus bretschneideri) BZR Gene PbBZR1 Act as a Transcriptional Repressor of Lignin Biosynthetic Genes in Fruits

BZR transcription factors play essential roles in plant growth and environmental stimuli, and they are also the positive regulators of Brassinosteroid (BR) signal transduction in diverse plants. In addition, BZR TFs, as crucial regulators of BR synthesis, may have multiple stress-resistance functions and their related regulatory mechanisms have been well illustrated in model plants. Here, we carried out a genome-wide identification of BZR members in Chinese pear (Pyrus bretschneideri) and identified 13 members. By comparative analysis in five Rosaceae genomes, BZR members in the pear genome may have undergone large-scale duplication events during evolution. Purifying selection played an important role in almost all of the orthologous and paralogous gene pairs. According to the expression analysis of the PbBZRs during fruit development, three PbBZRs were selected for detailed analysis. Transcriptional activation assays presented that PbBZR1 repressed the promoters of P. bretschneideri lignin biosynthetic genes, such as PbCES9, PbCOMT3, and PbHCT6. Our study traces the evolution of BZR gene family members in Rosaceae genomes and illustrates that the rates of gene loss and gain are far from equilibrium in different species. At the same time, our results suggest that PbBZR1 may be involved in the negative regulation of lignin biosynthesis.

Genome-Wide Analysis Characterization and Evolution of SBP Genes in Fragaria vesca, Pyrus bretschneideri, Prunus persica and Prunus mume

The SQUAMOSA promoter binding protein (SBP)-box proteins are plant-specific transcriptional factors in plants. SBP TFs are known to play important functions in a diverse development process and also related in the process of evolutionary novelties. SBP gene family has been characterized in several plant species, but little is known about molecular evolution, functional divergence and comprehensive study of SBP gene family in Rosacea. We carried out genome-wide investigations and identified 14, 32, 17, and 17 SBP genes from four Rosacea species (Fragaria vesca, Pyrus bretschneideri, Prunus persica and Prunus mume, respectively). According to phylogenetic analysis arranged the SBP protein sequences in seven groups. Localization of SBP genes presented an uneven distribution on corresponding chromosomes of Rosacea species. Our analyses designated that the SBP genes duplication events (segmental and tandem) and divergence. In addition, due to highly conserved structure pattern of SBP genes, recommended that highly conserved region of microsyneteny in the Rosacea species. Type I and II functional divergence was detected among various amino acids in SBP proteins, while there was no positive selection according to substitutional model analysis using PMAL software. These results recommended that the purifying selection might be leading force during the evolution process and dominate conservation of SBP genes in Rosacea species according to environmental selection pressure analysis. Our results will provide basic understanding and foundation for future research insights on the evolution of the SBP genes in Rosacea.

Genome-Wide Analysis Suggests the Relaxed Purifying Selection Affect the Evolution of WOX Genes in Pyrus bretschneideri, Prunus persica, Prunus mume, and Fragaria vesca

WUSCHEL-related homeobox (WOX) family is one of the largest group of transcription factors (TFs) specifically found in plant kingdom. WOX TFs play an important role in plant development processes and evolutionary novelties. Although the roles of WOXs in Arabidopsis and rice have been well-studied, however, little are known about the relationships among the main clades in the molecular evolution of these genes in Rosaceae. Here, we carried out a genome-wide analysis and identified 14, 10, 10, and 9 of WOX genes from four Rosaceae species (Fragaria vesca, Prunus persica, Prunus mume, and Pyrus bretschneideri, respectively). According to evolutionary analysis, as well as amino acid sequences of their homodomains, these genes were divided into three clades with nine subgroups. Furthermore, due to the conserved structural patterns among these WOX genes, it was proposed that there should exist some highly conserved regions of microsynteny in the four Rosaceae species. Moreover, most of WOX gene pairs were presented with the conserved orientation among syntenic genome regions. In addition, according to substitution models analysis using PMAL software, no significant positive selection was detected, but type I functional divergence was identified among certain amino acids in WOX protein. These results revealed that the relaxed purifying selection might be the main driving force during the evolution of WOX genes in the tested Rosaceae species. Our result will be useful for further precise research on evolution of the WOX genes in family Rosaceae.

Genome-Wide Association Mapping of Flowering and Ripening Periods in Apple

Deciphering the genetic control of flowering and ripening periods in apple is essential for breeding cultivars adapted to their growing environments. We implemented a large Genome-Wide Association Study (GWAS) at the European level using an association panel of 1,168 different apple genotypes distributed over six locations and phenotyped for these phenological traits. The panel was genotyped at a high-density of SNPs using the Axiom®Apple 480 K SNP array. We ran GWAS with a multi-locus mixed model (MLMM), which handles the putatively confounding effect of significant SNPs elsewhere on the genome. Genomic regions were further investigated to reveal candidate genes responsible for the phenotypic variation. At the whole population level, GWAS retained two SNPs as cofactors on chromosome 9 for flowering period, and six for ripening period (four on chromosome 3, one on chromosome 10 and one on chromosome 16) which, together accounted for 8.9 and 17.2% of the phenotypic variance, respectively. For both traits, SNPs in weak linkage disequilibrium were detected nearby, thus suggesting the existence of allelic heterogeneity. The geographic origins and relationships of apple cultivars accounted for large parts of the phenotypic variation. Variation in genotypic frequency of the SNPs associated with the two traits was connected to the geographic origin of the genotypes (grouped as North+East, West and South Europe), and indicated differential selection in different growing environments. Genes encoding transcription factors containing either NAC or MADS domains were identified as major candidates within the small confidence intervals computed for the associated genomic regions. A strong microsynteny between apple and peach was revealed in all the four confidence interval regions. This study shows how association genetics can unravel the genetic control of important horticultural traits in apple, as well as reduce the confidence intervals of the associated regions identified by linkage mapping approaches. Our findings can be used for the improvement of apple through marker-assisted breeding strategies that take advantage of the accumulating additive effects of the identified SNPs.

Comparative Genomic Analysis of the GRF Genes in Chinese Pear (Pyrus bretschneideri Rehd), Poplar (Populous), Grape (Vitis vinifera), Arabidopsis and Rice (Oryza sativa)

Growth-regulating factors (GRFs) are plant-specific transcription factors that have important functions in regulating plant growth and development. Previous studies on GRF family members focused either on a single or a small set of genes. Here, a comparative genomic analysis of the GRF gene family was performed in poplar (a model tree species), Arabidopsis (a model plant for annual herbaceous dicots), grape (one model plant for perennial dicots), rice (a model plant for monocots) and Chinese pear (one of the economical fruit crops). In total, 58 GRF genes were identified, 12 genes in rice (Oryza sativa), 8 genes in grape (Vitis vinifera), 9 genes in Arabidopsis thaliana, 19 genes in poplar (Populus trichocarpa) and 10 genes in Chinese pear (Pyrus bretschneideri). The GRF genes were divided into five subfamilies based on the phylogenetic analysis, which was supported by their structural analysis. Furthermore, microsynteny analysis indicated that highly conserved regions of microsynteny were identified in all of the five species tested. And Ka/Ks analysis revealed that purifying selection plays an important role in the maintenance of GRF genes. Our results provide basic information on GRF genes in five plant species and lay the foundation for future research on the functions of these genes.

Structural, Evolutionary, and Functional Analysis of the Class III Peroxidase Gene Family in Chinese Pear (Pyrus bretschneideri)

Peroxidases (PRXs) are widely existed in various organisms and could be divided into different types according to their structures and functions. Specifically, the Class III Peroxidase, a plant-specific multi-gene family, involves in many physiological processes, such as the metabolism of auxin, the extension and thickening of cell wall, as well as the formation of lignin. By searching the pear genome database, 94 non-redundant PRXs from Pyrus bretschneideri (PbPRXs) were identified. Subsequently, analysis of phylogenetic relationships, gene structures, conserved motifs, and microsynteny was performed. These PbPRXs were unevenly distributed among 17 chromosomes of pear. In addition, 26 segmental duplication events but only one tandem duplication were occurred in these PbPRXs, implying segmental duplication was the main contributor to the expansion of the PbPRX family. By the Ka/Ks analysis, 26 out of 27 duplicated PbPRXs has experienced purifying selection. Twenty motifs were identified in PbPRXs based on the MEME analysis, 11 of which were enriched in pear. A total of 41 expressed genes were identified from ESTs of pear fruit. According to qRT-PCR, the expression trends of five PbPRXs in subgroup C were consistent with the change of lignin content during pear fruit development. So we inferred that the five PbPRXs were candidate genes involved in the lignin synthesis pathway. These results provided useful information for further researches of PRX genes in pear.

Microstructure of a Brassica rapa genome segment homoeologous to the resistance gene cluster on Arabidopsis chromosome 4

Genome evolution is a continuous process and genomic rearrangement occurs both within and between species. With the sequencing of the Arabidopsis thaliana genome, comparative genetics and genomics offer new insights into plant biology. The genus Brassica offers excellent opportunities with which to compare genomic synteny so as to reveal genome evolution. During a previous genetic analysis of clubroot resistance in Brassica rapa, we identified a genetic region that is highly collinear with Arabidopsis chromosome 4. This region corresponds to a disease resistance gene cluster in the A. thaliana genome. Relying on synteny with Arabidopsis, we fine-mapped the region and found that the location and order of the markers showed good correspondence with those in Arabidopsis. Microsynteny on a physical map indicated an almost parallel correspondence, with a few rearrangements such as inversions and insertions. The results show that this genomic region of Brassica is conserved extensively with that of Arabidopsis and has potential as a disease resistance gene cluster, although the genera diverged 20 million years ago.

A BAC library of the SP80-3280 sugarcane variety (saccharum sp.) and its inferred microsynteny with the sorghum genome

BACKGROUND

Sugarcane breeding has significantly progressed in the last 30 years, but achieving additional yield gains has been difficult because of the constraints imposed by the complex ploidy of this crop. Sugarcane cultivars are interspecific hybrids between Saccharum officinarum and Saccharum spontaneum. S. officinarum is an octoploid with 2n = 80 chromosomes while S. spontaneum has 2n = 40 to 128 chromosomes and ploidy varying from 5 to 16. The hybrid genome is composed of 70-80% S. officinaram and 5-20% S. spontaneum chromosomes and a small proportion of recombinants. Sequencing the genome of this complex crop may help identify useful genes, either per se or through comparative genomics using closely related grasses. The construction and sequencing of a bacterial artificial chromosome (BAC) library of an elite commercial variety of sugarcane could help assembly the sugarcane genome.

RESULTS

A BAC library designated SS_SBa was constructed with DNA isolated from the commercial sugarcane variety SP80-3280. The library contains 36,864 clones with an average insert size of 125 Kb, 88% of which has inserts larger than 90 Kb. Based on the estimated genome size of 760-930 Mb, the library exhibits 5-6 times coverage the monoploid sugarcane genome. Bidirectional BAC end sequencing (BESs) from a random sample of 192 BAC clones sampled genes and repetitive elements of the sugarcane genome. Forty-five per cent of the total BES nucleotides represents repetitive elements, 83% of which belonging to LTR retrotransposons. Alignment of BESs corresponding to 42 BACs to the genome sequence of the 10 sorghum chromosomes revealed regions of microsynteny, with expansions and contractions of sorghum genome regions relative to the sugarcane BAC clones. In general, the sampled sorghum genome regions presented an average 29% expansion in relation to the sugarcane syntenic BACs.

CONCLUSION

The SS_SBa BAC library represents a new resource for sugarcane genome sequencing. An analysis of insert size, genome coverage and orthologous alignment with the sorghum genome revealed that the library presents whole genome coverage. The comparison of syntenic regions of the sorghum genome to 42 SS_SBa BES pairs revealed that the sorghum genome is expanded in relation to the sugarcane genome.

Local coexpression domains in the genome of rice show no microsynteny with Arabidopsis domains

Chromosomal coexpression domains are found in a number of different genomes under various developmental conditions. The size of these domains and the number of genes they contain vary. Here, we define local coexpression domains as adjacent genes where all possible pair-wise correlations of expression data are higher than 0.7. In rice, such local coexpression domains range from predominantly two genes, up to 4, and make up approximately 5% of the genomic neighboring genes, when examining different expression platforms from the public domain. The genes in local coexpression domains do not fall in the same ontology category significantly more than neighboring genes that are not coexpressed. Duplication, orientation or the distance between the genes does not solely explain coexpression. The regulation of coexpression is therefore thought to be regulated at the level of chromatin structure. The characteristics of the local coexpression domains in rice are strikingly similar to such domains in the Arabidopsis genome. Yet, no microsynteny between local coexpression domains in Arabidopsis and rice could be identified. Although the rice genome is not yet as extensively annotated as the Arabidopsis genome, the lack of conservation of local coexpression domains may indicate that such domains have not played a major role in the evolution of genome structure or in genome conservation.

An ancestral oomycete locus contains late blight avirulence gene Avr3a, encoding a protein that is recognized in the host cytoplasm

The oomycete Phytophthora infestans causes late blight, the potato disease that precipitated the Irish famines in 1846 and 1847. It represents a reemerging threat to potato production and is one of >70 species that are arguably the most devastating pathogens of dicotyledonous plants. Nevertheless, little is known about the molecular bases of pathogenicity in these algae-like organisms or of avirulence molecules that are perceived by host defenses. Disease resistance alleles, products of which recognize corresponding avirulence molecules in the pathogen, have been introgressed into the cultivated potato from a wild species, Solanum demissum, and R1 and R3a have been identified. We used association genetics to identify Avr3a and show that it encodes a protein that is recognized in the host cytoplasm, where it triggers R3a-dependent cell death. Avr3a resides in a region of the P. infestans genome that is colinear with the locus containing avirulence gene ATR1(NdWsB) in Hyaloperonospora parasitica, an oomycete pathogen of Arabidopsis. Remarkably, distances between conserved genes in these avirulence loci were often similar, despite intervening genomic variation. We suggest that Avr3a has undergone gene duplication and that an allele evading recognition by R3a arose under positive selection.