Using pseudogene database to identify lineage-specific genes and pseudogenes in humans and chimpanzees

Genome-scale evolution in local populations of wild chimpanzees

Analysis of genome-scale evolution has been difficult in large, endangered animals because opportunities to collect high-quality genetic samples are limited. There is a need for novel field-friendly, cost-effective genetic techniques. This study conducted an exome-wide analysis of a total of 42 chimpanzees (Pan troglodytes) across six African regions, providing insights into population discrimination techniques. Wild chimpanzee DNA was extracted noninvasively from collected fecal samples using the lysis-buffer storage method. To target genome-scale regions of host DNA, exome-capture sequencing was performed using cost-effective baits originally designed for humans (closely related to chimpanzees). Multivariate analysis effectively discriminated differences in local populations, aiding in the identification of samples' geographical origins. Exome-wide heterozygosity was negatively correlated significantly with genome-wide nonsynonymous-synonymous substitution ratios, suggesting that mutation loads exist at the local population level. Exome sequences revealed functional diversity and protein-coding gene divergence. Segregating pseudogenes were comprehensively annotated, with many being population-specific and others shared among populations. Focusing on multicopy chemosensory receptor genes, the segregating pseudogenes OR7D4 (an olfactory receptor) and TAS2R42 (a bitter taste receptor) were shared among western and eastern chimpanzees. Overall, our analytical framework offers ecological insights into chimpanzees and may be applicable to other organisms.

Identification, characterization and expression analysis of lineage-specific genes within mangrove species Aegiceras corniculatum

Lineage-specific genes (LSGs) are the genes that have no recognizable homology to any sequences in other species, which are important drivers for the generation of new functions, phenotypic changes, and facilitating species adaptation to environment. Aegiceras corniculatum is one of major mangrove plant species adapted to waterlogging and saline conditions, and the exploration of aegiceras-specific genes (ASGs) is important to reveal its adaptation to the harsh environment. Here, we performed a systematic analysis on ASGs, focusing on their sequence characterization, origination and expression patterns. Our results reveal that there are 4823 ASGs in the genome, approximately 11.84% of all protein-coding genes. High proportion (45.78%) of ASGs originate from gene duplication, and the time of gene duplication of ASGs is consistent with the timing of two genome-wide replication (WGD) events that occurred in A. corniculatum, and also coincides with a short period of global warming during the Paleocene-Eocene Maximum (PETM, 55.5 million years ago). Gene structure analysis showed that ASGs have shorter protein lengths, fewer exons, and higher isoelectric point. Expression patterns analysis showed that ASGs had low levels of expression and more tissue-specific expression. Weighted gene co-expression network analysis (WGCNA) revealed that 86 ASGs co-expressed gene modules were primarily involved in pathways related to adversity stress, including plant hormone signal transduction, phenylpropanoid biosynthesis, photosynthesis, peroxisome and pentose phosphate pathway. This study provides a comprehensive analysis of the characteristics and potential functions of ASGs and identifies key candidate genes, which will contribute to the subsequent further investigation of the adaptation of A. corniculatum to intertidal coastal wetland habitats.

Comprehensive elaboration of database resources utilized in next-generation sequencing-based tumor somatic mutation detection

The rapid evolution of next-generation sequencing (NGS)-based tumor genomic profile detection and the emergence of molecularly targeted therapies have enabled precision oncology. In NGS-based analysis, various types of databases have been developed to perform different functions. However, many problems still exist when using these public databases. Therefore, it is important to better understand the characteristics and limitations of each database and have them complement each other to provide useful clinical evidence for NGS testing. In this review, we elaborate on the important role of databases and their concrete applications in NGS-based somatic mutation detection. We introduce the typically used databases for sequence alignment, variant filtration, and variant interpretation, and compare the differences between the databases with similar functions. Subsequently, we determine the limitations of each database and provide the corresponding solutions. Furthermore, we present an overview diagram to clearly illustrate the database used in the entire NGS-based somatic mutation detection pipeline.

Major Histocompatibility Complex Genes Map to Two Chromosomes in an Evolutionarily Ancient Reptile, the Tuatara Sphenodon punctatus

Major histocompatibility complex (MHC) genes are a central component of the vertebrate immune system and usually exist in a single genomic region. However, considerable differences in MHC organization and size exist between different vertebrate lineages. Reptiles occupy a key evolutionary position for understanding how variation in MHC structure evolved in vertebrates, but information on the structure of the MHC region in reptiles is limited. In this study, we investigate the organization and cytogenetic location of MHC genes in the tuatara (Sphenodon punctatus), the sole extant representative of the early-diverging reptilian order Rhynchocephalia. Sequencing and mapping of 12 clones containing class I and II MHC genes from a bacterial artificial chromosome library indicated that the core MHC region is located on chromosome 13q. However, duplication and translocation of MHC genes outside of the core region was evident, because additional class I MHC genes were located on chromosome 4p. We found a total of seven class I sequences and 11 class II β sequences, with evidence for duplication and pseudogenization of genes within the tuatara lineage. The tuatara MHC is characterized by high repeat content and low gene density compared with other species and we found no antigen processing or MHC framework genes on the MHC gene-containing clones. Our findings indicate substantial differences in MHC organization in tuatara compared with mammalian and avian MHCs and highlight the dynamic nature of the MHC. Further sequencing and annotation of tuatara and other reptile MHCs will determine if the tuatara MHC is representative of nonavian reptiles in general.

Genetic variation in taste receptor pseudogenes provides evidence for a dynamic role in human evolution

BACKGROUND

Human bitter taste receptors are encoded by a gene family consisting of 25 functional TAS2R loci. In addition, humans carry 11 TAS2R pseudogenes, some of which display evidence for substantial diversification among species, showing lineage-specific loss of function. Since bitter taste is thought to help prevent the intake of toxic substances, diversity at TAS2R genes could reflect the action of natural selection on the ability to recognize some bitter compounds rather than others. Whether species-specific variation in TAS2R pseudogenes is solely the result of genetic drift or whether it may have been influenced by selection due to different feeding behaviors has been an open question.

RESULTS

In this study, we analyzed patterns of variation at human TAS2R pseudogenes in both African and non-African populations, and compared them to those observable in nonhuman primates and archaic human species. Our results showed a similar worldwide distribution of allelic variation for most of the pseudogenes, with the exception of the TAS2R6P and TAS2R18P loci, both of which presented an unexpected higher frequency of derived alleles outside Africa. At the TAS2R6P locus, two SNPs were found in strong linkage disequilibrium (r2 > 0.9) with variants in the functional TAS2R5 gene, which showed signatures of selection. The human TAS2R18P carried a species-specific stop-codon upstream of four polymorphic insertions in the reading frame. SNPs at this locus showed significant positive values in a number of neutrality statistics, and age estimates indicated that they arose after the homo-chimp divergence.

CONCLUSIONS

The similar distribution of variation of many human bitter receptor pseudogenes among human populations suggests that they arose from the ancestral forms by a unidirectional loss of function. However we explain the higher frequency of TAS2R6P derived alleles outside Africa as the effect of the balancing selection acting on the closely linked TAS2R5 gene. In contrast, TAS2R18P displayed a more complex history, suggesting an acquired function followed by a recent pseudogenization that predated the divergence of human modern and archaic species, which we hypothesize was associated with adaptions to dietary changes.