Identification and characterization of a rhamnosyltransferase involved in rutin biosynthesis in Fagopyrum esculentum (common buckwheat)

Rutin, a 3-rutinosyl quercetin, is a representative flavonoid distributed in many plant species, and is highlighted for its therapeutic potential. In this study, we purified uridine diphosphate-rhamnose: quercetin 3-O-glucoside 6″-O-rhamnosyltransferase and isolated the corresponding cDNA (FeF3G6″RhaT) from seedlings of common buckwheat (Fagopyrum esculentum). The recombinant FeF3G6″RhaT enzyme expressed in Escherichia coli exhibited 6″-O-rhamnosylation activity against flavonol 3-O-glucoside and flavonol 3-O-galactoside as substrates, but showed only faint activity against flavonoid 7-O-glucosides. Tobacco cells expressing FeF3G6″RhaT converted the administered quercetin into rutin, suggesting that FeF3G6″RhaT can function as a rhamnosyltransferase in planta. Quantitative PCR analysis on several organs of common buckwheat revealed that accumulation of FeF3G6″RhaT began during the early developmental stages of rutin-accumulating organs, such as flowers, leaves, and cotyledons. These results suggest that FeF3G6″RhaT is involved in rutin biosynthesis in common buckwheat.     

Chromosome‐level de novo genome assembly of Sarcophaga peregrina provides insights into the evolutionary adaptation of flesh flies

Sarcophaga peregrina is considered to be of great ecological, medical and forensic significance, and has unusual biological characteristics such as an ovoviviparous reproductive pattern and adaptation to feed on carrion. The availability of a high‐quality genome will help to further reveal the mechanisms underlying these charcateristics. Here we present a de novo‐assembled genome at chromosome scale for S. peregrina. The final assembled genome was 560.31 Mb with contig N50 of 3.84 Mb. Hi‐C scaffolding reliably anchored six pseudochromosomes, accounting for 97.76% of the assembled genome. Moreover, 45.70% of repeat elements were identified in the genome. A total of 14,476 protein‐coding genes were functionally annotated, accounting for 92.14% of all predicted genes. Phylogenetic analysis indicated that S. peregrina and S. bullata diverged ~ 7.14 million years ago. Comparative genomic analysis revealed expanded and positively selected genes related to biological features that aid in clarifying its ovoviviparous reproduction and carrion‐feeding adaptations, such as lipid metabolism, olfactory receptor activity, antioxidant enzymes, proteolysis and serine‐type endopeptidase activity. Protein‐coding genes associated with ovoviparity, such as yolk proteins, transferrin and acid sphingomyelinase, were identified. This study provides a valuable genomic resource for S. peregrina, and sheds insight into further revealing the underlying molecular mechanisms of adaptive evolution.     

Genomic analyses provide insights into the polyploidization-driven herbicide adaptation in Leptochloa weeds

Polyploidy confers a selective advantage under stress conditions; however, whether polyploidization mediates enhanced herbicide adaptation remains largely unknown. Tetraploid Leptochloa chinensis is a notorious weed in the rice ecosystem, causing severe yield loss in rice. In China, L. chinensis has only one sister species, the diploid L. panicea, whose damage is rarely reported. To gain insights into the effects of polyploidization on herbicide adaptation, we first assembled a high-quality genome of L. panicea and identified genome structure variations with L. chinensis. Moreover, we identified herbicide-resistance genes specifically expanded in L. chinensis, which may confer a greater herbicide adaptability in L. chinensis. Analysis of gene retention and loss showed that five herbicide target-site genes and several herbicide nontarget-site resistance gene families were retained during polyploidization. Notably, we identified three pairs of polyploidization-retained genes including LcABCC8, LcCYP76C1 and LcCYP76C4 that may enhance herbicide resistance. More importantly, we found that both copies of LcCYP76C4 were under herbicide selection during the spread of L. chinensis in China. Furthermore, we identified another gene potentially involved in herbicide resistance, LcCYP709B2, which is also retained during polyploidization and under selection. This study provides insights into the genomic basis of the enhanced herbicide adaptability of Leptochloa weeds during polyploidization and provides guidance for the precise and efficient control of polyploidy weeds.     

High-quality chromosome-level genome assembly of Litsea coreana L. provides insights into Magnoliids evolution and flavonoid biosynthesis

The magnoliid Litsea coreana has been the subject of a substantial amount of research owing to its production of many flavonoid metabolites, high food processing value, and a controversial phylogenetic position. For this study, we assembled a high-grade genome at the chromosome scale and annotation of L. coreana that was anchored to 12 chromosomes. The total genome was 1139.45 Mb, while the N50 scaffold was 97.18 Mb long. The analysis of phylogenetic trees constructed by different methods show that the phylogeny of Magnoliids is inconsistent, indicating that the differentiation process of monocots, eudicots, and Magnoliids still remains in dispute. An ancient whole-genome duplication (WGD) event was shown to have occurred before the Magnoliales and Laurels had differentiated. Subsequently, an independent WGD appeared in the Lauralean lineage. A total of 27 types of flavonoids were detected in all five tissues of L. coreana. Chalcone synthases (CHSs) that are responsible for production of flavonoids have been validated at the bioinformatics level. The retention of comparative genomic analyses of the CHS gene family showed that this family had contracted significantly in L. coreana. Our research further elaborated the evolution of Lauraceae and perfected the genetic basis of flavonoid biosynthesis in L. coreana.Significance statementProvides evidence that determines the evolutionary status of Magnoliids. The chalcone synthase gene family was significantly contracted in Litsea coreana.     

The Chimonanthus salicifolius genome provides insight into magnoliid evolution and flavonoid biosynthesis

Chimonanthus salicifolius, a member of the Calycanthaceae of magnoliids, is one of the most famous medicinal plants in Eastern China. Here, we report a chromosome‐level genome assembly of C. salicifolius, comprising 820.1 Mb of genomic sequence with a contig N50 of 2.3 Mb and containing 36 651 annotated protein‐coding genes. Phylogenetic analyses revealed that magnoliids were sister to the eudicots. Two rounds of ancient whole‐genome duplication were inferred in the C. salicifolious genome. One is shared by Calycanthaceae after its divergence with Lauraceae, and the other is in the ancestry of Magnoliales and Laurales. Notably, long genes with > 20 kb in length were much more prevalent in the magnoliid genomes compared with other angiosperms, which could be caused by the length expansion of introns inserted by transposon elements. Homologous genes within the flavonoid pathway for C. salicifolius were identified, and correlation of the gene expression and the contents of flavonoid metabolites revealed potential critical genes involved in flavonoids biosynthesis. This study not only provides an additional whole‐genome sequence from the magnoliids, but also opens the door to functional genomic research and molecular breeding of C. salicifolius.     

A near-complete genome assembly of Brassica rapa provides new insights into the evolution of centromeres

Brassica rapa comprises many important cultivated vegetables and oil crops. However, Chiifu v3.0, the current B. rapa reference genome, still contains hundreds of gaps. Here, we presented a near-complete genome assembly of B. rapa Chiifu v4.0, which was 424.59 Mb with only two gaps, using Oxford Nanopore Technology (ONT) ultralong-read sequencing and Hi-C technologies. The new assembly contains 12 contigs, with a contig N50 of 38.26 Mb. Eight of the ten chromosomes were entirely reconstructed in a single contig from telomere to telomere. We found that the centromeres were mainly invaded by ALE and CRM long terminal repeats (LTRs). Moreover, there is a high divergence of centromere length and sequence among B. rapa genomes. We further found that centromeres are enriched for Copia invaded at 0.14 MYA on average, while pericentromeres are enriched for Gypsy LTRs invaded at 0.51 MYA on average. These results indicated the different invasion mechanisms of LTRs between the two structures. In addition, a novel repetitive sequence PCR630 was identified in the pericentromeres of B. rapa. Overall, the near-complete genome assembly, B. rapa Chiifu v4.0, offers valuable tools for genomic and genetic studies of Brassica species and provides new insights into the evolution of centromeres.     

Chromosome-scale genome assembly and population genomics provide insights into the adaptation,domestication, and flavonoid metabolism of Chinese plum

Globally, commercialized plum cultivars are mostly diploid Chinese plums (Prunus salicina Lindl.), also known as Japanese plums, and are one of the most abundant and variable fruit tree species. To advance Prunus genomic research, we present a chromosome-scale P. salicina genome assembly, constructed using an integrated strategy that combines Illumina, Oxford Nanopore, and high-throughput chromosome conformation capture (Hi-C) sequencing. The high-quality genome assembly consists of a 318.6-Mb sequence (contig N50 length of 2.3 Mb) with eight pseudo-chromosomes. The expansion of the P. salicina genome is led by recent segmental duplications and a long terminal repeat burst of approximately 0.2 Mya. This resulted in a significant expansion of gene families associated with flavonoid metabolism and plant resistance, which impacted fruit flavor and increased species adaptability. Population structure and domestication history suggest that Chinese plum may have originated from South China and provides a domestication route with accompanying genomic variations. Selection sweep and genetic diversity analysis enabled the identification of several critical genes associated with flowering time, stress tolerance, and flavonoid metabolism, demonstrating the essential roles of related pathways during domestication. Furthermore, we reconstructed and exploited flavonoid–anthocyanin metabolism using multi-omics analysis in Chinese plum and proposed a complete metabolic pathway. Collectively, our results will facilitate further candidate gene discovery for important agronomic traits in Chinese plum and provide insights into future functional genomic studies and DNA-informed breeding.     

Disruption of specific flavonoid genes enhances the accumulation of flavonoid enzymes and end-products in Arabidopsis seedlings

Polyclonal antibodies were developed against the flavonoid biosynthetic enzymes, CHS, CHI, F3H, FLS, and LDOX from Arabidopsis thaliana. These antibodies were used to perform the first detailed analysis of coordinate expression of flavonoid metabolism at the protein level. The pattern of flavonoid enzyme expression over the course of seedling development was consistent with previous studies indicating that chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), and flavonol synthase (FLS) are encoded by early genes while leucoanthocyanidin dioxygenase (LDOX) is encoded by a late gene. This sequential expression may underlie the variations in flavonoid end-products produced during this developmental stage, as determined by HPLC analysis, which includes a shift in the ratio of the flavonols, quercetin and kaempferol. Moreover, immunoblot and HPLC analyses revealed that several transparent testa lines blocked at intermediate steps of the flavonoid pathway actually accumulated higher levels of specific flavonoid enzymes and end-products. These results suggest that specific intermediates may act as inducers of flavonoid metabolism.     

Horizontal gene transfer among microbial genomes: new insights from complete genome analysis

The determination and analysis of complete genome sequences has led to the suggestion that horizontal gene transfer may be much more extensive than previously appreciated. Many of these studies, however, rely on evidence that could be generated by forces other than gene transfer including selection, variable evolutionary rates, and biased sampling.     

A genome‐wide linkage map for the house sparrow (Passer domesticus) provides insights into the evolutionary history of the avian genome

The house sparrow is an important model species for studying physiological, ecological and evolutionary processes in wild populations. Here, we present a medium density, genome wide linkage map for house sparrow (Passer domesticus) that has aided the assembly of the house sparrow reference genome, and that will provide an important resource for ongoing mapping of genes controlling important traits in the ecology and evolution of this species. Using a custom house sparrow 10 K iSelect Illumina SNP chip we have assigned 6,498 SNPs to 29 autosomal linkage groups, based on a mean of 430 informative meioses per SNP. The map was constructed by combining the information from linkage with that of the physical position of SNPs within scaffold sequences in an iterative process. Averaged between the sexes; the linkage map had a total length of 2,004 cM, with a longer map for females (2,240 cM) than males (1,801 cM). Additionally, recombination rates also varied along the chromosomes. Comparison of the linkage map to the reference genomes of zebra finch, collared flycatcher and chicken, showed a chromosome fusion of the two avian chromosomes 8 and 4A in house sparrow. Lastly, information from the linkage map was utilized to conduct analysis of linkage disequilibrium (LD) in eight populations with different effective population sizes (N_e) in order to quantify the background level LD. Together, these results aid the design of future association studies, facilitate the development of new genomic tools and support the body of research that describes the evolution of the avian genome.     

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.     

Apolygus lucorum genome provides insights into omnivorousness and mesophyll feeding

Apolygus lucorum (Miridae) is an omnivorous pest that occurs worldwide and is notorious for the serious damage it causes to various crops and substantial economic losses. Although some studies have examined the biological characteristics of the mirid bug, no reference genome is available in Miridae, limiting in‐depth studies of this pest. Here, we present a chromosome‐scale reference genome of A. lucorum, the first sequenced Miridae species. The assembled genome size was 1.02 Gb with a contig N50 of 785 kb. With Hi‐C scaffolding, 1,016 Mb contig sequences were clustered, ordered and assembled into 17 large scaffolds with scaffold N50 length 68 Mb, each corresponding to a natural chromosome. Numerous transposable elements occur in this genome and contribute to the large genome size. Expansions of genes associated with omnivorousness and mesophyll feeding such as those related to digestion, chemosensory perception, and detoxification were observed in A. lucorum, suggesting that gene expansion contributed to its strong environmental adaptability and severe harm to crops. We clarified that a salivary enzyme polygalacturonase is unique in mirid bugs and has significantly expanded in A. lucorum, which may contribute to leaf damage from this pest. The reference genome of A. lucorum not only facilitates biological studies of Hemiptera as well as an understanding of the damage mechanism of mesophyll feeding, but also provides a basis on which to develop efficient control technologies for mirid bugs.     

Spatio‐temporal patterns of genome evolution in allotetraploid species of the genus Oryza

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 arisen <2 million years ago, and (v) there are no clear distinctions in the patterns of genome evolution in the diploid versus polyploid species.     

薏苡黑粉菌基因组测序揭示其物种分类及适应性进化

以薏苡黑粉菌病瘿为研究对象,采用固体和液体培养方法对薏苡黑粉菌进行分离纯化培养和形态观察;利用Illumina、Pacbio测序和Hi-C辅助基因组组装技术对黑粉菌进行测序、组装和注释,分析其基因组结构和组成特征。通过与其他6个黑粉菌物种基因组进行比较基因组学研究,分析薏苡黑粉菌的基因家族进化和系统发育关系。在pH 7的固体PDA和液体PDB培养基中实现了薏苡黑粉菌的纯化培养,固体培养菌落表现为白色、表面隆起、具褶皱,较湿润且不透明的形态特征。Illumina、Pacbio和Hi-C组装结果显示,薏苡黑粉菌具有20对染色体(2n=40),基因组大小20.093 417 Mb,GC含量53.81%;共预测和注释到7 476个蛋白编码基因和164个非编码基因,占整个基因组长度的61.02%和0.285 3%;查找到3 674个散在重复序列和8 139个串联重复序列,分别占整个基因组长度的2.754 4%和1.605 4%。同源基因家族聚类显示,7个黑粉菌物种Ustilago coicis、U. bromivora、U. maydis、U. hordei、Sporisorium reilianum、S. scitamineum和S. graminicola的编码基因被聚类为6 999个基因家族,其中共有基因家族5 379个。薏苡黑粉菌包含5 339个单拷贝基因、105个多拷贝基因、89个特异基因、752个其他基因和188个未分类基因。基因家族进化和系统发育聚类表明,薏苡黑粉菌U. coicis在~60.9 MYA分化出来,与雀麦黑粉菌U. bromivova、大麦坚黑粉菌U. hordei聚为一支,表现为更近的亲缘关系;而玉米黑粉菌U. maydis在~66.6 MYA分化出来与3个孢堆黑粉菌物种S. scitamineum、S. graminicola和S. reilianum聚为一支。     

Differential genome duplication and fish diversity

The duplication of genes and entire genomes arebelieved to be important mechanisms underlyingmorphological variation and functionalinnovation in the evolution of life andespecially for the broad diversity observed inthe speciation of fishes. How did these fishspecies and their genetic diversity arise? Theoccurrence of three rounds of genomeduplication during vertebrate evolution mightexplain why many gene families are typicallyabout half the size in land vertebrates as theyare in fishes. However, mechanisms of geneticdiversity in fish lineages need to be furtherexplained. Here we propose that differentialgenome duplication of from two to six roundsoccurred in different fish lines, offering newopportunities during the radiation of fishlineages. This model provides a fundamentalbasis for the understanding of theirspeciation, diversity and evolution.     

Purification,molecular cloning and functional characterization of flavonoid C‐glucosyltransferases from Fagopyrum esculentum M. (buckwheat) cotyledon

C‐Glycosides are characterized by their C–C bonds in which the anomeric carbon of the sugar moieties is directly bound to the carbon atom of aglycon. C‐Glycosides are remarkably stable, as their C–C bonds are resistant to glycosidase or acid hydrolysis. A variety of plant species are known to accumulate C‐glycosylflavonoids; however, the genes encoding for enzymes that catalyze C‐glycosylation of flavonoids have been identified only from Oryza sativa (rice) and Zea mays (maize), and have not been identified from dicot plants. In this study, we identified the C‐glucosyltransferase gene from the dicot plant Fagopyrum esculentum M. (buckwheat). We purified two isozymes from buckwheat seedlings that catalyze C‐glucosylation of 2‐hydroxyflavanones, which are expressed specifically in the cotyledon during seed germination. Following purification we isolated the cDNA corresponding to each isozyme [FeCGTa (UGT708C1) and FeCGTb (UGT708C2)]. When expressed in Escherichia coli, both proteins demonstrated C‐glucosylation activity towards 2‐hydroxyflavanones, dihydrochalcone, trihydroxyacetophenones and other related compounds with chemical structures similar to 2′,4′,6′‐trihydroxyacetophenone. Molecular phylogenetic analysis of plant glycosyltransferases shows that flavonoid C‐glycosyltransferases form a different clade with other functionally analyzed plant glycosyltransferases.