A genome-wide survey of segmental duplications that mediate common human genetic variation of chromosomal architecture

Genome-Wide Identification of the Peanut ASR Gene Family and Its Expression Analysis under Abiotic Stress

Peanut (Arachis hypogaea L.) is one of the most important oil and food legume crops worldwide. ASR (abscisic acid, stress, ripening) plays extremely important roles in plant growth and development, fruit ripening, pollen development, and stress. Here, six ASR genes were identified in peanut. Structural and conserved motif analyses were performed to identify common ABA/WDS structural domains. The vast majority of ASR genes encoded acidic proteins, all of which are hydrophilic proteins and localized on mitochondria and nucleus, respectively. The cis-element analysis revealed that some cis-regulatory elements were related to peanut growth and development, hormone, and stress response. Under normal growth conditions, AhASR4 and AhASR5 were expressed in all tissues of peanut plants. Quantitative real-time PCR (qRT-PCR) results indicated that peanut ASR genes exhibited complex expression patterns in response to abiotic stress. Notably, under drought and cadmium (Cd) stress, the expression levels of AhASR4 and AhASR5 were significantly upregulated, suggesting that these genes may play a crucial role in the peanut plant's resistance to such stressors. These results provide a theoretical basis for studying the evolution, expression, and function of the peanut ASR gene family and will provide valuable information in the identification and screening of genes for peanut stress tolerance breeding.

Genome-wide identification and expression analysis of the SPL gene family and its response to abiotic stress in barley (Hordeum vulgare L.)

BACKGROUND

Squamosa promoter-binding protein-like (SPL) is a plant-specific transcription factor that is widely involved in the regulation of plant growth and development, including flower and grain development, stress responses, and secondary metabolite synthesis. However, this gene family has not been comprehensively evaluated in barley, the most adaptable cereal crop with a high nutritional value.

RESULTS

In this study, a total of 15 HvSPL genes were identified based on the Hordeum vulgare genome. These genes were named HvSPL1 to HvSPL15 based on the chromosomal distribution of the HvSPL genes and were divided into seven groups (I, II, III, V, VI, VII, and VIII) based on the phylogenetic tree analysis. Chromosomal localization revealed one pair of tandem duplicated genes and one pair of segmental duplicated genes. The HvSPL genes exhibited the highest collinearity with the monocotyledonous plant, Zea mays (27 pairs), followed by Oryza sativa (18 pairs), Sorghum bicolor (16 pairs), and Arabidopsis thaliana (3 pairs), and the fewest homologous genes with Solanum lycopersicum (1 pair). The distribution of the HvSPL genes in the evolutionary tree was relatively scattered, and HvSPL proteins tended to cluster with SPL proteins from Z. mays and O. sativa, indicating a close relationship between HvSPL and SPL proteins from monocotyledonous plants. Finally, the spatial and temporal expression patterns of the 14 HvSPL genes from different subfamilies were determined using quantitative real-time polymerase chain reaction (qRT-PCR). Based on the results, the HvSPL gene family exhibited tissue-specific expression and played a regulatory role in grain development and abiotic stress. HvSPL genes are highly expressed in various tissues during seed development. The expression levels of HvSPL genes under the six abiotic stress conditions indicated that many genes responded to stress, especially HvSPL8, which exhibited high expression under multiple stress conditions, thereby warranting further attention.

CONCLUSION

In this study, 15 SPL gene family members were identified in the genome of Hordeum vulgare, and the phylogenetic relationships, gene structure, replication events, gene expression, and potential roles of these genes in millet development were studied. Our findings lay the foundation for exploring the HvSPL genes and performing molecular breeding of barley.

Genome-Wide Identification of NAC Family Genes and Their Expression Analyses in Response to Osmotic Stress in Cannabis sativa L

NAC (NAM, ATAF1/2, and CUC2) transcription factors are unique and essential for plant growth and development. Although the NAC gene family has been identified in a wide variety of plants, its chromosomal location and function in Cannabis sativa are still unknown. In this study, a total of 69 putative CsNACs were obtained, and chromosomal location analysis indicated that the CsNAC genes mapped unevenly to 10 chromosomes. Phylogenetic analyses showed that the 69 CsNACs could be divided into six subfamilies. Additionally, the CsNAC genes in group IV-a are specific to Cannabis sativa and contain a relatively large number of exons. Promoter analysis revealed that most CsNAC promoters contained cis-elements related to plant hormones, the light response, and abiotic stress. Furthermore, transcriptome expression profiling revealed that 24 CsNAC genes in two Cannabis sativa cultivars (YM1 and YM7) were significantly differentially expressed under osmotic stress, and these 12 genes presented differential expression patterns across different cultivars according to quantitative real-time PCR (RT-qPCR) analysis. Among these, the genes homologous to the CsNAC18, CsNAC24, and CsNAC61 genes have been proven to be involved in the response to abiotic stress and might be candidate genes for further exploration to determine their functions. The present study provides a comprehensive insight into the sequence characteristics, structural properties, evolutionary relationships, and expression patterns of NAC family genes under osmotic stress in Cannabis sativa and provides a basis for further functional characterization of CsNAC genes under osmotic stress to improve agricultural traits in Cannabis sativa.

Genome-Wide Analysis and Expression Profiling of Lectin Receptor-like Kinase Genes in Watermelon (Citrullus lanatus)

Watermelon is one of the most important edible plants worldwide. Owing to its special cultivation conditions, watermelon is exposed to many biological and abiotic stresses during its development. Lectin receptor-like kinases (LecRLKs) are plant-specific membrane proteins that play important roles in sensing and responding to environmental stimuli. Although the LecRLK gene family has been identified in a variety of plants, a comprehensive analysis has not yet been undertaken in watermelon. In this study, 61 putative LecRLK genes were identified in watermelon, consisting of 36 G-type, 24 L-type, and 1 C-type LecRLK genes. They were distributed in clusters on chromosomes, and members from the same subfamily were mostly clustered together. The analysis of the phylogenetic tree and conserved motif indicated that there were obvious differences among three ClaLecRLK subfamilies, and there was also rich diversity in the C-terminal within subfamilies. A collinear analysis revealed that the evolution of the ClaLecRLK gene family in different Cucurbitaceae crops was asynchronous. Furthermore, the analysis of the ClaLecRLK protein structure showed that not all proteins contained signal peptides and a single transmembrane domain. A subcellular localization assay confirmed that the number and position of transmembrane domains did not affect ClaLecRLK protein localization in cells. Transcriptome data revealed distinct expression patterns of LecRLK genes of watermelon in various tissues, and their responses to different fungi infection were also significantly different. Finally, the potential binding sites of the ClaLecRLK genes targeted by miRNA were predicted. This study enhances the understanding of the characteristics and functions of the LecRLK gene family in watermelon and opens up the possibility of exploring the roles that LecRLK genes may play in the life cycle of Cucurbitaceae plants.

Novel islands of GGC and GCC repeats coincide with human evolution

BACKGROUND

Because of high mutation rate, overrepresentation in genic regions, and link with various neurological, neurodegenerative, and movement disorders, GGC and GCC short tandem repeats (STRs) are prone to natural selection. Among a number of lacking data, the 3-repeats of these STRs remain widely unexplored.

RESULTS

In a genome-wide search in human, here we mapped GGC and GCC STRs of ≥3-repeats, and found novel islands of up to 45 of those STRs, populating spans of 1 to 2 kb of genomic DNA. RGPD4 and NOC4L harbored the densest (GGC)3 (probability 3.09061E-71) and (GCC)3 (probability 1.72376E-61) islands, respectively, and were human-specific. We also found prime instances of directional incremented density of STRs at specific loci in human versus other species, including the FOXK2 and SKI GGC islands. The genes containing those islands significantly diverged in expression in human versus other species, and the proteins encoded by those genes interact closely in a physical interaction network, consequence of which may be human-specific characteristics such as higher order brain functions.

CONCLUSION

We report novel islands of GGC and GCC STRs of evolutionary relevance to human. The density, and in some instances, periodicity of these islands support them as a novel genomic entity, which need to be further explored in evolutionary, mechanistic, and functional platforms.

Dyads of GGC and GCC form hotspot colonies that coincide with the evolution of human and other great apes

BACKGROUND

GGC and GCC short tandem repeats (STRs) are of various evolutionary, biological, and pathological implications. However, the fundamental two-repeats (dyads) of these STRs are widely unexplored.

RESULTS

On a genome-wide scale, we mapped (GGC)2 and (GCC)2 dyads in human, and found monumental colonies (distance between each dyad < 500 bp) of extraordinary density, and in some instances periodicity. The largest (GCC)2 and (GGC)2 colonies were intergenic, homogeneous, and human-specific, consisting of 219 (GCC)2 on chromosome 2 (probability < 1.545E-219) and 70 (GGC)2 on chromosome 9 (probability = 1.809E-148). We also found that several colonies were shared in other great apes, and directionally increased in density and complexity in human, such as a colony of 99 (GCC)2 on chromosome 20, that specifically expanded in great apes, and reached maximum complexity in human (probability 1.545E-220). Numerous other colonies of evolutionary relevance in human were detected in other largely overlooked regions of the genome, such as chromosome Y and pseudogenes. Several of the genes containing or nearest to those colonies were divergently expressed in human.

CONCLUSION

In conclusion, (GCC)2 and (GGC)2 form unprecedented genomic colonies that coincide with the evolution of human and other great apes. The extent of the genomic rearrangements leading to those colonies support overlooked recombination hotspots, shared across great apes. The identified colonies deserve to be studied in mechanistic, evolutionary, and functional platforms.

Genome-Wide Analysis of R2R3-MYB Genes and Functional Characterization of SmMYB75 in Eggplant Fruit Implications for Crop Improvement and Nutritional Enhancement

R2R3-MYB represents a substantial gene family that plays diverse roles in plant development. In this study, 102 SmR2R3-MYB genes were identified from eggplant fruit and classified into 31 subfamilies. Analysis indicated that segmental duplication events played a pivotal role in the expansion of the SmR2R3-MYB gene family. Furthermore, the prediction of miRNAs targeting SmR2R3-MYB genes revealed that 60 SmR2R3-MYBs are targeted by 57 miRNAs, with specific miRNAs displaying varying numbers of target genes, providing valuable insights into the regulatory functions of miRNAs in plant growth, development, and responses to stress conditions. Through expression profile analysis under various treatment conditions, including low temperature (4 °C), plant hormone (ABA, Abscisic acid), and drought stress (PEG, Polyethylene glycol), diverse and complex regulatory mechanisms governing SmR2R3-MYB gene expression were elucidated. Notably, EGP21875.1 and EGP21874.1 exhibited upregulation in expression under all treatment conditions. Transcriptome and metabolome analyses demonstrated that, apart from anthocyanins (delphinidin-3-O-glucoside, cyanidin-3-O-(6-O-p-coumaroyl)-glucoside, and malvidin-3-O-(6-O-p-coumaroyl)-glucoside), overexpression of SmMYB75 could also elevate the content of various beneficial compounds, such as flavonoids, phenolic acids, and terpenes, in eggplant pulp. This comprehensive study enhances our understanding of SmR2R3-MYB gene functions and provides a strong basis for further research on their roles in regulating anthocyanin synthesis and improving eggplant fruit quality.

Insights into ZmWAKL in maize kernel development: genome-wide investigation and GA-mediated transcription

BACKGROUND

The functional roles of the Wall Associated Kinase (WAK) and Wall Associated Kinase Like (WAKL) families in cellular expansion and developmental processes have been well-established. However, the molecular regulation of these kinases in maize development is limited due to the absence of comprehensive genome-wide studies.

RESULTS

Through an in-depth analysis, we identified 58 maize WAKL genes, and classified them into three distinct phylogenetic clusters. Moreover, structural prediction analysis showed functional conservation among WAKLs across maize. Promoter analysis uncovered the existence of cis-acting elements associated with the transcriptional regulation of ZmWAKL genes by Gibberellic acid (GA). To further elucidate the role of WAKL genes in maize kernels, we focused on three highly expressed genes, viz ZmWAKL38, ZmWAKL42 and ZmWAKL52. Co-expression analyses revealed that their expression patterns exhibited a remarkable correlation with GA-responsive transcription factors (TF) TF5, TF6, and TF8, which displayed preferential expression in kernels. RT-qPCR analysis validated the upregulation of ZmWAKL38, ZmWAKL42, ZmWAKL52, TF5, TF6, and TF8 following GA treatment. Additionally, ZmWAKL52 showed significant increase of transcription in the present of TF8, with ZmWAKL52 localizing in both the plasma membrane and cell wall. TF5 positively regulated ZmWAKL38, while TF6 positively regulated ZmWAKL42.

CONCLUSIONS

Collectively, these findings provide novel insights into the characterization and regulatory mechanisms of specific ZmWAKL genes involved in maize kernel development, offering prospects for their utilization in maize breeding programs.

Molecular evolution, diversification, and expression assessment of MADS gene family in Setaria italica, Setaria viridis, and Panicum virgatum

KEY MESSAGE

This paper sheds light on the evolution and expression patterns of MADS genes in Setaria and Panicum virgatum. SiMADS51 and SiMADS64 maybe involved in the ABA-dependent pathway of drought response. The MADS gene family is a key regulatory factor family that controls growth, reproduction, and response to abiotic stress in plants. However, the molecular evolution of this family is rarely reported. Here, a total of 265 MADS genes were identified in Setaria italica (foxtail millet), Setaria viridis (green millet), and Panicum virgatum (switchgrass) and analyzed by bioinformatics, including physicochemical characteristics, subcellular localization, chromosomal position and duplicate, motif distribution, genetic structure, genetic evolvement, and expression patterns. Phylogenetic analysis was used to categorize these genes into M and MIKC types. The distribution of motifs and gene structure were similar for the corresponding types. According to a collinearity study, the MADS genes have been mostly conserved during evolution. The principal cause of their expansion is segmental duplication. However, the MADS gene family tends to shrink in foxtail millet, green millet, and switchgrass. The MADS genes were subjected to purifying selection, but several positive selection sites were also identified in three species. And most of the promoters of MADS genes contain cis-elements related to stress and hormonal response. RNA-seq and quantitative Real-time PCR (qRT-PCR) analysis also were examined. SiMADS genes expression levels are considerably changed in reaction to various treatments, following qRT-PCR analysis. This sheds fresh light on the evolution and expansion of the MADS family in foxtail millet, green millet, and switchgrass, and lays the foundation for further research on their functions.

Gene Structural Specificity and Expression of MADS-Box Gene Family in Camellia chekiangoleosa

MADS-box genes encode transcription factors that affect plant growth and development. Camellia chekiangoleosa is an oil tree species with ornamental value, but there have been few molecular biological studies on the developmental regulation of this species. To explore their possible role in C. chekiangoleosa and lay a foundation for subsequent research, 89 MADS-box genes were identified across the whole genome of C. chekiangoleosa for the first time. These genes were present on all the chromosomes and were found to have expanded by tandem duplication and fragment duplication. Based on the results of a phylogenetic analysis, the 89 MADS-box genes could be divided into either type I (38) or type II (51). Both the number and proportion of the type II genes were significantly greater than those of Camellia sinensis and Arabidopsis thaliana, indicating that C. chekiangoleosa type II genes experienced a higher duplication rate or a lower loss rate. The results of both a sequence alignment and a conserved motif analysis suggest that the type II genes are more conserved, meaning that they may have originated and differentiated earlier than the type I genes did. At the same time, the presence of extra-long amino acid sequences may be an important feature of C. chekiangoleosa. Gene structure analysis revealed the number of introns of MADS-box genes: twenty-one type I genes had no introns, and 13 type I genes contained only 1~2 introns. The type II genes have far more introns and longer introns than the type I genes do. Some MIKC^C genes have super large introns (≥15 kb), which are rare in other species. The super large introns of these MIKC^C genes may indicate richer gene expression. Moreover, the results of a qPCR expression analysis of the roots, flowers, leaves and seeds of C. chekiangoleosa showed that the MADS-box genes were expressed in all those tissues. Overall, compared with that of the type I genes, the expression of the type II genes was significantly higher. The CchMADS31 and CchMADS58 genes (type II) were highly expressed specifically in the flowers, which may in turn regulate the size of the flower meristem and petals. CchMADS55 was expressed specifically in the seeds, which might affect seed development. This study provides additional information for the functional characterization of the MADS-box gene family and lays an important foundation for in-depth study of related genes, such as those involved in the development of the reproductive organs of C. chekiangoleosa.

Identification and characteristics of SnRK genes and cold stress-induced expression profiles in Liriodendron chinense

Background

The sucrose non-fermenting 1 (SNF1)-related protein kinases (SnRKs) play a vivid role in regulating plant metabolism and stress response, providing a pathway for regulation between metabolism and stress signals. Conducting identification and stress response studies on SnRKs in plants contributes to the development of strategies for tree species that are more tolerant to stress conditions.

Results

In the present study, a total of 30 LcSnRKs were identified in Liriodendron chinense (L. chinense) genome, which was distributed across 15 chromosomes and 4 scaffolds. It could be divided into three subfamilies: SnRK1, SnRK2, and SnRK3 based on phylogenetic analysis and domain types. The LcSnRK of the three subfamilies shared the same Ser/Thr kinase structure in gene structure and motif composition, while the functional domains, except for the kinase domain, showed significant differences. A total of 13 collinear gene pairs were detected in L. chinense and Arabidopsis thaliana (A. thaliana), and 18 pairs were detected in L. chinense and rice, suggesting that the LcSnRK family genes may be evolutionarily more closely related to rice. Cis-regulation element analysis showed that LcSnRKs were LTR and TC-rich, which could respond to different environmental stresses. Furthermore, the expression patterns of LcSnRKs are different at different times under low-temperature stress. LcSnRK1s expression tended to be down-regulated under low-temperature stress. The expression of LcSnRK2s tended to be up-regulated under low-temperature stress. The expression trend of LcSnRK3s under low-temperature stress was mainly up-or down-regulated.

Conclusion

The results of this study will provide valuable information for the functional identification of the LcSnRK gene in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08902-0.

Molecular characterization of SUT Gene Family in Solanaceae with emphasis on expression analysis of pepper genes during development and stresses

Sucrose, an essential carbohydrate, is transported from source to sink organs in the phloem and is involved in a variety of physiological and metabolic processes in plants. Sucrose transporter proteins (SUTs) may play significant parts in the phloem loading and unloading of sucrose. In our study, the SUT gene family was identified in four Solanaceae species (Capsicum annuum, Solanum lycopersicum, S. melongena, and S. tuberosum) and other 14 plant species ranged from lower and high plants. The comprehensive analysis was performed by integration of chromosomal distribution, gene structure, conserved motifs, evolutionary relationship and expression profiles during pepper growth under stresses. Chromosome mapping revealed that SUT genes in Solanaceae were distributed on chromosomes 4, 10 and 11. Gene structure analysis showed that the subgroup 1 members have the same number of introns and exons. All the SUTs had 12 transmembrane structural domains exception from CaSUT2 and SmSUT2, indicating that a structure variation might occurred among the Solanaceae SUT proteins. We also found a total of 20 conserved motifs, with over half of them shared by all SUT proteins, and the SUT proteins from the same subgroup shared common motifs. Phylogenetic analysis divided a total of 72 SUT genes in the plant species tested into three groups, and subgroup 1 might have diverged from a single common ancestor prior to the mono-dicot split. Finally, expression levels of CaSUTs were induced significantly under heat, cold, and salt treatments, indicating diverse functions of the CaSUTs to adapt to adverse environments.

Genome-wide identification, evolutionary and expression analyses of LEA gene family in peanut (Arachis hypogaea L.)

BACKGROUND

Late embryogenesis abundant (LEA) proteins are a group of highly hydrophilic glycine-rich proteins, which accumulate in the late stage of seed maturation and are associated with many abiotic stresses. However, few peanut LEA genes had been reported, and the research on the number, location, structure, molecular phylogeny and expression of AhLEAs was very limited.

RESULTS

In this study, 126 LEA genes were identified in the peanut genome through genome-wide analysis and were further divided into eight groups. Sequence analysis showed that most of the AhLEAs (85.7%) had no or only one intron. LEA genes were randomly distributed on 20 chromosomes. Compared with tandem duplication, segmental duplication played a more critical role in AhLEAs amplication, and 93 segmental duplication AhLEAs and 5 pairs of tandem duplication genes were identified. Synteny analysis showed that some AhLEAs genes come from a common ancestor, and genome rearrangement and translocation occurred among these genomes. Almost all promoters of LEAs contain ABRE, MYB recognition sites, MYC recognition sites, and ERE cis-acting elements, suggesting that the LEA genes were involved in stress response. Gene transcription analyses revealed that most of the LEAs were expressed in the late stages of peanut embryonic development. LEA3 (AH16G06810.1, AH06G03960.1), and Dehydrin (AH07G18700.1, AH17G19710.1) were highly expressed in roots, stems, leaves and flowers. Moreover, 100 AhLEAs were involved in response to drought, low-temperature, or Al stresses. Some LEAs that were regulated by different abiotic stresses were also regulated by hormones including ABA, brassinolide, ethylene and salicylic acid. Interestingly, AhLEAs that were up-regulated by ethylene and salicylic acid showed obvious subfamily preferences. Furthermore, three AhLEA genes, AhLEA1, AhLEA3-1, and AhLEA3-3, which were up-regulated by drought, low-temperature, or Al stresses was proved to enhance cold and Al tolerance in yeast, and AhLEA3-1 enhanced the drought tolerance in yeast.

CONCLUSIONS

AhLEAs are involved in abiotic stress response, and segmental duplication plays an important role in the evolution and amplification of AhLEAs. The genome-wide identification, classification, evolutionary and transcription analyses of the AhLEA gene family provide a foundation for further exploring the LEA genes' function in response to abiotic stress in peanuts.

Genomic Survey and Cold-Induced Expression Patterns of bHLH Transcription Factors in Liriodendron chinense (Hemsl) Sarg

bHLH transcription factors play an animated role in the plant kingdom during growth and development, and responses to various abiotic stress. In this current study, we conducted, the genome-wide survey of bHLH transcription factors in Liriodendron chinense (Hemsl) Sarg., 91 LcbHLH family members were identified. Identified LcbHLH gene family members were grouped into 19 different subfamilies based on the conserved motifs and phylogenetic analysis. Our results showed that LcbHLH genes clustered in the same subfamily exhibited a similar conservative exon-intron pattern. Hydrophilicity value analysis showed that all LcbHLH proteins were hydrophilic. The Molecular weight (Mw) of LcbHLH proteins ranged from 10.19 kD (LcbHLH15) to 88.40 kD (LcbHLH50). A greater proportion, ~63%, of LcbHLH proteins had a theoretical isoelectric point (pI) less than seven. Additional analysis on the collinear relationships within species and among dissimilar species illustrated that tandem and fragment duplication are the foremost factors of amplification of this family in the evolution process, and they are all purified and selected. RNA-seq and real-time quantitative PCR analysis of LcbHLH members showed that the expression of LcbHLH35, 55, and 86 are up-regulated, and the expression of LcbHLH9, 20, 39, 54, 56, and 69 is down-regulated during cold stress treatments while the expression of LcbHLH24 was up-regulated in the short term and then later down-regulated. From our results, we concluded that LcbHLH genes might participate in cold-responsive processes of L. chinense. These findings provide the basic information of bHLH gene in L. chinense and their regulatory roles in plant development and cold stress response.

Genome-wide identification of Tomato Golden 2-Like transcription factors and abiotic stress related members screening

BACKGROUND

Golden 2-Like (G2-like) transcription factors play an important role in plant development. However, the roles of these G2-like regulatory genes in response to abiotic stresses in tomato are not well understood.

RESULTS

In this study, we identified 66 putative G2-like genes in tomato (Solanum lycopersicum) and classified them into 5 groups (I to V) according to gene structure, motif composition and phylogenetic analysis. The G2-like genes were unevenly distributed across all 12 chromosomes. There were nine pairs of duplicated gene segments and four tandem duplicated SlGlk genes. Analysis of the cis-regulatory elements (CREs) showed that the promoter regions of SlGlks contain many kinds of stress- and hormone-related CREs. Based on RNA-seq, SlGlks were expressed in response to three abiotic stresses. Thirty-six differentially expressed SlGlks were identified; these genes have multiple functions according to Gene Ontology (GO) analysis and are enriched mainly in the zeatin biosynthesis pathway. Further studies exhibited that silencing SlGlk16 in tomato would reduce drought stress tolerance by earlier wilted, lower superoxide dismutase (SOD), peroxidase (POD) activities, less Pro contents and more MDA contents.

CONCLUSIONS

Overall, the results of this study provide comprehensive information on G2-like transcription factors and G2-like genes that may be expressed in response to abiotic stresses.

Genome-Wide Identification and Co-Expression Analysis of ARF and IAA Family Genes in Euscaphis konishii: Potential Regulators of Triterpenoids and Anthocyanin Biosynthesis

Euscaphis konishii is an evergreen plant that is widely planted as an industrial crop in Southern China. It produces red fruits with abundant secondary metabolites, giving E. konishii high medicinal and ornamental value. Auxin signaling mediated by members of the AUXIN RESPONSE FACTOR (ARF) and auxin/indole-3-acetic acid (Aux/IAA) protein families plays important roles during plant growth and development. Aux/IAA and ARF genes have been described in many plants but have not yet been described in E. konishii. In this study, we identified 34 EkIAA and 29 EkARF proteins encoded by the E. konishii genome through database searching using HMMER. We also performed a bioinformatic characterization of EkIAA and EkARF genes, including their phylogenetic relationships, gene structures, chromosomal distribution, and cis-element analysis, as well as conserved motifs in the proteins. Our results suggest that EkIAA and EkARF genes have been relatively conserved over evolutionary history. Furthermore, we conducted expression and co-expression analyses of EkIAA and EkARF genes in leaves, branches, and fruits, which identified a subset of seven EkARF genes as potential regulators of triterpenoids and anthocyanin biosynthesis. RT-qPCR, yeast one-hybrid, and transient expression analyses showed that EkARF5.1 can directly interact with auxin response elements and regulate downstream gene expression. Our results may pave the way to elucidating the function of EkIAA and EkARF gene families in E. konishii, laying a foundation for further research on high-yielding industrial products and E. konishii breeding.

Tandem Exon Duplications Expanding the Alternative Splicing Repertoire

Tandem exon duplications play an important role in the evolution of eukaryotic genes, providing a generic mechanism for adaptive regulation of protein function. In recent studies, tandem exon duplications have been linked to mutually exclusive exon choice, a pattern of alternative splicing in which one and only one exon from a group of tandemly arranged exons is included in the mature transcript. Here, we revisit the problem of identifying tandem exon duplications in eukaryotic genomes using bioinformatic methods and show that tandemly duplicated exons are abundant not only in the coding parts, but also in the untranslated regions. We present a number of remarkable examples of tandem exon duplications, identify unannotated duplicated exons, and provide statistical support for their expression using large panels of RNA-seq experiments.

Comparative genomic analysis of expansin superfamily gene members in zucchini and cucumber and their expression profiles under different abiotic stresses

Zucchini and cucumber belong to the Cucurbitaceae family, a group of economical and nutritious food plants that is consumed worldwide. Expansin superfamily proteins are generally localized in the cell wall of plants and are known to possess an effect on cell wall modification by causing the expansion of this region. Although the whole genome sequences of cucumber and zucchini plants have been resolved, the determination and characterization of expansin superfamily members in these plants using whole genomic data have not been implemented yet. In the current study, a genome-wide analysis of zucchini (Cucurbita pepo) and cucumber (Cucumis sativus) genomes was performed to determine the expansin superfamily genes. In total, 49 and 41 expansin genes were identified in zucchini and cucumber genomes, respectively. All expansin superfamily members were subjected to further bioinformatics analysis including gene and protein structure, ontology of the proteins, phylogenetic relations and conserved motifs, orthologous relations with other plants, targeting miRNAs of those genes and in silico gene expression profiles. In addition, various abiotic stress responses of zucchini and cucumber expansin genes were examined to determine their roles in stress tolerance. CsEXPB-04 and CsEXPA-11 from cucumber and CpEXPA-20 and CpEXPLA-14 from zucchini can be candidate genes for abiotic stress response and tolerance in addition to their roles in the normal developmental processes, which are supported by the gene expression analysis. This work can provide new perspectives for the roles of expansin superfamily genes and offers comprehensive knowledge for future studies investigating the modes of action of expansin proteins.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01108-w.

Genome-wide identification and expression analysis of the bHLH transcription factor family and its response to abiotic stress in foxtail millet (Setaria italica L.)

Background

Members of the basic helix-loop-helix (bHLH) transcription factor family perform indispensable functions in various biological processes, such as plant growth, seed maturation, and abiotic stress responses. However, the bHLH family in foxtail millet (Setaria italica), an important food and feed crop, has not been thoroughly studied.

Results

In this study, 187 bHLH genes of foxtail millet (SibHLHs) were identified and renamed according to the chromosomal distribution of the SibHLH genes. Based on the number of conserved domains and gene structure, the SibHLH genes were divided into 21 subfamilies and two orphan genes via phylogenetic tree analysis. According to the phylogenetic tree, the subfamilies 15 and 18 may have experienced stronger expansion in the process of evolution. Then, the motif compositions, gene structures, chromosomal spread, and gene duplication events were discussed in detail. A total of sixteen tandem repeat events and thirty-eight pairs of segment duplications were identified in bHLH family of foxtail millet. To further investigate the evolutionary relationship in the SibHLH family, we constructed the comparative syntenic maps of foxtail millet associated with representative monocotyledons and dicotyledons species. Finally, the gene expression response characteristics of 15 typical SibHLH genes in different tissues and fruit development stages, and eight different abiotic stresses were analysed. The results showed that there were significant differences in the transcription levels of some SibHLH members in different tissues and fruit development stages, and different abiotic stresses, implying that SibHLH members might have different physiological functions.

Conclusions

In this study, we identified 187 SibHLH genes in foxtail millet and further analysed the evolution and expression patterns of the encoded proteins. The findings provide a comprehensive understanding of the bHLH family in foxtail millet, which will inform further studies on the functional characteristics of SibHLH genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08095-y.

Genome-wide identification, classification, and expression analysis of the JmjC domain-containing histone demethylase gene family in birch

BACKGROUND

Histone methylation occurs primarily on lysine residues and requires a set of enzymes capable of reading, writing, and erasing to control its establishment and deletion, which is essential for maintaining chromatin structure and gene expression. Histone methylation and demethylation are contributed to plant growth and development, and are involved in adapting to environmental stresses. The JmjC domain-containing proteins are extensively studied for their function in histone lysine demethylation in plants, and play a critical role in sustaining histone methylation homeostasis.

RESULTS

In this study, a total of 21 JmjC domain-containing histone demethylase proteins (JHDMs) in birch were identified and classified into five subfamilies based on structural characteristics and phylogenetic relationships among Arabidopsis, rice, maize, and birch. Although the BpJMJ genes displayed significant schematic variation, their distribution on the chromosomes is relatively uniform. Additionally, the BpJMJ genes in birch have never experienced a tandem-duplication event proved by WGD analysis and were remaining underwent purifying selection (Ka/Ks < < 1). A typical JmjC domain was found in all BpJMJ genes, some of which have other essential domains for their functions. In the promoter regions of BpJMJ genes, cis-acting elements associated with hormone and abiotic stress responses were overrepresented. Under abiotic stresses, the transcriptome profile reveals two contrasting expression patterns within 21 BpJMJ genes. Furthermore, it was established that most BpJMJ genes had higher expression in young tissues under normal conditions, with BpJMJ06/16 having the highest expression in germinating seeds and participating in the regulation of BpGA3ox1/2 gene expression. Eventually, BpJMJ genes were found to directly interact with genes involved in the "intracellular membrane" in respond to cold stress.

CONCLUSIONS

The present study will provide a foundation for future experiments on histone demethylases in birch and a theoretical basis for epigenetic research on growth and development in response to abiotic stresses.

Genome-wide characterization and expression analysis of AP2/ERF genes in eggplant (Solanum melongena L.)

The AP2/ERF (APETALA2/Ethylene Response Factor) transcription factor superfamily plays crucial roles in a slew of physiological processes, such as plant growth and development, stress response, and secondary metabolites biosynthesis. Eggplant, especially the one rich with anthocyanins, is an economically important horticultural vegetable cultivated worldwide. In this study, we comprehensively analyzed the putative AP2/ERF gene family members and their response to abiotic stress in eggplant. As per the phylogenetic, conserved domains, and motif analysis, 178 AP2/ERF genes in this study belonged to five subfamilies. Chromosomal distributions analysis elucidated stochastic distribution of 178 putative SmAP2/ERF genes across the twelve chromosomes of eggplant. Expression profiles of sixteen selected AP2/ERF genes response to low temperature, drought, salt, abscisic acid, and ethylene treatments were analyzed, which revealed the involvement of SmAP2/ERF genes in diverse signaling pathways. In addition, we integrated RNA-Seq data on anthocyanin biosynthesis in eggplant with yeast one-hybrid and dual-luciferase assays and identified involvement of the SmAP2/ERF genes (Smechr0902114.1 and Smechr1102075.1) in the regulation of anthocyanin biosynthesis. This study will enable further functional characterization of AP2/ERF genes in eggplant and extend the current understanding of the role played by AP2/ERF genes in anthocyanin biosynthesis regulation.

Genome-wide identification and expression analysis of the bHLH transcription factor family and its response to abiotic stress in sorghum [Sorghum bicolor (L.) Moench]

BACKGROUND

Basic helix-loop-helix (bHLH) is a superfamily of transcription factors that is widely found in plants and animals, and is the second largest transcription factor family in eukaryotes after MYB. They have been shown to be important regulatory components in tissue development and many different biological processes. However, no systemic analysis of the bHLH transcription factor family has yet been reported in Sorghum bicolor.

RESULTS

We conducted the first genome-wide analysis of the bHLH transcription factor family of Sorghum bicolor and identified 174 SbbHLH genes. Phylogenetic analysis of SbbHLH proteins and 158 Arabidopsis thaliana bHLH proteins was performed to determine their homology. In addition, conserved motifs, gene structure, chromosomal spread, and gene duplication of SbbHLH genes were studied in depth. To further infer the phylogenetic mechanisms in the SbbHLH family, we constructed six comparative syntenic maps of S. bicolor associated with six representative species. Finally, we analyzed the gene-expression response and tissue-development characteristics of 12 typical SbbHLH genes in plants subjected to six different abiotic stresses. Gene expression during flower and fruit development was also examined.

CONCLUSIONS

This study is of great significance for functional identification and confirmation of the S. bicolor bHLH superfamily and for our understanding of the bHLH superfamily in higher plants.

Construction of a chromosome-level genome and variation map for the Pacific oyster Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is a widely distributed marine bivalve of great ecological and economic importance. In this study, we provide a high-quality chromosome-level genome assembled using Pacific Bioscience long reads and Hi-C-based and linkage-map-based scaffolding technologies and a high-resolution variation map constructed using large-scale resequencing analysis. The 586.8 Mb genome consists of 10 pseudochromosome sequences ranging from 38.6 to 78.9 Mb, containing 301 contigs with an N50 size of 3.1 Mb. A total of 30,078 protein-coding genes were predicted, of which 22,757 (75.7%) were high-reliability annotations supported by a homologous match to a curated protein in the SWISS-PROT database or transcript expression. Although a medium level of repeat components (57.2%) was detected, the genomic content of the segmental duplications reached 26.2%, which is the highest among the reported genomes. By whole genome resequencing analysis of 495 Pacific oysters, a comprehensive variation map was built, comprised of 4.78 million single nucleotide polymorphisms, 0.60 million short insertions and deletions, and 49,333 copy number variation regions. The structural variations can lead to an average interindividual genomic divergence of 0.21, indicating their crucial role in shaping the Pacific oyster genome diversity. The large amount of mosaic distributed repeat elements, small variations, and copy number variations indicate that the Pacific oyster is a diploid organism with an extremely high genomic complexity at the intra- and interindividual level. The genome and variation maps can improve our understanding of oyster genome diversity and enrich the resources for oyster molecular evolution, comparative genomics, and genetic research.

Genome-wide identification and expression pattern analysis of the ribonuclease T2 family in Eucommia ulmoides

The 2′,3′-cycling ribonuclease (RNase) genes are catalysts of RNA cleavage and include the RNase T2 gene family. RNase T2 genes perform important roles in plants and have been conserved in the genome of eukaryotic organisms. In this study we identified 21 EURNS genes in Eucommia ulmoides Oliver (E. ulmoides) and analyzed their structure, chromosomal location, phylogenetic tree, gene duplication, stress-related cis-elements, and expression patterns in different tissues. The length of 21 predicted EURNS proteins ranged from 143 to 374 amino acids (aa), their molecular weight (MW) ranged from 16.21 to 42.38 kDa, and their isoelectric point (PI) value ranged from 5.08 to 9.09. Two classifications (class I and class III) were obtained from the conserved domains analysis and phylogenetic tree. EURNS proteins contained a total of 15 motifs. Motif 1, motif 2, motif 3, and motif 7 were distributed in multiple sequences and were similar to the conserved domain of RNase T2. EURNS genes with similar structure and the predicted EURNS proteins with conserved motif compositions are in the same group in the phylogenetic tree. The results of RT-PCR and transcription data showed that EURNS genes have tissue-specific expression and exhibited obvious trends in different developmental stages. Gene duplication analysis results indicated that segment duplication may be the dominant duplication mode in this gene family. This study provides a theoretical basis for research on the RNase T2 gene family and lays a foundation for the further study of EURNS genes.

Genome-Wide Identification and Characterization of Lectin Receptor-Like Kinase Gene Family in Cucumber and Expression Profiling Analysis under Different Treatments

Lectin receptor-like kinases (LecRLKs) are a class of membrane proteins found in plants that are involved in diverse functions, including plant development and stress responses. Although LecRLK families have been identified in a variety of plants, a comprehensive analysis has not yet been undertaken in cucumber (Cucumis sativus L.). In this study, 46 putative LecRLK genes were identified in the cucumber genome, including 23 G-type and 22 L-type, and one C-type LecRLK gene. They were unequally distributed on all seven chromosomes, with a clustering tendency. Most of the genes in the cucumber LecRLK (CsLecRLK) gene family lacked introns. In addition, there were many regulatory elements associated with phytohormones and stress on these genes’ promoters. Transcriptome data demonstrated distinct expression patterns of CsLecRLK genes in various tissues. Furthermore, we found that each member of the CsLecRLK family had its own unique expression pattern under hormone and stress treatment by the quantitative real-time PCR (qRT-PCR) analysis. This study provides a better understanding of the character and function of the LecRLK gene family in cucumber and opens up the possibility to exploring the roles that LecRLKs might play in the life cycle of cucumber.

Comparative genome analysis of the SPL gene family reveals novel evolutionary features in maize

SPLs are plant-specific transcription factors that play important regulatory roles in plant growth and development. Systematic analysis of the SPL family has been performed in numerous plants, such as Arabidopsis, rice, and Populus. However, no comparative analysis has been performed across different species to examine evolutionary features. In this study, we present a comparative analysis of SPLs in different species. The results showed that 84 SPLs of the four species can be divided into six groups according to phylogeny. We found that most of the SPL-containing regions in maize showed extensive conservation with duplicated regions of rice and sorghum. A gene duplication analysis in maize indicated that ZmSPLs showed a significant excess of segmental duplication. The Ka/Ks analysis indicated that 9 out of 18 duplicated pairs in maize experienced positive selection, while SPL gene pairs of rice and sorghum mainly evolved under purifying selection, suggesting novel evolutionary features for ZmSPLs. The 31 ZmSPLs were further analyzed by describing their gene structure, phylogenetic relationships, chromosomal location, and expression, Among the ZmSPLs, 13 were predicated to be targeted by miR156s and involved in drought stress response. These results provide the foundation for future functional analyses of ZmSPLs.

Genome-wide analysis indicates diverse physiological roles of the turnip (Brassica rapa var. rapa) oligopeptide transporters gene family

Oligopeptide transporters (OPTs) encode integral membrane-localized proteins and have a broad range of substrate transport capabilities. Here, 28 BrrOPT genes were identified in the turnip. Phylogenetic analyses revealed two well-supported clades in the OPT family, containing 15 BrrOPTs and 13 BrrYSLs. The exon/intron structure of OPT clade are conserved but the yellow stripe-like (YSL) clade was different. The exon/intron of the YSL clade possesses structural differences, whereas the YSL class motifs structure are conserved. The OPT genes are distributed unevenly among the chromosomes of the turnip genome. Phylogenetic and chromosomal distribution analyses revealed that the expansion of the OPT gene family is mainly attributable to segmental duplication. For the expression profiles at different developmental stages, a comprehensive analysis provided insights into the possible functional divergence among members of the paralog OPT gene family. Different expression levels under a variety of ion deficiencies also indicated that the OPT family underwent functional divergence during long-term evolution. Furthermore, BrrOPT8.1, BrrYSL1.2, BrrYSL1.3, BrrYSL6 and BrrYSL9 responded to Fe(II) treatments and BrrYSL7 responded to calcium treatments, BrrYSL6 responded to multiple treatments in root, suggesting that turnip OPTs may be involved in mediating cross-talk among different ion deficiencies. Our data provide important information for further functional dissection of BrrOPTs, especially in transporting metal ions and nutrient deficiency stress adaptation.

Genome-Wide Identification and Characterization of the Potato bHLH Transcription Factor Family

Plant basic/helix–loop–helix (bHLH) transcription factors participate in a number of biological processes, such as growth, development and abiotic stress responses. The bHLH family has been identified in many plants, and several bHLH transcription factors have been functionally characterized in Arabidopsis. However, no systematic identification of bHLH family members has been reported in potato (Solanum tuberosum). Here, 124 StbHLH genes were identified and named according to their chromosomal locations. The intron numbers varied from zero to seven. Most StbHLH proteins had the highly conserved intron phase 0, which accounted for 86.2% of the introns. According to the Neighbor-joining phylogenetic tree, 259 bHLH proteins acquired from Arabidopsis and potato were divided into 15 groups. All of the StbHLH genes were randomly distributed on 12 chromosomes, and 20 tandem duplicated genes and four pairs of duplicated gene segments were detected in the StbHLH family. The gene ontology (GO) analysis revealed that StbHLH mainly function in protein and DNA binding. Through the RNA-seq and quantitative real time PCR (qRT-PCR) analyses, StbHLH were found to be expressed in various tissues and to respond to abiotic stresses, including salt, drought and heat. StbHLH1, 41 and 60 were highly expressed in flower tissues, and were predicted to be involved in flower development by GO annotation. StbHLH45 was highly expressed in salt, drought and heat stress, which suggested its important role in abiotic stress response. The results provide comprehensive information for further analyses of the molecular functions of the StbHLH gene family.

Genome-wide identification and comparative expression analysis of LEA genes in watermelon and melon genomes

Late embryogenesis abundant (LEA) proteins are large and diverse group of polypeptides which were first identified during seed dehydration and then in vegetative plant tissues during different stress responses. Now, gene family members of LEA proteins have been detected in various organisms. However, there is no report for this protein family in watermelon and melon until this study. A total of 73 LEA genes from watermelon (ClLEA) and 61 LEA genes from melon (CmLEA) were identified in this comprehensive study. They were classified into four and three distinct clusters in watermelon and melon, respectively. There was a correlation between gene structure and motif composition among each LEA groups. Segmental duplication played an important role for LEA gene expansion in watermelon. Maximum gene ontology of LEA genes was observed with poplar LEA genes. For evaluation of tissue specific expression patterns of ClLEA and CmLEA genes, publicly available RNA-seq data were analyzed. The expression analysis of selected LEA genes in root and leaf tissues of drought-stressed watermelon and melon were examined using qRT-PCR. Among them, ClLEA-12-17-46 genes were quickly induced after drought application. Therefore, they might be considered as early response genes for water limitation conditions in watermelon. In addition, CmLEA-42-43 genes were found to be up-regulated in both tissues of melon under drought stress. Our results can open up new frontiers about understanding of functions of these important family members under normal developmental stages and stress conditions by bioinformatics and transcriptomic approaches.

Genome-wide analysis of the fructose 1,6-bisphosphate aldolase (FBA) gene family and functional characterization of FBA7 in tomato

Fructose 1,6-bisphosphate aldolase (FBA) is a key enzyme in plants that is involved in glycolysis, gluconeogenesis, and the Calvin cycle. FBA genes play significant roles in biotic and abiotic stress responses and also regulate growth and development. Despite the importance of FBA genes, little is known about it in tomato. In this study, we identified 8 FBA genes in tomato and classified them into 2 subgroups based on a phylogenetic tree, gene structures, and conserved motifs. Five (SlFBA1, 2, 3, 4 and 5) and three (SlFBA6, 7, and 8) SlFBA proteins were predicted to be localized in chloroplasts and cytoplasm, respectively. The phylogenetic analysis of FBAs from tomato, Arabidopsis, rice, and other organisms suggested that SlFBA shared the highest protein homology with FBAs from other plants. Synteny analysis indicated that segmental duplication events contributed to the expansion of the tomato FBA family. The expression profiles revealed that all SlFBAs were involved in the response to low and high temperature stresses. SlFBA7 overexpression increased the expression and activities of other main enzymes in Calvin cycle, net photosynthetic rate (Pn), seed size and stem diameter. SlFBA7 overexpression enhanced tolerances in seed germination under suboptimal temperature stresses. Taken together, comprehensive analyses of SlFBAs would provide a basis for understanding of evolution and function of SlFBA family.

Genome-Wide Identification, Characterization and Expression Analysis of the Chalcone Synthase Family in Maize

Members of the chalcone synthase (CHS) family participate in the synthesis of a series of secondary metabolites in plants, fungi and bacteria. The metabolites play important roles in protecting land plants against various environmental stresses during the evolutionary process. Our research was conducted on comprehensive investigation of CHS genes in maize (Zea mays L.), including their phylogenetic relationships, gene structures, chromosomal locations and expression analysis. Fourteen CHS genes (ZmCHS01-14) were identified in the genome of maize, representing one of the largest numbers of CHS family members identified in one organism to date. The gene family was classified into four major classes (classes I-IV) based on their phylogenetic relationships. Most of them contained two exons and one intron. The 14 genes were unevenly located on six chromosomes. Two segmental duplication events were identified, which might contribute to the expansion of the maize CHS gene family to some extent. In addition, quantitative real-time PCR and microarray data analyses suggested that ZmCHS genes exhibited various expression patterns, indicating functional diversification of the ZmCHS genes. Our results will contribute to future studies of the complexity of the CHS gene family in maize and provide valuable information for the systematic analysis of the functions of the CHS gene family.

Genome-wide characterization and expression analysis of common bean bHLH transcription factors in response to excess salt concentration

Members of basic helix-loop-helix (bHLH) gene family found in all eukaryotes play crucial roles in response to stress. Though, most eukaryotes carry the proteins of this family, biological functions of the most bHLH family members are not deeply evaluated in plants. In this study, we conducted a comprehensive genome-wide analysis of bHLH transcription factors in salt tolerant common bean. We identified 155 bHLH protein-encoding genes (PvbHLH) by using in silico comparative genomics tools. Based on the phylogenetic tree, PvbHLH genes were classified into 8 main groups with 21 subfamilies. Exon-intron analysis indicated that proteins belonging to same main groups exhibited a closely related gene structure. While, the PvbHLH gene family has been mainly expanded through segmental duplications, a total of 11 tandem duplication were detected. Genome-wide expression analysis of bHLH genes showed that 63 PvbHLH genes were differentially expressed in at least one tissue. Three of them displayed higher expression values in both leaf and root tissues. The in silico micro-RNA target transcript analyses revealed that totally 100 PvHLH genes targeted by 86 plant miRNAs. The most abundant transcripts, which were targeted by all 18 plant miRNA, were belonging to PvHLH-22 and PvHLH-44 genes. The expression of 16 PvbHLH genes in the root and leaf tissues of salt-stressed common bean was evaluated using qRT-PCR. Among them, two of PvbHLHs, PvbHLH-54, PvbHLH-148, were found to be up-regulated in both tissues in correlation with RNA-seq measurements. The results of this study could help improve understanding of biological functions of common bean bHLH family under salt stress. Additionally, it may provide basic resources for analyzing bHLH protein function for improving economic, agronomic and ecological benefit in common bean and other species.

The Evolutionary History of R2R3-MYB Proteins Across 50 Eukaryotes: New Insights Into Subfamily Classification and Expansion

R2R3-MYB proteins (2R-MYBs) are one of the main transcription factor families in higher plants. Since the evolutionary history of this gene family across the eukaryotic kingdom remains unknown, we performed a comparative analysis of 2R-MYBs from 50 major eukaryotic lineages, with particular emphasis on land plants. A total of 1548 candidates were identified among diverse taxonomic groups, which allowed for an updated classification of 73 highly conserved subfamilies, including many newly identified subfamilies. Our results revealed that the protein architectures, intron patterns, and sequence characteristics were remarkably conserved in each subfamily. At least four subfamilies were derived from early land plants, 10 evolved from spermatophytes, and 19 from angiosperms, demonstrating the diversity and preferential expansion of this gene family in land plants. Moreover, we determined that their remarkable expansion was mainly attributed to whole genome and segmental duplication, where duplicates were preferentially retained within certain subfamilies that shared three homologous intron patterns (a, b, and c) even though up to 12 types of patterns existed. Through our integrated distributions, sequence characteristics, and phylogenetic tree analyses, we confirm that 2R-MYBs are old and postulate that 3R-MYBs may be evolutionarily derived from 2R-MYBs via intragenic domain duplication.

Identification and characterization of the RCI2 gene family in maize (Zea mays)

Rare-cold-inducible (RCI2) genes are structurally conserved members that encode small, highly hydrophobic proteins involved in response to various abiotic stresses. Phylogenetic and functional analyses of these genes have been conducted in Arabidopsis, but an extensive investigation of the RCI2 gene family has not yet been carried out in maize. In the present study, 10 RCI2 genes were identified in a fully sequenced maize genome. Structural characterization and expression pattern analysis of 10 ZmRCI2s (Zea mays RCI2 genes) were subsequently determined. Sequence and phylogenetic analyses indicated that ZmRCI2s are highly conserved, and most of them could be grouped with their orthologues from other organisms. Chromosomal location analysis indicated that ZmRCI2s were distributed unevenly on seven chromosomes with two segmental duplication events, suggesting that maize RCI2 gene family is an evolutionarily conserved family. Putative stress-responsive cis-elements were detected in the 2-kb promoter regions of the 10 ZmRCI2s. In addition, the 10 ZmRCI2s showed different expression patterns in maize development based on transcriptome analysis. Further, microarray and quantitative real-time PCR (qRT-PCR) analysis showed that each maize RCI2 genes were responsive to drought stress, suggesting their important roles in drought stress response. The results of this work provide a basis for future cloning and application studies of maize RCI2 genes.

A Populus TIR1 gene family survey reveals differential expression patterns and responses to 1-naphthaleneacetic acid and stress treatments

The plant hormone auxin is a central regulator of plant growth. TRANSPORT INHIBITOR RESPONSE 1/AUXIN SIGNALING F-BOX (TIR1/AFB) is a component of the E3 ubiquitin ligase complex SCF(TIR1/AFB) and acts as an auxin co-receptor for nuclear auxin signaling. The SCF(TIR1/AFB)-proteasome machinery plays a central regulatory role in development-related gene transcription. Populus trichocarpa, as a model tree, has a unique fast-growth trait to which auxin signaling may contribute. However, no systematic analyses of the genome organization, gene structure, and expression of TIR1-like genes have been undertaken in this woody model plant. In this study, we identified a total of eight TIR1 genes in the Populus genome that are phylogenetically clustered into four subgroups, PtrFBL1/PtrFBL2, PtrFBL3/PtrFBL4, PtrFBL5/PtrFBL6, and PtrFBL7/PtrFBL8, representing four paralogous pairs. In addition, the gene structure and motif composition were relatively conserved in each paralogous pair and all of the PtrFBL members were localized in the nucleus. Different sets of PtrFBLs were strongly expressed in the leaves, stems, roots, cambial zones, and immature xylem of Populus. Interestingly, PtrFBL1 and 7 were expressed mainly in vascular and cambial tissues, respectively, indicating their potential but different roles in wood formation. Furthermore, Populus FBLs responded differentially upon exposure to various stresses. Finally, over-expression studies indicated a role of FBL1 in poplar stem growth and response to drought stress. Collectively, these observations lay the foundation for further investigations into the potential roles of PtrFBL genes in tree growth and development.

Adolescent with tourette syndrome and bipolar disorder: a case report

Tourette syndrome consists of multiple motor tics and one or more vocal tics. Psychopathology occurs in approximately 90% of Tourette syndrome patients, with attention-deficit/hyperactivity, mood, and obsessive-compulsive disorders being common. Additionally, Tourette syndrome and bipolar disorder may be related in some individuals. However, it is unclear why bipolar disorder may be overrepresented in Tourette syndrome patients, and more research is needed. Herein, we report the case of a 15-year-old boy diagnosed with both Tourette syndrome and bipolar disorder, whose symptoms improved with aripiprazole, atomoxetine, and valproate. The patient was diagnosed with Tourette syndrome at 8 years of age when he developed tics and experienced his first depressive episode. The patient had a poor response to a variety of antidepressants and anti-tic medications. A combination of valproate and aripiprazole stabilized both the patient's tics and mood symptoms. It is important to assess individuals with Tourette syndrome for other disorders, including bipolar disorder. The treatment of children and adolescents with both Tourette syndrome and bipolar disorder is an important clinical issue.

Comprehensive analysis of the homeodomain-leucine zipper IV transcription factor family in Cucumis sativus

The class IV homeodomain-leucine zipper (HD-Zip IV) proteins are plant-specific transcriptional factors known to play crucial roles in plant growth and development. In this study, 11 cucumber (Cucumis sativus) HD-Zip IV genes were identified in the version 2 cucumber genome and found to be distributed unevenly across the chromosomes. The CsHDZIV (Cucumis sativus homeodomain-leucine zipper IV) gene family is smaller than in other studied species (except for rice) because of the absence of gene duplication events. Phylogenetic analysis showed that HD-Zip IV genes from cucumber, Arabidopsis, tomato, cotton, maize, and rice could be classified into five subgroups. All CsHDZIV genes appear to be derived from a basic module containing 11 exons in the coding region. Two conserved motifs of 21 and 19 nucleotides were found in the 3'-untranslated regions of six CsHDZIV genes, suggesting that post-transcriptional regulation may play a role in regulation of CsHDZIV genes. In addition, 6 of 11 CsHDZIV genes were found to undergo alternative splicing events. Reverse transcription PCR analysis showed that all CsHDZIV genes (except one) were expressed and showed preferential expression in reproductive organs.

Genome-wide analysis of the Populus Hsp90 gene family reveals differential expression patterns, localization, and heat stress responses

Background

Members of the heat shock protein 90 (Hsp90) class of proteins are evolutionarily conserved molecular chaperones. They are involved in protein folding, assembly, stabilization, activation, and degradation in many normal cellular processes and under stress conditions. Unlike many other well-characterized molecular chaperones, Hsp90s play key roles in signal transduction, cell-cycle control, genomic silencing, and protein trafficking. However, no systematic analysis of genome organization, gene structure, and expression compendium has been performed in the Populus model tree genus to date.

Results

We performed a comprehensive analysis of the Populus Hsp90 gene family and identified 10 Populus Hsp90 genes, which were phylogenetically clustered into two major groups. Gene structure and motif composition are relatively conserved in each group. In Populus trichocarpa, we identified three paralogous pairs, among which the PtHsp90-5a/PtHsp90-5b paralogous pair might be created by duplication of a genome segment. Subcellular localization analysis shows that PtHsp90 members are localized in different subcellular compartments. PtHsp90-3 is localized both in the nucleus and in the cytoplasm, PtHsp90-5a and PtHsp90-5b are in chloroplasts, and PtHsp90-7 is in the endoplasmic reticulum (ER). Furthermore, microarray and semi-quantitative real-time RT-PCR analyses show that a number of Populus Hsp90 genes are differentially expressed upon exposure to various stresses.

Conclusions

The gene structure and motif composition of PtHsp90s are highly conserved among group members, suggesting that members of the same group may also have conserved functions. Microarray and RT-PCR analyses show that most PtHsp90s were induced by various stresses, including heat stress. Collectively, these observations lay the foundation for future efforts to unravel the biological roles of PtHsp90 genes.

Genome-wide identification and characterization of R2R3MYB family in Cucumis sativus

BACKGROUND

The R2R3MYB proteins comprise one of the largest families of transcription factors in plants. Although genome-wide analysis of this family has been carried out in some species, little is known about R2R3MYB genes in cucumber (Cucumis sativus L.).

PRINCIPAL FINDINGS

This study has identified 55 R2R3MYB genes in the latest cucumber genome and the CsR2R3MYB family contained the smallest number of identified genes compared to other species that have been studied due to the absence of recent gene duplication events. These results were also supported by genome distribution and gene duplication analysis. Phylogenetic analysis showed that they could be classified into 11 subgroups. The evolutionary relationships and the intron-exon organizations that showed similarities with Arabidopsis, Vitis and Glycine R2R3MYB proteins were also analyzed and suggested strong gene conservation but also the expansions of particular functional genes during the evolution of the plant species. In addition, we found that 8 out of 55 (∼14.54%) cucumber R2R3MYB genes underwent alternative splicing events, producing a variety of transcripts from a single gene, which illustrated the extremely high complexity of transcriptome regulation. Tissue-specific expression profiles showed that 50 cucumber R2R3MYB genes were expressed in at least one of the tissues and the other 5 genes showed very low expression in all tissues tested, which suggested that cucumber R2R3MYB genes took part in many cellular processes. The transcript abundance level analysis during abiotic conditions (NaCl, ABA and low temperature treatments) identified a group of R2R3MYB genes that responded to one or more treatments.

CONCLUSIONS

This study has produced a comparative genomics analysis of the cucumber R2R3MYB gene family and has provided the first steps towards the selection of CsR2R3MYB genes for cloning and functional dissection that can be used in further studies to uncover their roles in cucumber growth and development.

Multiple pathways of duplication formation with and without recombination (RecA) in Salmonella enterica

Duplications are often attributed to "unequal recombination" between separated, directly repeated sequence elements (>100 bp), events that leave a recombinant element at the duplication junction. However, in the bacterial chromosome, duplications form at high rates (10(-3)-10(-5)/cell/division) even without recombination (RecA). Here we describe 1800 spontaneous lac duplications trapped nonselectively on the low-copy F'(128) plasmid, where lac is flanked by direct repeats of the transposable element IS3 (1258 bp) and by numerous quasipalindromic REP elements (30 bp). Duplications form at a high rate (10(-4)/cell/division) that is reduced only about 11-fold in the absence of RecA. With and without RecA, most duplications arise by recombination between IS3 elements (97%). Formation of these duplications is stimulated by IS3 transposase (Tnp) and plasmid transfer functions (TraI). Three duplication pathways are proposed. First, plasmid dimers form at a high rate stimulated by RecA and are then modified by deletions between IS3 elements (resolution) that leave a monomeric plasmid with an IS3-flanked lac duplication. Second, without RecA, duplications occur by single-strand annealing of DNA ends generated in different sister chromosomes after transposase nicks DNA near participating IS3 elements. The absence of RecA may stimulate annealing by allowing chromosome breaks to persist. Third, a minority of lac duplications (3%) have short (0-36 bp) junction sequences (SJ), some of which are located within REP elements. These duplication types form without RecA, Tnp, or Tra by a pathway in which the palindromic junctions of a tandem inversion duplication (TID) may stimulate deletions that leave the final duplication.

Systematic analysis of sequences and expression patterns of drought-responsive members of the HD-Zip gene family in maize

Background

Members of the homeodomain-leucine zipper (HD-Zip) gene family encode transcription factors that are unique to plants and have diverse functions in plant growth and development such as various stress responses, organ formation and vascular development. Although systematic characterization of this family has been carried out in Arabidopsis and rice, little is known about HD-Zip genes in maize (Zea mays L.).

Methods and Findings

In this study, we described the identification and structural characterization of HD-Zip genes in the maize genome. A complete set of 55 HD-Zip genes (Zmhdz1-55) were identified in the maize genome using Blast search tools and categorized into four classes (HD-Zip I-IV) based on phylogeny. Chromosomal location of these genes revealed that they are distributed unevenly across all 10 chromosomes. Segmental duplication contributed largely to the expansion of the maize HD-ZIP gene family, while tandem duplication was only responsible for the amplification of the HD-Zip II genes. Furthermore, most of the maize HD-Zip I genes were found to contain an overabundance of stress-related cis-elements in their promoter sequences. The expression levels of the 17 HD-Zip I genes under drought stress were also investigated by quantitative real-time PCR (qRT-PCR). All of the 17 maize HD-ZIP I genes were found to be regulated by drought stress, and the duplicated genes within a sister pair exhibited the similar expression patterns, suggesting their conserved functions during the process of evolution.

Conclusions

Our results reveal a comprehensive overview of the maize HD-Zip gene family and provide the first step towards the selection of Zmhdz genes for cloning and functional research to uncover their roles in maize growth and development.

Analyses of the oligopeptide transporter gene family in poplar and grape

Background

Oligopeptide transporters (OPTs) are a group of membrane-localized proteins that have a broad range of substrate transport capabilities and that are thought to contribute to many biological processes. The OPT proteins belong to a small gene family in plants, which includes about 25 members in Arabidopsis and rice. However, no comprehensive study incorporating phylogeny, chromosomal location, gene structure, expression profiling, functional divergence and selective pressure analysis has been reported thus far for Populus and Vitis.

Results

In the present study, a comprehensive analysis of the OPT gene family in Populus (P. trichocarpa) and Vitis (V. vinifera) was performed. A total of 20 and 18 full-length OPT genes have been identified in Populus and Vitis, respectively. Phylogenetic analyses indicate that these OPT genes consist of two classes that can be further subdivided into 11 groups. Gene structures are considerably conserved among the groups. The distribution of OPT genes was found to be non-random across chromosomes. A high proportion of the genes are preferentially clustered, indicating that tandem duplications may have contributed significantly to the expansion of the OPT gene family. Expression patterns based on our analyses of microarray data suggest that many OPT genes may be important in stress response and functional development of plants. Further analyses of functional divergence and adaptive evolution show that, while purifying selection may have been the main force driving the evolution of the OPTs, some of critical sites responsible for the functional divergence may have been under positive selection.

Conclusions

Overall, the data obtained from our investigation contribute to a better understanding of the complexity of the Populus and Vitis OPT gene family and of the function and evolution of the OPT gene family in higher plants.

Genome-wide identification and evolutionary analysis of Arabidopsis sm genes family

Sm proteins are members of a family of small proteins that are widespread in biosphere and found associated with RNA metabolism. To date, to our knowledge, only Arabidopsis SAD1 gene has been studied functionally in plant. In this study, 42 Sm genes are identified through comprehensive analysis in Arabidopsis. And a complete overview of this gene family is presented, including the gene structures, phylogeny, chromosome locations, selection pressure and expression. The results reveal that gene duplication contributes to the expansion of the Sm gene family in Arabidopsis genome, diverse expression patterns suggest their functional differentiation and divergence analysis indicates purifying selection as a key role in evolution. Our comparative genomics analysis of Sm genes will provide the first step towards the future experimental research on determining the functions of these genes.

Recurrent inversion events at 17q21.31 microdeletion locus are linked to the MAPT H2 haplotype

The chromosomal band 17q21.31, containing the microtubule-associated protein tau (MAPT) gene, is a hotspot for chromosomal rearrangements. It is known to contain a common inversion polymorphism of approximately 900 kb in populations with European ancestry. The inverted configuration is linked to a distinct MAPT haplotype, H2, which is relatively common in Europeans but nearly absent in Asian and African populations. Recent studies have demonstrated that the H2 haplotype is ancestral in hominoids, and under positive selection in Europeans. This haplotype is also linked to events leading to the 17q21.31 microdeletion syndrome, one of the most common causes of 'idiopathic' mental retardation in people of European descent. We performed direct analysis of the chromosome structure by fluorescence in situ hybridization and observed heterozygosity of the inversion status for the H2 chromosomes, but not for the H1 haplotype. Inversion heterozygosity was also observed in a mother homozygous for the H2 haplotype, who transmitted the chromosome with the deletion to a proband with 17q21.31 microdeletion syndrome. Our results highlight an allele-specific sensitivity to chromosome rearrangements and suggest that it is the heterozygosity of inversion status that predisposes to the 17q21.31 microdeletion syndrome.

Child and adolescent psychiatric genetics

The current status of child and adolescent psychiatric genetics appears promising in light of the initiation of genome-wide association studies (GWAS) for diverse polygenic disorders and the molecular elucidation of monogenic Rett syndrome, for which recent functional studies provide hope for pharmacological treatment strategies. Within the last 50 years, tremendous progress has been made in linking genetic variation to behavioral phenotypes and psychiatric disorders. We summarize the major findings of the Human Genome Project and dwell on largely unsuccessful candidate gene and linkage studies. GWAS for the first time offer the possibility to detect single nucleotide polymorphisms and copy number variants without a priori hypotheses as to their molecular etiology. At the same time it is becoming increasingly clear that very large sample sizes are required in order to enable genome wide significant findings, thus necessitating further large-scaled ascertainment schemes for the successful elucidation of the molecular genetics of childhood and adolescent psychiatric disorders. We conclude by reflecting on different scenarios for future research into the molecular basis of early onset psychiatric disorders. This review represents the introductory article of this special issue of the European Child and Adolescent Psychiatry.

Evolution in health and medicine Sackler colloquium: Genomic disorders: a window into human gene and genome evolution

Gene duplications alter the genetic constitution of organisms and can be a driving force of molecular evolution in humans and the great apes. In this context, the study of genomic disorders has uncovered the essential role played by the genomic architecture, especially low copy repeats (LCRs) or segmental duplications (SDs). In fact, regardless of the mechanism, LCRs can mediate or stimulate rearrangements, inciting genomic instability and generating dynamic and unstable regions prone to rapid molecular evolution. In humans, copy-number variation (CNV) has been implicated in common traits such as neuropathy, hypertension, color blindness, infertility, and behavioral traits including autism and schizophrenia, as well as disease susceptibility to HIV, lupus nephritis, and psoriasis among many other clinical phenotypes. The same mechanisms implicated in the origin of genomic disorders may also play a role in the emergence of segmental duplications and the evolution of new genes by means of genomic and gene duplication and triplication, exon shuffling, exon accretion, and fusion/fission events.