De novo genome assembly and Hi-C analysis reveal an association between chromatin architecture alterations and sex differentiation in the woody plant Jatropha curcas

Stage-specific dynamic reorganization of genome topology shapes transcriptional neighborhoods in developing human retinal organoids

We have generated a high-resolution Hi-C map of developing human retinal organoids to elucidate spatiotemporal dynamics of genomic architecture and its relationship with gene expression patterns. We demonstrate progressive stage-specific alterations in DNA topology and correlate these changes with transcription of cell-type-restricted gene markers during retinal differentiation. Temporal Hi-C reveals a shift toward A compartment for protein-coding genes and B compartment for non-coding RNAs, displaying high and low expression, respectively. Notably, retina-enriched genes are clustered near lost boundaries of topologically associated domains (TADs), and higher-order assemblages (i.e., TAD cliques) localize in active chromatin regions with binding sites for eye-field transcription factors. These genes gain chromatin contacts at their transcription start site as organoid differentiation proceeds. Our study provides a global view of chromatin architecture dynamics associated with diversification of cell types during retinal development and serves as a foundational resource for in-depth functional investigations of retinal developmental traits.

Emerging trends in genomic and epigenomic regulation of plant specialised metabolism

Regulation of specialised metabolism genes is multilayered and complex, influenced by an array of genomic, epigenetic and epigenomic mechanisms. Here, we review the most recent knowledge in this field, drawing from discoveries in several plant species. Our aim is to improve understanding of how plant genome structure and function influence specialised metabolism. We also highlight key areas for future exploration. Gene regulatory mechanisms influencing specialised metabolism include gene duplication and neo-functionalization, conservation of operon-like clusters of specialised metabolism genes, local chromatin modifications, and the organisation of higher order chromatin structures within the nucleus. Genomic and epigenomic research to-date in the discipline have focused on a relatively small number of plant species, primarily at whole organ or tissue level. This is largely due to the technical demands of the experimental methods needed. However, a high degree of cell-type specificity of function exists in specialised metabolism, driven by similarly specific gene regulation. In this review we focus on the genomic characteristics of genes that are found in different types of clusters within the genome. We propose that acquisition of cell-resolution epigenomic datasets in emerging models, such as the glandular trichomes of Cannabis sativa, will yield important advances. Data such as chromatin accessibility and histone modification profiles can pinpoint which regulatory sequences are active in individual cell types and at specific times in development. These could provide fundamental biological insight as well as novel targets for genetic engineering and crop improvement.

Genome-Wide Identification and Expression Analysis of Cytokinin Response Regulator (RR) Genes in the Woody Plant Jatropha curcas and Functional Analysis of JcRR12 in Arabidopsis

The cytokinin (CK) response regulator (RR) gene family plays a pivotal role in regulating the developmental and environmental responses of plants. Axillary bud outgrowth in the perennial woody plant Jatropha curcas is regulated by the crosstalk between CK and gibberellins (GA). In this study, we first analyzed the effects of gibberellin A3 (GA3), lovastatin (a CK synthesis inhibitor), decapitation, and their interaction, on the outgrowth of axillary buds. The results indicate that lovastatin completely inhibited GA-promoted axillary bud outgrowth and partially weakened the decapitation-promoted axillary bud outgrowth. To further characterize and understand the role of CK signaling in promoting the development of female flowers and branches, we performed bioinformatics and expression analyses to characterize the CK RR gene (JcRR) family in J. curcas. A total of 14 members of the JcRR family were identified; these genes were distributed on 10 chromosomes. Phylogenetic analysis indicated that the corresponding RR proteins are evolutionarily conserved across different plant species, and the Myb-like DNA-binding domain divides the 14 members of the JcRR family into type-A and type-B proteins. Further analysis of cis-acting elements in the promoter regions of JcRRs suggests that JcRRs are expressed in response to phytohormones, light, and abiotic stress factors; thus, JcRRs may be involved in some plant development processes. Genomic sequence comparison revealed that segmental duplication may have played crucial roles in the expansion of the JcRR gene family, and five pairs of duplicated genes were all subjected to purifying selection. By analyzing RNA sequencing (RNA-seq) and quantitative reverse transcription-polymerase chain reaction (qRT–PCR) data, we characterized that the temporospatial expression patterns of JcRRs during the development of various tissues and the response of these genes to phytohormones and abiotic stress. The JcRRs were mainly expressed in the roots, while they also exhibited differential expression patterns in other tissues. The expression levels of all six type-A and one type-B JcRRs increased in response to 6-benzylaminopurine (6-BA), while the four type-B JcRRs levels decreased. The expression levels of two type-B JcRRs increased in response to exogenous GA3 treatment, while those of three type-A and three type-B JcRRs decreased. We found that type-A JcRRs may play a positive role in the continuous growth of axillary buds, while the role of type-B JcRRs might be the opposite. In response to abiotic stress, the expression levels of two type-A and three type-B JcRRs strongly increased. The overexpression of JcRR12 in Arabidopsis thaliana slightly increased the numbers of rosette branches after decapitation, but not under normal conditions. In conclusion, our results provide detailed knowledge of JcRRs for further analysis of CK signaling and JcRR functions in J. curcas.

Genome-Wide Identification and Analysis of Lbd Transcription Factor Genes in Jatropha curcas and Related Species

Lateral organ boundaries domain (LBD) proteins are plant-specific transcription factors that play important roles in organ development and stress response. However, the function of LBD genes has not been reported in Euphorbiaceae. In this paper, we used Jatropha curcas as the main study object and added rubber tree (Hevea brasiliensis), cassava (Manihot esculenta Crantz) and castor (Ricinus communis L.) to take a phylogenetic analysis of LBD genes. Of LBD, 33, 58, 54 and 30 members were identified in J. curcas, rubber tree, cassava and castor, respectively. The phylogenetic analysis showed that LBD members of Euphorbiaceae could be classified into two major classes and seven subclasses (Ia-Ie,IIa-IIb), and LBD genes of Euphorbiaceae tended to cluster in the same branch. Further analysis showed that the LBD genes of Euphorbiaceae in the same clade usually had similar protein motifs and gene structures, and tissue expression patterns showed that they also have similar expression profiles. JcLBDs in class Ia and Ie are mainly expressed in male and female flowers, and there are multiple duplication genes with similar expression profiles in these clades. It was speculated that they are likely to play important regulatory roles in flower development. Our study provided a solid foundation for further investigation of the role of LBD genes in the sexual differentiaion of J. curcas.

An ortholog of the MADS-box gene SEPALLATA3 regulates stamen development in the woody plant Jatropha curcas

Overexpression of JcSEP3 causes defective stamen development in Jatropha curcas, in which brassinosteroid and gibberellin signaling pathways may be involved. SEPALLATAs (SEPs), the class E genes of the ABCE model, are required for floral organ determination. In this study, we investigated the role of the JcSEP3 gene in floral organ development in the woody plant Jatropha curcas. Transgenic Jatropha plants overexpressing JcSEP3 displayed abnormal phenotypes such as deficient anthers and pollen, as well as free stamen filaments, whereas JcSEP3-RNA interference (RNAi) transgenic plants had no obvious phenotypic changes, suggesting that JcSEP3 is redundant with other JcSEP genes in Jatropha. Moreover, we compared the transcriptomes of wild-type plants, JcSEP3-overexpressing, and JcSEP3-RNAi transgenic plants. In the JcSEP3-overexpressing transgenic plants, we discovered 25 upregulated genes involved in anther and pollen development, as well as 12 induced genes in brassinosteroid (BR) and gibberellin (GA) signaling pathways. These results suggest that JcSEP3 directly or indirectly regulates stamen development, concomitant with the regulation of BR and GA signaling pathways. Our findings help to understand the roles of SEP genes in stamen development in perennial woody plants.

Chromosome-Scale Assembly of the Dendrobium nobile Genome Provides Insights Into the Molecular Mechanism of the Biosynthesis of the Medicinal Active Ingredient of Dendrobium

Orchids constitute approximately 10% of flowering plant species. However, only about 10 orchid genomes have been published. Metabolites are the main way through which orchids respond to their environment. Dendrobium nobile, belonging to Dendrobium, the second largest genus in Orchidaceae, has high ornamental, medicinal, and ecological value. D. nobile is the source of many popular horticultural varieties. Among the Dendrobium species, D. nobile has the highest amount of dendrobine, which is regarded as one of the criteria for evaluating medicinal quality. Due to lack of data and analysis at the genomic level, the biosynthesis pathways of dendrobine and other related medicinal ingredients in D. nobile are unknown. In this paper, we report a chromosome-scale reference genome of D. nobile to facilitate the investigation of its genomic characteristics for comparison with other Dendrobium species. The assembled genome size of D. nobile was 1.19 Gb. Of the sequences, 99.45% were anchored to 19 chromosomes. Furthermore, we identified differences in gene number and gene expression patterns compared with two other Dendrobium species by integrating whole-genome sequencing and transcriptomic analysis [e.g., genes in the polysaccharide biosynthesis pathway and upstream of the alkaloid (dendrobine) biosynthesis pathway]. Differences in the TPS and CYP450 gene families were also found among orchid species. All the above differences might contribute to the species-specific medicinal ingredient biosynthesis pathways. The metabolic pathway-related analysis will provide further insight into orchid responses to the environment. Additionally, the reference genome will provide important insights for further molecular elucidation of the medicinal active ingredients of Dendrobium and enhance the understanding of orchid evolution.

Understanding 3D Genome Organization and Its Effect on Transcriptional Gene Regulation Under Environmental Stress in Plant: A Chromatin Perspective

The genome of a eukaryotic organism is comprised of a supra-molecular complex of chromatin fibers and intricately folded three-dimensional (3D) structures. Chromosomal interactions and topological changes in response to the developmental and/or environmental stimuli affect gene expression. Chromatin architecture plays important roles in DNA replication, gene expression, and genome integrity. Higher-order chromatin organizations like chromosome territories (CTs), A/B compartments, topologically associating domains (TADs), and chromatin loops vary among cells, tissues, and species depending on the developmental stage and/or environmental conditions (4D genomics). Every chromosome occupies a separate territory in the interphase nucleus and forms the top layer of hierarchical structure (CTs) in most of the eukaryotes. While the A and B compartments are associated with active (euchromatic) and inactive (heterochromatic) chromatin, respectively, having well-defined genomic/epigenomic features, TADs are the structural units of chromatin. Chromatin architecture like TADs as well as the local interactions between promoter and regulatory elements correlates with the chromatin activity, which alters during environmental stresses due to relocalization of the architectural proteins. Moreover, chromatin looping brings the gene and regulatory elements in close proximity for interactions. The intricate relationship between nucleotide sequence and chromatin architecture requires a more comprehensive understanding to unravel the genome organization and genetic plasticity. During the last decade, advances in chromatin conformation capture techniques for unravelling 3D genome organizations have improved our understanding of genome biology. However, the recent advances, such as Hi-C and ChIA-PET, have substantially increased the resolution, throughput as well our interest in analysing genome organizations. The present review provides an overview of the historical and contemporary perspectives of chromosome conformation capture technologies, their applications in functional genomics, and the constraints in predicting 3D genome organization. We also discuss the future perspectives of understanding high-order chromatin organizations in deciphering transcriptional regulation of gene expression under environmental stress (4D genomics). These might help design the climate-smart crop to meet the ever-growing demands of food, feed, and fodder.

High-quality genome assembly of an important biodiesel plant, Euphorbia lathyris L

Caper spurge, Euphorbia lathyris L., is an important energy crop and medicinal crop. Here, we generated a high-quality, chromosome-level genome assembly of caper spurge using Oxford Nanopore sequencing, Illumina sequencing, and Hi-C technology. The final genome assembly was ∼988.9 Mb in size, 99.8% of which could be grouped into 10 pseudochromosomes, with contig and scaffold N50 values of 32.6 and 95.7 Mb, respectively. A total of 651.4 Mb repetitive sequences and 36,342 protein-coding genes were predicted in the genome assembly. Comparative genomic analysis showed that caper spurge and castor bean clustered together. We found that no independent whole-genome duplication event had occurred in caper spurge after its split from castor bean, and recent substantial amplification of LTR retrotransposons (LTR-RTs) has contributed significantly to its genome expansion. Furthermore, based on gene homology searching, we identified a number of candidate genes involved in the biosynthesis of fatty acids and triacylglycerols. The reference genome presented here will be highly useful for the further study of the genetics, genomics, and breeding of this high-value crop, as well as for evolutionary studies of spurge family and angiosperms.

Sequencing an F1 hybrid of Silurus asotus and S. meridionalis enabled the assembly of high-quality parental genomes

Genome complexity such as heterozygosity may heavily influence its de novo assembly. Sequencing somatic cells of the F1 hybrids harboring two sets of genetic materials from both of the paternal and maternal species may avoid alleles discrimination during assembly. However, the feasibility of this strategy needs further assessments. We sequenced and assembled the genome of an F1 hybrid between Silurus asotus and S. meridionalis using the SequelII platform and Hi-C scaffolding technologies. More than 300 Gb raw data were generated, and the final assembly obtained 2344 scaffolds composed of 3017 contigs. The N50 length of scaffolds and contigs was 28.55 Mb and 7.49 Mb, respectively. Based on the mapping results of short reads generated for the paternal and maternal species, each of the 29 chromosomes originating from S. asotus and S. meridionalis was recognized. We recovered nearly 94% and 96% of the total length of S. asotus and S. meridionalis. BUSCO assessments and mapping analyses suggested that both genomes had high completeness and accuracy. Further analyses demonstrated the high collinearity between S. asotus, S. meridionalis, and the related Pelteobagrus fulvidraco. Comparison of the two genomes with that assembled only using the short reads from non-hybrid parental species detected a small portion of sequences that may be incorrectly assigned to the different species. We supposed that at least part of these situations may have resulted from mitotic recombination. The strategy of sequencing the F1 hybrid genome can recover the vast majority of the parental genomes and may improve the assembly of complex genomes.

Genome-wide characterization of bZIP transcription factors and their expression patterns in response to drought and salinity stress in Jatropha curcas

The basic leucine zipper (bZIP) family is one of the largest families of transcription factors (TFs) in plants and is responsible for various functions, including regulating development and responses to abiotic/biotic stresses. However, the roles of bZIPs in the regulation of responses to drought stress and salinity stress remain poorly understood in Jatropha curcas L., a biodiesel crop. In the present study, 50 JcbZIP genes were identified and classified into ten groups. Cis-element analysis indicated that JcbZIP genes are associated with abiotic stress. Gene expression patterns and quantitative real-time PCR (qRT-PCR) showed that four JcbZIP genes (JcbZIPs 34, 36, 49 and 50) are key resistance-related genes under both drought and salinity stress conditions. On the basis of the results of cis-element and phylogenetic analyses, JcbZIP49 and JcbZIP50 are likely involved in responses to drought and salinity stress; moreover, JcbZIP34 and JcbZIP36 might also play important roles in seed development and response to abiotic stress. These findings advance our understanding of the comprehensive characteristics of JcbZIP genes and provide new insights for functional validation in the further.

High-continuity genome assembly of the jellyfish Chrysaora quinquecirrha

The Atlantic sea nettle (Chrysaora quinquecirrha) has an important evolutionary position due to its high ecological value. However, due to limited sequencing technologies and complex jellyfish genomic sequences, the current C. quinquecirrha genome assembly is highly fragmented. Here, we used the most advanced high-throughput chromosome conformation capture (Hi-C) technology to obtain high-coverage sequencing data of the C. quinquecirrha genome. We then anchored these data to the previously published contig-level assembly to improve the genome. Finally, a high-continuity genome sequence of C. quinquecirrha was successfully assembled, which contained 1 882 scaffolds with a N50 length of 3.83 Mb. The N50 length of the genome assembly was 5.23 times longer than the previously released one, and additional analysis revealed that it had a high degree of genomic continuity and accuracy. Acquisition of the high-continuity genome sequence of C. quinquecirrha not only provides a basis for the study of jellyfish evolution through comparative genomics but also provides an important resource for studies on jellyfish growth and development.

Silencing of the Ortholog of DEFECTIVE IN ANTHER DEHISCENCE 1 Gene in the Woody Perennial Jatropha curcas Alters Flower and Fruit Development

DEFECTIVE IN ANTHER DEHISCENCE 1 (DAD1), a phospholipase A1, utilizes galactolipids (18:3) to generate α-linolenic acid (ALA) in the initial step of jasmonic acid (JA) biosynthesis in Arabidopsis thaliana. In this study, we isolated the JcDAD1 gene, an ortholog of Arabidopsis DAD1 in Jatropha curcas, and found that it is mainly expressed in the stems, roots, and male flowers of Jatropha. JcDAD1-RNAi transgenic plants with low endogenous jasmonate levels in inflorescences exhibited more and larger flowers, as well as a few abortive female flowers, although anther and pollen development were normal. In addition, fruit number was increased and the seed size, weight, and oil contents were reduced in the transgenic Jatropha plants. These results indicate that JcDAD1 regulates the development of flowers and fruits through the JA biosynthesis pathway, but does not alter androecium development in Jatropha. These findings strengthen our understanding of the roles of JA and DAD1 in the regulation of floral development in woody perennial plants.

De novo genome assembly and Hi-C analysis reveal an association between chromatin architecture alterations and sex differentiation in the woody plant Jatropha curcas

BACKGROUND

Chromatin architecture is an essential factor regulating gene transcription in different cell types and developmental phases. However, studies on chromatin architecture in perennial woody plants and on the function of chromatin organization in sex determination have not been reported.

RESULTS

Here, we produced a chromosome-scale de novo genome assembly of the woody plant Jatropha curcas with a total length of 379.5 Mb and a scaffold N50 of 30.7 Mb using Pacific Biosciences long reads combined with genome-wide chromosome conformation capture (Hi-C) technology. Based on this high-quality reference genome, we detected chromatin architecture differences between monoecious and gynoecious inflorescence buds of Jatropha. Differentially expressed genes were significantly enriched in the changed A/B compartments and topologically associated domain regions and occurred preferentially in differential contact regions between monoecious and gynoecious inflorescence buds. Twelve differentially expressed genes related to flower development or hormone synthesis displayed significantly different genomic interaction patterns in monoecious and gynoecious inflorescence buds. These results demonstrate that chromatin organization participates in the regulation of gene transcription during the process of sex differentiation in Jatropha.

CONCLUSIONS

We have revealed the features of chromatin architecture in perennial woody plants and investigated the possible function of chromatin organization in Jatropha sex differentiation. These findings will facilitate understanding of the regulatory mechanisms of sex determination in higher plants.