Draft genome sequence of the rubber tree Hevea brasiliensis

Pan-genome and phylogenomic analyses highlight Hevea species delineation and rubber trait evolution

The para rubber tree (Hevea brasiliensis) is the world's sole commercial source of natural rubber, a vital industrial raw material. However, the narrow genetic diversity of this crop poses challenges for rubber breeding. Here, we generate high-quality de novo genome assemblies for three H. brasiliensis cultivars, two H. brasiliensis wild accessions, and three other Hevea species (H. nitida, H. pauciflora, and H. benthamiana). Through analyzing genomes of 94 Hevea accessions, we identify five distinct lineages that do not align with their previous species delineations. We discover multiple accessions with hybrid origins between these lineages, indicating incomplete reproductive isolation between them. Only two out of four wild lineages have been introduced to commercial rubber cultivars. Furthermore, we reveal that the rubber production traits emerged following the development of a large REF/SRPP gene cluster and its functional specialization in rubber-producing laticifers within this genus. These findings would enhance rubber breeding and benefit research communities.

Plant Biotechnology-An Indispensable Tool for Crop Improvement

Traditional plant breeding has helped to increase food production dramatically over the past five decades, and many countries have managed to produce enough food for the growing population, particularly in the developing world [...].

Heterozygosity analysis of spontaneous 2n female gametes and centromere mapping of the diploid Hevea brasiliensis based on full-sib triploid populations

KEY MESSAGE

Unreduced megagametophytes via second-division restitution were confirmed through heterozygosity analysis, and four candidate physical centromeres of rubber were located for the first time. The evaluation of maternal heterozygosity restitution (MHR) is vital in identifying the mechanism of 2n gametogenesis and assessing the utilization value of 2n gametes. In this study, three full-sib triploid populations were employed to evaluate the MHR of 2n female gametes of rubber tree clone GT1 and to confirm their genetic derivation. The 2n female gametes of GT1 were derived from second-division restitution (SDR) and transmitted more than half of the parental heterozygosity. In addition, low recombination frequency markers were developed, and four candidate physical centromeres of rubber tree were located for the first time. The confirmation that 2n female gametes of rubber tree clone GT1 are derived from SDR provides insights into the molecular mechanisms of 2n gametogenesis. In addition, the identified centromere location will aid in the development of centromeric markers for the rapid identification of the 2n gametogenesis mechanism.

Advances in Genome Sequencing and Natural Rubber Biosynthesis in Rubber-Producing Plants

Natural rubber (cis-1,4-polyisoprene, NR) is an important raw material utilized widely in the manufacturing of medical, agricultural, and industrial products. Rubber tree (Hevea brasiliensis) and several alternative rubber-producing plants (Taraxacum kok-saghyz, Lactuca sativa, and Parthenium argentatum) have the capability to produce high-quality NR. With the progress of genome sequencing, similar rubber biosynthesis pathways have been discovered among different rubber-producing plant species. NR is synthesized and stored in rubber particles, which are specialized organelles comprising a hydrophobic NR core surrounded by a lipid monolayer and membrane-bound proteins. The rubber transferase complex is considered to be the pivotal enzyme involved in catalyzing NR biosynthesis. However, the exact compositions of the RT complex in rubber-producing plants remain elusive and poorly understood. Here, we review the progress of genome sequencing, natural rubber biosynthesis, and the components of the RT complex in rubber-producing plants. We emphasize that identifying the detailed components of the RT complex holds great significance for exploring the mechanism of NR biosynthesis and accelerating molecular breeding in rubber-producing plants.

Genomic insight into domestication of rubber tree

Understanding the genetic basis of rubber tree (Hevea brasiliensis) domestication is crucial for further improving natural rubber production to meet its increasing demand worldwide. Here we provide a high-quality H. brasiliensis genome assembly (1.58 Gb, contig N50 of 11.21 megabases), present a map of genome variations by resequencing 335 accessions and reveal domestication-related molecular signals and a major domestication trait, the higher number of laticifer rings. We further show that HbPSK5, encoding the small-peptide hormone phytosulfokine (PSK), is a key domestication gene and closely correlated with the major domestication trait. The transcriptional activation of HbPSK5 by myelocytomatosis (MYC) members links PSK signaling to jasmonates in regulating the laticifer differentiation in rubber tree. Heterologous overexpression of HbPSK5 in Russian dandelion (Taraxacum kok-saghyz) can increase rubber content by promoting laticifer formation. Our results provide an insight into target genes for improving rubber tree and accelerating the domestication of other rubber-producing plants.

Chromosome-level wild Hevea brasiliensis genome provides new tools for genomic-assisted breeding and valuable loci to elevate rubber yield

The rubber tree (Hevea brasiliensis) is grown in tropical regions and is the major source of natural rubber. Using traditional breeding approaches, the latex yield has increased by sixfold in the last century. However, the underlying genetic basis of rubber yield improvement is largely unknown. Here, we present a high-quality, chromosome-level genome sequence of the wild rubber tree, the first report on selection signatures and a genome-wide association study (GWAS) of its yield traits. Population genomic analysis revealed a moderate population divergence between the Wickham clones and wild accessions. Interestingly, it is suggestive that H. brasiliensis and six relatives of the Hevea genus might belong to the same species. The selective sweep analysis found 361 obvious signatures in the domesticated clones associated with 245 genes. In a 15-year field trial, GWAS identified 155 marker-trait associations with latex yield, in which 326 candidate genes were found. Notably, six genes related to sugar transport and metabolism, and four genes related to ethylene biosynthesis and signalling are associated with latex yield. The homozygote frequencies of the causal nonsynonymous SNPs have been greatly increased under selection, which may have contributed to the fast latex yield improvement during the short domestication history. Our study provides insights into the genetic basis of the latex yield trait and has implications for genomic-assisted breeding by offering valuable resources in this new domesticated crop.

The rubber tree kinome: Genome-wide characterization and insights into coexpression patterns associated with abiotic stress responses

The protein kinase (PK) superfamily constitutes one of the largest and most conserved protein families in eukaryotic genomes, comprising core components of signaling pathways in cell regulation. Despite its remarkable relevance, only a few kinase families have been studied in Hevea brasiliensis. A comprehensive characterization and global expression analysis of the PK superfamily, however, is currently lacking. In this study, with the aim of providing novel inferences about the mechanisms associated with the stress response developed by PKs and retained throughout evolution, we identified and characterized the entire set of PKs, also known as the kinome, present in the Hevea genome. Different RNA-sequencing datasets were employed to identify tissue-specific expression patterns and potential correspondences between different rubber tree genotypes. In addition, coexpression networks under several abiotic stress conditions, such as cold, drought and latex overexploitation, were employed to elucidate associations between families and tissues/stresses. A total of 1,809 PK genes were identified using the current reference genome assembly at the scaffold level, and 1,379 PK genes were identified using the latest chromosome-level assembly and combined into a single set of 2,842 PKs. These proteins were further classified into 20 different groups and 122 families, exhibiting high compositional similarities among family members and with two phylogenetically close species Manihot esculenta and Ricinus communis. Through the joint investigation of tandemly duplicated kinases, transposable elements, gene expression patterns, and coexpression events, we provided insights into the understanding of the cell regulation mechanisms in response to several conditions, which can often lead to a significant reduction in rubber yield.

Improving Deproteinization in Colombian Latex from Hevea brasiliensis: A Bibliometric Approximation

Natural Rubber Field Latex (NRFL) allergens restrict its use in some markets due to health-threatening allergic reactions. These molecules are proteins that are related to asymptomatic sensitization and hypersensitivity mediated by immunoglobulin E (IgE). Although NRFL allergens have been investigated since the 1980s, there are still gaps in knowledge regarding the development of deproteinized natural rubber (DPNR). Therefore, in this study, the deproteinization of NRFL from the lower basin of the Cauca River, Antioquia-Colombia was evaluated using eight systems. The highest removal value was 84.4% and was obtained from the treatment containing SDS (Sodium dodecyl sulfate), Urea, and Ethanol. It was also possible to determine that at high concentrations of SDS, removal percentages higher than 70% are reached. On the other hand, all deproteinizing systems decreased NRFL Zeta potentials without self-coagulation, suggesting enhanced colloidal stability in DPNR latex. On the other hand, the bibliometric analysis presented technological advances in DPRN through different parameters and bibliometric networks. The analysis presented makes an important contribution from the bibliometric approach that could be positive for the development of research on DPNR.

Characterization of the complete chloroplast genome of Hevea pauciflora (Euphorbiaceae), an important wild relative of the rubber tree

Hevea pauciflora belongs to the Euphorbiaceae family, an important wild relative of the rubber tree. This study sequenced, assembled, and annotated the complete chloroplast genome of H. pauciflora. The complete chloroplast genome is 161,123 bp with a canonical quadripartite structure containing a large single-copy (LSC) region (89,109 bp), a small single-copy (SSC) region (18,376 bp), and two inverted repeat regions (IRa and IRb) (26,819 bp, each). A total of 134 genes were annotated, including 86 protein-coding genes, four pseudogenes, 36 tRNA genes, and eight rRNA genes. The 134 genes include four major groups: ‘self-replication’, ‘photosynthesis’, ‘unknown function’, and ‘others’. A phylogenetic analysis clustered H. pauciflora, H. brasiliensis, H. camargoana, and H. benthamiana into one clade, consistent with traditional taxonomy. This study provides useful data for further studies of Hevea genus and the phylogenetic relationships of Euphorbiaceae species.

Genome-Wide Identification and Comparative Analysis of WOX Genes in Four Euphorbiaceae Species and Their Expression Patterns in Jatropha curcas

The WUSCHEL-related homeobox (WOX) proteins are widely distributed in plants and play important regulatory roles in growth and development processes such as embryonic development and organ development. Here, series of bioinformatics methods were utilized to unravel the structural basis and genetic hierarchy of WOX genes, followed by regulation of the WOX genes in four Euphorbiaceae species. A genome-wide survey identified 59 WOX genes in Hevea brasiliensis (H. brasiliensis: 20 genes), Jatropha curcas (J. curcas: 10 genes), Manihot esculenta (M. esculenta: 18 genes), and Ricinus communis (R. communis: 11 genes). The phylogenetic analysis revealed that these WOX members could be clustered into three close proximal clades, such as namely ancient, intermediate and modern/WUS clades. In addition, gene structures and conserved motif analyses further validated that the WOX genes were conserved within each phylogenetic clade. These results suggested the relationships among WOX members in the four Euphorbiaceae species. We found that WOX genes in H. brasiliensis and M. esculenta exhibit close genetic relationship with J. curcas and R. communis. Additionally, the presence of various cis-acting regulatory elements in the promoter of J. curcas WOX genes (JcWOXs) reflected distinct functions. These speculations were further validated with the differential expression profiles of various JcWOXs in seeds, reflecting the importance of two JcWOX genes (JcWOX6 and JcWOX13) during plant growth and development. Our quantitative real-time PCR (qRT-PCR) analysis demonstrated that the JcWOX11 gene plays an indispensable role in regulating plant callus. Taken together, the present study reports the comprehensive characteristics and relationships of WOX genes in four Euphorbiaceae species, providing new insights into their characterization.

Identification and expression assay of calcium-dependent protein kinase family genes in Hevea brasiliensis and determination of HbCDPK5 functions in disease resistance

Calcium (Ca2+) signaling is one of the earliest factors to coordinate plant adaptive responses. As direct sensors and activators of Ca2+ signals, calcium-dependent protein kinases (CDPKs) were reported to be widely involved in regulating different biotic and abiotic stress stimuli. In this study, 32 Hevea brasiliensis CDPK (HbCDPK) genes were predicted and classified into four subgroups. Among them, the full-length coding sequences of 28 HbCDPK genes were confirmed by RT-PCR and verified by sequencing. Putative cis-elements assay in the promoters of HbCDPKs showed that most of the HbCDPK genes contained gibberellic acid-responsive element (GARE), abscisic acid-responsive element (ABRE), salicylic acid-responsive element (SARE), defense and stress responsive element (TC-rich repeats) and low-temperature response element (LTR), which could be activated by different biotic and abiotic stresses. Real-time PCR analysis indicated that 28 HbCDPK genes respond to infection of pathogenic fungi and a variety of phytohormones. Subcellular localization was observed with most HbCDPKs located in cell membrane, cytoplasm or organelles. Some HbCDPKs were confirmed to cause reactive oxygen species (ROS) production and accumulation in rubber tree mesophyll protoplast directly. HbCDPK5 was strongly induced by the inoculation with Colletotrichum gloeosporioides and was chosen for further analysis. HbCDPK5 localized to the cell membrane and cytoplasm, and obviously regulated the accumulation of ROS in rubber tree mesophyll protoplast. Overexpression of HbCDPK5 in Arabidopsis enhanced the resistance to Botrytis cinerea. These results indicate that rubber tree CDPK genes play important roles in plant disease resistance.

Genomic insights into lineage-specific evolution of the oleosin family in Euphorbiaceae

BACKGROUND

Lipid droplets (LDs) present in land plants serve as an essential energy and carbon reserve for seed germination and seedling development. Oleosins, the most abundant structural proteins of LDs, comprise a small family involved in LD formation, stabilization and degradation. Despite their importance, our knowledge on oleosins is still poor in Euphorbiaceae, a large plant family that contains several important oil-bearing species.

RESULTS

To uncover lineage-specific evolution of oleosin genes in Euphorbiaceae, in this study, we performed a genome-wide identification and comprehensive comparison of the oleosin family in Euphorbiaceae species with available genome sequences, i.e. castor bean (Ricinus communis), physic nut (Jatropha curcas), tung tree (Vernicia fordii), Mercurialis annua, cassava (Manihot esculenta) and rubber tree (Hevea brasiliensis), and a number of five, five, five, five, eight and eight members were found, respectively. Synteny analysis revealed one-to-one collinear relationship of oleosin genes between the former four (i.e. castor bean, physic nut, tung tree and M. annua) as well as latter two species (i.e. cassava and rubber tree), whereas one-to-one and one-to-two collinear relationships were observed between physic nut and cassava, reflecting the occurrence of one recent whole-genome duplication (WGD) in the last common ancestor of cassava and rubber tree. The presence of five ortholog groups representing three previously defined clades (i.e. U, SL and SH) dates back at least to the Malpighiales ancestor, because they are also conserved in poplar (Populus trichocarpa), a tree having experienced one Salicaceae-specific recent WGD. As observed in poplar, WGD was shown to be the main driver for the family expansion in both cassava and rubber tree. Nevertheless, same retention patterns of WGD-derived duplicates observed in cassava and rubber tree are somewhat different from that of poplar, though certain homologous fragments are still present in rubber tree. Further transcriptional profiling revealed an apparent seed-predominant expression pattern of oleosin genes in physic nut, castor bean and rubber tree. Moreover, structure and expression divergence of paralogous pairs were also observed in both cassava and rubber tree.

CONCLUSION

Comparative genomics analysis of oleosin genes reported in this study improved our knowledge on lineage-specific family evolution in Euphorbiaceae, which also provides valuable information for further functional analysis and utilization of key members and their promoters.

Targeted designing functional markers revealed the role of retrotransposon derived miRNAs as mobile epigenetic regulators in adaptation responses of pistachio

We developed novel miRNA-based markers based on salt responsive miRNA sequences to detect polymorphisms in miRNA sequences and locations. The validation of 76 combined miRNA + miRNA and miRNA + ISSR markers in the three extreme pistachio populations led to the identification of three selected markers that could link salt tolerance phenotype to genotype and divided pistachio genotypes and Pistacia species into three clusters. This novel functional marker system, in addition to more efficient performance, has higher polymorphisms than previous miRNA-based marker systems. The functional importance of the target gene of five miRNAs in the structure of the three selected markers in regulation of different genes such as ECA2, ALA10, PFK, PHT1;4, PTR3, KUP2, GRAS, TCP, bHLH, PHD finger, PLATZ and genes involved in developmental, signaling and biosynthetic processes shows that the polymorphism associated with these selected miRNAs can make a significant phenotypic difference between salt sensitive and tolerant pistachio genotypes. The sequencing results of selected bands showed the presence of conserved miRNAs in the structure of the mitochondrial genome. Further notable findings of this study are that the sequences of PCR products of two selected markers were annotated as Gypsy and Copia retrotransposable elements. The transposition of retrotransposons with related miRNAs by increasing the number of miRNA copies and changing their location between nuclear and organellar genomes can affect the regulatory activity of these molecules. These findings show the crucial role of retrotransposon-derived miRNAs as mobile epigenetic regulators between intracellular genomes in regulating salt stress responses as well as creating new and tolerant phenotypes for adaptation to environmental conditions.

Identification of Candidate Gene-Based Markers for Girth Growth in Rubber Trees

Girth growth is an important factor in both latex and timber production of the rubber tree. In this study, we performed candidate gene association mapping for girth growth in rubber trees using intron length polymorphism markers (ILP) in identifying the candidate genes responsible for girth growth. The COBL064_1 marker developed from the candidate gene (COBL4) regulating cellulose deposition and oriented cell expansion in the plant cell wall showed the strongest association with girth growth across two seasons in the Amazonian population and was validated in the breeding lines. We then applied single molecule real-time (SMRT) circular consensus sequencing (CCS) to analyze a wider gene region of the COBL4 to pinpoint the single nucleotide polymorphism (SNP) that best explains the association with the traits. A SNP in the 3’ UTR showing linkage disequilibrium with the COBL064_1 most associated with girth growth. This study showed that the cost-effective method of ILP gene-based markers can assist in identification of SNPs in the candidate gene associated with girth growth. The SNP markers identified in this study added useful markers for the improvement of girth growth in rubber tree breeding programs.

K-mer-based machine learning method to classify LTR-retrotransposons in plant genomes

Every day more plant genomes are available in public databases and additional massive sequencing projects (i.e., that aim to sequence thousands of individuals) are formulated and released. Nevertheless, there are not enough automatic tools to analyze this large amount of genomic information. LTR retrotransposons are the most frequent repetitive sequences in plant genomes; however, their detection and classification are commonly performed using semi-automatic and time-consuming programs. Despite the availability of several bioinformatic tools that follow different approaches to detect and classify them, none of these tools can individually obtain accurate results. Here, we used Machine Learning algorithms based on k-mer counts to classify LTR retrotransposons from other genomic sequences and into lineages/families with an F1-Score of 95%, contributing to develop a free-alignment and automatic method to analyze these sequences.

InpactorDB: A Classified Lineage-Level Plant LTR Retrotransposon Reference Library for Free-Alignment Methods Based on Machine Learning

Long terminal repeat (LTR) retrotransposons are mobile elements that constitute the major fraction of most plant genomes. The identification and annotation of these elements via bioinformatics approaches represent a major challenge in the era of massive plant genome sequencing. In addition to their involvement in genome size variation, LTR retrotransposons are also associated with the function and structure of different chromosomal regions and can alter the function of coding regions, among others. Several sequence databases of plant LTR retrotransposons are available for public access, such as PGSB and RepetDB, or restricted access such as Repbase. Although these databases are useful to identify LTR-RTs in new genomes by similarity, the elements of these databases are not fully classified to the lineage (also called family) level. Here, we present InpactorDB, a semi-curated dataset composed of 130,439 elements from 195 plant genomes (belonging to 108 plant species) classified to the lineage level. This dataset has been used to train two deep neural networks (i.e., one fully connected and one convolutional) for the rapid classification of these elements. In lineage-level classification approaches, we obtain up to 98% performance, indicated by the F1-score, precision and recall scores.

Hevea brasiliensis latex proteomics: a review of analytical methods and the way forward

Natural rubber or latex from the Hevea brasiliensis is an important commodity in various economic sectors in today's modern society. Proteins have been detected in latex since the early twentieth century, and they are known to regulate various biological pathways within the H. brasiliensis trees such as the natural rubber biosynthesis, defence against pathogens, wound healing, and stress tolerance. However, the exact mechanisms of the pathways are still not clear. Proteomic analyses on latex have found various proteins and revealed how they fit into the mechanisms of the biological pathways. In the past three decades, there has been rapid latex protein identification due to the improvement of latex protein extraction methods, as well as the emergence of two-dimensional gel electrophoresis (2-DE) and mass spectrometry (MS). In this manuscript, we reviewed the methods of latex protein extraction that keeps on improving over the past three decades as well as the results of numerous latex protein identification and quantitation.

Supernumerary B chromosomes of Aegilops speltoides undergo precise elimination in roots early in embryo development

Not necessarily all cells of an organism contain the same genome. Some eukaryotes exhibit dramatic differences between cells of different organs, resulting from programmed elimination of chromosomes or their fragments. Here, we present a detailed analysis of programmed B chromosome elimination in plants. Using goatgrass Aegilops speltoides as a model, we demonstrate that the elimination of B chromosomes is a strictly controlled and highly efficient root-specific process. At the onset of embryo differentiation B chromosomes undergo elimination in proto-root cells. Independent of centromere activity, B chromosomes demonstrate nondisjunction of chromatids and lagging in anaphase, leading to micronucleation. Chromatin structure and DNA replication differ between micronuclei and primary nuclei and degradation of micronucleated DNA is the final step of B chromosome elimination. This process might allow root tissues to survive the detrimental expression, or overexpression of B chromosome-located root-specific genes with paralogs located on standard chromosomes.

B chromosomes are supernumerary chromosomes exhibiting dramatic differences between different organs in same species. Here, the authors show programmed B chromosome elimination in goatgrass starts at the onset of embryo differentiation by nondisjunction of chromatids, anaphase lagging, and ends with the degradation of micronucleated DNA.

Complete chloroplast genome of Hevea benthamiana, a SALB-resistant relative wild species of rubber tree

Hevea benthamiana is a SALB-resistant wild species of H. brasiliensis, the only source of mass production of high quality natural rubber. This study sequenced and analyzed the chloroplast genome of H. benthamiana. The chloroplast genome of H. benthamiana contains 161,124 bp and consists of 51,495 bp of A (31.96%), 52,022 bp of T (32.29%), 28,915 bp of G (17.95%), and 28,692 bp of C (17.81%). The ring-shaped genome includes four regions: a large single-copy region (LSC, 89,110 bp), a small single copy (SSC, 18,376 bp) region, and two inverted repeat regions (IRs, 26,819 bp). A total of 134 genes were annotated, of which 86 encode proteins; four are pseudogenes; 36 are tRNA genes, and eight are rRNA genes. Phylogenetic analyses showed that H. benthamiana is very closely related to H. Brasiliensis, this result indicates that H. benthamiana is highly valuable for the breeding of SALB-resistant varieties of rubber trees.

Molecular Mechanisms of Natural Rubber Biosynthesis

Natural rubber (NR), principally comprising cis-1,4-polyisoprene, is an industrially important natural hydrocarbon polymer because of its unique physical properties, which render it suitable for manufacturing items such as tires. Presently, industrial NR production depends solely on latex obtained from the Pará rubber tree, Hevea brasiliensis. In latex, NR is enclosed in rubber particles, which are specialized organelles comprising a hydrophobic NR core surrounded by a lipid monolayer and membrane-bound proteins. The similarity of the basic carbon skeleton structure between NR and dolichols and polyprenols, which are found in most organisms, suggests that the NR biosynthetic pathway is related to the polyisoprenoid biosynthetic pathway and that rubber transferase, which is the key enzyme in NR biosynthesis, belongs to the cis-prenyltransferase family. Here, we review recent progress in the elucidation of molecular mechanisms underlying NR biosynthesis through the identification of the enzymes that are responsible for the formation of the NR backbone structure.

Genome sequencing of Rigidoporus microporus provides insights on genes important for wood decay, latex tolerance and interspecific fungal interactions

Fungal plant pathogens remain a serious threat to the sustainable agriculture and forestry, despite the extensive efforts undertaken to control their spread. White root rot disease is threatening rubber tree (Hevea brasiliensis) plantations throughout South and Southeast Asia and Western Africa, causing tree mortality and severe yield losses. Here, we report the complete genome sequence of the basidiomycete fungus Rigidoporus microporus, a causative agent of the disease. Our phylogenetic analysis confirmed the position of R. microporus among the members of Hymenochaetales, an understudied group of basidiomycetes. Our analysis further identified pathogen's genes with a predicted role in the decay of plant cell wall polymers, in the utilization of latex components and in interspecific interactions between the pathogen and other fungi. We also detected putative horizontal gene transfer events in the genome of R. microporus. The reported first genome sequence of a tropical rubber tree pathogen R. microporus should contribute to the better understanding of how the fungus is able to facilitate wood decay and nutrient cycling as well as tolerate latex and utilize resinous extractives.

Complete Genome Insights into Lactococcus petauri CF11 Isolated from a Healthy Human Gut Using Second- and Third-Generation Sequencing

Lactococcus petauri CF11 was originally isolated from the gut of healthy humans. To determine the underlying molecular and genetic mechanisms of the probiotic potential of CF11, we performed complete genome sequencing, annotation, and comparative genome analysis. The complete genome of L. petauri CF11 comprised of 1,997,720 bp, with a DNA G+C content of 38.21 mol% containing 1982 protein coding genes and 16 rRNA operons. We found that 1206 genes (56.05%) were assigned a putative function using the gene ontology (GO) resource. The gene products of CF11 were primarily concentrated in molecular function and biological processes, such as catalysis, binding, metabolism, and cellular processes. Furthermore, 1,365 (68.87%) genes were assigned an illative function using COGs. CF11 proteins were associated with carbohydrate transport and metabolism, and amino acid transport and metabolism. This indicates that CF11 bacteria can perform active energy exchange. We classified 1,111 (56.05%) genes into six KEGG functional categories; fructose-bisphosphate aldolase and the phosphoenol pyruvate:phosphotransferase system (PTS), which are necessary in producing short-chain fatty acids (SCFAs), were excited in the carbohydrate metabolic pathway. This suggests that L. petauri CF11 produces SCFAs via glycolysis. The genomic island revealed that some regions contain fragments of antibiotic resistance and bacteriostatic genes. In addition, ANI analysis showed that L. petauri CF11 had the closest relationship with L. petauri 159469^T, with an average nucleotide consistency of 98.03%. Taken together, the present study offers further insights into the functional and potential role of L. petauri CF11 in health care.

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.

The complete chloroplast genome of Hevea camargoana

Hevea camargoana is a natural latex producing tropical plant and a close relative of H. brasiliensis, the primary commercial source of natural rubber. This study sequenced and analyzed the chloroplast genome of H. camargoana. The circular chloroplast genome of H. camargoana contains 161,291 bp with a GC content of 35.72%. This region contains two inverted repeat regions (26,819 bp), a large single-copy region (89,281 bp), and a small single-copy (18,372 bp) region in the complete chloroplast genome. A total of 134 genes were annotated, including 86 protein-coding genes, 36 transfer RNA genes, 8 ribosomal RNA genes, and 4 pseudogenes. The results showed that H. camargoana and H. brasiliensis were closely related, suggesting that H. camargoana may be used for the future variety improvement of rubber trees.

A systematic review of the application of machine learning in the detection and classification of transposable elements

Background

Transposable elements (TEs) constitute the most common repeated sequences in eukaryotic genomes. Recent studies demonstrated their deep impact on species diversity, adaptation to the environment and diseases. Although there are many conventional bioinformatics algorithms for detecting and classifying TEs, none have achieved reliable results on different types of TEs. Machine learning (ML) techniques can automatically extract hidden patterns and novel information from labeled or non-labeled data and have been applied to solving several scientific problems.

Methodology

We followed the Systematic Literature Review (SLR) process, applying the six stages of the review protocol from it, but added a previous stage, which aims to detect the need for a review. Then search equations were formulated and executed in several literature databases. Relevant publications were scanned and used to extract evidence to answer research questions.

Results

Several ML approaches have already been tested on other bioinformatics problems with promising results, yet there are few algorithms and architectures available in literature focused specifically on TEs, despite representing the majority of the nuclear DNA of many organisms. Only 35 articles were found and categorized as relevant in TE or related fields.

Conclusions

ML is a powerful tool that can be used to address many problems. Although ML techniques have been used widely in other biological tasks, their utilization in TE analyses is still limited. Following the SLR, it was possible to notice that the use of ML for TE analyses (detection and classification) is an open problem, and this new field of research is growing in interest.

Genome-Wide Identification and Characterization of Heat-Shock Transcription Factors in Rubber Tree

Heat-shock transcription factors (Hsfs) play a pivotal role in the response of plants to various stresses. The present study aimed to characterize the Hsf genes in the rubber tree, a primary global source of natural rubber. In this study, 30 Hsf genes were identified in the rubber tree using genome-wide analysis. They possessed a structurally conserved DNA-binding domain and an oligomerization domain. On the basis of the length of the insert region between HR-A and HR-B in the oligomerization domain, the 30 members were clustered into three classes, Classes A (18), B (10), and C (2). Members within the same class shared highly conserved gene structures and protein motifs. The background expression levels of 11 genes in cold-tolerant rubber-tree clone 93-14 were significantly higher than those in cold-sensitive rubber-tree clone Reken501, while four genes exhibited inverse expression patterns. Upon cold stress, 20 genes were significantly upregulated in 93-114. Of the upregulated genes, HbHsfA2b, HbHsfA3a, and HbHsfA7a were also significantly upregulated in three other cold-tolerant rubber-tree clones at one or more time intervals upon cold stress. Their nuclear localization was verified, and the protein–protein interaction network was predicted. This study provides a basis for dissecting Hsf function in the enhanced cold tolerance of the rubber tree.

Hybrid de novo transcriptome assembly of poinsettia (Euphorbia pulcherrima Willd. Ex Klotsch) bracts

BACKGROUND

Poinsettia is a popular and important ornamental crop, mostly during the Christmas season. Its bract coloration ranges from pink/red to creamy/white shades. Despite its ornamental value, there is a lack of knowledge about the genetics and molecular biology of poinsettia, especially on the mechanisms of color formation. We performed an RNA-Seq analysis in order to shed light on the transcriptome of poinsettia bracts. Moreover, we analyzed the transcriptome differences of red- and white-bracted poinsettia varieties during bract development and coloration. For the assembly of a bract transcriptome, two paired-end cDNA libraries from a red and white poinsettia pair were sequenced with the Illumina technology, and one library from a red-bracted variety was used for PacBio sequencing. Both short and long reads were assembled using a hybrid de novo strategy. Samples of red- and white-bracted poinsettias were sequenced and comparatively analyzed in three color developmental stages in order to understand the mechanisms of color formation and accumulation in the species.

RESULTS

The final transcriptome contains 288,524 contigs, with 33% showing confident protein annotation against the TAIR10 database. The BUSCO pipeline, which is based on near-universal orthologous gene groups, was applied to assess the transcriptome completeness. From a total of 1440 BUSCO groups searched, 77% were categorized as complete (41% as single-copy and 36% as duplicated), 10% as fragmented and 13% as missing BUSCOs. The gene expression comparison between red and white varieties of poinsettia showed a differential regulation of the flavonoid biosynthesis pathway only at particular stages of bract development. An initial impairment of the flavonoid pathway early in the color accumulation process for the white poinsettia variety was observed, but these differences were no longer present in the subsequent stages of bract development. Nonetheless, GSTF11 and UGT79B10 showed a lower expression in the last stage of bract development for the white variety and, therefore, are potential candidates for further studies on poinsettia coloration.

CONCLUSIONS

In summary, this transcriptome analysis provides a valuable foundation for further studies on poinsettia, such as plant breeding and genetics, and highlights crucial information on the molecular mechanism of color formation.

Single-nucleotide polymorphism markers within MVA and MEP pathways among Hevea brasiliensis clones through transcriptomic analysis

In this study, RNA sequencing of several Hevea brasiliensis clones grown in Malaysia with different annual rubber production yields and disease resistance was performed on the Illumina platform. A total of 29,862,548 reads were generated, resulting in 101,269 assembled transcripts that were used as the reference transcripts. A similarity search against the non-redundant (nr) protein databases presented 83,771 (83%) positive BLASTx hits. The transcriptome was annotated using gene ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG) and the Pfam database. A search for putative molecular markers was performed to identify single-nucleotide polymorphisms (SNPs). Overall, 3,210,629 SNPs were detected and a total of 1314 SNPs associated with the genes involved in MVA and MEP pathways were identified. A total of 176 SNP primer pairs were designed from sequences that were related to the MVA and MEP pathways. The transcriptome of RRIM 3001 and RRIM 712 were subjected to pairwise comparison and the results revealed that there were 1262 significantly differentially expressed genes unique to RRIM 3001, 1499 significantly differentially expressed genes unique to RRIM 712 and several genes related to the MVA and MEP pathways such as AACT, HMGS, PMK, MVD, DXS and HDS were included. The results will facilitate the characterization of H. brasiliensis transcriptomes and the development of a new set of molecular markers in the form of SNPs from transcriptome assembly for the genotype identification of various rubber varieties with superior traits in Malaysia.

Proteomic Landscape Has Revealed Small Rubber Particles Are Crucial Rubber Biosynthetic Machines for Ethylene-Stimulation in Natural Rubber Production

Rubber particles are a specific organelle for natural rubber biosynthesis (NRB) and storage. Ethylene can significantly improve rubber latex production by increasing the generation of small rubber particles (SRPs), regulating protein accumulation, and activating many enzyme activities. We conducted a quantitative proteomics study of different SRPs upon ethylene stimulation by differential in-gel electrophoresis (DIGE) and using isobaric tags for relative and absolute quantification (iTRAQ) methods. In DIGE, 79 differentially accumulated proteins (DAPs) were determined as ethylene responsive proteins. Our results show that the abundance of many NRB-related proteins has been sharply induced upon ethylene stimulation. Among them, 23 proteins were identified as rubber elongation factor (REF) and small rubber particle protein (SRPP) family members, including 16 REF and 7 SRPP isoforms. Then, 138 unique phosphorylated peptides, containing 129 phosphorylated amino acids from the 64 REF/SRPP family members, were identified, and most serine and threonine were phosphorylated. Furthermore, we identified 226 DAPs from more than 2000 SRP proteins by iTRAQ. Integrative analysis revealed that almost all NRB-related proteins can be detected in SRPs, and many proteins are positively responsive to ethylene stimulation. These results indicate that ethylene may stimulate latex production by regulating the accumulation of some key proteins. The phosphorylation modification of REF and SRPP isoforms might be crucial for NRB, and SRP may act as a complex natural rubber biosynthetic machine.

Identification and characterization of jasmonic acid- and linolenic acid-mediated transcriptional regulation of secondary laticifer differentiation in Hevea brasiliensis

Hevea brasiliensis remains the primary crop commercially exploited to obtain latex, which is produced from the articulated secondary laticifer. Here, we described the transcriptional events related to jasmonic acid (JA)- and linolenic acid (LA)-induced secondary laticifer differentiation (SLD) in H. brasiliensis clone RRIM 600 based on RNA-seq approach. Histochemical approach proved that JA- and LA-treated samples resulted in SLD in H. brasiliensis when compared to ethephon and untreated control. RNA-seq data resulted in 86,614 unigenes, of which 2,664 genes were differentially expressed in JA and LA-induced secondary laticifer harvested from H. brasiliensis bark samples. Among these, 450 genes were unique to JA and LA as they were not differentially expressed in ethephon-treated samples compared with the untreated samples. Most transcription factors from the JA- and LA-specific dataset were classified under MYB, APETALA2/ethylene response factor (AP2/ERF), and basic-helix-loop-helix (bHLH) gene families that were involved in tissue developmental pathways, and we proposed that Bel5-GA2 oxidase 1-KNOTTED-like homeobox complex are likely involved in JA- and LA-induced SLD in H. brasiliensis. We also discovered alternative spliced transcripts, putative novel transcripts, and cis-natural antisense transcript pairs related to SLD event. This study has advanced understanding on the transcriptional regulatory network of SLD in H. brasiliensis.

A draft genome and transcriptome of common milkweed (Asclepias syriaca) as resources for evolutionary, ecological, and molecular studies in milkweeds and Apocynaceae

Milkweeds (Asclepias) are used in wide-ranging studies including floral development, pollination biology, plant-insect interactions and co-evolution, secondary metabolite chemistry, and rapid diversification. We present a transcriptome and draft nuclear genome assembly of the common milkweed, Asclepias syriaca. This reconstruction of the nuclear genome is augmented by linkage group information, adding to existing chloroplast and mitochondrial genomic resources for this member of the Apocynaceae subfamily Asclepiadoideae. The genome was sequenced to 80.4× depth and the draft assembly contains 54,266 scaffolds ≥1 kbp, with N50 = 3,415 bp, representing 37% (156.6 Mbp) of the estimated 420 Mbp genome. A total of 14,474 protein-coding genes were identified based on transcript evidence, closely related proteins, and ab initio models, and 95% of genes were annotated. A large proportion of gene space is represented in the assembly, with 96.7% of Asclepias transcripts, 88.4% of transcripts from the related genus Calotropis, and 90.6% of proteins from Coffea mapping to the assembly. Scaffolds covering 75 Mbp of the Asclepias assembly formed 11 linkage groups. Comparisons of these groups with pseudochromosomes in Coffea found that six chromosomes show consistent stability in gene content, while one may have a long history of fragmentation and rearrangement. The progesterone 5β-reductase gene family, a key component of cardenolide production, is likely reduced in Asclepias relative to other Apocynaceae. The genome and transcriptome of common milkweed provide a rich resource for future studies of the ecology and evolution of a charismatic plant family.

Identification and characterization of the MADS-box genes highly expressed in the laticifer cells of Hevea brasiliensis

MADS-box transcription factors possess many functions in plant reproduction and development. However, few MADS-box genes related to secondary metabolites regulation have been identified. In Hevea brasiliensis, natural rubber is a representative cis-polyisoprenoids in secondary metabolism which occurs in the rubber laticifer cells, the molecular regulation basis of natural rubber biosynthesis is not clear. Here, a total of 24 MADS-box genes including 4 type I MADS-box genes and 20 type II MADS-box genes were identified in the transcriptome of rubber tree latex. The phylogenetic analysis was performed to clarify the evolutionary relationships of all the 24 rubber tree MADS-box proteins with MADS-box transcription factors from Arabidopsis thaliana and Oryza sativa. Four type I MADS-box genes were subdivided into Mα (3 genes) and Mβ (1 gene). Twenty type II MADS-box genes were subclassified into MIKC* (8 genes) and MIKC^c (12 genes). Eight MADS-box genes (HblMADS3, 5, 6, 7, 9, 13, 23, 24) were predominant expression in laticifers. ABA up-regulated the expression of HblMADS9, and the expression of HblMADS3, HblMADS5, HblMADS24 were up-regulated by MeJA. The function of HblMADS24 was elucidated. HblMADS24 bound HbFPS1 promoter in yeast and HblMADS24 activated HbFPS1 promoter in tobacco plants. Moreover, we proposed that HblMADS24 is a transcription activator of HbFPS1 which taking part in natural rubber biosynthesis.

Identification of histone methylation modifiers and their expression patterns during somatic embryogenesis in Hevea brasiliensis

Histone methylation plays a crucial role in various biological processes, from heterochromatin formation to transcriptional regulation. Currently, no information is available regarding histone methylation modifiers in the important rubber-producing plant Hevea brasiliensis. Here, we identified 47 histone methyltransferase (HMT) genes and 25 histone demethylase (HDM) genes as possible members of the histone methylation modifiers in the rubber tree genome. According to the structural features of HMT and HDM, the HbHMTs were classified into two groups (HbPRMs and HbSDGs), the HbHDMs have two groups (HbLSDs and HbJMJs). Expression patterns were analyzed in five different tissues and at different phases of somatic embryogenesis. HbSDG10, 21, 25, 33, HbJMJ2, 18, 20 were with high expression at different phases of somatic embryogenesis. HbSDG10,14, 20, 21, 33 and HbPRMT4 were expressed highly in anther, HbSDG14, 20, 21, 22, 23, 33, 35 and HbPRMT1 HbJMJ7 and HbLSD1, 2, 3, 4 showed high expression levels in callus. HbSDG1, 7, 10, 13, 14, 18, 19, 21, 22, 23, 35, HbPRMT1, 8, HbJMJ5, 7, 11, 16, 20 and HbLSD2, 3, 4 were expressed highly in somatic embryo. HbSDG10, 21, 25, 33, HbLSD2, 3 were expressed highly in bud of regenerated plant. The analyses reveal that HbHMTs and HbHDMs exhibit different expression patterns at different phases during somatic embryogenesis, implying that some HbHMTs and HbHDMs play important roles during somatic embryogenesis. This study provide fundamental information for further studies on histone methylation in Hevea brasiliensis.

Retrotransposons in Plant Genomes: Structure, Identification, and Classification through Bioinformatics and Machine Learning

Transposable elements (TEs) are genomic units able to move within the genome of virtually all organisms. Due to their natural repetitive numbers and their high structural diversity, the identification and classification of TEs remain a challenge in sequenced genomes. Although TEs were initially regarded as "junk DNA", it has been demonstrated that they play key roles in chromosome structures, gene expression, and regulation, as well as adaptation and evolution. A highly reliable annotation of these elements is, therefore, crucial to better understand genome functions and their evolution. To date, much bioinformatics software has been developed to address TE detection and classification processes, but many problematic aspects remain, such as the reliability, precision, and speed of the analyses. Machine learning and deep learning are algorithms that can make automatic predictions and decisions in a wide variety of scientific applications. They have been tested in bioinformatics and, more specifically for TEs, classification with encouraging results. In this review, we will discuss important aspects of TEs, such as their structure, importance in the evolution and architecture of the host, and their current classifications and nomenclatures. We will also address current methods and their limitations in identifying and classifying TEs.

Early Sex-Chromosome Evolution in the Diploid Dioecious Plant Mercurialis annua

Suppressed recombination allows divergence between homologous sex chromosomes and the functionality of their genes. Here, we reveal patterns of the earliest stages of sex-chromosome evolution in the diploid dioecious herb Mercurialis annua on the basis of cytological analysis, de novo genome assembly and annotation, genetic mapping, exome resequencing of natural populations, and transcriptome analysis. The genome assembly contained 34,105 expressed genes, of which 10,076 were assigned to linkage groups. Genetic mapping and exome resequencing of individuals across the species range both identified the largest linkage group, LG1, as the sex chromosome. Although the sex chromosomes of M. annua are karyotypically homomorphic, we estimate that about one-third of the Y chromosome, containing 568 transcripts and spanning 22.3 cM in the corresponding female map, has ceased recombining. Nevertheless, we found limited evidence for Y-chromosome degeneration in terms of gene loss and pseudogenization, and most X- and Y-linked genes appear to have diverged in the period subsequent to speciation between M. annua and its sister species M. huetii, which shares the same sex-determining region. Taken together, our results suggest that the M. annua Y chromosome has at least two evolutionary strata: a small old stratum shared with M. huetii, and a more recent larger stratum that is probably unique to M. annua and that stopped recombining ∼1 MYA. Patterns of gene expression within the nonrecombining region are consistent with the idea that sexually antagonistic selection may have played a role in favoring suppressed recombination.

Genomic technologies for Hevea breeding

The commercial production of high quality natural rubber (NR) solely depends on Hevea brasiliensis Muell. Arg, (Para rubber tree) and accounts for >98% of total production worldwide. NR with its unique properties is an essential commodity for the automobile industry and its synthetic counterparts are in no way substitute to it. The rubber tree genome is very complex and plays an important role in delivering the unique properties of Hevea. But a lack of knowledge on the molecular mechanisms of rubber biosynthesis, disease resistance, etc., in elite clones of rubber still persists. Marker-assisted selection and transgenic techniques were proved to be advantageous in improving the breeding efficiency for latex yield, disease resistance, etc. The suppression subtractive hybridization (SSH), in the form of subtracted cDNA libraries and microarrays, can assist in searching the functions of expressed genes (candidate gene approach). Expressed sequence tags (ESTs) related to various metabolic aspects are well utilized to create EST banks that broadly represent the genes expressed in one tissue, such as latex cells, that assists in the study of gene function and regulation. Transcriptome analysis and gene mapping have been accomplished in Hevea at various stages. However, a selection criterion to delineate high yielding genotypes at the juvenile stage has not been accomplished so far. This is the main pit fall for rubber breeding apart from stock-scion interactions leading to yield differences among a clonally multiplied population. At least four draft genome sequences have been published on Hevea rubber, and all give different genome size and contig lengths-a comprehensive and acceptable genomic map remains unfulfilled. The progress made in molecular markers, latex biosynthesis genes, transcriptome analysis, chloroplast and mitochondrial DNA diversity, paternity identification through Breeding without Breeding (BwB), stimulated latex production and its molecular intricacies, molecular biology of tapping panel dryness, genomics for changed climates and genome mapping are discussed in this review. These information can be utilized to improvise the molecular breeding programs of Hevea in future.

Genome-Wide Identification and Characterization of the JAZ Gene Family in Rubber Tree (Hevea brasiliensis)

Jasmonate signaling plays a vital role in the regulation of secondary laticifer differentiation and natural rubber biosynthesis in Hevea brasiliensis. Jasmonate ZIM-domain (JAZ) proteins are the master regulators of jasmonate signaling. Although several JAZs have been reported in the laticifer cells of H. brasiliensis, the genome-wide screening of HbJAZ members has not yet been explored. In the present study, 18 HbJAZs were identified based on the recent H. brasiliensis genome. Phylogenetic construction revealed that the HbJAZs were clustered into five subgroups and that members within the same subgroup shared highly conserved gene structures and protein motifs. Cis-element analysis of HbJAZ promoters suggested the presence of hormone, stress and development-related cis-elements. HbJAZ1.0, HbJAZ2.0, and HbJAZ5.0 interacted with CORONATINE INSENSITIVE1 (COI1) in the presence of coronatine (COR, a JA mimic). HbJAZ1.0, HbJAZ2.0, HbJAZ5.0, and HbJAZ12.0 could also interact with each other. Of the 18 HbJAZs, transcripts of 15 HbJAZs were present in the vascular cambium region except for that of HbJAZ7.0, HbJAZ8.0d, and HbJAZ13.0. Fourteen of the 15 HbJAZs were significantly up-regulated upon COR treatment. The transcripts of three genes that were absent from vascular cambium region were also absent from the latex. Among the 15 HbJAZs in the latex, the expression patterns of 13 HbJAZs were different between the tapping and ethrel treatments. Eight of the 14 COR-up-regulated HbJAZs in the vascular cambium region were also activated by tapping in latex. Of the eight tapping-activated HbJAZs, 5 HbJAZs were repressed by ethrel application. Based on the computational analyses and gene expression patterns described in this study, the HbJAZ5.0 and HbJAZ10.0b may be associated with laticifer differentiation while the HbJAZ8.0b is a negative regulator for natural rubber biosynthesis in H. brasiliensis.

Identification of the targets of HbEIN3/EILs in genomic wide in Hevea brasiliensis

EIN3/EILs are key regulators in ET signaling pathway. In this work, 4 members of EIN3/EILs of Hevea brasiliensis (HbEIN3/EILs) showed interaction with two F box proteins, HbEBF1 and HbEBF2. HbEIN3 located in nucleus and exhibited strong transcriptional activity. HbEIN3 was induced by ET treatment in C-serum, but not in B-serum of latex. HbEIN3/EILs bound to G-box cis-element. To globally search the potential targets of HbEIN3/EILs, genomic sequences of H. brasiliensis was re-annotated and an HCES (Hevea Cis-Elements Scanning) program was developed ( www.h-brasiliensis.com ). HCES scanning results showed that ET- and JA- responsive cis-elements distribute overlapping in gene promoters. 3146 genes containing G-box in promoters are potential targets of HbEIN3, including 41 genes involved in biosynthesis and drainage of latex, of which 7 rate-limiting genes of latex production were regulated by both ET and JA, suggesting that ET and JA signaling pathways coordinated the latex biosynthesis and drainage in H. brasiliensis. Abbreviations: ABRE: ABA responsive elements; bHLH: basic helix-loop-helix; COG: Orthologous Groups; DRE: dehydration response element; ERE: ethylene responsive element; ET: Ethylene; GO: Gene Ontology; HCES: Hevea Cis-Elements Scanning; JA: jasmonates; JRE: Jasmonate-responsive element; KEGG: Kyoto Encyclopedia of Genes and Genomes; NR: non-redundant database; PLACE: Plant Cis-acting Regulatory DNA Elements; qRT-PCR: quantitative real-time RT-PCR.

Multi-omics approaches for strategic improvement of stress tolerance in underutilized crop species: A climate change perspective

For several decades, researchers are working toward improving the "major" crops for better adaptability and tolerance to environmental stresses. However, little or no research attention is given toward neglected and underutilized crop species (NUCS) which hold the potential to ensure food and nutritional security among the ever-growing global population. NUCS are predominantly climate resilient, but their yield and quality are compromised due to selective breeding. In this context, the importance of omics technologies namely genomics, transcriptomics, proteomics, phenomics and ionomics in delineating the complex molecular machinery governing growth, development and stress responses of NUCS is underlined. However, gaining insights through individual omics approaches will not be sufficient to address the research questions, whereas integrating these technologies could be an effective strategy to decipher the gene function, genome structures, biological pathways, metabolic and regulatory networks underlying complex traits. Given this, the chapter enlists the importance of NUCS in food and nutritional security and provides an overview of deploying omics approaches to study the NUCS. Also, the chapter enumerates the status of crop improvement programs in NUCS and suggests implementing "integrating omics" for gaining a better understanding of crops' response to abiotic and biotic stresses.

Biosynthesis of Natural Rubber: Current State and Perspectives

Natural rubber is a kind of indispensable biopolymers with great use and strategic importance in human society. However, its production relies almost exclusively on rubber-producing plants Hevea brasiliensis, which have high requirements for growth conditions, and the mechanism of natural rubber biosynthesis remains largely unknown. In the past two decades, details of the rubber chain polymerization and proteins involved in natural rubber biosynthesis have been investigated intensively. Meanwhile, omics and other advanced biotechnologies bring new insight into rubber production and development of new rubber-producing plants. This review summarizes the achievements of the past two decades in understanding the biosynthesis of natural rubber, especially the massive information obtained from the omics analyses. Possibilities of natural rubber biosynthesis in vitro or in genetically engineered microorganisms are also discussed.

Chromosome level high-density integrated genetic maps improve the Pyrus bretschneideri 'DangshanSuli' v1.0 genome

Background

Chromosomal level reference genomes provide a crucial foundation for genomics research such as genome-wide association studies (GWAS) and whole genome selection. The chromosomal-level sequences of both the European (Pyrus communis) and Chinese (P. bretschneideri) pear genomes have not been published in public databases so far.

Results

To anchor the scaffolds of P. bretschneideri ‘DangshanSuli’ (DS) v1.0 genome into pseudo-chromosomes, two genetic maps (MH and YM maps) were constructed using half sibling populations of Chinese pear crosses, ‘Mantianhong’ (MTH) × ‘Hongxiangsu’ (HXS) and ‘Yuluxiang’ (YLX) × MTH, from 345 and 162 seedlings, respectively, which were prepared for SNP discovery using genotyping-by-sequencing (GBS) technology. The MH and YM maps, each with 17 linkage groups (LGs), were constructed from 2606 and 2489 SNP markers and spanned 1847 and 1668 cM, respectively, with average marker intervals of 0.7. The two maps were further merged with a previously published genetic map (BD) based on the cross ‘Bayuehong’ (BYH) × ‘Dangshansuli’ (DS) to build a new integrated MH-YM-BD map. By using 7757 markers located on the integrated MH-YM-BD map, 898 scaffolds (400.57 Mb) of the DS v1.0 assembly were successfully anchored into 17 pseudo-chromosomes, accounting for 78.8% of the assembled genome size. About 88.31% of them (793 scaffolds) were directionally anchored with two or more markers on the pseudo-chromosomes. Furthermore, the errors in each pseudo-chromosome (especially 1, 5, 7 and 11) were manually corrected and pseudo-chromosomes 1, 5 and 7 were extended by adding 19, 12 and 14 scaffolds respectively in the newly constructed DS v1.1 genome. Synteny analyses revealed that the DS v1.1 genome had high collinearity with the apple genome, and the homologous fragments between pseudo-chromosomes were similar to those found in previous studies. Moreover, the red-skin trait of Asian pear was mapped to an identical locus as identified previously.

Conclusions

The accuracy of DS v1.1 genome was improved by using larger mapping populations and merged genetic map. With more than 400 MB anchored to 17 pseudo-chromosomes, the new DS v1.1 genome provides a critical tool that is essential for studies of pear genetics, genomics and molecular breeding.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5224-6) contains supplementary material, which is available to authorized users.

A multiscale study of fungal endophyte communities of the foliar endosphere of native rubber trees in Eastern Amazon

Hevea brasiliensis is a native hyperdiverse tree species in the Amazon basin with great economic importance since it produces the highest quality natural rubber. H. brasiliensis, in its natural habitat, may harbor fungal endophytes that help defend against phytopathogenic fungi. In this work, we investigated the fungal endophytic communities in two pristine areas in Eastern Amazon (Anavilhanas National Park - ANP and Caxiuanã National Forest - CNF) at different spatial scales: regional, local, individual (tree), and intra-individual (leaflet). Using a culture-based approach, 210 fungal endophytes were isolated from 240 sampling units and assigned to 46 distinct MOTUs based on sequencing of the nrITS DNA. The community compositions of the endophytomes are different at both regional and local scales, dominated by very few taxa and highly skewed toward rare taxa, with many endophytes infrequently isolated across hosts in sampled space. Colletotrichum sp. 1, a probably latent pathogen, was the most abundant endophytic putative species and was obtained from all individual host trees in both study areas. Although the second most abundant putative species differed between the two collection sites, Clonostachys sp. 1 and Trichoderma sp. 1, they are phylogenetically related (Hypocreales) mycoparasites. Thus, they probably exhibit the same ecological function in the foliar endosphere of rubber tree as antagonists of its fungal pathogens.

Use of a draft genome of coffee (Coffea arabica) to identify SNPs associated with caffeine content

Arabica coffee (Coffea arabica) has a small gene pool limiting genetic improvement. Selection for caffeine content within this gene pool would be assisted by identification of the genes controlling this important trait. Sequencing of DNA bulks from 18 genotypes with extreme high- or low-caffeine content from a population of 232 genotypes was used to identify linked polymorphisms. To obtain a reference genome, a whole genome assembly of arabica coffee (variety K7) was achieved by sequencing using short read (Illumina) and long-read (PacBio) technology. Assembly was performed using a range of assembly tools resulting in 76 409 scaffolds with a scaffold N50 of 54 544 bp and a total scaffold length of 1448 Mb. Validation of the genome assembly using different tools showed high completeness of the genome. More than 99% of transcriptome sequences mapped to the C. arabica draft genome, and 89% of BUSCOs were present. The assembled genome annotated using AUGUSTUS yielded 99 829 gene models. Using the draft arabica genome as reference in mapping and variant calling allowed the detection of 1444 nonsynonymous single nucleotide polymorphisms (SNPs) associated with caffeine content. Based on Kyoto Encyclopaedia of Genes and Genomes pathway-based analysis, 65 caffeine-associated SNPs were discovered, among which 11 SNPs were associated with genes encoding enzymes involved in the conversion of substrates, which participate in the caffeine biosynthesis pathways. This analysis demonstrated the complex genetic control of this key trait in coffee.

Overview of Next-generation Sequencing Platforms Used in Published Draft Plant Genomes in Light of Genotypization of Immortelle Plant (Helichrysium Arenarium)

INTRODUCTION

Major advancements in DNA sequencing methods introduced in the first decade of the new millennium initiated a rapid expansion of sequencing studies, which yielded a tremendous amount of DNA sequence data, including whole sequenced genomes of various species, including plants. A set of novel sequencing platforms, often collectively named as "next-generation sequencing" (NGS) completely transformed the life sciences, by allowing extensive throughput, while greatly reducing the necessary time, labor and cost of any sequencing endeavor.

PURPOSE

of this paper is to present an overview NGS platforms used to produce the current compendium of published draft genomes of various plants, namely the Roche/454, ABI/SOLiD, and Solexa/Illumina, and to determine the most frequently used platform for the whole genome sequencing of plants in light of genotypization of immortelle plant.

MATERIALS AND METHODS

45 papers were selected (with 47 presented plant genome draft sequences), and utilized sequencing techniques and NGS platforms (Roche/454, ABI/SOLiD and Illumina/Solexa) in selected papers were determined. Subsequently, frequency of usage of each platform or combination of platforms was calculated.

RESULTS

Illumina/Solexa platforms are by used either as sole sequencing tool in 40.42% of published genomes, or in combination with other platforms - additional 48.94% of published genomes, followed by Roche/454 platforms, used in combination with traditional Sanger sequencing method (10.64%), and never as a sole tool. ABI/SOLiD was only used in combination with Illumina/Solexa and Roche/454 in 4.25% of publications.

CONCLUSIONS

Illumina/Solexa platforms are by far most preferred by researchers, most probably due to most affordable sequencing costs. Taking into consideration the current economic situation in the Balkans region, Illumina Solexa is the best (if not the only) platform choice if the sequencing of immortelle plant (Helichrysium arenarium) is to be performed by the researchers in this region.