The genome of Lespedeza potaninii reveals biased subgenome evolution and drought adaptation

Lespedeza potaninii, a xerophytic subshrub belonging to the legume family, is native to the Tengger Desert and is highly adapted to drought. It has important ecological value due to its drought adaptability, but the underlying molecular mechanisms remain largely unknown. Here, we report a 1.24 Gb chromosome-scale assembly of the L. potaninii genome (contig N50=15.75 Mb). Our results indicate that L. potaninii underwent an allopolyploid event with two subgenomes, A and B, presenting asymmetric evolution and B subgenome dominance. We estimate that the two diploid progenitors of L. potaninii diverged around 3.6 MYA and merged around 1.0 MYA. We revealed that the expansion of hub genes associated with drought responses, such as the binding partner 1 of accelerated cell death 11 (ACD11) (BPA1), facilitated environmental adaptations of L. potaninii to desert habitats. We found a novel function of the BPA1 family in abiotic stress tolerance in addition to the known role in regulating the plant immune response, which could improve drought tolerance by positively regulating reactive oxygen species homeostasis in plants. We revealed that bZIP transcription factors could bind to the BPA1 promoter and activate its transcription. Our work fills the genomic data gap in the Lespedeza genus and the tribe Desmodieae, which should provide both theoretical support in the study of drought tolerance and in the molecular breeding of legume crops.

Expansion and diversity of caspases in Mytilus coruscus contribute to larval metamorphosis and environmental adaptation

BACKGROUND

Apoptosis is involved (directly and indirectly) in several physiological processes including tissue remodeling during the development, the turnover of immune cells, and a defense against harmful stimuli. The disordered apoptotic process participates in the pathogenesis of various diseases, such as neoplasms, and chronic inflammatory or systemic autoimmune diseases, which are associated with its inadequate regulation. Caspases are vital components of the apoptotic pathway that are involved in developmental and immune processes. However, genome-wide identification and functional analysis of caspase have not been conducted in Mytilus coruscus, which is an economically important bivalve.

RESULTS

Here, 47 caspase genes were identified from the genomes of M. coruscus, and the expansion of caspase-2/9 and caspase-3/6/7 genes were observed. Tandem duplication acts as an essential driver of gene expansion. The expanded caspase genes were highly diverse in terms of sequence, domain structure, and spatiotemporal expression profiles, suggesting their functional differentiation. The high expression of the expanded caspase genes at the pediveliger larvae stage and the result of apoptosis location in the velum suggest that the apoptosis mediated by them plays a critical role in the metamorphosis of M. coruscus larvae. In gill, caspase genes respond differently to the challenge of different strains, and most caspase-2/9 and caspase-3/6/7 genes were induced by copper stress, whereas caspase-8/10 genes were suppressed. Additionally, most caspase genes were upregulated in the mantle under ocean acidification which could weaken the biomineralization capacity of the mantle tissue.

CONCLUSIONS

These results provide a comprehensive overview of the evolution and function of the caspase family and enhanced the understanding of the biological function of caspases in M. coruscus larval development and response to biotic and abiotic challenges.

Do cord blood CD34 cells stay the same or change after expansion: A comprehensive analysis of gene expression results before and after seven days of cell culture

INTRODUCTION

Hematopoietic stem cell transplantation is used to treat hematopoietic malignancies and bone marrow failure syndromes. Due to the difficulties of using these cells isolated from the bone marrow, an additional source for receiving essential hematopoietic cells is umbilical cord blood. But the main problem with using the umbilical cord is its insufficient blood volume. The ex-vivo reproductive system of hematopoietic stem cells can overcome this problem that has the ability for Transplantation and hematopoiesis in the long term. This study aimed to evaluate the expression profile of CD133 + umbilical cord blood microRNAs in different stages of hematopoietic stem cells before and after ex-vivo proliferation.

MATERIALS AND METHODS

The expression profile of 1034 types of microRNA of CD34 + hematopoietic stem cells of the umbilical cord was analyzed before and after ex vivo proliferation. After isolating CD34 + cells from the umbilical cord blood of normal-born babies using MACS (Magnetic-activated cell sorting) column, these cells were cultured in a Stem span culture medium containing SCF, FLT3, and TPO cytokines in 24-well plates. The expression profile of microRNAs was investigated on days 0 and 7 days after cultivation by the Real-Time PCR method.

RESULTS

The results showed that the production of two-thirds of the micro-RNAs was reduced during the proliferation process. The micro-RNA expression of hematopoietic stem cell proliferation was also lower. At the same time, micro-RNA expression related to differentiation was higher.

CONCLUSION

The observed reduction in miRNA expression may be attributed to enhanced differentiation through proliferation. Therefore, miRNAs appear to be a viable option for regulating the proliferation processes of Hematopoietic stem cell transplantation.

Genome-wide identification and expression analyses of Toll-like receptors provide new insights on adaptation to intertidal benthic environments in Urechis unicinctus (Annelida, Echiura)

Toll-like receptors (TLR) are an important class of pattern recognition receptors involved in innate immunity that recognize pathogen-associated and damage-associated molecular patterns. Although the role of TLRs in immunity has been extensively studied, a systematic investigation of their function in environmental adaptation is still in its infancy. In this study, a genome-wide search was conducted to systematically investigate TLR family members of Urechis unicinctus, a typical benthic organism in intertidal mudflats. A total of 28 TLR genes were identified in the U. unicinctus genome, and their fundamental physiological and biochemical properties were characterized. Gene copy number analysis among species in different habitats indicated that TLR family gene expansion may be probably related with benthic environmental adaptation. To further investigate the expression patterns of TLR members under environmental stress, transcriptome data was analyzed from different developmental stages and the hindgut under sulfide stress. Transcriptome analysis of different developmental stages showed that most TLR genes were highly expressed during a key period of benthic environment adaptation (worm-shaped larva). Transcriptome analysis of the hindgut under sulfide stress showed that the expression of 12 TLR members was significantly induced under sulfide stress. These results indicate that the regulation of TLR gene expression may be probably involved in the adaptation of U. unicinctus to the benthic intertidal zone environment. Taken together, this study may lay the foundation for future functional analysis of the specific role of TLRs in host immune responses against sulfide exposure and benthic environmental stress in annelid.

Massive expansion of P-selectin genes in two Venerida species, Sinonovacula constricta and Mercenaria mercenaria: evidence from comparative genomics of Bivalvia

Background

P-selectin is a molecule participating in the inflammatory response through mediating cellular adhesion and essential for wound repair. However, studies regarding P-selectin in Bivalvia are rare. This study identified 90 P-selectin genes among nine bivalve genomes and classified them into 4 subfamilies according to phylogenetic analysis.

Results

Notable P-selectin gene expansion was observed in two Venerida species, Sinonovacula constricta and Mercenaria mercenaria. The synteny analysis revealed that P-selectin gene expansion was mostly caused by tandem duplication. In addition, the expression profiles of P-selectin genes in S. constricta showed that many P-selectins were specifically highly expressed in the gills, and the P-selectin expression patterns changed dramatically under low salt stress and ammonia nitrogen stress.

Conclusions

The massive expansion of P-selectins may facilitate the tolerance to environmental stresses. This study sheds light on the characterizations and expression profiles of P-selectin genes in Bivalvia and provides an integrated framework for further investigation of the role of P-selectins in the environmental tolerance of bivalves.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08861-6.

Genome-wide identification, characterization and expression analysis of glutaredoxin gene family (Grxs) in Phaseolus vulgaris

Glutaredoxins (Grxs) are ubiquitous oxidoreductase proteins implicated in development and abiotic stress response mainly through maintaining redox homoeostasis. Here, we conducted the first systematic analysis of the Grx gene family (PvGrx) in the most popular legume Phaseolus vulgaris (common bean). A total of 50 PvGrx genes were identified, and divided into four classes (CC-type, CGFS-type, CPYC-type and Grl-type) based on the phylogenetic analysis. The different classes have different introns-exons structures and conserved motifs, indicating functional divergence in the PvGrx family. Both tandem and segmental duplications were found to be involved in the expansion of PvGrx family that underwent a purifying selection by excluding the deleterious loss-of-function mutations. Cis-acting regulatory elements and gene ontology analyses predicted their role of distinctive members in abiotic stress response and hormonal signalling. RNA-seq based expression analysis revealed their differential expression pattern during plant development. On the other hand, RT q-PCR analysis revealed that target PvGrx isoforms were associated with nodule organogenesis and symbiosis based on their expression profiles. In addition, a battery of PvGrx candidates were markedly upregulated by different abiotic stressors suggesting their broad spectrum of functions. These findings serve as a reference for functional analysis and genetic improvement in P. vulgaris and related legume species.

Identification and expression analysis of phosphate transporter (PHT) gene family in Lupinus albus cluster root under phosphorus stress

According to global estimation, 5.7 billion hectares of agricultural land contain limited phosphorus (P) availability leading to insufficient plant growth and productivity. Internal phosphate transporters play an essential role in mediating P mobilization and uptake from the soil. White lupin (Lupinus albus) is a cluster root (CR) forming crop with great potential to survive under P limited soil. However, it is imperative to identify and characterize the phosphate transporter (PHT) gene family in plants to validate their involvement in solving P deficiency problems. The recent availability of white lupin high-quality genome allowed us an exhaustive searches in the whole genome and identified five phosphates transporters subfamilies, including 35 putative genes that are unevenly distributed on 16 chromosomes. The LaPHT1 subfamily contained eight genes, LaPHT2 subfamily have three, LaPHT3 subfamily have eight, LaPHT4 subfamily have nine, and LaPHO subfamily has seven. Gene structure and duplication were also examined in detail. Syntenic analysis revealed that white lupin PHT family members had maximum the collinear relationship with those in L. angustifolius followed by Phaseolus vulgaris but showed the least collinear relationship with those in Arabidopsis. Gene ontology (GO) analysis revealed that the in white lupin PHT genes were enriched in functions regulated P uptake, transport, and recycling mechanisms. RT-qPCR was performed to evaluate the transcript levels of LaPHT genes in different parts of CR under P deficient hydroponic culture. Our study would provide better understanding the genetic evolution and expression phosphate of phosphate transporters in L. albus CR under P deficiency. It will also be helpful for further functional-based studies to solve P deficiency-related issues and mitigate P stress responses.

Gene duplication and subsequent functional diversification of maltase in fig wasp (Chalcidoidea, Hymenoptera)

Maltase can catalyze the hydrolysis of α-1,4-glucosidic linkages and release α-d-glucoses that are used as a source of energy by insects. Maltase has been extensively studied in Lepidoptera and Diptera, while the characterization and evolutionary history of maltase are largely unknown in Hymenoptera. Here, we undertook a bioinformatics study and identified 105 maltase genes in 12 fig wasp species. Together with the maltase genes of Nasonia vitripennis and Apis mellifera, phylogenetic analysis showed that all the maltase genes were clustered into three clades. Clade I and III included maltase genes from all the fig wasp species, while clade II contained the maltase genes from non-pollinating fig wasps (NPFWs) only. Interestingly, the maltase genes located in clade II were intronless. Fig pollinators and NPFWs had lineage-specific gene expansion in clade I and II respectively, which were mainly derived from tandem duplications. The three clades displayed distinct gene structures. Furthermore, maltase showed significant functional divergence among the three clades and the critical amino acid sites were detected. These sites could be responsible for the ligand-binding preference and hydrolytic specificity. Overall, our results demonstrated that maltase might contribute to the discrepancy of life histories and feeding regimes between fig pollinators and NPFWs.

Toxin- and species-dependent regulation of ATP-binding cassette (ABC) transporters in scallops after exposure to paralytic shellfish toxin-producing dinoflagellates

ATP-binding cassette (ABC) transporters are membrane-bound proteins involved in exporting various xenobiotic compounds from living cells. Bivalve mollusks can accumulate large amounts of paralytic shellfish toxins (PSTs) from marine dinoflagellates. For aquatic invertebrates, the importance of ABC proteins in multi-xenobiotic resistance has been demonstrated, however, the systematic identification of ABC transporters is very limited. In this study, 64 and 67 ABC genes containing all eight described subfamilies (A to H) were identified in Yesso scallop (Patinopecten yessoensis) and Zhikong scallop (Chlamys farreri), respectively, with massive gene expansion being observed in the ABCC and ABCG subfamilies. The kidney harbored more specifically expressed ABC genes than other organs/tissues, most of which belonged to ABCB, ABCC, and ABCG subfamilies. After feeding the scallops with PST-producing dinoflagellates, the expression of scallop ABC genes in the kidney was regulated in toxin- and species-dependent manners. In total, 20 and 24 ABC genes in Zhikong scallop (CfABCs) were induced after exposure to Alexandrium minutum and A. catenella, with the up-regulated members from both ABCC and ABCG subfamilies mainly showing acute and chronic induction by A. minutum and A. catenella, respectively, while the up-regulated CfABCBs mainly showing chronic induction by both dinoflagellates. In Yesso scallop, only eight ABC genes (PyABCs) were regulated after A. catenella exposure, and all the five up-regulated PyABCs were acutely induced. Our findings imply the functional diversity of scallop ABC genes in coping with PST accumulation, which may contribute to the lineage-specific adaptation of scallops for dealing with algal toxins challenge.

Molecular evolution of melatonin receptor genes (mtnr) in vertebrates and its shedding light on mtnr1c

Melatonin receptors (MTNRs) play important roles in regulation of circadian rhythms and seasonal reproduction. However, their origin and evolution in vertebrates have not been investigated. Here, we performed a comprehensive examination by comparative genome mining of MTNRs in vertebrates. We successfully extracted 164 putative encoding sequences for MTNRs (including 57 mtnr1a, 59 mtnr1b and 48 mtnr1c) from 45 high-quality representative genomes. Interestingly, the putative expansions of mtnr1a and mtnr1b in zebrafish were also identified in other Cyprinifomes, but not in other orders of teleost. Using phylogenetic interference, we observed this expansion to be clustered into a primitive position of the Actinopterygii, which may be resulted from teleost-specific genome duplication. The C-terminal extension of MTNR1C, predicted to be proteoglycan 4 (PRG4), originated after the speciation of Monotremata or Marsupialia. Our present genomics survey provides novel insights into the evolution of MTNRs in vertebrates and updates our understanding of these proteins.

Distinct evolutionary trajectories of V1R clades across mouse species

BACKGROUND

Many animals rely heavily on olfaction to navigate their environment. Among rodents, olfaction is crucial for a wide range of social behaviors. The vomeronasal olfactory system in particular plays an important role in mediating social communication, including the detection of pheromones and recognition signals. In this study we examine patterns of vomeronasal type-1 receptor (V1R) evolution in the house mouse and related species within the genus Mus. We report the extent of gene repertoire turnover and conservation among species and clades, as well as the prevalence of positive selection on gene sequences across the V1R tree. By exploring the evolution of these receptors, we provide insight into the functional roles of receptor subtypes as well as the dynamics of gene family evolution.

RESULTS

We generated transcriptomes from the vomeronasal organs of 5 Mus species, and produced high quality V1R repertoires for each species. We find that V1R clades in the house mouse and relatives exhibit distinct evolutionary trajectories. We identify putative species-specific gene expansions, including a large clade D expansion in the house mouse. While gene gains are abundant, we detect very few gene losses. We describe a novel V1R clade and highlight candidate receptors for future study. We find evidence for distinct evolutionary processes across different clades, from largescale turnover to highly conserved repertoires. Patterns of positive selection are similarly variable, as some clades exhibit abundant positive selection while others display high gene sequence conservation. Based on clade-level evolutionary patterns, we identify receptor families that are strong candidates for detecting social signals and predator cues. Our results reveal clades with receptors detecting female reproductive status are among the most conserved across species, suggesting an important role in V1R chemosensation.

CONCLUSION

Analysis of clade-level evolution is critical for understanding species' chemosensory adaptations. This study provides clear evidence that V1R clades are characterized by distinct evolutionary trajectories. As receptor evolution is shaped by ligand identity, these results provide a framework for examining the functional roles of receptors.

Genome of the webworm Hyphantria cunea unveils genetic adaptations supporting its rapid invasion and spread

BACKGROUND

The fall webworm Hyphantria cunea is an invasive and polyphagous defoliator pest that feeds on nearly any type of deciduous tree worldwide. The silk web of H. cunea aids its aggregating behavior, provides thermal regulation and is regarded as one of causes for its rapid spread. In addition, both chemosensory and detoxification genes are vital for host adaptation in insects.

RESULTS

Here, a high-quality genome of H. cunea was obtained. Silk-web-related genes were identified from the genome, and successful silencing of the silk protein gene HcunFib-H resulted in a significant decrease in silk web shelter production. The CAFE analysis showed that some chemosensory and detoxification gene families, such as CSPs, CCEs, GSTs and UGTs, were expanded. A transcriptome analysis using the newly sequenced H. cunea genome showed that most chemosensory genes were specifically expressed in the antennae, while most detoxification genes were highly expressed during the feeding peak. Moreover, we found that many nutrient-related genes and one detoxification gene, HcunP450 (CYP306A1), were under significant positive selection, suggesting a crucial role of these genes in host adaptation in H. cunea. At the metagenomic level, several microbial communities in H. cunea gut and their metabolic pathways might be beneficial to H. cunea for nutrient metabolism and detoxification, and might also contribute to its host adaptation.

CONCLUSIONS

These findings explain the host and environmental adaptations of H. cunea at the genetic level and provide partial evidence for the cause of its rapid invasion and potential gene targets for innovative pest management strategies.

Expansion of sweet taste receptor genes in grass carp (Ctenopharyngodon idellus) coincided with vegetarian adaptation

BACKGROUND

Taste is fundamental to diet selection in vertebrates. Genetic basis of sweet taste receptor in the shaping of food habits has been extensively studied in mammals and birds, but scarcely studied in fishes. Grass carp is an excellent model for studying vegetarian adaptation, as it exhibits food habit transition from carnivory to herbivory.

RESULTS

We identified six sweet taste receptors (gcT1R2A-F) in grass carp. The four gcT1R2s (gcT1R2C-F) have been suggested to be evolved from and paralogous to the two original gcT1R2s (gcT1R2A and gcT1R2B). All gcT1R2s were expressed in taste organs and mediated glucose-, fructose- or arginine-induced intracellular calcium signaling, revealing they were functional. In addition, grass carp was performed to prefer fructose to glucose under a behavioral experiment. Parallelly, compared with gcT1R2A-F/gcT1R3 co-transfected cells, gcT1R2C-F/gcT1R3 co-transfected cells showed a higher response to plant-specific fructose. Moreover, food habit transition from carnivory to herbivory in grass carp was accompanied by increased gene expression of certain gcT1R2s.

CONCLUSIONS

We suggested that the gene expansion of T1R2s in grass carp was an adaptive strategy to accommodate the change in food environment. Moreover, the selected gene expression of gcT1R2s might drive the food habit transition from carnivory to herbivory in grass carp. This study provided some evolutional and physiological clues for the formation of herbivory in grass carp.

Massive expansion and diversity of nicotinic acetylcholine receptors in lophotrochozoans

Background

Nicotinic acetylcholine receptors (nAChRs) are among the oldest and most conserved transmembrane receptors involved in signal transduction. Despite the prevalence and significance of cholinergic signaling, the diversity and evolution of nAChRs are not fully understood.

Result

By comparative genomic analysis, we found massive expansions of nAChR genes in molluscs and some other lophotrochozoans. The expansion is particularly pronounced in stationary bivalve molluscs with simple nervous systems, with the number of nAChR genes ranging from 99 to 217 in five bivalves, compared with 10 to 29 in five ecdysozoans and vertebrates. The expanded molluscan nAChR genes tend to be intronless and in tandem arrays due to retroposition followed by tandem duplication. Phylogenetic analysis revealed diverse nAChR families in the common ancestor of bilaterians, which subsequently experienced lineage-specific expansions or contractions. The expanded molluscan nAChR genes are highly diverse in sequence, domain structure, temporal and spatial expression profiles, implying diversified functions. Some molluscan nAChR genes are expressed in early development before the development of the nervous system, while others are involved in immune and stress responses.

Conclusion

The massive expansion and diversification of nAChR genes in bivalve molluscs may be a compensation for reduced nervous systems as part of adaptation to stationary life under dynamic environments, while in vertebrates a subset of specialized nAChRs are retained to work with advanced nervous systems. The unprecedented diversity identified in molluscs broadens our view on the evolution and function of nAChRs that are critical to animal physiology and human health.

Diverse expression regulation of Hsp70 genes in scallops after exposure to toxic Alexandrium dinoflagellates

Heat shock proteins 70KD (Hsp70s) are highly conserved molecular chaperones with essential roles against biotic and abiotic stressors. Marine bivalves inhabit highly complex environments and could accumulate paralytic shellfish toxins (PSTs), the well-noted neurotoxins generated during harmful algal blooms. Here, we systematically analyzed Hsp70 genes (CfHsp70s) in Zhikong scallop (Chlamys farreri), an important aquaculture mollusk in China. Sixty-five CfHsp70s from eight sub-families were identified, and 47 of these genes showed expansion in the Hspa12 sub-family. After exposure to different PST-producing dinoflagellates, Alexandrium minutum and Alexandrium catenella, diverse CfHsp70s regulation presented in scallop hepatopancreas, mainly accumulating incoming PSTs, and kidneys, transforming PSTs into higher toxic analogs. All the up-regulated CfHsp70s were from CfHsp70B2, CfHspa12, and CfHspa5 sub-families. CfHsp70B2 sub-family was mainly induced in the hepatopancreas, and CfHspa12 sub-family was highly induced in the kidneys. CfHsp70s up-regulation under two dinoflagellates exposure was stronger in the kidneys (log₂FC: 19.5 and 18.6) than that in hepatopancreas (log₂FC: 4.3 and 6.1). Exposure to different Alexandrium species had varying effects, that in hepatopancreas, CfHsp70B2s were chronically induced only after A. catenella exposure, whereas in kidney, CfHspa12s were more acutely induced after exposure of A. minutum than A. caenella. Moreover, in Yesso scallops (Patinopecten yessoensis), only Hspa12s were up-regulated in hepatopancreas after A. catenella exposure, and all the Hsp70B2s were down-regulated. These organ-, toxin-, and species-dependent Hsp70 regulation suggested the functional diversity of duplicated Hsp70s in response to the stress by PST-producing algae. Our findings provide insights into the evolution and functional characteristics of Hsp70s in scallops.

Genome-wide analysis of the Tetrahymena thermophila glutathione S-transferase gene superfamily

The ciliate Tetrahymena thermophila has a rapid response to detoxify xenobiotics, which presents opportunity to study the diversification of Glutathione S-Transferase superfamily. In-silico identification of putative GST genes were resulted with 70 GST genes; 49 TtGSTmu, 7 TtGSTomega, 5 TtGSTtheta, 2 TtGSTzeta, 4 TtMAPEG and 3 TtEF1G. TtGST superfamily has short intron carrying or intronless genes. The most expressed mRNAs of TtGST are limited to 4 members at all life stages. TtGST genes are widely distributed to all five micronuclear chromosomes with the highest diversified members from different classes in chromosome 4. The clustering and the orientation of some TtGSTs in the T. thermophila genome give clues about the recent gene duplication. Analysis of GSH affinity-purified GST proteins with Western blot and activity assay showed GST activity carrying purified TtGST populations. In conclusion, the enhanced genome capacity of TtGST superfamily may have evolved through improved GST enzymatic activity.

Genome-wide identification and characterization of the Dof gene family in cassava (Manihot esculenta)

Dof proteins comprise a plant-specific transcription factor family involved in regulation of many plant processes. Despite their diverse roles, little information is available in cassava (Manihot esculenta, Euphorbiaceae), the sixth major staple crop in the world. This study presents the first genome-wide identification of Dof family genes in cassava, resulting in 45 family members that are unevenly distributed across 17 out of the 18 chromosomes. Phylogenetic analysis assigned these genes into eight groups representing four subfamilies. Synteny analysis and BRH (Best Reciprocal Hit) method were also adopted to study evolutionary relationships. Among 21 MeDof duplicates identified, the majority of them were derived from the recent whole-genome duplication (WGD), whereas only one from local duplication. Transcriptional profiling revealed diverse expression patterns among various tissues, and further comparison of conserved motifs and expression profiles revealed divergence of MeDof duplicates. These findings provide valuable information for functional analysis of Dof genes in cassava.

Expansion and stress responses of the AP2/EREBP superfamily in cotton

Background

The allotetraploid cotton originated from one hybridization event between an extant progenitor of Gosssypium herbaceum (A₁) or G. arboreum (A₂) and another progenitor, G. raimondii Ulbrich (D₅) 1–1.5 million years ago (Mya). The APETALA2/ethylene-responsive element binding protein (AP2/EREBP) transcription factors constitute one of the largest and most conserved gene families in plants. They are characterized by their AP2 domain, which comprises 60–70 amino acids, and are classified into four main subfamilies: the APETALA2 (AP2), Related to ABI3/VP1 (RAV), Dehydration-Responsive Element Binding protein (DREB) and Ethylene-Responsive Factor (ERF) subfamilies. The AP2/EREBP genes play crucial roles in plant growth, development and biotic and abiotic stress responses. Hence, understanding the molecular characteristics of cotton stress tolerance and gene family expansion would undoubtedly facilitate cotton resistance breeding and evolution research.

Results

A total of 269 AP2/EREBP genes were identified in the G. raimondii (D5) cotton genome. The protein domain architecture and intron/exon structure are simple and relatively conserved within each subfamily. They are distributed throughout all chromosomes but are clustered on various chromosomes due to genomic tandem duplication. We identified 73 tandem duplicated genes and 221 segmental duplicated gene pairs which contributed to the expansion of AP2/EREBP superfamily. Of them, tandem duplication was the most important force of the expansion of the B3 group. Transcriptome analysis showed that 504 AP2/EREBP genes were expressed in at least one tested G. hirsutum TM-1 tissues. In G. hirsutum, 151 non-repeated genes of the DREB and ERF subfamily genes were responsive to different stresses: 132 genes were induced by cold, 63 genes by drought and 94 genes by heat. qRT-PCR confirmed that 13 GhDREB and 15 GhERF genes were induced by cold and/or drought. No transcripts detected for 53 of the 111 tandem duplicated genes in TM-1. In addition, some homoeologous genes showed biased expression toward either A-or D-subgenome.

Conclusions

The AP2/EREBP genes were obviously expanded in Gossypium. The GhDREB and GhERF genes play crucial roles in cotton stress responses. Our genome-wide analysis of AP2/EREBP genes in cotton provides valuable information for characterizing the molecular functions of AP2/EREBP genes and reveals insights into their evolution in polyploid plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3517-9) contains supplementary material, which is available to authorized users.

Hsp70 gene expansions in the scallop Patinopecten yessoensis and their expression regulation after exposure to the toxic dinoflagellate Alexandrium catenella

Heat shock protein 70 (Hsp70s) family members are present in virtually all living organisms and perform a fundamental role against different types of environmental stressors and pathogenic organisms. Marine bivalves live in highly dynamic environments and may accumulate paralytic shellfish toxins (PSTs), a class of well-known neurotoxins closely associated with harmful algal blooms (HABs). Here, we provide a systematic analysis of Hsp70 genes (PyHsp70s) in the genome of Yesso scallop (Patinopecten yessoensis), an important aquaculture species in China, through in silico analysis using transcriptome and genome databases. Phylogenetic analyses indicated extensive expansion of Hsp70 genes from the Hspa12 sub-family in the Yesso scallop and also the bivalve lineages, with gene duplication events before or after the split between the Yesso scallop and the Pacific oyster. In addition, we determined the expression patterns of PyHsp70s after exposure to Alexandrium catenella, the dinoflagellate producing PSTs. Our results confirmed the inducible expression patterns of PyHsp70s under PSTs stress, and the responses to the toxic stress may have arisen through the adaptive recruitment of tandem duplication of Hsp70 genes. These findings provide a thorough overview of the evolution and modification of the Hsp70 family, which will gain insights into the functional characteristics of scallop Hsp70 genes in response to different stresses.

Genomic analyses and expression evaluation of thaumatin-like gene family in the cacao fungal pathogen Moniliophthora perniciosa

Thaumatin-like proteins (TLPs) are found in diverse eukaryotes. Plant TLPs, known as Pathogenicity Related Protein (PR-5), are considered fungal inhibitors. However, genes encoding TLPs are frequently found in fungal genomes. In this work, we have identified that Moniliophthora perniciosa, a basidiomycete pathogen that causes the Witches' Broom Disease (WBD) of cacao, presents thirteen putative TLPs from which four are expressed during WBD progression. One of them is similar to small TLPs, which are present in phytopathogenic basidiomycete, such as wheat stem rust fungus Puccinia graminis. Fungi genomes annotation and phylogenetic data revealed a larger number of TLPs in basidiomycetes when comparing with ascomycetes, suggesting that these proteins could be involved in specific traits of mushroom-forming species. Based on the present data, we discuss the contribution of TLPs in the combat against fungal competitors and hypothesize a role of these proteins in M. perniciosa pathogenicity.