Feature co-localization landscape of the human genome

Forward-reverse mutation cycles in cancer cell lines under chemical treatments

Spontaneous forward-reverse mutations were reported by us earlier in clinical samples from various types of cancers and in HeLa cells under normal culture conditions. To investigate the effects of chemical stimulations on such mutation cycles, the present study examined single nucleotide variations (SNVs) and copy number variations (CNVs) in HeLa and A549 cells exposed to wogonin-containing or acidic medium. In wogonin, both cell lines showed a mutation cycle during days 16-18. In acidic medium, both cell lines displayed multiple mutation cycles of different magnitudes. Genomic feature colocalization analysis suggests that CNVs tend to occur in expanded and unstable regions, and near promoters, histones, and non-coding transcription sites. Moreover, phenotypic variations in cell morphology occurred during the forward-reverse mutation cycles under both types of chemical treatments. In conclusion, chemical stresses imposed by wogonin or acidity promoted cyclic forward-reverse mutations in both HeLa and A549 cells to different extents.

Copy number variation profile-based genomic typing of premenstrual dysphoric disorder in Chinese

Premenstrual dysphoric disorder (PMDD) affects nearly 5% women of reproductive age. Symptomatic heterogeneity, together with largely unknown genetics, have greatly hindered its effective treatment. In the present study, analysis of genomic sequencing-based copy-number-variations (CNVs) called from 100-kb white blood cell DNA sequence windows by means of semi-supervised clustering led to the segregation of patient genomes into the D and V groups, which correlated with the depression and invasion clinical types respectively with 89.0% consistency. Application of diagnostic CNV features selected using the correlation-based machine-learning method enabled the classification of the CNVs obtained into the D group, V group, total-patient group and control group with an average accuracy of 83.0%. The power of the diagnostic CNV features was 0.98 on average, suggesting that these CNV features could be employed for the molecular diagnosis of the major clinical types of PMDD. This demonstrated concordnce between the CNV profiles and clinical types of PMDD supported the validity of symptom-based diagnosis of PMDD for differentiating between its two major clinical types, as well as the predominanly genetic nature of PMDD with a host of overlaps between multiple susceptibility genes/pathways and the diagnostic CNV features as indicators of involvement in PMDD etiology.

Genetic-variant hotspots and hotspot clusters in the human genome facilitating adaptation while increasing instability

Background

Genetic variants, underlining phenotypic diversity, are known to distribute unevenly in the human genome. A comprehensive understanding of the distributions of different genetic variants is important for insights into genetic functions and disorders.

Methods

Herein, a sliding-window scan of regional densities of eight kinds of germline genetic variants, including single-nucleotide-polymorphisms (SNPs) and four size-classes of copy-number-variations (CNVs) in the human genome has been performed.

Results

The study has identified 44,379 hotspots with high genetic-variant densities, and 1135 hotspot clusters comprising more than one type of hotspots, accounting for 3.1% and 0.2% of the genome respectively. The hotspots and clusters are found to co-localize with different functional genomic features, as exemplified by the associations of hotspots of middle-size CNVs with histone-modification sites, work with balancing and positive selections to meet the need for diversity in immune proteins, and facilitate the development of sensory-perception and neuroactive ligand-receptor interaction pathways in the function-sparse late-replicating genomic sequences. Genetic variants of different lengths co-localize with retrotransposons of different ages on a “long-with-young” and “short-with-all” basis. Hotspots and clusters are highly associated with tumor suppressor genes and oncogenes (p < 10⁻¹⁰), and enriched with somatic tumor CNVs and the trait- and disease-associated SNPs identified by genome-wise association studies, exceeding tenfold enrichment in clusters comprising SNPs and extra-long CNVs.

Conclusions

In conclusion, the genetic-variant hotspots and clusters represent two-edged swords that spearhead both positive and negative genomic changes. Their strong associations with complex traits and diseases also open up a potential “Common Disease-Hotspot Variant” approach to the missing heritability problem.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-021-00318-3.

Highly Recurrent Copy Number Variations in GABRB2 Associated With Schizophrenia and Premenstrual Dysphoric Disorder

OBJECTIVE

Although single-nucleotide polymorphisms in GABRB2, the gene encoding for GABAA receptors β2 subunit, have been associated with schizophrenia (SCZ), it is unknown whether there is any association of copy number variations (CNVs) in this gene with either SCZ or premenstrual dysphoric disorder (PMDD).

METHODS

In this study, the occurrences of the recurrent CNVs esv2730987 in Intron 6 and nsv1177513 in Exon 11 of GABRB2 in Chinese and German SCZ, and Chinese PMDD patients were compared to controls of same ethnicity and gender by quantitative PCR (qPCR).

RESULTS

The results demonstrated that copy-number-gains were enriched in both SCZ and PMDD patients with significant odds ratios (OR). For combined-gender SCZ patients versus controls, about two-fold increases were observed in both ethnic groups at both esv2730987 (OR = 2.15, p = 5.32E-4 in Chinese group; OR = 2.79, p = 8.84E-3 in German group) and nsv1177513 (OR = 3.29, p = 1.28E-11 in Chinese group; OR = 2.44, p = 6.17E-5 in German group). The most significant copy-number-gains were observed in Chinese females at nsv1177513 (OR = 3.41), and German females at esv2730987 (OR=3.96). Copy-number-gains were also enriched in Chinese PMDD patients versus controls at esv2730987 (OR = 10.53, p = 4.34E-26) and nsv1177513 (OR = 2.39, p = 3.19E-5).

CONCLUSION

These findings established for the first time the association of highly recurrent CNVs with SCZ and PMDD, suggesting the presence of an overlapping genetic basis with shared biomarkers for these two common psychiatric disorders.

Genomic subtyping of liver cancers with prognostic application

BACKGROUND

Cancer subtyping has mainly relied on pathological and molecular means. Massively parallel sequencing-enabled subtyping requires genomic markers to be developed based on global features rather than individual mutations for effective implementation.

METHODS

In the present study, the whole genome sequences (WGS) of 110 liver cancers of Japanese patients published with different pathologies were analyzed with respect to their single nucleotide variations (SNVs) comprising both gain-of-heterozygosity (GOH) and loss-of-heterozygosity (LOH) mutations, the signatures of combined GOH and LOH mutations, along with recurrent copy number variations (CNVs).

RESULTS

The results, obtained based on the WGS sequences as well as the Exome subset within the WGSs that covered ~ 2.0% of the WGS and the AluScan-subset within the WGSs that were amplifiable by Alu element-consensus primers and covered ~ 2.1% of the WGS, indicated that the WGS samples could be employed with the mutational parameters of SNV load, LOH%, the Signature α%, and survival-associated recurrent CNVs (srCNVs) as genomic markers for subtyping to stratify liver cancer patients prognostically into the long and short survival subgroups. The usage of the AluScan-subset data, which could be implemented with sub-micrograms of DNA samples and vastly reduced sequencing analysis task, outperformed the usage of WGS data when LOH% was employed as stratifying criterion.

CONCLUSIONS

Thus genomic subtyping performed with novel genomic markers identified in this study was effective in predicting patient-survival duration, with cohorts of hepatocellular carcinomas alone and those including intrahepatic cholangiocarcinomas. Such relatively heterogeneity-insensitive genomic subtyping merits further studies with a broader spectrum of cancers.

A nonparametric test for association with multiple loci in the retrospective case-control study

The genome-wide association studies aim at identifying common or rare variants associated with common diseases and explaining more heritability. It is well known that common diseases are influenced by multiple single nucleotide polymorphisms (SNPs) that are usually correlated in location or function. In order to powerfully detect association signals, it is highly desirable to take account of correlations or linkage disequilibrium (LD) information among multiple SNPs in testing for association. In this article, we propose a test SLIDE that depicts the difference of the average multi-locus genotypes between cases and controls and derive its variance-covariance matrix in the retrospective design. This matrix is composed of the pairwise LD between SNPs. Thus SLIDE can borrow the strength from an external database in the population of interest with a few thousands to hundreds of thousands individuals to improve the power for detecting association. Extensive simulations show that SLIDE has apparent superiority over the existing methods, especially in the situation involving both common and rare variants, both protective and deleterious variants. Furthermore, the efficiency of the proposed method is demonstrated in the application to the data from the Wellcome Trust Case Control Consortium.

Forward and reverse mutations in stages of cancer development

Background

Massive occurrences of interstitial loss of heterozygosity (LOH) likely resulting from gene conversions were found by us in different cancers as a type of single-nucleotide variations (SNVs), comparable in abundance to the commonly investigated gain of heterozygosity (GOH) type of SNVs, raising the question of the relationships between these two opposing types of cancer mutations.

Methods

In the present study, SNVs in 12 tetra sample and 17 trio sample sets from four cancer types along with copy number variations (CNVs) were analyzed by AluScan sequencing, comparing tumor with white blood cells as well as tissues vicinal to the tumor. Four published “nontumor”-tumor metastasis trios and 246 pan-cancer pairs analyzed by whole-genome sequencing (WGS) and 67 trios by whole-exome sequencing (WES) were also examined.

Results

Widespread GOHs enriched with CG-to-TG changes and associated with nearby CNVs and LOHs enriched with TG-to-CG changes were observed. Occurrences of GOH were 1.9-fold higher than LOH in “nontumor” tissues more than 2 cm away from the tumors, and a majority of these GOHs and LOHs were reversed in “paratumor” tissues within 2 cm of the tumors, forming forward-reverse mutation cycles where the revertant LOHs displayed strong lineage effects that pointed to a sequential instead of parallel development from “nontumor” to “paratumor” and onto tumor cells, which was also supported by the relative frequencies of 26 distinct classes of CNVs between these three types of cell populations.

Conclusions

These findings suggest that developing cancer cells undergo sequential changes that enable the “nontumor” cells to acquire a wide range of forward mutations including ones that are essential for oncogenicity, followed by revertant mutations in the “paratumor” cells to avoid growth retardation by excessive mutation load. Such utilization of forward-reverse mutation cycles as an adaptive mechanism was also observed in cultured HeLa cells upon successive replatings. An understanding of forward-reverse mutation cycles in cancer development could provide a genomic basis for improved early diagnosis, staging, and treatment of cancers.

Electronic supplementary material

The online version of this article (10.1186/s40246-018-0170-6) contains supplementary material, which is available to authorized users.

Genome-wide characterization of the aldehyde dehydrogenase gene superfamily in soybean and its potential role in drought stress response

BACKGROUND

Aldehyde dehydrogenases (ALDHs) represent a group of enzymes that detoxify aldehydes by facilitating their oxidation to carboxylic acids, and have been shown to play roles in plant response to abiotic stresses. However, the comprehensive analysis of ALDH superfamily in soybean (Glycine max) has been limited.

RESULTS

In present study, a total of 53 GmALDHs were identified in soybean, and grouped into 10 ALDH families according to the ALDH Gene Nomenclature Committee and phylogenetic analysis. These groupings were supported by their gene structures and conserved motifs. Soybean ALDH superfamily expanded mainly by whole genome duplication/segmental duplications. Gene network analysis identified 1146 putative co-functional genes of 51 GmALDHs. Gene Ontology (GO) enrichment analysis suggested the co-functional genes of these 51 GmALDHs were enriched (FDR < 1e-3) in the process of lipid metabolism, photosynthesis, proline catabolism, and small molecule catabolism. In addition, 22 co-functional genes of GmALDHs are related to plant response to water deprivation/water transport. GmALDHs exhibited various expression patterns in different soybean tissues. The expression levels of 13 GmALDHs were significantly up-regulated and 14 down-regulated in response to water deficit. The occurrence frequencies of three drought-responsive cis-elements (ABRE, CRT/DRE, and GTGCnTGC/G) were compared in GmALDH genes that were up-, down-, or non-regulated by water deficit. Higher frequency of these three cis-elements was observed for the group of up-regulated GmALDH genes as compared to the group of down- or non- regulated GmALDHs by drought stress, implying their potential roles in the regulation of soybean response to drought stress.

CONCLUSIONS

A total of 53 ALDH genes were identified in soybean genome and their phylogenetic relationships and duplication patterns were analyzed. The potential functions of GmALDHs were predicted by analyses of their co-functional gene networks, gene expression profiles, and cis-regulatory elements. Three GmALDH genes, including GmALDH3H2, GmALDH12A2 and GmALDH18B3, were highly induced by drought stress in soybean leaves. Our study provides a foundation for future investigations of GmALDH gene function in soybean.

Tripartite genome of all species

Neutral theory has dominated the molecular evolution field for more than half a century, but it has been severely challenged by the recently emerged Maximum Genetic Diversity (MGD) theory. However, based on our recent work of tripartite human genome architecture, we found that MGD theory may have overlooked the regulatory but variable genomic regions that increase with species complexity. Here we propose a new molecular evolution theory named Increasing Functional Variation (IFV) hypothesis. According to the IFV hypothesis, the genome of all species is divided into three regions that are ‘functional and invariable’, ‘functional and variable’ and ‘non-functional and variable’. While the ‘non-functional and variable’ region decreases as species become more complex, the other two regions increase.