Conserved noncoding sequences correlate with distant gene contacts in Arabidopsis and Brassica.

Physical contact between genes distant on chromosomes is a potentially important way for genes to coordinate their expressions. To investigate the potential importance of distant contacts, we performed high-throughput chromatin conformation capture (Hi-C) experiments on leaf nuclei isolated from Brassica rapa and Brassica oleracea. We then combined our results with published Hi-C data from Arabidopsis thaliana. We found that distant genes come into physical contact and do so preferentially between the proximal promoter of one gene and the downstream region of another gene. Genes with higher numbers of conserved noncoding sequences (CNSs) nearby were more likely to have contact with distant genes. With more CNSs came higher numbers of transcription factor binding sites and more histone modifications associated with the activity. In addition, for the genes we studied, distant contacting genes with CNSs were more likely to be transcriptionally coordinated. These observations suggest that CNSs may enrich active histone modifications and recruit transcription factors, correlating with distant contacts to ensure coordinated expression. This study advances our knowledge of gene contacts and provides insights into the relationship between CNSs and distant gene contacts in plants.

URA3 as a Selectable Marker for Disruption and Functional Assessment of PacC Gene in the Entomopathogenic Fungus Isaria javanica.

An effective selection marker is necessary for genetic engineering and functional genomics research in the post-genomic era. Isaria javanica is an important entomopathogenic fungus with a broad host range and prospective biocontrol potentials. Given that no antibiotic marker is available currently in this fungus, developing an effective selection marker is necessary. In this study, by applying overlap PCR and split-marker deletion strategy, combining PEG-mediated protoplasm transformation method, the uridine auxotrophy gene (ura3) in the I. javanica genome was knocked out. Then, using this transformation system, the pH response transcription factor gene (IjpacC) was disrupted successfully. Loss of IjpacC gene results in an obvious decrease in conidial production, but little impact on mycelial growth. The virulence of the ΔIjpacC mutant on caterpillars is similar to that of the wild-type strain. RT-qPCR detection shows that expression level of an acidic-expressed S53 gene (IF1G_06234) in ΔIjpacC mutant is more significantly upregulated than in the wild-type strain during the fungal infection on caterpillars. Our results indicate that a markerless transformation system based upon complementation of uridine auxotrophy is successfully developed in I. javanica, which is useful for exploring gene function and for genetic engineering to enhance biological control potential of the fungus.

Chromosomal-level genome assembly of potato tuberworm, Phthorimaea operculella: a pest of solanaceous crops.

The potato tuberworm, Phthorimaea operculella Zeller, is an oligophagous pest feeding on crops mainly belonging to the family Solanaceae. It is one of the most destructive pests of potato worldwide and attacks foliage and tubers in the field and in storage. However, the lack of a high-quality reference genome has hindered the association of phenotypic traits with their genetic basis. Here, we report on the genome assembly of P. operculella at the chromosomal level. Using Illumina, Nanopore and Hi-C sequencing, a 648.2 Mb genome was generated from 665 contigs, with an N50 length of 3.2 Mb, and 92.0% (596/648.2 Mb) of the assembly was anchored to 29 chromosomes. In total, 16619 genes were annotated, and 92.4% of BUSCO genes were fully represented. The chromosome-level genome of P. operculella will provide a significant resource for understanding the genetic basis for the biological study of this insect, and for promoting the integrative management of this pest in future.

Investigation of Brassica and its relative genomes in the post-genomics era.

The Brassicaceae family includes many economically important crop species, as well as cosmopolitan agricultural weed species. In addition, Arabidopsis thaliana, a member of this family, is used as a molecular model plant species. The genus Brassica is mesopolyploid, and the genus comprises comparatively recently originated tetrapolyploid species. With these characteristics, Brassicas have achieved the commonly accepted status of model organisms for genomic studies. This paper reviews the rapid research progress in the Brassicaceae family from diverse omics studies, including genomics, transcriptomics, epigenomics, and three-dimensional (3D) genomics, with a focus on cultivated crops. The morphological plasticity of Brassicaceae crops is largely due to their highly variable genomes. The origin of several important Brassicaceae crops has been established. Genes or loci domesticated or contributing to important traits are summarized. Epigenetic alterations and 3D structures have been found to play roles in subgenome dominance, either in tetraploid Brassica species or their diploid ancestors. Based on this progress, we propose future directions and prospects for the genomic investigation of Brassicaceae crops.

Secreted Glycosyltransferase RsIA_GT of Rhizoctonia solani AG-1 IA Inhibits Defense Responses in Nicotiana benthamiana.

Anastomosis group AG-1 IA of Rhizoctonia solani Khün has a wide host range and threatens crop production. Various glycosyltransferases secreted by phytopathogenic fungi play an essential role in pathogenicity. Previously, we identified a glycosyltransferase RsIA_GT (AG11A_09161) as a secreted protein-encoding gene of R. solani AG-1 IA, whose expression levels increased during infection in rice. In this study, we further characterized the virulence function of RsIA_GT. It is conserved not only in Basidiomycota, including multiple anastomosis groups of R. solani, but also in other primary fungal taxonomic categories. RsIA_GT possesses a signal peptide (SP) for protein secretion, and its functionality was proven using yeast and Nicotiana benthamiana. The SP-truncated form of RsIA_GT (RsIA_GT(ΔS)) expressed in Escherichia coli-induced lesion-like phenotype in rice leaves when applied to punched leaves. However, Agrobacterium-mediated transient expressions of both the full-length RsIA_GT and RsIA_GT(ΔS) did not induce cell death in N. benthamiana leaves. Instead, only RsIA_GT(ΔS) suppressed the cell death induced by two reference cell death factors BAX and INF1 in N.benthamiana. RsIA_GT(ΔS)R154A D168A D170A, a mutant RsIA_GT(ΔS) for the glycosyltransferase catalytic domain, still suppressed the BAX- or INF1-induced cell death, suggesting that the cell death suppression activity of RsIA_GT(ΔS) would be independent from its enzymatic activity. RsIA_GT(ΔS) also suppressed the H2O2 production and callose deposition and showed an effect on the induction of defense genes associated with the expression of BAX and INF1. The transient expression of RsIA_GT(ΔS) in N. benthamiana enhanced the lesion area caused by R. solani AG-1 IA. The secreted glycosyltransferase, RsIA_GT, of R. solani AG-1 IA is likely to have a dual role in virulence inside and outside of host cells.

Roles of Species-Specific Legumains in Pathogenicity of the Pinewood Nematode Bursaphelenchus xylophilus.

Peptidases are very important to parasites, which have central roles in parasite biology and pathogenesis. In this study, by comparative genome analysis, genome-wide peptidase diversities among plant-parasitic nematodes are estimated. We find that genes encoding cysteine peptidases in family C13 (legumain) are significantly abundant in pine wood nematodes Bursaphelenchus genomes, compared to those in other plant-parasitic nematodes. By phylogenetic analysis, a clade of B. xylophilus-specific legumain is identified. RT-qPCR detection shows that these genes are highly expressed at early stage during the nematode infection process. Utilizing transgene technology, cDNAs of three species-specific legumain were introduced into the Arabidopsis γvpe mutant. Functional complementation assay shows that these B. xylophilus legumains can fully complement the activity of Arabidopsis γVPE to mediate plant cell death triggered by the fungal toxin FB1. Secretory activities of these legumains are experimentally validated. By comparative transcriptome analysis, genes involved in plant cell death mediated by legumains are identified, which enrich in GO terms related to ubiquitin protein transferase activity in category molecular function, and response to stimuli in category biological process. Our results suggest that B. xylophilu-specific legumains have potential as effectors to be involved in nematode-plant interaction and can be related to host cell death.

Improved Reference Genome Annotation of Brassica rapa by Pacific Biosciences RNA Sequencing.

The species Brassica rapa includes several important vegetable crops. The draft reference genome of B. rapa ssp. pekinensis was completed in 2011, and it has since been updated twice. The pangenome with structural variations of 18 B. rapa accessions was published in 2021. Although extensive genomic analysis has been conducted on B. rapa, a comprehensive genome annotation including gene structure, alternative splicing (AS) events, and non-coding genes is still lacking. Therefore, we used the Pacific Biosciences (PacBio) single-molecular long-read technology to improve gene models and produced the annotated genome version 3.5. In total, we obtained 753,041 full-length non-chimeric (FLNC) reads and collapsed these into 92,810 non-redundant consensus isoforms, capturing 48% of the genes annotated in the B. rapa reference genome annotation v3.1. Based on the isoform data, we identified 830 novel protein-coding genes that were missed in previous genome annotations, defined the untranslated regions (UTRs) of 20,340 annotated genes and corrected 886 wrongly spliced genes. We also identified 28,564 AS events and 1,480 long non-coding RNAs (lncRNAs). We produced a relatively complete and high-quality reference transcriptome for B. rapa that can facilitate further functional genomic research.

Transposable element insertion: a hidden major source of domesticated phenotypic variation in Brassica rapa.

Transposable element (TE) is prevalent in plant genomes. However, studies on their impact on phenotypic evolution in crop plants are relatively rare, because systematically identifying TE insertions within a species has been a challenge. Here, we present a novel approach for uncovering TE insertion polymorphisms (TIPs) using pan-genome analysis combined with population-scale resequencing, and we adopt this pipeline to retrieve TIPs in a Brassica rapa germplasm collection. We found that 23% of genes within the reference Chiifu-401-42 genome harbored TIPs. TIPs tended to have large transcriptional effects, including modifying gene expression levels and altering gene structure by introducing new introns. Among 524 diverse accessions, TIPs broadly influenced genes related to traits and acted a crucial role in the domestication of B. rapa morphotypes. As examples, four specific TIP-containing genes were found to be candidates that potentially involved in various climatic conditions, promoting the formation of diverse vegetable crops in B. rapa. Our work reveals the hitherto hidden TIPs implicated in agronomic traits and highlights their widespread utility in studies of crop domestication.

Expanding the genetic variation of Brassica juncea by introgression of the Brassica rapa Genome.

Brassica juncea is an important vegetable and oil crop cultivated worldwide. To increase its genetic variation, we introgressed the A-genome of Brassica rapa into B. juncea. We used three each of heading and semi-heading B. juncea accessions as recipient parents and a B. rapa line, B9008, as the donor parent. We obtained 101 BC1S1 lines in total with expanded phenotypic variations such as leafy head shapes. We developed 132 SNP markers, which could distinguish the A genome of B. juncea from B. rapa genome, and tracked the introgression of B. rapa segments in the new B. juncea germplasm. On average, 59.2% of the B. juncea A genome in the B. juncea introgression lines was covered by the donor segments. We also identified three markers whose donor genotype frequencies were significantly lower than the theoretical value, suggesting strong selection of the recipient genotype during the introgression process. We provide an effective strategy to evaluate the diversity of the new germplasm based on the combination of parental re-sequencing data and marker genotyping results. Further genetic analysis of 1642 SNPs showed that the genetic diversity of the new B. juncea germplasm with the introgressed B. rapa genome was significantly increased. This study illustrated the potential for expanding the genetic diversity of B. juncea through the introgression of the B. rapa genome.

BRAD V3.0: an upgraded Brassicaceae database.

The Brassicaceae Database (BRAD version 3.0, BRAD V3.0; http://brassicadb.cn) has evolved from the former Brassica Database (BRAD V2.0), and represents an important community portal hosting genome information for multiple Brassica and related Brassicaceae plant species. Since the last update in 2015, the complex genomes of numerous Brassicaceae species have been decoded, accompanied by many omics datasets. To provide an up-to-date service, we report here a major upgrade of the portal. The Model-View-ViewModel (MVVM) framework of BRAD has been re-engineered to enable easy and sustainable maintenance of the database. The collection of genomes has been increased to 26 species, along with optimization of the user interface. Features of the previous version have been retained, with additional new tools for exploring syntenic genes, gene expression and variation data. In the 'Syntenic Gene @ Subgenome' module, we added features to view the sequence alignment and phylogenetic relationships of syntenic genes. New modules include 'MicroSynteny' for viewing synteny of selected fragment pairs, and 'Polymorph' for retrieval of variation data. The updated BRAD provides a substantial expansion of genomic data and a comprehensive improvement of the service available to the Brassicaceae research community.

Impacts of allopolyploidization and structural variation on intraspecific diversification in Brassica rapa.

BACKGROUND: Despite the prevalence and recurrence of polyploidization in the speciation of flowering plants, its impacts on crop intraspecific genome diversification are largely unknown. Brassica rapa is a mesopolyploid species that is domesticated into many subspecies with distinctive morphotypes. RESULTS: Herein, we report the consequences of the whole-genome triplication (WGT) on intraspecific diversification using a pan-genome analysis of 16 de novo assembled and two reported genomes. Among the genes that derive from WGT, 13.42% of polyploidy-derived genes accumulate more transposable elements and non-synonymous mutations than other genes during individual genome evolution. We denote such genes as being "flexible." We construct the Brassica rapa ancestral genome and observe the continuing influence of the dominant subgenome on intraspecific diversification in B. rapa. The gene flexibility is biased to the more fractionated subgenomes (MFs), in contrast to the more intact gene content of the dominant LF (least fractionated) subgenome. Furthermore, polyploidy-derived flexible syntenic genes are implicated in the response to stimulus and the phytohormone auxin; this may reflect adaptation to the environment. Using an integrated graph-based genome, we investigate the structural variation (SV) landscapes in 524 B. rapa genomes. We observe that SVs track morphotype domestication. Four out of 266 candidate genes for Chinese cabbage domestication are speculated to be involved in the leafy head formation. CONCLUSIONS: This pan-genome uncovers the possible contributions of allopolyploidization on intraspecific diversification and the possible and underexplored role of SVs in favorable trait domestication. Collectively, our work serves as a rich resource for genome-based B. rapa improvement.

Divergence of three BRX homoeologs in Brassica rapa and its effect on leaf morphology.

The leafy head characteristic is a special phenotype of Chinese cabbage resulting from artificial selection during domestication and breeding. BREVIS RADIX (BRX) has been suggested to control root elongation, shoot growth, and tiller angle in Arabidopsis and rice. In Brassica rapa, three BrBRX homoeologs have been identified, but only BrBRX.1 and BrBRX.2 were found to be under selection in leaf-heading accessions, indicating their functional diversification in leafy head formation. Here, we show that these three BrBRX genes belong to a plant-specific BRX gene family but that they have significantly diverged from other BRX-like members on the basis of different phylogenetic classifications, motif compositions and expression patterns. Moreover, although the expression of these three BrBRX genes differed, compared with BrBRX.3, BrBRX.1, and BrBRX.2 displayed similar expression patterns. Arabidopsis mutant complementation studies showed that only BrBRX.1 could rescue the brx root phenotype, whereas BrBRX.2 and BrBRX.3 could not. However, overexpression of each of the three BrBRX genes in Arabidopsis resulted in similar pleiotropic leaf phenotypes, including epinastic leaf morphology, with an increase in leaf number and leaf petiole length and a reduction in leaf angle. These leaf traits are associated with leafy head formation. Further testing of a SNP (T/C) in BrBRX.2 confirmed that this allele in the heading accessions was strongly associated with the leaf-heading trait of B. rapa. Our results revealed that all three BrBRX genes may be involved in the leaf-heading trait, but they may have functionally diverged on the basis of their differential expression.

Comparative Mitogenomic Analysis and the Evolution of Rhizoctonia solani Anastomosis Groups.

Mitochondria are the major energy source for cell functions. However, for the plant fungal pathogens, mitogenome variations and their roles during the host infection processes remain largely unknown. Rhizoctonia solani, an important soil-borne pathogen, forms different anastomosis groups (AGs) and adapts to a broad range of hosts in nature. Here, we reported three complete mitogenomes of AG1-IA RSIA1, AG1-IB RSIB1, and AG1-IC, and performed a comparative analysis with nine published Rhizoctonia mitogenomes (AG1-IA XN, AG1-IB 7/3/14, AG3, AG4, and five Rhizoctonia sp. mitogenomes). These mitogenomes encoded 15 typical proteins (cox1-3, cob, atp6, atp8-9, nad1-6, nad4L, and rps3) and several LAGLIDADG/GIY-YIG endonucleases with sizes ranging from 109,017 bp (Rhizoctonia sp. SM) to 235,849 bp (AG3). We found that their large sizes were mainly contributed by repeat sequences and genes encoding endonucleases. We identified the complete sequence of the rps3 gene in 10 Rhizoctonia mitogenomes, which contained 14 positively selected sites. Moreover, we inferred a robust maximum-likelihood phylogeny of 32 Basidiomycota mitogenomes, representing that seven R. solani and other five Rhizoctonia sp. lineages formed two parallel branches in Agaricomycotina. The comparative analysis showed that mitogenomes of Basidiomycota pathogens had high GC content and mitogenomes of R. solani had high repeat content. Compared to other strains, the AG1-IC strain had low substitution rates, which may affect its mitochondrial phylogenetic placement in the R. solani clade. Additionally, with the published RNA-seq data, we investigated gene expression patterns from different AGs during host infection stages. The expressed genes from AG1-IA (host: rice) and AG3 (host: potato) mainly formed four groups by k-mean partitioning analysis. However, conserved genes represented varied expression patterns, and only the patterns of rps3-nad2 and nad1-m3g18/mag28 (an LAGLIDADG endonuclease) were conserved in AG1-IA and AG3 as shown by the correlation coefficient analysis, suggesting regulation of gene repertoires adapting to infect varied hosts. The results of variations in mitogenome characteristics and the gene substitution rates and expression patterns may provide insights into the evolution of R. solani mitogenomes.

Series-Spatial Transcriptome Profiling of Leafy Head Reveals the Key Transition Leaves for Head Formation in Chinese Cabbage.

Chinese cabbage is an important leaf heading vegetable crop. At the heading stage, its leaves across inner to outer show significant morphological differentiation. However, the genetic control of this complex leaf morphological differentiation remains unclear. Here, we reported the transcriptome profiling of Chinese cabbage plant at the heading stage using 24 spatially dissected tissues representing different regions of the inner to outer leaves. Genome-wide transcriptome analysis clearly separated the inner leaf tissues from the outer leaf tissues. In particular, we identified the key transition leaf by the spatial expression analysis of key genes for leaf development and sugar metabolism. We observed that the key transition leaves were the first inwardly curved ones. Surprisingly, most of the heading candidate genes identified by domestication selection analysis obviously showed a corresponding expression transition, supporting that key transition leaves are related to leafy head formation. The key transition leaves were controlled by a complex signal network, including not only internal hormones and protein kinases but also external light and other stimuli. Our findings provide new insights and the rich resource to unravel the genetic control of heading traits.

Improved Brassica oleracea JZS assembly reveals significant changing of LTR-RT dynamics in different morphotypes.

Brassica oleracea is an important vegetable crop that has provided ancestor genomes of the two most important Brassica oil crops, Brassica napus and Brassica carinata. The current B. oleracea reference genome (JZS, also named 02-12) displays problems of large mis-assemblies, low sequence continuity, and low assembly integrity, thus limiting genomic analysis. We reported an updated assembly of the B. oleracea reference genome (JZS v2) obtained through single-molecule sequencing and chromosome conformation capture technologies. We assembled an additional 83.16 Mb of genomic sequences, and the updated genome features a contig N50 size of 2.37 Mb, representing an ~ 88-fold improvement. We detected a new round of long terminal repeat retrotransposon (LTR-RT) burst in the new assembly. Comparative analysis with the reported genome sequences of two other genomes of B. oleracea (TO1000 and HDEM) identified extensive gene order and gene structural variation. In addition, we found that the genome-specific amplification of Gypsy-like LTR-RTs occurred around 0-1 million years ago (MYA). In particular, the athila, tat, and Del families were extensively amplified in JZS around 0-1 MYA. Moreover, we identified that the syntenic genes were modified due to the insertion of genome-specific LTR-RTs. These results indicated that the genome-specific LTR-RT dynamics were associated with genome diversification in B. oleracea.

Integrative transcriptome analysis discloses the molecular basis of a heterogeneous fungal phytopathogen complex, Rhizoctonia solani AG-1 subgroups.

Rhizoctonia solani is a fungal species complex that causes necrotrophic crop diseases. It comprises several anastomosis groups, some of which include intra-subgroups, such as AG-1 IA and AG-1 IB, exhibiting varying pathogenicity. Owing to its heterozygous and multinucleate features, genomic analyses of R. solani are still challenging, and understanding of its genetic diversity and genic components is limited. In this study, in order to elucidate the molecular basis of this phytopathogen complex, an integrated transcriptome analysis was undertaken for three subgroups of AG-1, i.e. AG-1 IA, AG-1 IB, and AG-1 IC. Sequence variations suggested substantial evolutionary distances within AG-1. Transcript simple sequence repeats showed comparable characteristics among AG-1, but contained polymorphic sites. Intra-subgroup polymorphisms suggested varying genic heterozygosity within AG-1, suggesting their independent evolutionary trajectory. Sequences of pathogenic factors, phytotoxin biosynthesis pathway enzymes, secreted lignocellulosic enzymes, secreted reactive oxygen species detoxification enzymes, apoplastic/cytoplasmic effector candidates, were conserved among those subgroups. dN/dS ratios of a secretome subset suggested core secreted proteins in AG-1 and distinct evolution of Cys-rich small secreted proteins after differentiation of AG-1 subgroups. Identification of likely pathogenic factors including allergen protein homologues, oxidative phosphorylation and ethylene biosynthesis pathways, and diversification of polysaccharide monooxygenases provides molecular insight into key genomic components that play a role in R. solani pathogenesis.

Genome sequence of Isaria javanica and comparative genome analysis insights into family S53 peptidase evolution in fungal entomopathogens.

The fungus Isaria javanica is an important entomopathogen that parasitizes various insects and is effective for pest control. In this study, we sequenced and assembled the genomes (IJ1G and IJ2G) of two I. javanica strains isolated from different insects. The genomes were approximately 35 Mb in size with 11,441 and 11,143 protein-coding genes, respectively. Using a phylogenomic approach, we evaluated genome evolution across five entomopathogenic fungi in Cordycipitaceae. By comparative genome analysis, it was found that family S53 serine peptidases were expanded in Cordycipitaceae entomopathogens, particularly in I. javanica. Gene duplication events were identified based on phylogenetic relationships inferred from 82 S53 peptidases within six entomopathogenic fungal genomes. Moreover, we found that carbohydrate-active enzymes and proteinases were the largest secretory protein groups encoded in the I. javanica genome, especially chitinases (GH18), serine and aspartic peptidases (S53, S08, S10, A01). Pathogenesis-related genes and genes for bacterial-like toxins and secondary metabolites were also identified. By comparative transcriptome analysis, differentially expressed genes in response to insect nutrients (in vitro) were identified. Moreover, most S53 peptidases were detected to be significantly upregulated during the initial fungal infection process in insects (in vivo) by RT-qPCR. Our results provide new clues about understanding evolution of pathogenic proteases and may suggest that abundant S53 peptidases in the I. javanica genome may contribute to its effective parasitism on various insects.

Antibacterial Radicicol Analogues from Pochonia chlamydosporia and Their Biosynthetic Gene Cluster.

Chemical investigation of fungus Pochonia chlamydosporia strain 170, derived from rice fermentation sediment samples, afforded seven radicicol analogues, including two new compounds, monocillin VI (1) and monocillin VII (2), and five known compounds, monocillin II (3), monorden D (4), monocillin IV (5), monocillin V (6), and pochonin M (7). The structures of compounds 1-7 were established primarily by analysis of nuclear magnetic resonance data, and the absolute configurations of the secondary alcohol in compounds 1 and 2 were assigned by the modified Mosher method. All seven compounds have modest antibacterial activities, with a minimal inhibitory concentration (MIC) of 25.6 µg/mL for compounds 1 and 3-7 and 51.2 µg/mL for compound 2, on inhibition of the growth of the plant pathogen Xanthomonas campestris (the positive control ampicillin showed a MIC value of 12.8 µg/mL), indicating that the fungus has the potential to control bacterial disease. The biosynthetic gene cluster and putative biosynthetic pathways of these radicicol analogues in the P. chlamydosporia genome were proposed. These findings increase our knowledge of the chemical potential of P. chlamydosporia and may allow us to better utilize the fungus as a biological control agent.