Meliaceae genomes provide insights into wood development and limonoids biosynthesis.

Meliaceae is a useful plant family owing to its high-quality timber and its many limonoids that have pharmacological and biological activities. Although some genomes of Meliaceae species have been reported, many questions regarding their unique family features, namely wood quality and natural products, have not been answered. In this study, we provide the whole-genome sequence of Melia azedarach comprising 237.16 Mb with a contig N50 of 8.07 Mb, and an improved genome sequence of Azadirachta indica comprising 223.66 Mb with a contig N50 of 8.91 Mb. Moreover, genome skimming data, transcriptomes and other published genomes were comprehensively analysed to determine the genes and proteins that produce superior wood and valuable limonoids. Phylogenetic analysis of chloroplast genomes, single-copy gene families and single-nucleotide polymorphisms revealed that Meliaceae should be classified into two subfamilies: Cedreloideae and Melioideae. Although the Meliaceae species did not undergo additional whole-genome duplication events, the secondary wall biosynthetic genes of the woody Cedreloideae species, Toona sinensis, expanded significantly compared to those of A. indica and M. azedarach, especially in downstream transcription factors and cellulose/hemicellulose biosynthesis-related genes. Moreover, expanded special oxidosqualene cyclase catalogues can help diversify Sapindales skeletons, and the clustered genes that regulate terpene chain elongation, cyclization and modification would support their roles in limonoid biosynthesis. The expanded clans of terpene synthase, O-methyltransferase and cytochrome P450, which are mainly derived from tandem duplication, are responsible for the different limonoid classes among the species. These results are beneficial for further investigations of wood development and limonoid biosynthesis.

Linolenic Acid-Derived Oxylipins Inhibit Aflatoxin Biosynthesis in Aspergillus flavus through Activation of Imizoquin Biosynthesis.

Oxylipins play important signaling roles in aflatoxin (AF) biosynthesis in Aspergillus flavus. We previously showed that exogenous supply of autoxidated linolenic acid (AL) inhibited AF biosynthesis in A. flavus via oxylipins, but the molecular mechanism is still unknown. Here, we performed multiomics analyses of A. flavus grown in media with or without AL. Targeted metabolite analyses and quantitative reverse transcription (qRT)-polymerase chain reaction (PCR) showed that the imizoquin (IMQ) biosynthetic pathway was distinctly upregulated in the presence of AL. 13C-glucose labeling confirmed in parallel that the tricarboxylic acid cycle was also enhanced by AL, consistent with observed increases in mycelial growth. Moreover, we integrated thermal proteome profiling and molecular dynamics simulations to identify a potential receptor of AL; AL was found to interact with a transporter (ImqJ) located in the IMQ gene cluster, primarily through hydrophobic interactions. Further analyses of strains with an IMQ pathway transcription factor overexpressed or knocked out confirmed that this pathway was critical for AL-mediated inhibition of AF biosynthesis. Comparison of 22 assembled A. flavus and Aspergillus oryzae genomes showed that genes involved in the IMQ pathway were positively selected in A. oryzae. Taken together, the results of our study provide novel insights into oxylipin-mediated regulation of AF biosynthesis and suggest potential methods for preventing AF contamination of crops.

Identification and Functional Analysis of Long Non-Coding RNA (lncRNA) in Response to Seed Aging in Rice.

Many lncRNAs have been shown to play a vital role in aging processes. However, how lncRNAs regulate seed aging remains unknown. In this study, we performed whole transcriptome strand-specific RNA sequencing of samples from rice embryos, analyzed the differences in expression of rice seed lncRNAs before and after artificial aging treatment (AAT), and systematically screened 6002 rice lncRNAs. During the AAT period, the expression levels of most lncRNAs (454) were downregulated and only four were upregulated among the 458 differentially expressed lncRNAs (DELs). Cis- or trans-regulated target genes of the four upregulated lncRNAs were mainly related to base repair, while 454 downregulated lncRNAs were related to plant-pathogen interaction, plant hormones, energy metabolism, and secondary metabolism. The pathways of DEL target genes were similar with those of differentially expressed mRNAs (DEGs). A competing endogenous RNA (ceRNA) network composed of 34 lncRNAs, 24 microRNAs (miRNA), and 161 mRNAs was obtained. The cDNA sequence of lncRNA LNC_037529 was obtained by rapid amplification of cDNA ends (RACE) cloning with a total length of 1325 bp, a conserved 5' end, and a non-conserved 3' end. Together, our findings indicate that genome-wide selection for lncRNA downregulation was an important mechanism for rice seed aging. LncRNAs can be used as markers of seed aging in rice. These findings provide a future path to decipher the underlying mechanism associated with lncRNAs in seed aging.

Genome-Wide Analysis of CCT Transcript Factors to Identify Genes Contributing to Photoperiodic Flowering in Oryza rufipogon.

Photoperiod sensitivity is a dominant determinant for the phase transition in cereal crops. CCT (CONSTANS, CO-like, and TOC1) transcription factors (TFs) are involved in many physiological functions including the regulation of the photoperiodic flowering. However, the functional roles of CCT TFs have not been elucidated in the wild progenitors of crops. In this study, we identified 41 CCT TFs, including 19 CMF, 17 COL, and five PRR TFs in Oryza rufipogon, the presumed wild ancestor of Asian cultivated rice. There are thirty-eight orthologous CCT genes in Oryza sativa, of which ten pairs of duplicated CCT TFs are shared with O. rufipogon. We investigated daily expression patterns, showing that 36 OrCCT genes exhibited circadian rhythmic expression. A total of thirteen OrCCT genes were identified as putative flowering suppressors in O. rufipogon based on rhythmic and developmental expression patterns and transgenic phenotypes. We propose that OrCCT08, OrCCT24, and OrCCT26 are the strong functional alleles of rice DTH2, Ghd7, and OsPRR37, respectively. The SD treatment at 80 DAG stimulated flowering of the LD-grown O. rufipogon plants. Our results further showed that the nine OrCCT genes were significantly downregulated under the treatment. Our findings would provide valuable information for the construction of photoperiodic flowering regulatory network and functional characterization of the CCT TFs in both O. rufipogon and O. sativa.

Draft genome of the herbaceous bamboo Raddia distichophylla.

Bamboos are important nontimber forest plants widely distributed in the tropical and subtropical regions of Asia, Africa, America, and Pacific islands. They comprise the Bambusoideae in the grass family (Poaceae), including approximately 1700 described species in 127 genera. In spite of the widespread uses of bamboo for food, construction, and bioenergy, the gene repertoire of bamboo still remains largely unexplored. Raddia distichophylla (Schrad. ex Nees) Chase, belonging to the tribe Olyreae (Bambusoideae, Poaceae), is a diploid herbaceous bamboo with only slightly lignified stems. In this study, we report a draft genome assembly of the ∼589 Mb whole-genome sequence of R. distichophylla with a contig N50 length of 86.36 Kb. Repeat sequences account for ∼49.08% of the genome assembly, of which LTR retrotransposons occupy ∼35.99% of the whole genome. A total of 30,763 protein-coding genes were annotated in the R. distichophylla genome with an average transcript size of 2887 bp. Access to this herbaceous bamboo genome sequence will provide novel insights into biochemistry, molecular marker-assisted breeding programs, and germplasm conservation for bamboo species worldwide.

Comparative Multi-Omics of Tender Shoots from a Novel Evergrowing Tea Cultivar Provide Insight into the Winter Adaptation Mechanism.

Tea (Camellia sinensis [L.] O. Kuntze) tree is a perennial plant in which winter dormancy is an important biological adaptation to environmental changes. We discovered and reported a novel tea tree cultivar that can generate tender shoots in winter several years ago, but the molecular mechanism for this unique phenotype remains unknown . Here, we conducted comparative transcriptomics, proteomics and metabolomics along with phytohormone quantitation between the winter and spring tender shoots to investigate the physiological basis and putative regulatory mechanisms of its evergrowing character during winter. Our multi-omics study has led to the following findings. Gibberellin (GA) levels and key enzymes for GA biosynthesis and the signal transduction pathway were increased in the winter shoots, causing the ABA/GA content ratio to decrease, which might play a key regulatory role in maintaining normal growth during winter. The abundance of proteins, genes and metabolites involved in energy metabolism was all increased in winter shoots, indicating that energy is critical for continuous growth under the relatively weak-light and low-temperature environment. Abiotic resistance-related proteins and free amino acids were also increased in abundance in the winter shoots, which possibly represents an adaptation response to winter conditions. These results allowed us to hypothesize a novel molecular mechanism of adaptation for this unique tender shoot evergrowing in winter.

Improved hybrid de novo genome assembly and annotation of African wild rice, Oryza longistaminata, from Illumina and PacBio sequencing reads.

African wild rice Oryza longistaminata, one of the eight AA- genome species in the genus Oryza, possesses highly valued traits, such as the rhizomatousness for perennial rice breeding, strong tolerance to biotic and abiotic stresses, and high biomass production on poor soils. To obtain the high-quality reference genome for O. longistaminata we employed a hybrid assembly approach through incorporating Illumina and PacBio sequencing datasets. The final genome assembly comprised only 107 scaffolds and was approximately ∼363.5 Mb, representing ∼92.7% of the estimated African wild rice genome (∼392 Mb). The N50 lengths of the assembled contigs and scaffolds were ∼46.49 Kb and ∼6.83 Mb, indicating ∼3.72-fold and ∼18.8-fold improvement in length compared to the earlier released assembly (∼12.5 Kb and 364 Kb, respectively). Aided with Hi-C data and syntenic relationship with O. sativa, these assembled scaffolds were anchored into 12 pseudo-chromosomes. Genome annotation and comparative genomic analysis reveal that lineage-specific expansion of gene families that respond to biotic- and abiotic stresses are of great potential for mining novel alleles to overcome major diseases and abiotic adaptation in rice breeding programs. This reference genome of African wild rice will greatly enlarge the existing database of rice genome resources and unquestionably form a solid base to understand genomic basis underlying highly valued phenotypic traits and search for novel gene sources in O. longistaminata for the future rice breeding programs.

Draft genomes of two outcrossing wild rice, Oryza rufipogon and O. longistaminata, reveal genomic features associated with mating-system evolution.

Oryza rufipogon and O. longistaminata are important wild relatives of cultivated rice, harboring a promising source of novel genes for rice breeding programs. Here, we present de novo assembled draft genomes and annotation of O. rufipogon and O. longistaminata. Our analysis reveals a considerable number of lineage-specific gene families associated with the self-incompatibility (SI) and formation of reproductive separation. We show how lineage-specific expansion or contraction of gene families with functional enrichment of the recognition of pollen, thus enlightening their reproductive diversification. We documented a large number of lineage-specific gene families enriched in salt stress, antifungal response, and disease resistance. Our comparative analysis further shows a genome-wide expansion of genes encoding NBS-LRR proteins in these two outcrossing wild species in contrast to six other selfing rice species. Conserved noncoding sequences (CNSs) in the two wild rice genomes rapidly evolve relative to selfing rice species, resulting in the reduction of genomic variation owing to shifts of mating systems. We find that numerous genes related to these rapidly evolving CNSs are enriched in reproductive structure development, flower development, and postembryonic development, which may associate with SI in O. rufipogon and O. longistaminata.

The draft genome sequence of an upland wild rice species, Oryza granulata.

Exploiting novel gene sources from wild relatives has proven to be an efficient approach to advance crop genetic breeding efforts. Oryza granulata, with the GG genome type, occupies the basal position of the Oryza phylogeny and has the second largest genome (~882 Mb). As an upland wild rice species, it possesses renowned traits that distinguish it from other Oryza species, such as tolerance to shade and drought, immunity to bacterial blight and resistance to the brown planthopper. Here, we generated a 736.66-Mb genome assembly of O. granulata with 40,131 predicted protein-coding genes. With Hi-C data, for the first time, we anchored ~98.2% of the genome assembly to the twelve pseudo-chromosomes. This chromosome-length genome assembly of O. granulata will provide novel insights into rice genome evolution, enhance our efforts to search for new genes for future rice breeding programmes and facilitate the conservation of germplasm of this endangered wild rice species.

SMRT sequencing of the Oryza rufipogon genome reveals the genomic basis of rice adaptation.

Asian cultivated rice is believed to have been domesticated from a wild progenitor, Oryza rufipogon, offering promising sources of alleles for world rice improvement. Here we first present a high-quality chromosome-scale genome of the typical O. rufipogon. Comparative genomic analyses of O. sativa and its two wild progenitors, O. nivara and O. rufipogon, identified many dispensable genes functionally enriched in the reproductive process. We detected millions of genomic variants, of which large-effect mutations could affect agronomically relevant traits. We demonstrate how lineage-specific expansion of gene families may have contributed to the formation of reproduction isolation. We document thousands of genes with signatures of positive selection that are mainly involved in the reproduction and response to biotic- and abiotic stresses. We show that selection pressures may serve as forces to govern substantial genomic alterations that form the genetic basis of rapid evolution of mating and reproductive systems under diverse habitats.

The Chromosome-Based Rubber Tree Genome Provides New Insights into Spurge Genome Evolution and Rubber Biosynthesis.

The rubber tree, Hevea brasiliensis, produces natural rubber that serves as an essential industrial raw material. Here, we present a high-quality reference genome for a rubber tree cultivar GT1 using single-molecule real-time sequencing (SMRT) and Hi-C technologies to anchor the ∼1.47-Gb genome assembly into 18 pseudochromosomes. The chromosome-based genome analysis enabled us to establish a model of spurge chromosome evolution, since the common paleopolyploid event occurred before the split of Hevea and Manihot. We show recent and rapid bursts of the three Hevea-specific LTR-retrotransposon families during the last 10 million years, leading to the massive expansion by ∼65.88% (∼970 Mbp) of the whole rubber tree genome since the divergence from Manihot. We identify large-scale expansion of genes associated with whole rubber biosynthesis processes, such as basal metabolic processes, ethylene biosynthesis, and the activation of polysaccharide and glycoprotein lectin, which are important properties for latex production. A map of genomic variation between the cultivated and wild rubber trees was obtained, which contains ∼15.7 million high-quality single-nucleotide polymorphisms. We identified hundreds of candidate domestication genes with drastically lowered genomic diversity in the cultivated but not wild rubber trees despite a relatively short domestication history of rubber tree, some of which are involved in rubber biosynthesis. This genome assembly represents key resources for future rubber tree research and breeding, providing novel targets for improving plant biotic and abiotic tolerance and rubber production.

Development and Genome Sequencing of a Laboratory-Inbred Miniature Pig Facilitates Study of Human Diabetic Disease.

Pig has been proved to be a valuable large animal model used for research on diabetic disease. However, their translational value is limited given their distinct anatomy and physiology. For the last 30 years, we have been developing a laboratory Asian miniature pig inbred line (Bama miniature pig [BM]) from the primitive Bama xiang pig via long-term selective inbreeding. Here, we assembled a BM reference genome at full chromosome-scale resolution with a total length of 2.49 Gb. Comparative and evolutionary genomic analyses identified numerous variations between the BM and commercial pig (Duroc), particularly those in the genetic loci associated with the features advantageous to diabetes studies. Resequencing analyses revealed many differentiated gene loci associated with inbreeding and other selective forces. These together with transcriptome analyses of diabetic pig models provide a comprehensive genetic basis for resistance to diabetogenic environment, especially related to energy metabolism.

DNA methylome analysis provides evidence that the expansion of the tea genome is linked to TE bursts.

DNA methylation is essential for gene regulation, imprinting and silencing of transposable elements (TEs). Although bursts of transposable elements are common in many plant lineages, how plant DNA methylation is related to transposon bursts remains unclear. Here we explore the landscape of DNA methylation of tea, a species thought to have experienced a recent transposon burst event. This species possesses more transposable elements than any other sequenced asterids (potato, tomato, coffee, pepper and tobacco). The overall average DNA methylation levels were found to differ among the tea, potato and tomato genomes, and methylation at CHG sequence sites was found to be significantly higher in tea than that in potato or tomato. Moreover, the abundant TEs resulting from burst events not only resulted in tea developing a very large genome size, but also affected many genes involved in importantly biological processes, including caffeine, theanine and flavonoid metabolic pathway genes. In addition, recently transposed TEs were more heavily methylated than ancient ones, implying that DNA methylation is proportionate to the degree of TE silencing, especially on recent active ones. Taken together, our results show that DNA methylation regulates transposon silencing and may play a role in genome size expansion.

The Tea Tree Genome Provides Insights into Tea Flavor and Independent Evolution of Caffeine Biosynthesis.

Tea is the world's oldest and most popular caffeine-containing beverage with immense economic, medicinal, and cultural importance. Here, we present the first high-quality nucleotide sequence of the repeat-rich (80.9%), 3.02-Gb genome of the cultivated tea tree Camellia sinensis. We show that an extraordinarily large genome size of tea tree is resulted from the slow, steady, and long-term amplification of a few LTR retrotransposon families. In addition to a recent whole-genome duplication event, lineage-specific expansions of genes associated with flavonoid metabolic biosynthesis were discovered, which enhance catechin production, terpene enzyme activation, and stress tolerance, important features for tea flavor and adaptation. We demonstrate an independent and rapid evolution of the tea caffeine synthesis pathway relative to cacao and coffee. A comparative study among 25 Camellia species revealed that higher expression levels of most flavonoid- and caffeine- but not theanine-related genes contribute to the increased production of catechins and caffeine and thus enhance tea-processing suitability and tea quality. These novel findings pave the way for further metabolomic and functional genomic refinement of characteristic biosynthesis pathways and will help develop a more diversified set of tea flavors that would eventually satisfy and attract more tea drinkers worldwide.

The caterpillar fungus, Ophiocordyceps sinensis, genome provides insights into highland adaptation of fungal pathogenicity.

To understand the potential genetic basis of highland adaptation of fungal pathogenicity, we present here the ~116 Mb de novo assembled high-quality genome of Ophiocordyceps sinensis endemic to the Qinghai-Tibetan Plateau. Compared with other plain-dwelling fungi, we find about 3.4-fold inflation of the O. sinensis genome due to a rapid amplification of long terminal repeat retrotransposons that occurred ~38 million years ago in concert with the uplift of the plateau. We also observe massive removal of thousands of genes related to the transport process and energy metabolism. O. sinensis displays considerable lineage-specific expansion of gene families functionally enriched in the adaptability of low-temperature of cold tolerance, fungal pathogenicity and specialized host infection. We detect signals of positive selection for genes involved in peroxidase and hypoxia to enable its highland adaptation. Resequencing and analyzing 31 whole genomes of O. sinensis, representing nearly all of its geographic range, exhibits latitude-based population divergence and nature selection for population inhabitation towards higher altitudes on the Qinghai-Tibetan Plateau.

Rapid and Recent Evolution of LTR Retrotransposons Drives Rice Genome Evolution During the Speciation of AA-Genome Oryza Species.

The dynamics of long terminal repeat (LTR) retrotransposons and their contribution to genome evolution during plant speciation have remained largely unanswered. Here, we perform a genome-wide comparison of all eight Oryza AA-genome species, and identify 3911 intact LTR retrotransposons classified into 790 families. The top 44 most abundant LTR retrotransposon families show patterns of rapid and distinct diversification since the species split over the last ∼4.8 MY (million years). Phylogenetic and read depth analyses of 11 representative retrotransposon families further provide a comprehensive evolutionary landscape of these changes. Compared with Ty1-copia, independent bursts of Ty3-gypsy retrotransposon expansions have occurred with the three largest showing signatures of lineage-specific evolution. The estimated insertion times of 2213 complete retrotransposons from the top 23 most abundant families reveal divergent life histories marked by speedy accumulation, decline, and extinction that differed radically between species. We hypothesize that this rapid evolution of LTR retrotransposons not only divergently shaped the architecture of rice genomes but also contributed to the process of speciation and diversification of rice.

Overexpressed Proteins in Hypervirulent Clade 8 and Clade 6 Strains of Escherichia coli O157:H7 Compared to E. coli O157:H7 EDL933 Clade 3 Strain.

Escherichia coli O157:H7 is responsible for severe diarrhea and hemolytic uremic syndrome (HUS), and predominantly affects children under 5 years. The major virulence traits are Shiga toxins, necessary to develop HUS and the Type III Secretion System (T3SS) through which bacteria translocate effector proteins directly into the host cell. By SNPs typing, E. coli O157:H7 was separated into nine different clades. Clade 8 and clade 6 strains were more frequently associated with severe disease and HUS. In this study, we aimed to identify differentially expressed proteins in two strains of E. coli O157:H7 (clade 8 and clade 6), obtained from cattle and compared them with the well characterized reference EDL933 strain (clade 3). Clade 8 and clade 6 strains show enhanced pathogenicity in a mouse model and virulence-related properties. Proteins were extracted and analyzed using the TMT-6plex labeling strategy associated with two dimensional liquid chromatography and mass spectrometry in tandem. We detected 2241 proteins in the cell extract and 1787 proteins in the culture supernatants. Attention was focused on the proteins related to virulence, overexpressed in clade 6 and 8 strains compared to EDL933 strain. The proteins relevant overexpressed in clade 8 strain were the curli protein CsgC, a transcriptional activator (PchE), phage proteins, Stx2, FlgM and FlgD, a dienelactone hydrolase, CheW and CheY, and the SPATE protease EspP. For clade 6 strain, a high overexpression of phage proteins was detected, mostly from Stx2 encoding phage, including Stx2, flagellin and the protease TagA, EDL933_p0016, dienelactone hydrolase, and Haemolysin A, amongst others with unknown function. Some of these proteins were analyzed by RT-qPCR to corroborate the proteomic data. Clade 6 and clade 8 strains showed enhanced transcription of 10 out of 12 genes compared to EDL933. These results may provide new insights in E. coli O157:H7 mechanisms of pathogenesis.

The complete chloroplast genome sequence of desert poplar (Populus euphratica).

The complete chloroplast sequence of the desert poplar (Populus euphratica), a plant well-adapted to salt stress, was determined in this study. The genome consists of 156,766 bp containing a pair of inverted repeats (IRs) of 16,591 bp separated by a large single-copy region and a small single-copy region of 84,888 bp and 27,646 bp, respectively. The chloroplast genome contains 130 known genes, including 89 protein-coding genes, 8 ribosomal RNA genes, and 37 tRNA genes; 18 of these are located in the inverted repeat region.