Reductive genome evolution from the mother of Rickettsia

The Rickettsia genus is a group of obligate intracellular α-proteobacteria representing a paradigm of reductive evolution. Here, we investigate the evolutionary processes that shaped the genomes of the genus. The reconstruction of ancestral genomes indicates that their last common ancestor contained more genes, but already possessed most traits associated with cellular parasitism. The differences in gene repertoires across modern Rickettsia are mainly the result of differential gene losses from the ancestor. We demonstrate using computer simulation that the propensity of loss was variable across genes during this process. We also analyzed the ratio of nonsynonymous to synonymous changes (Ka/Ks) calculated as an average over large sets of genes to assay the strength of selection acting on the genomes of Rickettsia, Anaplasmataceae, and free-living γ-proteobacteria. As a general trend, Ka/Ks were found to decrease with increasing divergence between genomes. The high Ka/Ks for closely related genomes are probably due to a lag in the removal of slightly deleterious nonsynonymous mutations by natural selection. Interestingly, we also observed a decrease of the rate of gene loss with increasing divergence, suggesting a similar lag in the removal of slightly deleterious pseudogene alleles. For larger divergence (Ks > 0.2), Ka/Ks converge toward similar values indicating that the levels of selection are roughly equivalent between intracellular α-proteobacteria and their free-living relatives. This contrasts with the view that obligate endocellular microorganisms tend to evolve faster as a consequence of reduced effectiveness of selection, and suggests a major role of enhanced background mutation rates on the fast protein divergence in the obligate intracellular α-proteobacteria.     

Large-scale profiling and identification of potential regulatory mechanisms for allelic gene expression in colorectal cancer cells

Allelic variation in gene expression is common in humans and this variation is associated with phenotypic variation. In this study, we employed high-density single nucleotide polymorphism (SNP) chips containing 13,900 exonic SNPs to identify genes with allelic gene expression in cells from colorectal cancer cell lines. We found 2 monoallelically expressed genes (ERAP2 and MYLK4), 32 genes with an allelic imbalance in their expression, and 13 genes showing allele substitution by RNA editing. Among a total of 34 allelically expressed genes in colorectal cancer cells, 15 genes (44.1%) were associated with cis-acting eQTL, indicating that large portions of allelically expressed genes are regulated by cis-acting mechanisms of gene expression. In addition, potential regulatory variants present in the proximal promoter regions of genes showing either monoallelic expression or allelic imbalance were not tightly linked with coding SNPs, which were detected with allelic gene expression. These results suggest that multiple rare variants could be involved in the cis-acting regulatory mechanism of allelic gene expression. In the comparison with allelic gene expression data from Centre d'Etude du Polymorphisme Humain (CEPH) family B cells, 12 genes showed B-cell specific allelic imbalance and 1 noncoding SNP showed colorectal cancer cell-specific allelic imbalance. In addition, different patterns of allele substitution were observed between B cells and colorectal cancer cells. Overall, our study not only indicates that allelic gene expression is common in colorectal cancer cells, but our study also provides a better understanding of allele-specific gene expression in colorectal cancer cells.     

百合属5种植物的比较转录组分析

百合属植物是多年生的宿根草本植物,具有重要的观赏价值,有的种可食用也可药用。本研究对百合属4种植物(细叶百合、川百合、卷丹和百合)花的转录组进行测序,并结合前期的野百合的转录组数据进行比较转录组分析,鉴别花发育相关基因。并筛选出5种植物的直系同源基因,进而检测正选择基因,分析它们对环境的适应性。通过组装分别在细叶百合、川百合、卷丹以及百合花中获得了44565、51413、41638和44716个unigene。为了了解unigene功能信息,将其在公共数据库中进行比对和注释,发现了13个花发育相关基因。另外,在5个物种间共获得8247对直系同源基因,并计算其Ka、Ks及Ka/Ks值。在5个物种中有33对直系同源基因受到了强烈的正选择作用,随后对其进行GO、KEGG功能富集分析及基因功能注释,发现这些基因主要是与抗性以及生物合成相关的基因。     

Patterns of divergence among conifer ESTs and polymorphism in Pinus sylvestris identify putative selective sweeps

Finding genes that are under positive selection is a difficult task, especially in non-model organisms. Here, we have analyzed expressed sequence tag (EST) data from 4 species (Pinus pinaster, Pinus taeda, Picea glauca, and Pseudotsuga menziesii) to investigate selection patterns during their evolution and to identify genes likely to be under positive selection. To confirm selection, population samples of these genes have been sequenced in Pinus sylvestris, a species that was not included in the EST data set. The estimates of branch-specific Ka/Ks (nonsynonymous/synonymous substitution rates) across all genes in the EST data set were similar or smaller than estimates from other higher plant species. There was no evidence for the traditional indication of positive selection, Ka/Ks above 1. However, several lines of evidence based on polymorphism patterns suggest that genes with high Ka/Ks (0.20-0.52) in the EST data set are in fact more affected by positive selection in P. sylvestris than genes with low Ka/Ks (0.01-0.04). The high Ka/Ks genes have a lower level of polymorphism and more negative Tajima's D than the low Ka/Ks genes. Further, in the high Ka/Ks group, the Hudson-Kreitman-Aguade test is significant. This suggests that the EST data set is a good starting point for finding genes under positive selection in conifers and that even moderate Ka/Ks values could be indicative of selection. A group of 5 genes with high Ka/Ks collectively show evidence for positive selection within P. sylvestris.     

Pericentromeric effects shape the patterns of divergence, retention, and expression of duplicated genes in the paleopolyploid soybean

The evolutionary forces that govern the divergence and retention of duplicated genes in polyploids are poorly understood. In this study, we first investigated the rates of nonsynonymous substitution (Ka) and the rates of synonymous substitution (Ks) for a nearly complete set of genes in the paleopolyploid soybean (Glycine max) by comparing the orthologs between soybean and its progenitor species Glycine soja and then compared the patterns of gene divergence and expression between pericentromeric regions and chromosomal arms in different gene categories. Our results reveal strong associations between duplication status and Ka and gene expression levels and overall low Ks and low levels of gene expression in pericentromeric regions. It is theorized that deleterious mutations can easily accumulate in recombination-suppressed regions, because of Hill-Robertson effects. Intriguingly, the genes in pericentromeric regions-the cold spots for meiotic recombination in soybean-showed significantly lower Ka and higher levels of expression than their homoeologs in chromosomal arms. This asymmetric evolution of two members of individual whole genome duplication (WGD)-derived gene pairs, echoing the biased accumulation of singletons in pericentromeric regions, suggests that distinct genomic features between the two distinct chromatin types are important determinants shaping the patterns of divergence and retention of WGD-derived genes.     

Allelic variation in gene expression identified through computational analysis of the dbEST database

Differential expression between the two alleles of an individual and between people with different genotypes has been commonly observed. Quantitative differences in gene expression between people may provide the genetic basis for the phenotypic difference between individuals and may be the primary cause of complex diseases. In this paper, we developed a computational method to identify genes that displayed allelic variation in gene expression in human EST libraries. To model allele-specific gene expression, we first identified EST libraries in which both A and B alleles were expressed and then identified allelic variation in gene expression based on the EST counts for each allele using a binomial test. Among 1107 SNPs that had a sufficient number of ESTs for the analysis, 524 (47%) displayed allelic variation in at least one cDNA library. We verified experimentally the allelic variation in gene expression for 6 of these SNPs. The frequency of allelic variation observed in EST libraries was similar to the previous studies using the SNP chip and primer extension method. We found that genes that displayed allelic variation were distributed throughout the human genome and were enriched in certain chromosome regions. The SNPs and genes identified in this study will provide a rich source for evaluating the effects of those SNPs and associated haplotypes in human health and diseases.     

Genome-wide identification,characterization and expression analysis of the LIM transcription factor family in quinoa

Lim family members play an important role in the regulation of plant cell development and stress response. However, there are few studies on LIM family in quinoa. In this study, we identified nine LIMS (named cqlim01-cqlim09) from quinoa, which were divided into three subfamilies (α Lim1, γ lim2 and δ lim2) according to phylogeny. The differences in gene structure and motif composition among different subfamilies have been observed. In addition, we studied the repetitive events of the members of the family. The Ka/Ks (non synchronous substitution rate / synchronous substitution rate) ratio analysis showed that the repetitive CqLIMs probably experienced purifying selection pressure. Promoter analysis showed that the family genes contained a variety of hormones, stress and tissue-specific expression elements, and protein interactions showed that these genes had actin stabilizing effect. In addition, QRT PCR results showed that all CqLIM genes were positively regulated under three stresses (low temperature, salt and ABA). These results provide a theoretical basis of further study of LIM gene in quinoa. Supplementary informationThe online version contains supplementary material available at 10.1007/s12298-021-00988-2.     

SNP-based large-scale identification of allele-specific gene expression in human B cells

Polymorphism and variations in gene expression provide the genetic basis for human variation. Allelic variation of gene expression, in particular, may play a crucial role in phenotypic variation and disease susceptibility. To identify genes with allelic expression in human cells, we genotyped genomic DNA and cDNA isolated from 31 immortalized B cell lines from three Centre d'Etude du Polymorphisme Humain (CEPH) families using high-density single-nucleotide polymorphism (SNP) chips containing 13,900 exonic SNPs. We identified seven SNPs in five genes with monoallelic expression, 146 SNPs in 125 genes with allelic imbalance in expression with preferentially higher expression of one allele in a heterozygous individual. The monoallelically expressed genes (ERAP2, MDGA1, LOC644422, SDCCAG3P1 and CLTCL1) were regulated by cis-acting, non-imprinted differential allelic control. In addition, all monoallelic gene expression patterns and allelic imbalances in gene expression in B cells were transmitted from parents to offspring in the pedigree, indicating genetic transmission of allelic gene expression. Furthermore, frequent allele substitution, probably due to RNA editing, was also observed in 21 genes in 23 SNPs as well as in 48 SNPs located in regions containing no known genes. In this study, we demonstrated that allelic gene expression is frequently observed in human B cells, and SNP chips are very useful tools for detecting allelic gene expression. Overall, our data provide a valuable framework for better understanding allelic gene expression in human B cells.     

Expression diversity and evolutionary dynamics of rice duplicate genes

Duplicate genes are believed to be a major source of new gene functions over evolutionary time. In order to evaluate the evolutionary dynamics of rice duplicate genes, formed principally by paleoployploidization prior to the speciation of the Poaceae family, we have employed a public microarray dataset including 155 gene expression omnibus sample plates and bioinformatics tools. At least 57.4% of old ~70 million years ago (MYA) duplicate gene pairs exhibit divergences in expression over the given experimental set, whereas at least 50.9% of young ~7.7-MYA duplicate gene pairs were shown to be divergent. When grouping the rice duplicate genes according to functional categories, we noted a striking and significant enrichment of divergent duplicate metabolism-associated genes, as compared to that observed in non-divergent duplicate genes. While both non-synonymous substitution (Ka) and synonymous substitution (Ks) values between non- and divergent duplicate gene pairs evidenced significant differences, the Ka/Ks values between them exhibited no significant differences. Interestingly, the average numbers of conserved motifs of the duplicate gene pairs revealed a pattern of decline along with an increase in expression diversity, partially supporting the subfunctionalization model with degenerative complementation in regulatory motifs. Duplicate gene pairs with high local similarity (HLS) segments, which might be formed via conversion between rice paleologs, evidenced higher expression correlations than were observed in the gene pairs without the HLS segments; this probably resulted in an increased likelihood of gene conversion in promoters of the gene pairs harboring HLS segments. More than 60% of the rice gene families exhibited similar high expression diversity between members as compared to that of randomly selected gene pairs. These findings are likely reflective of the evolutionary dynamics of rice duplicate genes for gene retention. An erratum to this article can be found at     

Novel and functional regulatory SNPs in the promoter region of <Emphasis Type="Italic">FOXP3</Emphasis> gene in a Gabonese population

Parasites exert a selection pressure on their hosts and are accountable for driving diversity within gene families and immune gene polymorphisms in a host population. The overwhelming response of regulatory T cells during infectious challenges directs the host immune system to lose the ability to mount parasite specific T cell responses. The underlying idea of this study is that regulatory single nucleotide polymorphism (SNPs) can cause significant changes in gene expression in functional immune genes. We identified and investigated regulatory SNPs in the promoter region of the FOXP3 gene in a group of Gabonese individuals exposed to a variety of parasitic infections. We identified two novel and one promoter variants in 40 individual subjects. We further validated these promoter variants for their allelic gene expression using transient transfection assays. Two promoter variants, −794 (C/G) and the other variant −734/−540 (C/T) revealed a significant higher expression of the reporter gene at basal level in comparison to the major allele. The identification of SNPs that modify function in the promoter regions could provide a basis for studying parasite susceptibility in association studies.     

Genome- Wide Analysis and Characterization of the TRX Gene Family in Upland Cotton

Thioredoxins (TRX) are small molecules of proteins that are present in all organisms. TRXs play an important role in diverse functions of plant growth and development. In this study, we performed genome-wide, characterization and expression levels of TRX gene family in cotton. A total of 150 GhTRX proteins were identified in upland cotton and classified into five subfamilies based on their domain compositions. Phylogenetic tree analysis divided TRX genes into seven subgroups. GhTRX genes covered all upland cotton chromosomes, with duplicated gene events. Ka/Ks ratio of three gene pairs was less than 1, suggesting purifying selection. The functions of GhTRX genes were studied using gene ontology, protein localization, and promoter analysis. Furthermore, six GhTRX genes were randomly selected to examine their expression level in cotton development and under various exogenous treatments. The genes showed high expressions in various tissues and at different stages of leaf senescence, also showed high expression under abscisic acid, ethylene, drought, and salinity. This study reveals the first report of TRX family genes in upland cotton. However further studies are needed to elucidate their specific functions in cotton plant.