Characterization and acceleration of genome shuffling and ploidy reduction in synthetic allopolyploids by genome sequencing and editing.

Polyploidy and the subsequent ploidy reduction and genome shuffling are the major driving forces of genome evolution. Here, we revealed short-term allopolyploid genome evolution by sequencing a synthetic intergeneric hybrid (Raphanobrassica, RRCC). In this allotetraploid, the genome deletion was quick, while rearrangement was slow. The core and high-frequency genes tended to be retained while the specific and low-frequency genes tended to be deleted in the hybrid. The large-fragment deletions were enriched in the heterochromatin region and probably derived from chromosome breaks. The intergeneric translocations were primarily of short fragments dependent on homoeology, indicating a gene conversion origin. To accelerate genome shuffling, we developed an efficient genome editing platform for Raphanobrassica. By editing Fanconi Anemia Complementation Group M (FANCM) genes, homoeologous recombination, chromosome deletion and secondary meiosis with additional ploidy reduction were accelerated. FANCM was shown to be a checkpoint of meiosis and controller of ploidy stability. By simultaneously editing FLIP genes, gene conversion was precisely introduced, and mosaic genes were produced around the target site. This intergeneric hybrid and genome editing platform not only provides models that facilitate experimental evolution research by speeding up genome shuffling and conversion but also accelerates plant breeding by enhancing intergeneric genetic exchange and creating new genes.

Proteome and Transcriptome Analysis of Gonads Reveals Intersex in Gigantidas haimaensis.

Sex has proven to be one of the most intriguing areas of research across evolution, development, and ecology. Intersex or sex change occurs frequently in molluscs. The deep-sea mussel Gigantidas haimaensis often dominates within Haima cold seep ecosystems, but details of their reproduction remain unknown. Herein, we conducted a combined proteomic and transcriptomic analysis of G. haimaensis gonads to provide a systematic understanding of sexual development in deep-sea bivalves. A total of 2,452 out of 42,238 genes (5.81%) and 288 out of 7,089 proteins (4.06%) were significantly differentially expressed between ovaries and testes with a false discovery rate (FDR) <0.05. Candidate genes involved in sexual development were identified; among 12 differentially expressed genes between sexes, four ovary-biased genes (ß-catenin, fem-1, forkhead box L2 and membrane progestin receptor α) were expressed significantly higher in males than females. Combining histological characteristics, we speculate that the males maybe intersex undergoing sex change, and implied that these genes may be involved in the process of male testis converting into female gonads in G. haimaensis. The results suggest that this adaptation may be based on local environmental factors, sedentary lifestyles, and patchy distribution, and sex change may facilitate adaptation to a changing environment and expansion of the population. The findings provide a valuable genetic resource to better understand the mechanisms of sex change and survival strategies in deep-sea bivalves.

Molecular and Functional Characterization of a Novel Kunitz-Type Toxin-like Peptide in the Giant Triton Snail Charonia tritonis.

It has been reported that the giant triton snail (Charonia tritonis) inserts its large proboscis and then injects venom or acid saliva from its salivary gland into its prey, the crown-of-thorns starfish Acanthaster planci (COTS), paralyzing it. A full-length cDNA sequence of the C. tritonis Ct-kunitzin gene was obtained by RACE PCR based on a transcriptomic database constructed by our laboratory (data not published), which contains an open reading frame (ORF) sequence with a length of 384 bp including a 1-32aa Kunitz domain. The Ct-kunitzin peptide was synthesized by solid-phase polypeptide methods according to its conserved amino acid sequence, with a molecular weight of 3746.0 as well as two disulfide bonds. Renatured Ct-kunitzin was injected into mice ventricles to evaluate its potential function. Compared with the normal control group (physiological saline), the spontaneous locomotor activity of the Ct-kunitzin group decreased significantly. There was a significant effect on Ct-kunitzin on mice grip strength in the grip strength test. In addition, Ct-kunitzin exhibited remarkable biological activity in suppressing pain in the pain thresholds test. There were no significant differences between the Ct-kunitzin group and the normal control group in terms of various hematological indexes and histopathological observations. When tested in COTS, the most significant histological change was the destruction, disorganization, and significant reduction in the amount of COTS tube feet tissues. Altogether, the potential paralyzing effect on mice suggests that Ct-kunitzin is a possible agent for novel drug development.

Molecular and Functional Characterization of a Short-Type Peptidoglycan Recognition Protein, Ct-PGRP-S1 in the Giant Triton Snail Charonia tritonis.

Peptidoglycan recognition proteins (PGRPs) are a family of pattern recognition receptors (PRRs) involved in host antibacterial responses, and their functions have been characterized in most invertebrate and vertebrate animals. However, little information is available regarding the potential function of PGRPs in the giant triton snail Charonia tritonis. In this study, a short-type PGRP gene (termed Ct-PGRP-S1) was identified in C. tritonis. Ct-PGRP-S1 was predicted to contain several structural features known in PGRPs, including a typical PGRP domain (Amidase_2) and Src homology-3 (SH3) domain. The Ct-PGRP-S1 gene was constitutively expressed in all tissues examined except in proboscis, with the highest expression level observed in the liver. As a typical PRR, Ct-PGRP-S1 has an ability to degrade peptidoglycan (PGN) and was proven to have non-Zn2+-dependent amidase activity and antibacterial activity against Vibrioalginolyticus and Staphylococcus aureus. It is the first report to reveal the peptidoglycan recognition protein in C. tritonis, and these results suggest that peptidoglycan recognition protein Ct-PGRP-S1 is an important effector of C. tritonis that modulates bacterial infection resistance of V. alginolyticus and S. aureus, and this study may provide crucial basic data for the understanding of an innate immunity system of C. tritonis.

Genome resequencing reveals demographic history and genetic architecture of seed salinity tolerance in Populus euphratica.

Populus euphratica is a dominant tree species in desert riparian forests and possesses extraordinary adaptation to salinity stress. Exploration of its genomic variation and molecular underpinning of salinity tolerance is important for elucidating population evolution and identifying stress-related genes. Here, we identify approximately 3.15 million single nucleotide polymorphisms using whole-genome resequencing. The natural populations of P. euphratica in northwest China are divided into four distinct clades that exhibit strong geographical distribution patterns. Pleistocene climatic fluctuations and tectonic deformation jointly shaped the extant genetic patterns. A seed germination rate-based salinity tolerance index was used to evaluate seed salinity tolerance of P. euphratica and a genome-wide association study was implemented. A total of 38 single nucleotide polymorphisms were associated with seed salinity tolerance and were located within or near 82 genes. Expression profiles showed that most of these genes were regulated under salt stress, revealing the genetic complexity of seed salinity tolerance. Furthermore, DEAD-box ATP-dependent RNA helicase 57 and one undescribed gene (CCG029559) were demonstrated to improve the seed salinity tolerance in transgenic Arabidopsis. These results provide new insights into the demographic history and genetic architecture of seed salinity tolerance in desert poplar.

Comparative metabolomics analysis reveals different metabolic responses to drought in tolerant and susceptible poplar species.

Drought is one of the critical factors limiting tree growth and survival. Clarifying the adaptation to drought will facilitate the cultivation of drought-tolerant varieties. Metabolites, as direct signatures of biochemical functions, can uncover the biochemical pathways involved in drought responses. Here, we investigated the physiological and metabolic responses of drought-tolerant Populus simonii and drought-susceptible Populus deltoides cv. 'Danhong' to drought. Under drought conditions, P. simonii grew better and had a higher photosynthetic rate than P. deltoides cv. 'Danhong'. Global untargeted metabolite profiling was analyzed using gas chromatography time-of-flight mass spectrometry system. A total of 69 and 53 differentially accumulated metabolites were identified in drought-stressed P. simonii and P. deltoides cv. 'Danhong', respectively. The metabolisms of carbohydrate, amino acid, lipid and energy were involved in the drought responses common to both poplar species. The citric acid cycle was significantly inhibited to conserve energy, whereas multiple carbohydrates acting as osmolytes and osmoprotectants were induced to alleviate the adverse effects of drought stress. Unlike P. deltoides cv. 'Danhong', P. simonii underwent a specific metabolic reprogramming that enhanced non-enzymatic antioxidants, coordinated the cellular carbon/nitrogen balance and regulated wax biosynthesis. These results provide a reference for characterizing the mechanisms involved in poplar response to drought and for enhancing the drought tolerance of forest trees.

Genome-wide transcriptomic analysis of a desert willow, Salix psammophila, reveals the function of hub genes SpMDP1 and SpWRKY33 in drought tolerance.

BACKGROUND: Drought is a major environmental constraint to plant growth, development and productivity. Compared with most willows that are generally susceptible to drought, the desert willow Salix psammophila has extraordinary adaptation to drought stress. However, its molecular basis of drought tolerance is still largely unknown. RESULTS: During polyethylene glycol 6000 (PEG 6000)-simulated drought stress, we found that the osmotic adjustment substances were accumulated and the antioxidant enzyme activities were enhanced in S. psammophila roots. A total of 8172 differentially expressed genes were identified in roots of S. psammophila through RNA-Sequencing. Based on K-means clustering, their expression patterns were classified into nine clusters, which were enriched in several stress-related processes including transcriptional regulation, response to various stresses, cell death, etc. Moreover, 672 transcription factors from 45 gene families were differentially expressed under drought stress. Furthermore, a weighted gene co-expression network was constructed, and eight genes were identified as hub genes. We demonstrated the function of two hub genes, magnesium-dependent phosphatase 1 (SpMDP1) and SpWRKY33, through overexpression in Arabidopsis thaliana. Overexpression of the two hub genes enhanced the drought tolerance in transgenic plants, suggesting that the identification of candidate drought tolerance genes in this study was highly efficient and credible. CONCLUSIONS: Our study analyzed the physiological and molecular responses to drought stress in S. psammophila, and these results contribute to dissect the mechanism of drought tolerance of S. psammophila and facilitate identification of critical genes involved in drought tolerance for willow breeding.

A Stress-Associated Protein, PtSAP13, From Populus trichocarpa Provides Tolerance to Salt Stress.

The growth and production of poplars are usually affected by unfavorable environmental conditions such as soil salinization. Thus, enhancing salt tolerance of poplars will promote their better adaptation to environmental stresses and improve their biomass production. Stress-associated proteins (SAPs) are a novel class of A20/AN1 zinc finger proteins that have been shown to confer plants' tolerance to multiple abiotic stresses. However, the precise functions of SAP genes in poplars are still largely unknown. Here, the expression profiles of Populus trichocarpa SAPs in response to salt stress revealed that PtSAP13 with two AN1 domains was up-regulated dramatically during salt treatment. The ß-glucuronidase (GUS) staining showed that PtSAP13 was accumulated dominantly in leaf and root, and the GUS signal was increased under salt condition. The Arabidopsis transgenic plants overexpressing PtSAP13 exhibited higher seed germination and better growth than wild-type (WT) plants under salt stress, demonstrating that overexpression of PtSAP13 increased salt tolerance. Higher activities of antioxidant enzymes were found in PtSAP13-overexpressing plants than in WT plants under salt stress. Transcriptome analysis revealed that some stress-related genes, including Glutathione peroxidase 8, NADP-malic enzyme 2, Response to ABA and Salt 1, WRKYs, Glutathione S-Transferase, and MYBs, were induced by salt in transgenic plants. Moreover, the pathways of flavonoid biosynthesis and metabolic processes, regulation of response to stress, response to ethylene, dioxygenase activity, glucosyltransferase activity, monooxygenase activity, and oxidoreductase activity were specially enriched in transgenic plants under salt condition. Taken together, our results demonstrate that PtSAP13 enhances salt tolerance through up-regulating the expression of stress-related genes and mediating multiple biological pathways.

Deciphering Genetic Architecture of Adventitious Root and Related Shoot Traits in Populus Using QTL Mapping and RNA-Seq Data.

Understanding the genetic architecture of adventitious root and related shoot traits will facilitate the cultivation of superior genotypes. In this study, we measured 12 adventitious root and related shoot traits of 434 F1 genotypes originating from Populus deltoides 'Danhong' × Populus simonii 'Tongliao1' and conducted an integrative analysis of quantitative trait locus (QTL) mapping and RNA-Seq data to dissect their genetic architecture and regulatory genes. Extensive segregation, high repeatability, and significant correlation relationship were detected for the investigated traits. A total of 150 QTLs were associated with adventitious root traits, explaining 3.1-6.1% of phenotypic variation (PVE); while 83 QTLs were associated with shoot traits, explaining 3.1-19.8% of PVE. Twenty-five QTL clusters and 40 QTL hotspots were identified for the investigated traits. Ten QTL clusters were overlapped in both adventitious root traits and related shoot traits. Transcriptome analysis identified 10,172 differentially expressed genes (DEGs) among two parents, three fine rooting and three poor-rooting genotypes, 143 of which were physically located within the QTL intervals. K-means cluster and weighted gene co-expression network analysis showed that PtAAAP19 (Potri.004G111400) encoding amino acid transport protein was tightly associated with adventitious roots and highly expressed in fine-rooting genotypes. Compare with 'Danhong', 153 bp deletion in the coding sequence of PtAAAP19 in 'Tongliao1' gave rise to lack one transmembrane domain, which might cause the variation of adventitious roots. Taken together, this study deciphered the genetic basis of adventitious root and related shoot traits and provided potential function genes for genetic improvement of poplar breeding.

Genome-Wide Characterization of the sHsp Gene Family in Salix suchowensis Reveals Its Functions under Different Abiotic Stresses.

Small heat shock proteins (sHsps) function mainly as molecular chaperones that play vital roles in response to diverse stresses, especially high temperature. However, little is known about the molecular characteristics and evolutionary history of the sHsp family in Salix suchowensis, an important bioenergy woody plant. In this study, 35 non-redundant sHsp genes were identified in S. suchowensis, and they were divided into four subfamilies (C, CP, PX, and MT) based on their phylogenetic relationships and predicted subcellular localization. Though the gene structure and conserved motif were relatively conserved, the sequences of the Hsp20 domain were diversified. Eight paralogous pairs were identified in the Ssu-sHsp family, in which five pairs were generated by tandem duplication events. Ka/Ks analysis indicated that Ssu-sHsps had undergone purifying selection. The expression profiles analysis showed Ssu-Hsps tissue-specific expression patterns, and they were induced by at least one abiotic stress. The expression correlation between two paralogous pairs (Ssu-sHsp22.2-CV/23.0-CV and 23.8-MT/25.6-MT) were less than 0.6, indicating that they were divergent during the evolution. Various cis-acting elements related to stress responses, hormone or development, were detected in the promoter of Ssu-sHsps. Furthermore, the co-expression network revealed the potential mechanism of Ssu-sHsps under stress tolerance and development. These results provide a foundation for further functional research on the Ssu-sHsp gene family in S. suchowensis.

The Populus trichocarpa PtHSP17.8 involved in heat and salt stress tolerances.

KEY MESSAGE: PtHSP17.8 was regulated by various abiotic stresses. Overexpression of PtHSP17.8 enhanced the tolerance to heat and salt stresses through maintain ROS homeostasis and cooperate with stress-related genes in Arabidopsis. Small heat shock proteins (sHSPs) play important roles in response to diverse biotic and abiotic stresses, especially in heat tolerance. However, limited information is available on the stress tolerance roles of sHSPs in woody species. To explore the function of sHSPs in poplar, we isolated and characterized PtHSP17.8 from Populus trichocarpa. Phylogenetic analysis and subcellular localization revealed that PtHSP17.8 was a cytosolic class I sHSP. The gene expression profile of PtHSP17.8 in various tissues showed that it was significantly accumulated in stem and root, which was consistent with the GUS expression pattern driven by promoter of PtHSP17.8. The expression of PtHSP17.8 could be induced by various abiotic stresses and significantly activated by heat stress. Overexpression of PtHSP17.8 enhanced the tolerance to heat and salt stresses in Arabidopsis. The seedling survival rate, root length, relative water content, antioxidative enzyme activities, proline, and soluble sugar content were increased in transgenic Arabidopsis under heat and salt stresses, but not in normal condition. The co-expression network of PtHSP17.8 were constructed and demonstrated many stress responsive genes included. The stress-related genes in the co-expression network were up-regulated in the PtHSP17.8 overexpression seedlings. These results suggest that PtHSP17.8 confers heat and salt tolerances in plants.

Comparative Proteomic Analysis of Mature Pollen in Triploid and Diploid Populus deltoides.

Ploidy affects plant growth vigor and cell size, but the relative effects of pollen fertility and allergenicity between triploid and diploid have not been systematically examined. Here we performed comparative analyses of fertility, proteome, and abundances of putative allergenic proteins of pollen in triploid poplar 'ZhongHuai1' ('ZH1', triploid) and 'ZhongHuai2' ('ZH2', diploid) generated from the same parents. The mature pollen was sterile in triploid poplar 'ZH1'. By applying two-dimensional gel electrophoresis (2-DE), a total of 72 differentially expressed protein spots (DEPs) were detected in triploid poplar pollen. Among them, 24 upregulated and 43 downregulated proteins were identified in triploid poplar pollen using matrix-assisted laser desorption/ionisation coupled with time of-flight tandem mass spectrometer analysis (MALDI-TOF/TOF MS/MS). The main functions of these DEPs were related with "S-adenosylmethionine metabolism", "actin cytoskeleton organization", or "translational elongation". The infertility of triploid poplar pollen might be related to its abnormal cytoskeletal system. In addition, the abundances of previously identified 28 putative allergenic proteins were compared among three poplar varieties ('ZH1', 'ZH2', and '2KEN8'). Most putative allergenic proteins were downregulated in triploid poplar pollen. This work provides an insight into understanding the protein regulation mechanism of pollen infertility and low allergenicity in triploid poplar, and gives a clue to improving poplar polyploidy breeding and decreasing the pollen allergenicity.

[Genetic background of lily germplasm resources based on SSR markers].

To study the genetic background of lily (Lilium spp.) germplasm resources, and accurately evaluate and select excellent germplasm for genetic improvement of lily, we analyzed the genetic background of 62 lily germplasm accessions from 11 provinces of China by using simple sequence repeat (SSR) molecular markers. The results showed that 15 out of 83 pairs of lily SSR primers were polymorphic. A total of 157 allelic loci were amplified, with the number of alleles per locus ranging from 5 to 19 and the average number of effective alleles per locus being 4.162 8. The average observed heterozygosity and expected heterozygosity were 0.228 2 and 0.694 1, respectively. The average polymorphic information content was 0.678 8. The average Nei's diversity index and Shannon's information index were 0.694 1 and 1.594 9, respectively, indicating that the tested lily germplasm had high genetic diversity. The 62 germplasm accessions were classified into 5 groups by the unweighted pair group method with arithmetic mean (UPGMA) and into 3 groups by the principal component analysis. The two analyses revealed a geographic correlation among different groups. The majority of lily germplasm accessions from the same source tended to cluster together. The population structure analysis classified the lily accessions into 4 populations and 1 mixed population. The above results provide a theoretical basis and genetic resources for the precise identification and breeding of lily germplasm resources.

[Genome-wide identification and characterization of the WOX gene family in Brassica juncea].

The WUSCHEL related-homeobox (WOX) family is one of the plant-specific transcription factor families, playing important roles in plant growth and development. In this study, 51 WOX gene family members were identified from the genome data of Brassica juncea by searching and screening with HUMMER, Smart and other software. Their protein molecular weight, amino acids numbers, and isoelectric point were analyzed by using Expasy online software. Furthermore, bioinformatics software was used to systematically analyze the evolutionary relationship, conservative region, and gene structure of the WOX gene family. The mustard WOX gene family was divided into three subfamilies: ancient clade, intermediate clade, and WUS clade/modern clade. Structural analysis showed that the type, organization form and gene structure of the conservative domain of WOX transcription factor family members in the same subfamily were highly consistent, while there was a certain diversity among different subfamilies. 51 WOX genes are distributed unevenly on 18 chromosomes of mustard. Most of the promoters of these genes contain cis acting elements related to light, hormone and abiotic stress. Using transcriptome data and real-time fluorescence quantitative PCR (qRT-PCR) analysis, it was found that the expression of mustard WOX gene was spatio-temporal specific, among which BjuWOX25, BjuWOX33, and BjuWOX49 might play an important role in the development of silique, and BjuWOX10, BjuWOX32, and BjuWOX11, BjuWOX23 respectively might play an important role in the response to drought and high temperature stresses. The above results may facilitate the functional study of mustard WOX gene family.

Large-scale population structure and genetic architecture of agronomic traits of garlic.

Garlic, an asexually propagated crop, is the second important bulb crop after the onion and is used as a vegetable and medicinal plant. Abundant and diverse garlic resources have been formed over thousands of years of cultivation. However, genome variation, population structure and genetic architecture of garlic agronomic traits were still not well elucidated. Here, 1 100 258 single nucleotide polymorphisms (SNPs) were identified using genotyping-by-sequencing in 606 garlic accessions collected from 43 countries. Population structure, principal component and phylogenetic analysis showed that these accessions were divided into five subpopulations. Twenty agronomic traits, including above-ground growth traits, bulb-related and bolt-related traits in two consecutive years were implemented in a genome-wide association study. In total, 542 SNPs were associated with these agronomic traits, among which 188 SNPs were repeatedly associated with more than two traits. One SNP (chr6: 1896135972) was repeatedly associated with ten traits. These associated SNPs were located within or near 858 genes, 56 of which were transcription factors. Interestingly, one non-synonymous SNP (Chr4: 166524085) in ribosomal protein S5 was repeatedly associated with above-ground growth and bulb-related traits. Additionally, gene ontology enrichment analysis of candidate genes for genomic selection regions between complete-bolting and non-bolting accessions showed that these genes were significantly enriched in 'vegetative to reproductive phase transition of meristem', 'shoot system development', 'reproductive process', etc. These results provide valuable information for the reliable and efficient selection of candidate genes to achieve garlic genetic improvement and superior varieties.

Combined widely targeted metabolomics and transcriptomics analysis reveals differentially accumulated metabolites and the underlying molecular bases in fleshy taproots of distinct radish genotypes.

Radish is an important taproot crop with medicinal and edible uses that is cultivated worldwide. However, the differences in metabolites and the underlying molecular bases among different radish types remain largely unknown. In the present study, a combined analysis of liquid chromatography-electrospray ionization-tandem mass spectrometry (LC-ESI-MS/MS) and RNA-Seq data was conducted to uncover important differentially accumulated metabolites (DAMs) among radish accessions with green, white and red taproot flesh colours. A total of 657 metabolites were identified and 138 DAMs were commonly present in red vs. green and red vs. white accessions. Red accessions were rich in phenolic compounds, while green and white accessions had more amino acids. Additionally, 41 metabolites and 98 genes encoding 37 enzymes were enriched in the shikimate and phenolic biosynthesis pathways. CHS is the rate-limiting enzyme determining flavonoid differences among accessions. A total of 119 candidate genes might contribute to red accession-specific accumulated metabolites. Specifically, one gene cluster consisting of 16 genes, including one RsMYB1, RsMYBL2, RsTT8, RsDFR, RsANS, Rs4CL3, RsSCPL10, Rs3AT1 and RsSAP2 gene, two RsTT19 and RsWRKY44 genes and three RsUGT genes, might be involved in anthocyanin accumulation in red radish fleshy taproots. More importantly, an InDel marker was developed based on an RsMYB1 promoter mutation, and the accuracy reached 95.9% when it was used to select red-fleshed radishes. This study provides comprehensive insights into the metabolite differences and underlying molecular mechanisms in fleshy taproots among different radish genotypes and will be beneficial for the genetic improvement of radish nutritional quality.

Genome resequencing reveals the evolutionary history of garlic reproduction traits.

The propagation of cultivated garlic relies on vegetative cloves, thus flowers become non-essential for reproduction in this species, driving the evolution of reproductive feature-derived traits. To obtain insights into the evolutionary alteration of reproductive traits in the clonally propagated garlic, the evolutionary histories of two main reproduction-related traits, bolting and flower differentiation, were explored by genome analyses using 134 accessions displaying wide diversity in these two traits. Resequencing identified 272.8 million variations in the garlic genome, 198.0 million of which represent novel variants. Population analysis identified five garlic groups that have evolved into two clades. Gene expression, single-cell transcriptome sequencing, and genome-wide trait association analyses have identified numerous candidates that correlate with reproductive transition and flower development, some of which display distinct selection signatures. Selective forces acting on the B-box zinc finger protein-encoding Asa2G00291.1, the global transcription factor group E protein-encoding Asa5G01527.1, and VERNALIZATION INSENSITIVE 3-like Asa3G03399.1 appear to be representative of the evolution of garlic bolting. Plenty of novel genomic variations and trait-related candidates represent valuable resources for biological studies of garlic. Numerous selective signatures from genes associated with the two chosen reproductive traits provide important insights into the evolutionary history of reproduction in this clonally propagated crop.

Effects of scale worm parasitism on interactions between the symbiotic gill microbiome and gene regulation in deep sea mussel hosts.

Diverse adaptations to the challenging deep sea environment are expected to be found across all deep sea organisms. Scale worms Branchipolynoe pettiboneae are believed to adapt to the deep sea environment by parasitizing deep sea mussels; this biotic interaction is one of most known in the deep sea chemosynthetic ecosystem. However, the mechanisms underlying the effects of scale worm parasitism on hosts are unclear. Previous studies have revealed that the microbiota plays an important role in host adaptability. Here, we compared gill-microbiota, gene expression and host-microorganism interactions in a group of deep sea mussels (Gigantidas haimaensis) parasitized by scale worm (PA group) and a no parasitic control group (NPA group). The symbiotic microorganism diversity of the PA group significantly decreased than NPA group, while the relative abundance of chemoautotrophic symbiotic bacteria that provide the host with organic carbon compounds significantly increased in PA. Interestingly, RNA-seq revealed that G. haimaensis hosts responded to B. pettiboneaei parasitism through significant upregulation of protein and lipid anabolism related genes, and that this parasitism may enhance host mussel nutrient anabolism but inhibit the host's ability to absorb nutrients, thus potentially helping the parasite obtain nutrients from the host. In an integrated analysis of the interactions between changes in the microbiota and host gene dysregulation, we found an agreement between the microbiota and transcriptomic responses to B. pettiboneaei parasitism. Together, our findings provide new insights into the effects of parasite scale worms on changes in symbiotic bacteria and gene expression in deep sea mussel hosts. We explored the potential role of host-microorganism interactions between scale worms and deep sea mussels, and revealed the mechanisms through which scale worm parasitism affects hosts in deep sea chemosynthetic ecosystem.

Corrigendum: Effects of scale worm parasitism on interactions between the symbiotic gill microbiome and gene regulation in deep sea mussel hosts.

[This corrects the article DOI: 10.3389/fmicb.2022.940766.].

[Identification and analysis of the TALE transcription factor family in radish].

Three-amino acid loop extension (TALE) transcription factors play important roles in plant growth and cell differentiation. There are plenty of studies on TALE transcription factors in several model plants, but not in radish (Raphanus sativas). A genome-wide bioinformatics analysis identified 33 TALE family genes in the Xiang-Ya-Bai (XYB) radish, These genes, are distributed on nine chromosomes and all contain 4-6 exons. The 33 TALE genes in radish showed a co-linearity relationship with the 17 homologous genes in Arabidopsis thaliana. Moreover, a large number of stress response cis-elements were found in the promoter regions of these genes. Expression analysis showed that four genes in the BELL subfamily were highly expressed in roots, and two genes in the KNOX subfamily were highly expressed in shoots of bolting plants and callus. All radish TALE genes contain sequences encoding the conserved HOX domain, except for the gene RSA10037940, which is homologous to Arabidopsis KNATM. The deduced 3D structures of the TALE proteins irrespective of subtypes are highly similar. All the encoded proteins were weakly acidic and hydrophilic. The radish TALE gene family is relatively evolutionarily conserved, which was consistent with results from Arabidopsis, but quite different from that of rice. This study provides important clues for studying the biological functions of TALE transcription factors in radish.