Evolutionary implications of intron-exon distribution and the properties and sequences of the RPL10A gene in eukaryotes

GhLCYε-3 characterized as a lycopene cyclase gene responding to drought stress in cotton

Drought stress significantly affects crop productivity. Carotenoids are essential photosynthetic pigment for plants, bacteria, and algae, with signaling and antioxidant functions. Lutein is a crucial branch product in the carotenoid synthesis pathway, which effectively improves the stress tolerance of higher plants. lycopene cyclase, a central enzyme for lutein synthesis, holds great significance in regulating lutein production. This research establishes a correlation between lutein content and stress resistance by measuring the drought resistance and lutein content of various cotton materials. To identify which crucial genes are associated with lutein, the lycopene cyclase family (LCYs) was analyzed. The research found that LCYs form a highly conserved family divided into two subfamilies, LCY-ε (lycopene ε-cyclase) and LCY-β (lycopene β-cyclase). Most members of the LCY family contain photoresponsive elements and abscisic acid elements. qRT-PCR demonstrates showed that most genes responded positively to drought stress, and GhLCYε-3 was expressed significantly differently under drought stress. Virus-induced gene silencing (VIGS) assay showed that the content of GhLCYε-3 was significantly increased with MDA and PRO, and the contents of chlorophyll and lutein were significantly decreased in pYL156 plants. The decrease in GhLCYε-3 expression is speculated to lead to reduced lutein content in vivo, resulting in the accumulation of reactive oxygen species (ROS) and decreased drought tolerance. This research enriched the understanding of LCY gene family and lutein function, and provided a new reference for cotton planting in arid areas.

Synopsis

Lycopene cyclase plays an important role in enhancing the ability of scavenging ROS and drought resistance of plants.

Genome-wide identification, phylogenetic and expression pattern analysis of HSF family genes in the Rye (Secale cereale L.)

BACKGROUND

Heat shock factor (HSF), a typical class of transcription factors in plants, has played an essential role in plant growth and developmental stages, signal transduction, and response to biotic and abiotic stresses. The HSF genes families has been identified and characterized in many species through leveraging whole genome sequencing (WGS). However, the identification and systematic analysis of HSF family genes in Rye is limited.

RESULTS

In this study, 31 HSF genes were identified in Rye, which were unevenly distributed on seven chromosomes. Based on the homology of A. thaliana, we analyzed the number of conserved domains and gene structures of ScHSF genes that were classified into seven subfamilies. To better understand the developmental mechanisms of ScHSF family during evolution, we selected one monocotyledon (Arabidopsis thaliana) and five (Triticum aestivum L., Hordeum vulgare L., Oryza sativa L., Zea mays L., and Aegilops tauschii Coss.) specific representative dicotyledons associated with Rye for comparative homology mapping. The results showed that fragment replication events modulated the expansion of the ScHSF genes family. In addition, interactions between ScHSF proteins and promoters containing hormone- and stress-responsive cis-acting elements suggest that the regulation of ScHSF expression was complex. A total of 15 representative genes were targeted from seven subfamilies to characterize their gene expression responses in different tissues, fruit developmental stages, three hormones, and six different abiotic stresses.

CONCLUSIONS

This study demonstrated that ScHSF genes, especially ScHSF1 and ScHSF3, played a key role in Rye development and its response to various hormones and abiotic stresses. These results provided new insights into the evolution of HSF genes in Rye, which could help the success of molecular breeding in Rye.

Identification of genomic biomarkers and their pathway crosstalks for deciphering mechanistic links in glioblastoma

Glioblastoma is a grade IV pernicious neoplasm occurring in the supratentorial region of brain. As its causes are largely unknown, it is essential to understand its dynamics at the molecular level. This necessitates the identification of better diagnostic and prognostic molecular candidates. Blood-based liquid biopsies are emerging as a novel tool for cancer biomarker discovery, guiding the treatment and improving its early detection based on their tumour origin. There exist previous studies focusing on the identification of tumour-based biomarkers for glioblastoma. However, these biomarkers inadequately represent the underlying pathological state and incompletely illustrate the tumour because of non-recursive nature of this approach to monitor the disease. Also, contrary to the tumour biopsies, liquid biopsies are non-invasive and can be performed at any interval during the disease span to surveil the disease. Therefore, in this study, a unique dataset of blood-based liquid biopsies obtained primarily from tumour-educated blood platelets (TEP) is utilised. This RNA-seq data from ArrayExpress is acquired comprising human cohort with 39 glioblastoma subjects and 43 healthy subjects. Canonical and machine learning approaches are applied for identification of the genomic biomarkers for glioblastoma and their crosstalks. In our study, 97 genes appeared enriched in 7 oncogenic pathways (RAF-MAPK, P53, PRC2-EZH2, YAP conserved, MEK-MAPK, ErbB2 and STK33 signalling pathways) using GSEA, out of which 17 have been identified participating actively in crosstalks. Using PCA, 42 genes are found enriched in 7 pathways (cytoplasmic ribosomal proteins, translation factors, electron transport chain, ribosome, Huntington's disease, primary immunodeficiency pathways, and interferon type I signalling pathway) harbouring tumour when altered, out of which 25 actively participate in crosstalks. All the 14 pathways foster well-known cancer hallmarks and the identified DEGs can serve as genomic biomarkers, not only for the diagnosis and prognosis of Glioblastoma but also in providing a molecular foothold for oncogenic decision making in order to fathom the disease dynamics. Moreover, SNP analysis for the identified DEGs is performed to investigate their roles in disease dynamics in an elaborated manner. These results suggest that TEPs are capable of providing disease insights just like tumour cells with an advantage of being extracted anytime during the course of disease in order to monitor it.

Genome-wide identification, molecular evolution and expression analysis of the non-specific lipid transfer protein (nsLTP) family in Setaria italica

BACKGROUND

Foxtail millet (Setaria italica L.) is a millet species with high tolerance to stressful environments. Plant non-specific lipid transfer proteins (nsLTPs) are a kind of small, basic proteins involved in many biological processes. So far, the genome of S. italica has been fully sequenced, and a comprehensive understanding of the evolution and expression of the nsLTP family is still lacking in foxtail millet.

RESULTS

Forty-five nsLTP genes were identified in S. italica and clustered into 5 subfamilies except three single genes (SinsLTP38, SinsLTP7, and SinsLTP44). The proportion of SinsLTPs was different in each subfamily, and members within the same subgroup shared conserved exon-intron structures. Besides, 5 SinsLTP duplication events were investigated. Both tandem and segmental duplication contributed to nsLTP expansion in S. italica, and the duplicated SinsLTPs had mainly undergone purifying selection pressure, which suggested that the function of the duplicated SinsLTPs might not diverge much. Moreover, we identified the nsLTP members in 5 other monocots, and 41, 13, 10, 4, and 1 orthologous gene pairs were identified between S. italica and S. viridis, S. bicolor, Z. mays, O. sativa, and B. distachyon, respectively. The functional divergence within the nsLTP orthologous genes might be limited. In addition, the tissue-specific expression patterns of the SinsLTPs were investigated, and the expression profiles of the SinsLTPs in response to abiotic stress were analyzed, all the 10 selected SinsLTPs were responsive to drought, salt, and cold stress. Among the selected SinsLTPs, 2 paired duplicated genes shared almost equivalent expression profiles, suggesting that these duplicated genes might retain some essential functions during subsequent evolution.

CONCLUSIONS

The present study provided the first systematic analysis for the phylogenetic classification, conserved domain and gene structure, expansion pattern, and expression profile of the nsLTP family in S. italica. These findings could pave a way for further comparative genomic and evolution analysis of nsLTP family in foxtail millet and related monocots, and lay the foundation for the functional analysis of the nsLTPs in S. italica.

Genome-wide identification and expression analysis of the SPL transcription factor family and its response to abiotic stress in Quinoa (Chenopodium quinoa)

BACKGROUND

Squamous promoter binding protein-like (SPL) proteins are a class of transcription factors that play essential roles in plant growth and development, signal transduction, and responses to biotic and abiotic stresses. The rapid development of whole genome sequencing has enabled the identification and characterization of SPL gene families in many plant species, but to date this has not been performed in quinoa (Chenopodium quinoa).

RESULTS

This study identified 23 SPL genes in quinoa, which were unevenly distributed on 18 quinoa chromosomes. Quinoa SPL genes were then classified into eight subfamilies based on homology to Arabidopsis thaliana SPL genes. We selected three dicotyledonous and monocotyledonous representative species, each associated with C. quinoa, for comparative sympatric mapping to better understand the evolution of the developmental mechanisms of the CqSPL family. Furthermore, we also used 15 representative genes from eight subfamilies to characterize CqSPLs gene expression in different tissues and at different fruit developmental stages under six different abiotic stress conditions.

CONCLUSIONS

This study, the first to identify and characterize SPL genes in quinoa, reported that CqSPL genes, especially CqSPL1, play a critical role in quinoa development and in its response to various abiotic stresses.

Identification and Expression Analysis of bZIP Members under Abiotic Stress in Mung Bean (Vigna radiata)

The main aim of this study was to identify the bZIP family members in mung bean and explore their expression patterns under several abiotic stresses, with the overarching goal of elucidating their biological functions. Results identified 75 bZIP members in mung bean, which were unevenly distributed in the chromosomes (1-11), and all had a highly conserved bZIP domain. Phylogenetic analysis divided the members into 10 subgroups, with members in the same subgroup having similar structure and motif. The cis-acting elements in the promoter region revealed that most of the bZIP members might have the connection with abscisic acid, ethylene, and stress responsive elements. The transcriptome data demonstrated that bZIP members could respond to salt stress at different degrees in leaves, but the expression patterns could vary at different time points under stress. Differentially expressed genes (DEGs), such as VrbZIP12, VrbZIP37, and VrZIP45, were annotated into the plant hormone signal transduction pathway, which might be regulated the expression of abiotic stress-related gene (ABF). Quantitative real-time polymerase chain reaction (qRT-PCR) was applied to determine the expression of bZIP members in roots and leaves under drought, alkali, and low-temperature stress. Results showed that bZIP members respond differently to diverse stresses, and their expression was tissue-specific, which suggests that they may have different regulatory mechanism in different tissues. Overall, this study will provide a reference for further research on the functions of bZIP members in mung bean.

Identification, evolutionary analysis and functional diversification of RAV gene family in cotton (G. hirsutum L.)

Genome wide analysis, expression pattern analysis, and functional characterization of RAV genes highlight their roles in roots, stem development and hormonal response. RAV (Related to ABI3 and VP1) gene family members have been involved in tissues/organs growth and hormone signaling in various plant species. Here, we identified 247 RAVs from 12 different species with 33 RAV genes from G. hirsutum. Phylogenetic analysis classified RAV genes into four distinct groups. Analysis of gene structure showed that most GhRAVs lack introns. Motif distribution pattern and protein sequence logos indicated that GhRAV genes were highly conserved during the process of evolution. Promotor cis-acting elements revealed that promotor regions of GhRAV genes encode numerous elements related to plant growth, abiotic stresses and phytohormones. Chromosomal location information showed uneven distribution of 33 GhRAV genes on different chromosomes. Collinearity analysis identified 628 and 52 orthologous/ paralogous gene pairs in G. hirsutum and G. barbadense, respectively. Ka/Ks values indicated that GhRAV and GbRAV genes underwent strong purifying selection pressure. Selecton model and codon model selection revealed that GhRAV amino acids were under purifying selection and adaptive evolution exists among GhRAV proteins. Three dimensional structure of GhRAVs indicated the presence of numerous alpha helix and beta-barrels. Expression level revealed that some GhRAV genes exhibited high expression in roots (GhRAV3, GhRAV4, GhRAV11, GhRAV18, GhRAV20 and GhRAV30) and stem (GhRAV3 and GhRAV18), indicating their potential role in roots and stem development. GhRAV genes can be regulated by phytohormonal stresses (BL, JA and IAA). Our study provides a reference for future studies related to the functional analysis of GhRAVs in cotton.

Identification of Distant Regulatory Elements Using Expression Quantitative Trait Loci Mapping for Heat-Responsive Genes in Oysters

Many marine ectotherms, especially those inhabiting highly variable intertidal zones, develop high phenotypic plasticity in response to rapid climate change by modulating gene expression levels. Herein, we examined the regulatory architecture of heat-responsive gene expression plasticity in oysters using expression quantitative trait loci (eQTL) analysis. Using a backcross family of Crassostrea gigas and its sister species Crassostrea angulata under acute stress, 56 distant regulatory regions accounting for 6–26.6% of the gene expression variation were identified for 19 heat-responsive genes. In total, 831 genes and 164 single nucleotide polymorphisms (SNPs) that could potentially regulate expression of the target genes were screened in the eQTL region. The association between three SNPs and the corresponding target genes was verified in an independent family. Specifically, Marker13973 was identified for heat shock protein (HSP) family A member 9 (HspA9). Ribosomal protein L10a (RPL10A) was detected approximately 2 kb downstream of the distant regulatory SNP. Further, Marker14346-48 and Marker14346-85 were in complete linkage disequilibrium and identified for autophagy-related gene 7 (ATG7). Nuclear respiratory factor 1 (NRF1) was detected approximately 3 kb upstream of the two SNPs. These results suggested regulatory relationships between RPL10A and HSPA9 and between NRF1 and ATG7. Our findings indicate that distant regulatory mutations play an important role in the regulation of gene expression plasticity by altering upstream regulatory factors in response to heat stress. The identified eQTLs provide candidate biomarkers for predicting the persistence of oysters under future climate change scenarios.

Genome-Wide Identification and Expression Analysis of the β-Amylase Gene Family in Chenopodium quinoa

β-Amylase (BAM) is an important starch hydrolase, playing a role in a variety of plant growth and development processes. In this study, 22 BAM gene family members (GFMs) were identified in quinoa (Chenopodium quinoa), an ancient crop gaining modern consumer acceptance because of its nutritional qualities. The genetic structure, phylogenetic and evolutionary relationships, and expression patterns of CqBAM GFMs in different tissues, were analyzed. Phylogenetic analyses assigned the CqBAMs, AtBAMs, and OsBAMs into four clades. The CqBAM gene family had expanded due to segmental duplication. RNA-seq analysis revealed expression of the duplicated pairs to be similar, with the expression of CqBAM GFM pairs showing a degree of tissue specificity that was confirmed by reverse transcription quantitative PCR (RT-qPCR). Several CqBAM GFMs were also responsive to abiotic stresses in shoots and/or roots. In conclusion, the BAM gene family in quinoa was identified and systematically analyzed using bioinformatics and experimental methods. These results will help to elucidate the evolutionary relationship and biological functions of the BAM gene family in quinoa.

Ribosomal Protein L10: From Function to Dysfunction

Eukaryotic cytoplasmic ribosomes are highly structured macromolecular complexes made up of four different ribosomal RNAs (rRNAs) and 80 ribosomal proteins (RPs), which play a central role in the decoding of genetic code for the synthesis of new proteins. Over the past 25 years, studies on yeast and human models have made it possible to identify RPL10 (ribosomal protein L10 gene), which is a constituent of the large subunit of the ribosome, as an important player in the final stages of ribosome biogenesis and in ribosome function. Here, we reviewed the literature to give an overview of the role of RPL10 in physiologic and pathologic processes, including inherited disease and cancer.

Genome-wide survey, characterization, and expression analysis of bZIP transcription factors in Chenopodium quinoa

BACKGROUND

Chenopodium quinoa Willd. (quinoa) is a pseudocereal crop of the Amaranthaceae family and represents a promising species with the nutritional content and high tolerance to stressful environments, such as soils affected by high salinity. The basic leucine zipper (bZIP) transcription factor represents exclusively in eukaryotes and can be related to many biological processes. So far, the genomes of quinoa and 3 other Amaranthaceae crops (Spinacia oleracea, Beta vulgaris, and Amaranthus hypochondriacus) have been fully sequenced. However, information about the bZIPs in these Amaranthaceae species is limited, and genome-wide analysis of the bZIP family is lacking in quinoa.

RESULTS

We identified 94 bZIPs in quinoa (named as CqbZIP1-CqbZIP94). All the CqbZIPs were phylogenetically splitted into 12 distinct subfamilies. The proportion of CqbZIPs was different in each subfamily, and members within the same subgroup shared conserved exon-intron structures and protein motifs. Besides, 32 duplicated CqbZIP gene pairs were investigated, and the duplicated CqbZIPs had mainly undergone purifying selection pressure, which suggested that the functions of the duplicated CqbZIPs might not diverge much. Moreover, we identified the bZIP members in 3 other Amaranthaceae species, and 41, 32, and 16 orthologous gene pairs were identified between quinoa and S. oleracea, B. vulgaris, and A. hypochondriacus, respectively. Among them, most were a single copy being present in S. oleracea, B. vulgaris, and A. hypochondriacus, and two copies being present in allotetraploid quinoa. The function divergence within the bZIP orthologous genes might be limited. Additionally, 11 selected CqbZIPs had specific spatial expression patterns, and 6 of 11 CqbZIPs were up-regulated in response to salt stress. Among the selected CqbZIPs, 3 of 4 duplicated gene pairs shared similar expression patterns, suggesting that these duplicated genes might retain some essential functions during subsequent evolution.

CONCLUSIONS

The present study provided the first systematic analysis for the phylogenetic classification, motif and gene structure, expansion pattern, and expression profile of the bZIP family in quinoa. Our results would lay an important foundation for functional and evolutionary analysis of CqbZIPs, and provide promising candidate genes for further investigation in tissue specificity and their functional involvement in quinoa's resistance to salt stress.

Genome-Wide Characterization and Analysis of CIPK Gene Family in Two Cultivated Allopolyploid Cotton Species: Sequence Variation, Association with Seed Oil Content, and the Role of GhCIPK6

Calcineurin B-like protein-interacting protein kinases (CIPKs), as key regulators, play an important role in plant growth and development and the response to various stresses. In the present study, we identified 80 and 78 CIPK genes in the Gossypium hirsutum and G. barbadense, respectively. The phylogenetic and gene structure analysis divided the cotton CIPK genes into five groups which were classified into an exon-rich clade and an exon-poor clade. A synteny analysis showed that segmental duplication contributed to the expansion of Gossypium CIPK gene family, and purifying selection played a major role in the evolution of the gene family in cotton. Analyses of expression profiles showed that GhCIPK genes had temporal and spatial specificity and could be induced by various abiotic stresses. Fourteen GhCIPK genes were found to contain 17 non-synonymous single nucleotide polymorphisms (SNPs) and co-localized with oil or protein content quantitative trait loci (QTLs). Additionally, five SNPs from four GhCIPKs were found to be significantly associated with oil content in one of the three field tests. Although most GhCIPK genes were not associated with natural variations in cotton oil content, the overexpression of the GhCIPK6 gene reduced the oil content and increased C18:1 and C18:1+C18:1d6 in transgenic cotton as compared to wild-type plants. In addition, we predicted the potential molecular regulatory mechanisms of the GhCIPK genes. In brief, these results enhance our understanding of the roles of CIPK genes in oil synthesis and stress responses.

Genome-Wide Identification and Characterization of the Vacuolar H+-ATPase Subunit H Gene Family in Crop Plants

The vacuolar H⁺-ATPase (V-ATPase) plays many important roles in cell growth and in response to stresses in plants. The V-ATPase subunit H (VHA-H) is required to form a stable and active V-ATPase. Genome-wide analyses of VHA-H genes in crops contribute significantly to a systematic understanding of their functions. A total of 22 VHA-H genes were identified from 11 plants representing major crops including cotton, rice, millet, sorghum, rapeseed, maize, wheat, soybean, barley, potato, and beet. All of these VHA-H genes shared exon-intron structures similar to those of Arabidopsis thaliana. The C-terminal domain of VHA-H was shorter and more conserved than the N-terminal domain. The VHA-H gene was effectively used as a genetic marker to infer the phylogenetic relationships among plants, which were congruent with currently accepted taxonomic groupings. The VHA-H genes from six species of crops (Gossypium raimondii, Brassica napus, Glycine max, Solanum tuberosum, Triticum aestivum, and Zea mays) showed high gene structural diversity. This resulted from the gains and losses of introns. Seven VHA-H genes in six species of crops (Gossypium raimondii, Hordeum vulgare, Solanum tuberosum, Setaria italica, Triticum aestivum, and Zea mays) contained multiple transcript isoforms arising from alternative splicing. The study of cis-acting elements of gene promoters and RNA-seq gene expression patterns confirms the role of VHA-H genes as eco-enzymes. The gene structural diversity and proteomic diversity of VHA-H genes in our crop sampling facilitate understanding of their functional diversity, including stress responses and traits important for crop improvement.

Genome-Wide Identification, Characterization, and Expression Analysis of the NAC Transcription Factor in Chenopodium quinoa

The NAC (NAM, ATAF, and CUC) family is one of the largest families of plant-specific transcription factors. It is involved in many plant growth and development processes, as well as abiotic/biotic stress responses. So far, little is known about the NAC family in Chenopodium quinoa. In the present study, a total of 90 NACs were identified in quinoa (named as CqNAC1-CqNAC90) and phylogenetically divided into 14 distinct subfamilies. Different subfamilies showed diversities in gene proportions, exon-intron structures, and motif compositions. In addition, 28 CqNAC duplication events were investigated, and a strong subfamily preference was found during the NAC expansion in quinoa, indicating that the duplication event was not random across NAC subfamilies during quinoa evolution. Moreover, the analysis of Ka/Ks (non-synonymous substitution rate/synonymous substitution rate) ratios suggested that the duplicated CqNACs might have mainly experienced purifying selection pressure with limited functional divergence. Additionally, 11 selected CqNACs showed significant tissue-specific expression patterns, and all the CqNACs were positively regulated in response to salt stress. The result provided evidence for selecting candidate genes for further characterization in tissue/organ specificity and their functional involvement in quinoa's strong salinity tolerance.

Dynamic evolution of mitochondrial genomes in Trebouxiophyceae, including the first completely assembled mtDNA from a lichen-symbiont microalga (Trebouxia sp. TR9)

Trebouxiophyceae (Chlorophyta) is a species-rich class of green algae with a remarkable morphological and ecological diversity. Currently, there are a few completely sequenced mitochondrial genomes (mtDNA) from diverse Trebouxiophyceae but none from lichen symbionts. Here, we report the mitochondrial genome sequence of Trebouxia sp. TR9 as the first complete mtDNA sequence available for a lichen-symbiont microalga. A comparative study of the mitochondrial genome of Trebouxia sp. TR9 with other chlorophytes showed important organizational changes, even between closely related taxa. The most remarkable change is the enlargement of the genome in certain Trebouxiophyceae, which is principally due to larger intergenic spacers and seems to be related to a high number of large tandem repeats. Another noticeable change is the presence of a relatively large number of group II introns interrupting a variety of tRNA genes in a single group of Trebouxiophyceae, which includes Trebouxiales and Prasiolales. In addition, a fairly well-resolved phylogeny of Trebouxiophyceae, along with other Chlorophyta lineages, was obtained based on a set of seven well-conserved mitochondrial genes.

Molecular Evolution and Stress and Phytohormone Responsiveness of SUT Genes in Gossypium hirsutum

Sucrose transporters (SUTs) play key roles in allocating the translocation of assimilates from source to sink tissues. Although the characteristics and biological roles of SUTs have been intensively investigated in higher plants, this gene family has not been functionally characterized in cotton. In this study, we performed a comprehensive analysis of SUT genes in the tetraploid cotton Gossypium hirsutum. A total of 18 G. hirsutum SUT genes were identified and classified into three groups based on their evolutionary relationships. Up to eight SUT genes in G. hirsutum were placed in the dicot-specific SUT1 group, while four and six SUT genes were, respectively, clustered into SUT4 and SUT2 groups together with members from both dicot and monocot species. The G. hirsutum SUT genes within the same group displayed similar exon/intron characteristics, and homologous genes in G. hirsutum At and Dt subgenomes, G. arboreum, and G. raimondii exhibited one-to-one relationships. Additionally, the duplicated genes in the diploid and polyploid cotton species have evolved through purifying selection, suggesting the strong conservation of SUT loci in these species. Expression analysis in different tissues indicated that SUT genes might play significant roles in cotton fiber elongation. Moreover, analyses of cis-acting regulatory elements in promoter regions and expression profiling under different abiotic stress and exogenous phytohormone treatments implied that SUT genes, especially GhSUT6A/D, might participate in plant responses to diverse abiotic stresses and phytohormones. Our findings provide valuable information for future studies on the evolution and function of SUT genes in cotton.

Genomic Identification and Comparative Expansion Analysis of the Non-Specific Lipid Transfer Protein Gene Family in Gossypium

Plant non-specific lipid transfer proteins (nsLTPs) are involved in many biological processes. In this study, 51, 47 and 91 nsLTPs were identified in Gossypium arboreum, G. raimondii and their descendant allotetraploid G. hirsutum, respectively. All the nsLTPs were phylogenetically divided into 8 distinct subfamilies. Besides, the recent duplication, which is considered cotton-specific whole genome duplication, may have led to nsLTP expansion in Gossypium. Both tandem and segmental duplication contributed to nsLTP expansion in G. arboreum and G. hirsutum, while tandem duplication was the dominant pattern in G. raimondii. Additionally, the interspecific orthologous gene pairs in Gossypium were identified. Some GaLTPs and GrLTPs lost their orthologs in the At and Dt subgenomes, respectively, of G. hirsutum. The distribution of these GrLTPs and GaLTPs within each subfamily was complementary, suggesting that the loss and retention of nsLTPs in G. hirsutum might not be random. Moreover, the nsLTPs in the At and Dt subgenomes might have evolved symmetrically. Furthermore, both intraspecific and interspecific orthologous genes showed considerable expression variation, suggesting that their functions were strongly differentiated. Our results lay an important foundation for expansion and evolutionary analysis of the nsLTP family in Gossypium, and advance nsLTP studies in other plants, especially polyploid plants.

Whole genome/proteome based phylogeny reconstruction for prokaryotes using higher order Markov model and chaos game representation

Traditional methods for sequence comparison and phylogeny reconstruction rely on pair wise and multiple sequence alignments. But alignment could not be directly applied to whole genome/proteome comparison and phylogenomic studies due to their high computational complexity. Hence alignment-free methods became popular in recent years. Here we propose a fast alignment-free method for whole genome/proteome comparison and phylogeny reconstruction using higher order Markov model and chaos game representation. In the present method, we use the transition matrices of higher order Markov models to characterize amino acid or DNA sequences for their comparison. The order of the Markov model is uniquely identified by maximizing the average Shannon entropy of conditional probability distributions. Using one-dimensional chaos game representation and linked list, this method can reduce large memory and time consumption which is due to the large-scale conditional probability distributions. To illustrate the effectiveness of our method, we employ it for fast phylogeny reconstruction based on genome/proteome sequences of two species data sets used in previous published papers. Our results demonstrate that the present method is useful and efficient.

AVAILABILITY AND IMPLEMENTATION

The source codes for our algorithm to get the distance matrix and genome/proteome sequences can be downloaded from ftp://121.199.20.25/. The software Phylip and EvolView we used to construct phylogenetic trees can be referred from their websites.