An evolving view of copy number variants

Selection for robust metabolism in domesticated yeasts is driven by adaptation to Hsp90 stress

Protein folding both promotes and constrains adaptive evolution. We uncover this surprising duality in the role of the protein-folding chaperone heat shock protein 90 (Hsp90) in maintaining the integrity of yeast metabolism amid proteotoxic stressors within industrial domestication niches. Ethanol disrupts critical Hsp90-dependent metabolic pathways and exerts strong selective pressure for redundant duplications of key genes within these pathways, yielding the classical genomic signatures of beer and bread domestication. This work demonstrates a mechanism of adaptive canalization in an ecology of major economic importance and highlights Hsp90-dependent variation as an important source of phantom heritability in complex traits.

Analysis of copy number variants detected by sequencing in spontaneous abortion

BACKGROUND

The incidence of spontaneous abortion (SA), which affects approximately 15-20% of pregnancies, is the most common complication of early pregnancy. Pathogenic copy number variations (CNVs) are recognized as potential genetic causes of SA. However, CNVs of variants of uncertain significance (VOUS) have been identified in products of conceptions (POCs), and their correlation with SA remains uncertain.

RESULTS

Of 189 spontaneous abortion cases, trisomy 16 was the most common numerical chromosome abnormality, followed by monosomy X. CNVs most often occurred on chromosomes 4 and 8. Gene Ontology and signaling pathway analysis revealed significant enrichment of genes related to nervous system development, transmembrane transport, cell adhesion, and structural components of chromatin. Furthermore, genes within the VOUS CNVs were screened by integrating human placental expression profiles, PhyloP scores, and Residual Variance Intolerance Score (RVIS) percentiles to identify potential candidate genes associated with spontaneous abortion. Fourteen potential candidate genes (LZTR1, TSHZ1, AMIGO2, H1-4, H2BC4, H2AC7, H3C8, H4C3, H3C6, PHKG2, PRR14, RNF40, SRCAP, ZNF629) were identified. Variations in LZTR1, TSHZ1, and H4C3 may contribute to embryonic lethality.

CONCLUSIONS

CNV sequencing (CNV-seq) analysis is an effective technique for detecting chromosomal abnormalities in POCs and identifying potential candidate genes for SA.

Template switching between the leading and lagging strands at replication forks generates inverted copy number variants through hairpin-capped extrachromosomal DNA

Inherited and germ-line de novo copy number variants (CNVs) are increasingly found to be correlated with human developmental and cancerous phenotypes. Several models for template switching during replication have been proposed to explain the generation of these gross chromosomal rearrangements. We proposed a model of template switching (ODIRA-origin dependent inverted repeat amplification) in which simultaneous ligation of the leading and lagging strands at diverging replication forks could generate segmental inverted triplications through an extrachromosomal inverted circular intermediate. Here, we created a genetic assay using split-ura3 cassettes to trap the proposed inverted intermediate. However, instead of recovering circular inverted intermediates, we found inverted linear chromosomal fragments ending in native telomeres-suggesting that a template switch had occurred at the centromere-proximal fork of a replication bubble. As telomeric inverted hairpin fragments can also be created through double strand breaks we tested whether replication errors or repair of double stranded DNA breaks were the most likely initiating event. The results from CRISPR/Cas9 cleavage experiments and growth in the replication inhibitor hydroxyurea indicate that it is a replication error, not a double stranded break that creates the inverted junctions. Since inverted amplicons of the SUL1 gene occur during long-term growth in sulfate-limited chemostats, we sequenced evolved populations to look for evidence of linear intermediates formed by an error in replication. All of the data are compatible with a two-step version of the ODIRA model in which sequential template switching at short inverted repeats between the leading and lagging strands at a replication fork, followed by integration via homologous recombination, generates inverted interstitial triplications.

Genetic analysis of chorionic villus tissues in early missed abortions

Chromosomal abnormalities are the most common etiology of early spontaneous miscarriage. However, traditional karyotyping of chorionic villus samples (CVSs) is limited by cell culture and its low resolution. The objective of our study was to investigate the efficiency of molecular karyotyping technology for genetic diagnosis of early missed abortion tissues. Chromosome analysis of 1191 abortion CVSs in early pregnancy was conducted from August 2016 to June 2021; 463 cases were conducted via copy-number variations sequencing (CNV-seq)/quantitative fluorescent-polymerase chain reaction (QF-PCR) and 728 cases were conducted using SNP array. Clinically significant CNVs of CVSs were identified to clarify the cause of miscarriage and to guide the couples' subsequent pregnancies. Among these, 31 cases with significant maternal cell contamination were removed from the study. Among the remaining 1160 samples, 751 cases (64.7%) with genetic abnormalities were identified, of which, 531 (45.8%) were single aneuploidies, 31 (2.7%) were multiple aneuploidies, 50 (4.3%) were polyploidies, 54 (4.7%) were partial aneuploidies, 77 (6.6%) had submicroscopic CNVs (including 25 with clinically significant CNVs and 52 had variants of uncertain significance), and 8 cases (0.7%) were uniparental disomies. Our study suggests that both SNP array and CNV-seq/QF-PCR are reliable, robust, and high-resolution technologies for genetic diagnosis of miscarriage.

Copy number variation of the restorer Rf4 underlies human selection of three-line hybrid rice breeding

Cytoplasmic male sterility (CMS) lines are important for breeding hybrid crops, and utilization of CMS lines requires strong fertility restorer (Rf) genes. Rf4, a major Rf for Wild-Abortive CMS (CMS-WA), has been cloned in rice. However, the Rf4 evolution and formation of CMS-WA/Rf system remain elusive. Here, we show that the Rf4 locus emerges earlier than the CMS-WA gene WA352 in wild rice, and 69 haplotypes of the Rf4 locus are generated in the Oryza genus through the copy number and sequence variations. Eight of these haplotypes of the Rf4 locus are enriched in modern rice cultivars during natural and human selections, whereas non-functional rf4i is preferentially selected for breeding current CMS-WA lines. We further verify that varieties carrying two-copy Rf4 haplotype have stronger fertility restoration ability and are widely used in three-line hybrid rice breeding. Our findings increase our understanding of CMS/Rf systems and will likely benefit crop breeding.

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast

Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene overexpression, revealing that the fitness costs of copy-number variation can vary substantially with genetic background in a common-garden environment. But the interplay between copy-number variation tolerance and environment remains unexplored on a genomic scale. Here, we measured the tolerance to gene overexpression in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride stress. Overexpressed genes that are commonly deleterious during sodium chloride stress recapitulated those commonly deleterious under standard conditions. However, sodium chloride stress uncovered novel differences in strain responses to gene overexpression. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to sodium chloride stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to overexpression of specific genes. Although most genes were deleterious, hundreds were beneficial when overexpressed-remarkably, most of these effects were strain specific. Few beneficial genes were shared between the sodium chloride-sensitive isolates, implicating mechanistic differences behind their sodium chloride sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural copy-number variation of a sodium export pump that likely contributes to strain-specific responses to overexpression of other genes. Our results reveal extensive strain-by-environment interactions in the response to gene copy-number variation, raising important implications for the accessibility of copy-number variation-dependent evolutionary routes under times of stress.

Evolutionary Dynamics of FLC-like MADS-Box Genes in Brassicaceae

MADS-box genes encode transcription factors that play important roles in the development and evolution of plants. There are more than a dozen clades of MADS-box genes in angiosperms, of which those with functions in the specification of floral organ identity are especially well-known. From what has been elucidated in the model plant Arabidopsis thaliana, the clade of FLC-like MADS-box genes, comprising FLC-like genes sensu strictu and MAF-like genes, are somewhat special among the MADS-box genes of plants since FLC-like genes, especially MAF-like genes, show unusual evolutionary dynamics, in that they generate clusters of tandemly duplicated genes. Here, we make use of the latest genomic data of Brassicaceae to study this remarkable feature of the FLC-like genes in a phylogenetic context. We have identified all FLC-like genes in the genomes of 29 species of Brassicaceae and reconstructed the phylogeny of these genes employing a Maximum Likelihood method. In addition, we conducted selection analyses using PAML. Our results reveal that there are three major clades of FLC-like genes in Brassicaceae that all evolve under purifying selection but with remarkably different strengths. We confirm that the tandem arrangement of MAF-like genes in the genomes of Brassicaceae resulted in a high rate of duplications and losses. Interestingly, MAF-like genes also seem to be prone to transposition. Considering the role of FLC-like genes sensu lato (s.l.) in the timing of floral transition, we hypothesize that this rapid evolution of the MAF-like genes was a main contributor to the successful adaptation of Brassicaceae to different environments.

Multi-Omics Analysis of NCI-60 Cell Line Data Reveals Novel Metabolic Processes Linked with Resistance to Alkylating Anti-Cancer Agents

This study aimed to elucidate the molecular determinants influencing the response of cancer cells to alkylating agents, a major class of chemotherapeutic drugs used in cancer treatment. The study utilized data from the National Cancer Institute (NCI)-60 cell line screening program and employed a comprehensive multi-omics approach integrating transcriptomic, proteomic, metabolomic, and SNP data. Through integrated pathway analysis, the study identified key metabolic pathways, such as cysteine and methionine metabolism, starch and sucrose metabolism, pyrimidine metabolism, and purine metabolism, that differentiate drug-sensitive and drug-resistant cancer cells. The analysis also revealed potential druggable targets within these pathways. Furthermore, copy number variant (CNV) analysis, derived from SNP data, between sensitive and resistant cells identified notable differences in genes associated with metabolic changes (WWOX, CNTN5, DDAH1, PGR), protein trafficking (ARL17B, VAT1L), and miRNAs (MIR1302-2, MIR3163, MIR1244-3, MIR1302-9). The findings of this study provide a holistic view of the molecular landscape and dysregulated pathways underlying the response of cancer cells to alkylating agents. The insights gained from this research can contribute to the development of more effective therapeutic strategies and personalized treatment approaches, ultimately improving patient outcomes in cancer treatment.

Copy number variation alters local and global mutational tolerance

Copy number variants (CNVs), duplications and deletions of genomic sequences, contribute to evolutionary adaptation but can also confer deleterious effects and cause disease. Whereas the effects of amplifying individual genes or whole chromosomes (i.e., aneuploidy) have been studied extensively, much less is known about the genetic and functional effects of CNVs of differing sizes and structures. Here, we investigated Saccharomyces cerevisiae (yeast) strains that acquired adaptive CNVs of variable structures and copy numbers following experimental evolution in glutamine-limited chemostats. Although beneficial in the selective environment, CNVs result in decreased fitness compared with the euploid ancestor in rich media. We used transposon mutagenesis to investigate mutational tolerance and genome-wide genetic interactions in CNV strains. We find that CNVs increase mutational target size, confer increased mutational tolerance in amplified essential genes, and result in novel genetic interactions with unlinked genes. We validated a novel genetic interaction between different CNVs and BMH1 that was common to multiple strains. We also analyzed global gene expression and found that transcriptional dosage compensation does not affect most genes amplified by CNVs, although gene-specific transcriptional dosage compensation does occur for ∼12% of amplified genes. Furthermore, we find that CNV strains do not show previously described transcriptional signatures of aneuploidy. Our study reveals the extent to which local and global mutational tolerance is modified by CNVs with implications for genome evolution and CNV-associated diseases, such as cancer.

Regulatory role of Chitinase 3-like 1 gene in papillary thyroid carcinoma proved by integration analyses of single-cell sequencing with cohort and experimental validations

Papillary thyroid carcinoma (PTC) is one of the most common thyroid carcinomas. The gross extrathyroidal extension and extensive metastases of PTC lead to high rates of recurrence and poor clinical outcomes. However, the mechanisms underlying PTC development are poorly understood. In this study, using single-cell RNA sequencing, the transcriptome profiles of two PTC patients were addressed, including PTC1 with low malignancy and good prognosis and PTC2 with high malignancy and poor prognosis. We found that epithelial subcluster Epi02 was the most associated with the malignant development of PTC cells, with which the fold change of Chitinase 3-like 1 (CHI3L1) is on the top of the differentially expressed genes between PTC1 and PTC2 (P < 0.001). However CHI3L1 is rarely investigated in PTC as far. We then studied its role in PTC with a series of experiments. Firstly, qRT-PCR analysis of 14 PTC patients showed that the expression of CHI3L1 was positively correlated with malignancy. In addition, overexpression or silencing of CHI3L1 in TPC-1 cells, a PTC cell line, cultured in vitro showed that the proliferation, invasion, and metastasis of the cells were promoted or alleviated by CHI3L1. Further, immunohistochemistry analysis of 110 PTC cases revealed a significant relationship between CHI3L1 protein expression and PTC progression, especially the T (P < 0.001), N (P < 0.001), M stages (P = 0.007) and gross ETE (P < 0.001). Together, our results prove that CHI3L1 is a positive regulator of malignant development of PTC, and it promotes proliferation, invasion, and metastasis of PTC cells. Our study improves understanding of the molecular mechanisms underlying the progression of PTC and provides new insights for the clinical diagnosis and treatment of PTC.

Ethanol Drives Evolution of Hsp90-Dependent Robustness by Redundancy in Yeast Domestication

Protein folding promotes and constrains adaptive evolution. We uncover this surprising duality in the role the protein-folding chaperone Hsp90 plays in mediating the interplay between proteome and the genome which acts to maintain the integrity of yeast metabolism in the face of proteotoxic stressors in anthropic niches. Of great industrial relevance, ethanol concentrations generated by fermentation in the making of beer and bread disrupt critical Hsp90-dependent nodes of metabolism and exert strong selective pressure for increased copy number of key genes encoding components of these nodes, yielding the classical genetic signatures of beer and bread domestication. This work establishes a mechanism of adaptive canalization in an ecology of major economic significance and highlights Hsp90-contingent variation as an important source of phantom heritability in complex traits.

Nanopore sequencing of PCR products enables multicopy gene family reconstruction

The importance of gene amplifications in evolution is more and more recognized. Yet, tools to study multi-copy gene families are still scarce, and many such families are overlooked using common sequencing methods. Haplotype reconstruction is even harder for polymorphic multi-copy gene families. Here, we show that all variants (or haplotypes) of a multi-copy gene family present in a single genome, can be obtained using Oxford Nanopore Technologies sequencing of PCR products, followed by steps of mapping, SNP calling and haplotyping. As a proof of concept, we acquired the sequences of highly similar variants of the cidA and cidB genes present in the genome of the Wolbachia wPip, a bacterium infecting Culex pipiens mosquitoes. Our method relies on a wide database of cid genes, previously acquired by cloning and Sanger sequencing. We addressed problems commonly faced when using mapping approaches for multi-copy gene families with highly similar variants. In addition, we confirmed that PCR amplification causes frequent chimeras which have to be carefully considered when working on families of recombinant genes. We tested the robustness of the method using a combination of bioinformatics (read simulations) and molecular biology approaches (sequence acquisitions through cloning and Sanger sequencing, specific PCRs and digital droplet PCR). When different haplotypes present within a single genome cannot be reconstructed from short reads sequencing, this pipeline confers a high throughput acquisition, gives reliable results as well as insights of the relative copy numbers of the different variants.

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast

Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation (CNV) provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene over-expression (OE), revealing that the fitness costs of CNV can vary substantially with genetic background in a common-garden environment. But the interplay between CNV tolerance and environment remains unexplored on a genomic scale. Here we measured the tolerance to gene OE in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride (NaCl) stress. OE genes that are commonly deleterious during NaCl stress recapitulated those commonly deleterious under standard conditions. However, NaCl stress uncovered novel differences in strain responses to gene OE. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to NaCl stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to gene OE. Although most genes were deleterious, hundreds were beneficial when overexpressed - remarkably, most of these effects were strain specific. Few beneficial genes were shared between the NaCl-sensitive isolates, implicating mechanistic differences behind their NaCl sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural CNV of a sodium export pump that likely contributes to strain-specific responses to OE of other genes. Our results reveal extensive strain-by-environment interaction in the response to gene CNV, raising important implications for the accessibility of CNV-dependent evolutionary routes under times of stress.

Goat Pleomorphic Adenoma Gene 1 (PLAG1): mRNA Expression, CNV Detection and Associations with Growth Traits

The pleomorphic adenoma gene 1 (PLAG1) gene, as the major gene responsible for growth, plays a vital role in myogenesis. Meanwhile, the relationship between copy number variation (CNV) of this gene and growth traits in goats remains unclear. Therefore, this study investigated four aspects: bioinformatics analysis, mRNA expression (n = 6), CNV detection (n = 224), and association analysis. The findings indicated that the gene had a large number of conserved motifs, and the gene expression level was higher in fetal goats than in adult goats. Three CNV loci were selected from the database, among which CNV1 was located in the bidirectional promoter region and was associated with goat growth traits. CNV analysis showed that CNV2 and CNV3 of the PLAG1 gene were associated with growth traits such as body weight, heart girth, height at hip cross, and hip width (p < 0.05), with CNV1 loss genotype being the superior genotype, and CNV2 and CNV3 median and gain genotypes of being superior genotypes. This finding further confirms that the PLAG1 gene is the dominant gene for growth traits, which will serve as theoretical guidance for goat breeding.

Rare CNVs and phenome-wide profiling highlight brain structural divergence and phenotypical convergence

Copy number variations (CNVs) are rare genomic deletions and duplications that can affect brain and behaviour. Previous reports of CNV pleiotropy imply that they converge on shared mechanisms at some level of pathway cascades, from genes to large-scale neural circuits to the phenome. However, existing studies have primarily examined single CNV loci in small clinical cohorts. It remains unknown, for example, how distinct CNVs escalate vulnerability for the same developmental and psychiatric disorders. Here we quantitatively dissect the associations between brain organization and behavioural differentiation across 8 key CNVs. In 534 CNV carriers, we explored CNV-specific brain morphology patterns. CNVs were characteristic of disparate morphological changes involving multiple large-scale networks. We extensively annotated these CNV-associated patterns with ~1,000 lifestyle indicators through the UK Biobank resource. The resulting phenotypic profiles largely overlap and have body-wide implications, including the cardiovascular, endocrine, skeletal and nervous systems. Our population-level investigation established brain structural divergences and phenotypical convergences of CNVs, with direct relevance to major brain disorders.

Targeted nanopore sequencing enables complete characterisation of structural deletions initially identified using exon-based short-read sequencing strategies

BACKGROUND

The widespread adoption of exome sequencing has greatly increased the rate of genetic diagnosis for inherited conditions. However, the detection and validation of large deletions remains challenging. While numerous bioinformatics approaches have been developed to detect deletions from whole - exome sequencing and targeted panels, further work is typically required to define the physical breakpoints or integration sites. Accurate characterisation requires either expensive follow - up whole - genome sequencing or the time - consuming, laborious process of PCR walking, both of which are challenging when dealing with the repeat sequences which frequently intersect deletion breakpoints. The aim of this study was to develop a cost-effective, long-range sequencing method to characterise deletions.

METHODS

Genomic DNA was amplified with primers spanning the deletion using long-range PCR and the products purified. Sequencing was performed on MinION flongle flowcells. The resulting fast5 files were basecalled using Guppy, trimmed using Porechop and aligned using Minimap2. Filtering was performed using NanoFilt. Nanopore sequencing results were verified by Sanger sequencing.

RESULTS

Four cases with deletions detected following comparative read-depth analysis of targeted short-read sequencing were analysed. Nanopore sequencing defined breakpoints at the molecular level in all cases including homozygous breakpoints in EYS, CNGA1 and CNGB1 and a heterozygous deletion in PRPF31. All breakpoints were verified by Sanger sequencing.

CONCLUSIONS

In this study, a quick, accurate and cost - effective method is described to characterise deletions identified from exome, and similar data, using nanopore sequencing.

Multi-omics profiling and digital image analysis reveal the potential prognostic and immunotherapeutic properties of CD93 in stomach adenocarcinoma

Background

Recent evidence highlights the fact that immunotherapy has significantly improved patient outcomes. CD93, as a type I transmembrane glycoprotein, was correlated with tumor-associated angiogenesis; however, how CD93 correlates with immunotherapy in stomach adenocarcinoma (STAD) remains unclear.

Methods

TCGA, GTEx, GEO, TIMER2.0, HPA, TISIDB, TCIA, cBioPortal, LinkedOmics, and ImmuCellAI public databases were used to elucidate CD93 in STAD. Visualization and statistical analysis of data were performed by R (Version 4.1.3), GraphPad (Version 8.0.1), and QuPath (Version 0.3.2).

Results

CD93 was highly expressed in STAD compared with adjacent normal tissues. The overexpression of CD93 was significantly correlated with a poor prognosis in STAD. There was a negative correlation between CD93 expression levels with CD93 mutation and methylation in STAD. Our results revealed that CD93 expression was positively associated with most immunosuppressive genes (including PD-1, PD-L1, CTLA-4, and LAG3), immunostimulatory genes, HLA, chemokine, and chemokine receptor proteins in STAD. Furthermore, in STAD, CD93 was noticeably associated with the abundance of multiple immune cell infiltration levels. Functional HALLMARK and KEGG term enhancement analysis of CD93 through Gene Set Enrichment Analysis was correlated with the process of the angiogenesis pathway. Subsequently, digital image analysis results by QuPath revealed that the properties of CD93⁺ cells were statistically significant in different regions of stomach cancer and normal stomach tissue. Finally, we utilized external databases, including GEO, TISIDB, ImmuCellAI, and TCIA, to validate that CD93 plays a key role in the immunotherapy of STAD.

Conclusion

Our study reveals that CD93 is a potential oncogene and is an indicative biomarker of a worse prognosis and exerts its immunomodulatory properties and potential possibilities for immunotherapy in STAD.

Analysis of eEF1A2 gene expression and copy number in cervical carcinoma

OBJECTIVE

To explore and analyze the expression of eukaryotic translation elongation factor 1 alpha 2 (eEF1A2) gene in cervical cancer tissues, its relationship with patient survival, gene mutations, and changes in copy number in cervical cancer and chronic cervicitis tissues.

METHODS

The expression of the eEF1A2 gene in cervical cancer and its relationship with patient survival were analyzed using gene expression profile interactive analysis. Changes in eEF1A2 expression in cervical cancer tissues were analyzed using cBioPortal, a portal for cancer genomics analysis. The eEF1A2 copy number in cervical cancer tissues and chronic cervicitis tissues was determined by real-time fluorescence quantitative polymerase chain reaction. The relationship between the expression of eEF1A2 protein and the clinical stage, pathological grade, and patient survival of cervical cancer was analyzed by the database: The Human Protein Atlas, an integrated repository portal for tumor-immune system interactions.

RESULTS

Gene expression profile interactive analysis database analysis showed no significant differences in the expression of eEF1A2 between cervical cancer and normal cervical tissues (P > .05). The eEF1A2 gene expression level was not correlated with the survival of cervical cancer patients (P > .05). Analysis of the cBioPortal database showed that 18 of 297 cervical cancer patients had eEF1A2 gene changes, including missense mutation, splice mutation, amplification, and messenger RNA increase. There was no significant difference in eEF1A2 gene copy number between cervical cancer and chronic cervicitis (P > .05). The Human Protein Atlas and an integrated repository portal for tumor-immune system interactions database analysis of immunohistochemical data showed that eEF1A2 protein expression was no significant difference in clinical stage, pathological grade and patient survival of cervical cancer (P > .05).

CONCLUSION

The eEF1A2 gene was mutated in cervical cancer tissues. The eEF1A2 gene copy number was not associated with changes in the expression of the eEF1A2 gene in cervical cancer tissues.

Inferencing Bulk Tumor and Single-Cell Multi-Omics Regulatory Networks for Discovery of Biomarkers and Therapeutic Targets

There are insufficient accurate biomarkers and effective therapeutic targets in current cancer treatment. Multi-omics regulatory networks in patient bulk tumors and single cells can shed light on molecular disease mechanisms. Integration of multi-omics data with large-scale patient electronic medical records (EMRs) can lead to the discovery of biomarkers and therapeutic targets. In this review, multi-omics data harmonization methods were introduced, and common approaches to molecular network inference were summarized. Our Prediction Logic Boolean Implication Networks (PLBINs) have advantages over other methods in constructing genome-scale multi-omics networks in bulk tumors and single cells in terms of computational efficiency, scalability, and accuracy. Based on the constructed multi-modal regulatory networks, graph theory network centrality metrics can be used in the prioritization of candidates for discovering biomarkers and therapeutic targets. Our approach to integrating multi-omics profiles in a patient cohort with large-scale patient EMRs such as the SEER-Medicare cancer registry combined with extensive external validation can identify potential biomarkers applicable in large patient populations. These methodologies form a conceptually innovative framework to analyze various available information from research laboratories and healthcare systems, accelerating the discovery of biomarkers and therapeutic targets to ultimately improve cancer patient survival outcomes.

The Evolution of tRNA Copy Number and Repertoire in Cellular Life

tRNAs are universal decoders that bridge the gap between transcriptome and proteome. They can also be processed into small RNA fragments with regulatory functions. In this work, we show that tRNA copy number is largely controlled by genome size in all cellular organisms, in contrast to what is observed for protein-coding genes that stop expanding between ~20,000 and ~35,000 loci per haploid genome in eukaryotes, regardless of genome size. Our analyses indicate that after the bacteria/archaea split, the tRNA gene pool experienced the evolution of increased anticodon diversity in the archaeal lineage, along with a tRNA gene size increase and mature tRNA size decrease. The evolution and diversification of eukaryotes from archaeal ancestors involved further expansion of the tRNA anticodon repertoire, additional increase in tRNA gene size and decrease in mature tRNA length, along with an explosion of the tRNA gene copy number that emerged coupled with accelerated genome size expansion. Our findings support the notion that macroscopic eukaryotes with a high diversity of cell types, such as land plants and vertebrates, independently evolved a high diversity of tRNA anticodons along with high gene redundancy caused by the expansion of the tRNA copy number. The results presented here suggest that the evolution of tRNA genes played important roles in the early split between bacteria and archaea, and in eukaryogenesis and the later emergence of complex eukaryotes, with potential implications in protein translation and gene regulation through tRNA-derived RNA fragments.

Analysis of Copy Number Variation in the Whole Genome of Normal-Haired and Long-Haired Tianzhu White Yaks

Long-haired individuals in the Tianzhu white yak population are a unique genetic resource, and have important landscape value. Copy number variation (CNV) is an important source of phenotypic variation in mammals. In this study, we used resequencing technology to detect the whole genome of 10 long-haired Tianzhu white yaks (LTWY) and 10 normal-haired Tianzhu white yaks (NTWY), and analyzed the differences of CNV in the genome of LTWYs and NTWYs. A total of 110268 CNVs were identified, 2006 CNVRs were defined, and the distribution map of these CNVRs on chromosomes was constructed. The comparison of LTWYs and NTWYs identified 80 differential CNVR-harbored genes, which were enriched in lipid metabolism, cell migration and other functions. Notably, some differential genes were identified as associated with hair growth and hair-follicle development (e.g., ASTN2, ATM, COL22A1, GK5, SLIT3, PM20D1, and SGCZ). In general, we present the first genome-wide analysis of CNV in LTWYs and NTWYs. Our results can provide new insights into the phenotypic variation of different hair lengths in Tianzhu white yaks.

Immunosuppressive role of SPP1-CD44 in the tumor microenvironment of intrahepatic cholangiocarcinoma assessed by single-cell RNA sequencing

PURPOSE

To demonstrate the biological function of Secreted Phosphoprotein 1(SPP1) and its immune suppressive role in the progression intrahepatic cholangiocarcinoma (ICC).

METHODS

We collected 62,770 cells' published transcriptome data of nine patients whose paired adjacent liver and tumor tissues were both available. We applied differential gene expression analysis to screen potential ICC marker genes, survival analysis to verify the prognostic value of SPP1, and correlation analysis to decipher factors that are related to SPP1 expression. The CellChat was used to distinguish interactions between cancer and T cells. CytoSig was applied to query cytokines that modulate CD44. Further, we established a proliferation score and correlated the score with inhibitory signals to determine the proliferation-suppressive function of SPP1-CD44.

RESULTS

SPP1 expression is significantly upregulated in tumoral epitheliums, and patients with higher SPP1 expression have worse survival (P < 0.05). Tumor cells communicate with T cells via SPP1-CD44 interactions. The average expression of SPP1 in malignant cells (SPP1m) and CD44 in T cells (CD44t) is moderately negatively correlated with T cell proliferation score. Immunosuppressive cytokine TGFβ-3 identified as an inducer of CD44 and was significantly negatively correlated with proliferation score (R = - 0.88, P < 0.01), and the negative correlation was aggravated in samples with high CD44 expression.

CONCLUSION

SPP1 is a prognostic marker of ICC and is associated with the genome heterogeneity. SPP1-CD44 hinders sustained proliferation of T cells, but immunosuppressive T cells in the tumor microenvironment may evade this inhibition by reducing CD44 expression.

Exome sequencing as a first-tier test for copy number variant detection: retrospective evaluation and prospective screening in 2418 cases

BACKGROUND

Despite the availability of whole exome (WES) and genome sequencing (WGS), chromosomal microarray (CMA) remains the first-line diagnostic test in most rare disorders diagnostic workup, looking for copy number variations (CNVs), with a diagnostic yield of 10%-20%. The question of the equivalence of CMA and WES in CNV calling is an organisational and economic question, especially when ordering a WGS after a negative CMA and/or WES.

METHODS

This study measures the equivalence between CMA and GATK4 exome sequencing depth of coverage method in detecting coding CNVs on a retrospective cohort of 615 unrelated individuals. A prospective detection of WES-CNV on a cohort of 2418 unrelated individuals, including the 615 individuals from the validation cohort, was performed.

RESULTS

On the retrospective validation cohort, every CNV detectable by the method (ie, a CNV with at least one exon not in a dark zone) was accurately called (64/64 events). In the prospective cohort, 32 diagnoses were performed among the 2418 individuals with CNVs ranging from 704 bp to aneuploidy. An incidental finding was reported. The overall increase in diagnostic yield was of 1.7%, varying from 1.2% in individuals with multiple congenital anomalies to 1.9% in individuals with chronic kidney failure.

CONCLUSION

Combining single-nucleotide variant (SNV) and CNV detection increases the suitability of exome sequencing as a first-tier diagnostic test for suspected rare Mendelian disorders. Before considering the prescription of a WGS after a negative WES, a careful reanalysis with updated CNV calling and SNV annotation should be considered.

Recent advances and current challenges in population genomics of structural variation in animals and plants

The field of population genomics has seen a surge of studies on genomic structural variation over the past two decades. These studies witnessed that structural variation is taxonomically ubiquitous and represent a dominant form of genetic variation within species. Recent advances in technology, especially the development of long-read sequencing platforms, have enabled the discovery of structural variants (SVs) in previously inaccessible genomic regions which unlocked additional structural variation for population studies and revealed that more SVs contribute to evolution than previously perceived. An increasing number of studies suggest that SVs of all types and sizes may have a large effect on phenotype and consequently major impact on rapid adaptation, population divergence, and speciation. However, the functional effect of the vast majority of SVs is unknown and the field generally lacks evidence on the phenotypic consequences of most SVs that are suggested to have adaptive potential. Non-human genomes are heavily under-represented in population-scale studies of SVs. We argue that more research on other species is needed to objectively estimate the contribution of SVs to evolution. We discuss technical challenges associated with SV detection and outline the most recent advances towards more representative reference genomes, which opens a new era in population-scale studies of structural variation.

Complex Genomic Rearrangements following Selection in a Glutamine-Limited Medium over Hundreds of Generations

Large-scale genomic changes, including copy number variations (CNVs), are frequently observed in long-term evolution experiments (LTEEs). We have previously reported the detection of recurrent CNVs in Saccharomyces cerevisiae populations adapting to glutamine-limited conditions over hundreds of generations. Here, we present the whole-genome sequencing (WGS) assemblies of 7 LTEE strains and their ancestor.

A comprehensive pan-cancer analysis of necroptosis molecules in four gynecologic cancers

Background

In recent years, it has been proved that necroptosis plays an important role in the occurrence, development, invasion, metastasis and drug resistance of malignant tumors. Hence, further evaluation and targeting of necroptosis may be of clinical benefit for gynecologic cancers (GCs).

Methods

To compare consistency and difference, we explored the expression pattern and prognostic value of necroptosis-related genes (NRGs) in pan-GC analysis through Linear regression and Empirical Bayesian, Univariate Cox analysis, and public databases from TCGA and Genotype-Tissue Expression (GTEx), including CESC, OV, UCEC, and UCS. We explored the copy number variation (CNV), methylation level and enrichment pathways of NRGs in the four GCs. Based on LASSO Cox regression analysis or principal component analysis, we established the prognostic NRG-signature or necroptosis-score for the four GCs. In addition, we predicted and compared functional pathways, tumor mutational burden (TMB), somatic mutation features, immunity status, immunotherapy, chemotherapeutic drug sensitivity of the NRG-signature based on NRGs. We also examined the expression level of several NRGs in OV samples that we collected using Quantitative Real-time PCR.

Results

We confirmed the presence of NRGs in expression, prognosis, CNV, and methylation for four GCs, thus comparing the consistency and difference among the four GCs. The prognosis and independent prognostic value of the risk signatures based on NRGs were determined. Through the results of subclass mapping, we found that GC patients with lower risk score may be more sensitive to PDL1 response and more sensitive to immune checkpoint blockade therapy. Drug susceptibility analysis showed that, 51, 45, 64, and 29 drugs with differences between risk groups were yielded in CESC, OV, UCEC, and UCS respectively. For OV, the expression differences of several NRGs in the tissues we collected were similar to that in TCGA.

Conclusion

Our comprehensive analysis of NRGs and NRG-signature demonstrated their similarity and difference, as well as their potential roles in prognosis and could guide therapeutic strategies, thus improving the outcome of GC patients.

Advances of multiplex ligation-dependent probe amplification technology in molecular diagnostics

Multiplex ligation-dependent probe amplification (MLPA) is a multiplex copy number analysis tool which is routinely used to detect large mutations in genetic diseases. With continuous modifications, MLPA has been extended for the detection of DNA methylation variation, single nucleotide polymorphisms, chromosome abnormalities and other forms of genomic variation. The combination with other techniques has even enlarged the application of MLPA in molecular diagnostics of various human diseases. In this review, the principle of MLPA-based techniques as well as their main and latest applications in clinical detection are described. It is believed that with improved automation, increased multiplexing, lower cost and the combination with other technologies, MLPA will play an increasingly important role in molecular diagnosis of human disease.

Systematic Analysis of Copy Number Variations in the Pathogenic Yeast Candida parapsilosis Identifies a Gene Amplification in RTA3 That is Associated with Drug Resistance

We analyzed the genomes of 170 C. parapsilosis isolates and identified multiple copy number variations (CNVs). We identified two genes, RTA3 (CPAR2_104610) and ARR3 (CPAR2_601050), each of which was the target of multiple independent amplification events. Phylogenetic analysis shows that most of these amplifications originated only once. For ARR3, which encodes a putative arsenate transporter, 8 distinct CNVs were identified, ranging in size from 2.3 kb to 10.5 kb with 3 to 23 copies. For RTA3, 16 distinct CNVs were identified, ranging in size from 0.3 kb to 4.5 kb with 2 to ~50 copies. One unusual amplification resulted in a DUP-TRP/INV-DUP structure similar to some human CNVs. RTA3 encodes a putative phosphatidylcholine (PC) floppase which is known to regulate the inward translocation of PC in Candida albicans. We found that an increased copy number of RTA3 correlated with resistance to miltefosine, an alkylphosphocholine drug that affects PC metabolism. Additionally, we conducted an adaptive laboratory evolution experiment in which two C. parapsilosis isolates were cultured in increasing concentrations of miltefosine. Two genes, CPAR2_303950 and CPAR2_102700, coding for putative PC flippases homologous to S. cerevisiae DNF1 gained homozygous protein-disrupting mutations in the evolved strains. Overall, our results show that C. parapsilosis can gain resistance to miltefosine, a drug that has recently been granted orphan drug designation approval by the United States Food and Drug Administration for the treatment of invasive candidiasis, through both CNVs or loss-of-function alleles in one of the flippase genes. IMPORTANCE Copy number variations (CNVs) are an important source of genomic diversity that have been associated with drug resistance. We identify two unusual CNVs in the human fungal pathogen Candida parapsilosis. Both target a single gene (RTA3 or ARR3), and they have occurred multiple times in multiple isolates. The copy number of RTA3, a putative floppase that controls the inward translocation of lipids in the cell membrane, correlates with resistance to miltefosine, a derivative of phosphatidylcholine (PC) that was originally developed as an anticancer drug. In 2021, miltefosine was designated an orphan drug by the United States Food and Drug Administration for the treatment of invasive candidiasis. Importantly, we find that resistance to miltefosine is also caused by mutations in flippases, which control the outward movement of lipids, and that many C. parapsilosis isolates are prone to easily acquiring an increased resistance to miltefosine.

Identification of pyroptosis-related lncRNA signature and AC005253.1 as a pyroptosis-related oncogene in prostate cancer

Background

Pyroptosis and prostate cancer (PCa) are closely related. The role of pyroptosis-related long non-coding RNAs (lncRNAs) (PRLs) in PCa remains elusive. This study aimed to explore the relationship between PRL and PCa prognosis.

Methods

Gene expression and clinical signatures were obtained from The Cancer Genome Atlas and Gene Expression Omnibus databases. A PRL risk prediction model was established by survival random forest analysis and least absolute shrinkage and selection operator regression. Functional enrichment, immune status, immune checkpoints, genetic mutations, and drug susceptibility analyses related to risk scores were performed by the single-sample gene set enrichment analysis, gene set variation analysis, and copy number variation analysis. PRL expression was verified in PCa cells. Cell Counting Kit-8, 5-ethynyl-2′-deoxyuridine, wound healing, transwell, and Western blotting assay were used to detect the proliferation, migration, invasion, and pyroptosis of PCa cells, respectively.

Results

Prognostic features based on six PRL (AC129507.1, AC005253.1, AC127502.2, AC068580.3, LIMD1-AS1, and LINC01852) were constructed, and patients in the high-score group had a worse prognosis than those in the low-score group. This feature was determined to be independent by Cox regression analysis, and the area under the curve of the 1-, 3-, and 5-year receiver operating characteristic curves in the testing cohort was 1, 0.93, and 0.92, respectively. Moreover, the external cohort validation confirmed the robustness of the PRL risk prediction model. There was a clear distinction between the immune status of the two groups. The expression of multiple immune checkpoints was also reduced in the high-score group. Gene mutation proportion in the high-score group increased, and the sensitivity to drugs increased significantly. Six PRLs were upregulated in PCa cells. Silencing of AC005253.1 inhibited cell proliferation, migration, and invasion in DU145 and PC-3 cells. Moreover, silencing of AC005253.1 promoted pyroptosis and inflammasome AIM2 expression.

Conclusions

Overall, we constructed a prognostic model of PCa with six PRLs and identified their expression in PCa cells. The experimental verification showed that AC005253.1 could affect the proliferation, migration, and invasion abilities of PCa cells. Meanwhile, AC005253.1 may play an important role in PCa by affecting pyroptosis through the AIM2 inflammasome. This result requires further research for verification.

A Novel hepatocellular carcinoma specific hypoxic related signature for predicting prognosis and therapeutic responses

Hypoxia is an important feature of the tumor microenvironment(TME) and is closely associated with cancer metastasis, immune evasion, and drug resistance. However, the precise role of hypoxia in hepatocellular carcinoma(HCC), as well as its influence on the TME, and drug sensitivity remains unclear. We found the excellent survival prediction value of Hypoxia_DEGs_Score model. In hypoxic HCC, somatic mutation, copy number variation, and DNA methylation were closely related to hypoxic changes and affected tumorigenesis, progression, metastasis, and drug resistance. In HCC, aggravated hypoxic stress was found to be accompanied by an immune exclusion phenotype and increased infiltration of immunosuppressive cells. In the validation cohort, patients with high Hypoxia_DEGs_Score were found to have worse immunotherapeutic outcomes and prognoses, and may benefit from drugs against cell cycle signaling pathways rather than those inhibiting the PI3K/mTOR pathway. Hypoxia_DEGs_Score has an excellent predictive capability of changes in the TME, the efficacy of immunotherapy, and the response of drugs. Therefore, Hypoxia_DEGs_Score can help develop personalized immunotherapy regimens and improve the prognosis of HCC patients.

Convergent copy number increase of genes associated with freshwater colonization in fishes

Copy number variation (CNV) can cause phenotypic changes. However, in contrast to amino acid substitutions and cis-regulatory changes, little is known about the functional categories of genes in which CNV is important for adaptation to novel environments. It is also unclear whether the same genes repeatedly change the copy numbers for adapting to similar environments. Here, we investigate CNV associated with freshwater colonization in fishes, which was observed multiple times across different lineages. Using 48 ray-finned fishes across diverse orders, we identified 23 genes whose copy number increases were associated with freshwater colonization. These genes showed enrichment for peptide receptor activity, hexosyltransferase activity and unsaturated fatty acid metabolism. We further revealed that three of the genes showed copy number increases in freshwater populations compared to marine ancestral populations of the stickleback genus Gasterosteus. These results indicate that copy number increases of genes involved in fatty acid metabolism (FADS2), immune function (PSMB8a) and thyroid hormone metabolism (UGT2) may be important for freshwater colonization at both the inter-order macroevolutionary scale and at the intra-genus microevolutionary scale. Further analysis across diverse taxa will help to understand the role of CNV in the adaptation to novel environments. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Chromosomal Copy Number Variation Analysis in Pregnancy Products from Recurrent and Sporadic Miscarriage Using Next-Generation Sequencing

Chromosomal abnormality is one of the causes of fetal miscarriage. The potential differences of fetal chromosomal abnormalities in sporadic miscarriage (SM) and recurrent miscarriage (RM) remain unclear. The purpose of this study was to investigate copy number variations (CNVs) in SM and RM to provide useful genetic guidance for pregnancy and prenatal diagnosis. Four hundred eight samples of aborted fetuses were analyzed by CNV sequencing, and further functional enrichment analysis was performed. Chromosomal abnormalities were identified in 218 (53.4%) fetuses. There were 62 cases (15.2%) with structural chromosomal abnormalities, including 41 with VUS CNVs, 8 with pathogenic CNVs (pCNVs), and 5 with likely pCNVs. Duplications or deletions of 7p22, 8p22, 8p23, and Xp22.31 were significantly more common in RM cases and therefore believed to be related to RM. A total of 289 genes were identified, and 29 different functions were enriched as potential RM candidate genes and functions, which were mainly concentrated in 4 functional categories: chemokines and chemotaxis, protease activity and protein modification, defense response to bacterial and fungal infections, and immune response. The results of this study may improve our understanding of the etiology of RM and contribute to the establishment of a population-based genetic marker information for RM.

Application of Droplet Digital PCR Technology in Muscular Dystrophies Research

Although they are considered rare disorders, muscular dystrophies have a strong impact on people’s health. Increased disease severity with age, frequently accompanied by the loss of ability to walk in some people, and the lack of treatment, have directed the researchers towards the development of more effective therapeutic strategies aimed to improve the quality of life and life expectancy, slow down the progression, and delay the onset or convert a severe phenotype into a milder one. Improved understanding of the complex pathology of these diseases together with the tremendous advances in molecular biology technologies has led to personalized therapeutic procedures. Different approaches that are currently under extensive investigation require more efficient, sensitive, and less invasive methods. Due to its remarkable analytical sensitivity, droplet digital PCR has become a promising tool for accurate measurement of biomarkers that monitor disease progression and quantification of various therapeutic efficiency and can be considered a tool for non-invasive prenatal diagnosis and newborn screening. Here, we summarize the recent applications of droplet digital PCR in muscular dystrophy research and discuss the factors that should be considered to get the best performance with this technology.

Selection of lncRNAs That Influence the Prognosis of Osteosarcoma Based on Copy Number Variation Data

Osteosarcoma is the most common primary malignancy in the musculoskeletal system. It is reported that copy number variation- (CNV-) derived lncRNAs contribute to the progression of osteosarcoma. However, whether CNV-derived lncRNAs affect the prognosis of osteosarcoma remains unclear. Here, we obtained osteosarcoma-related CNV data and gene expression profiles from The Cancer Genome Atlas (TCGA) database. CNV landscape analysis indicated that copy number amplification of lncRNAs was more frequent than deletion in osteosarcoma samples. Thirty-four CNV-lncRNAs with DNA-CNV frequencies greater than 30% and their corresponding 294 mRNAs were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway enrichment analyses revealed that these mRNAs were mainly enriched in olfaction, olfactory receptor activity, and olfactory transduction processes. Furthermore, we predicted that a total of 23 genes were cis-regulated by 16 CNV-lncRNAs, while 30 transcription factors (TFs) were trans-regulated by 5 CNV-lncRNAs. Through $t$ -tests, univariate Cox regression analysis, and the least absolute shrinkage and selection operator (LASSO), we constructed a CNV-related risk model including 3 lncRNAs (AC129492.1, PSMB1, and AC037459.4). The Kaplan-Meier (K-M) curves indicated that patients with high-risk scores showed poor prognoses. The areas under the receiver operating characteristic (ROC) curves (AUC) for predicting 3-, 5-, and 7-year overall survival (OS) were greater than 0.7, showing a satisfactory predictive efficiency. Gene set enrichment analysis (GSEA) revealed that the prognostic signature was intimately linked to skeletal system development, immune regulation, and inflammatory response. Collectively, our study developed a novel 3-CNV-lncRNA prognostic signature that would provide theoretical guidance for the clinical prognostic management of osteosarcoma.

Genome-wide analysis of copy number variants and normal facial variation in a large cohort of Bantu Africans

Similarity in facial characteristics between relatives suggests a strong genetic component underlies facial variation. While there have been numerous studies of the genetics of facial abnormalities and, more recently, single nucleotide polymorphism (SNP) genome-wide association studies (GWASs) of normal facial variation, little is known about the role of genetic structural variation in determining facial shape. In a sample of Bantu African children, we found that only 9% of common copy number variants (CNVs) and 10-kb CNV analysis windows are well tagged by SNPs (r² ≥ 0.8), indicating that associations with our internally called CNVs were not captured by previous SNP-based GWASs. Here, we present a GWAS and gene set analysis of the relationship between normal facial variation and CNVs in a sample of Bantu African children. We report the top five regions, which had p values ≤ 9.35 × 10⁻⁶ and find nominal evidence of independent CNV association (p < 0.05) in three regions previously identified in SNP-based GWASs. The CNV region with strongest association (p = 1.16 × 10⁻⁶, 55 losses and seven gains) contains NFATC1, which has been linked to facial morphogenesis and Cherubism, a syndrome involving abnormal lower facial development. Genomic loss in the region is associated with smaller average lower facial depth. Importantly, new loci identified here were not identified in a SNP-based GWAS, suggesting that CNVs are likely involved in determining facial shape variation. Given the plethora of SNP-based GWASs, calling CNVs from existing data may be a relatively inexpensive way to aid in the study of complex traits.

Establishment of a Dihydrofolate Reductase Gene Knock-In Zebrafish Strain to Aid Preliminary Analysis of Congenital Heart Disease Mechanisms

Background: The dihydrofolate reductase (DHFR) gene is imperative in development, therefore it is essential to explore its effects on heart development. Thus, here a dhfr zebrafish knock-in (KI) strain was constructed.

Methods: CRISPR/Cas9 technology was used to establish the dhfr KI zebrafish strain. This strain was hybridized with TgG fluorescent strain zebrafish to observe the phenotypes of heart shape, size, and circularization direction. Wild-type (WT) and KI zebrafish were then dissected and histologically stained to observe pathological changes. Western blot analysis was used to verify the increased expressions of zebrafish genes after KI. Hybridization experiments were used to confirm the presence of abnormal gonadal dysplasia.

Results: The zebrafish dhfr KI strain was successfully constructed through CRISPR/Cas9 technology. At 6 days post fertilization (dpf), microscopic examinations of KI (homozygous) specimens revealed pericardial effusions, heart compressions, and curled tails. Compared with WT, the Hematoxylin and Eosin (H&E) tissue sections of KI-homozygous zebrafish showed defects such as reduced atria and ventricles. Western blot analysis indicated that the expression of the DHFR protein increased in both heterozygotes and homozygotes of dhfr KI zebrafish. Hybridization experiments revealed that dhfr KI may affect gonadal function.

Conclusion: The DHFR gene plays an important regulatory role in the process of heart development, and copy number variations (CNVs) of this gene may constitute a new pathogenic mechanism of congenital heart disease (CHD).

Effects of copy number variations on brain structure and risk for psychiatric illness: Large-scale studies from the ENIGMA working groups on CNVs

The Enhancing NeuroImaging Genetics through Meta-Analysis copy number variant (ENIGMA-CNV) and 22q11.2 Deletion Syndrome Working Groups (22q-ENIGMA WGs) were created to gain insight into the involvement of genetic factors in human brain development and related cognitive, psychiatric and behavioral manifestations. To that end, the ENIGMA-CNV WG has collated CNV and magnetic resonance imaging (MRI) data from ~49,000 individuals across 38 global research sites, yielding one of the largest studies to date on the effects of CNVs on brain structures in the general population. The 22q-ENIGMA WG includes 12 international research centers that assessed over 533 individuals with a confirmed 22q11.2 deletion syndrome, 40 with 22q11.2 duplications, and 333 typically developing controls, creating the largest-ever 22q11.2 CNV neuroimaging data set. In this review, we outline the ENIGMA infrastructure and procedures for multi-site analysis of CNVs and MRI data. So far, ENIGMA has identified effects of the 22q11.2, 16p11.2 distal, 15q11.2, and 1q21.1 distal CNVs on subcortical and cortical brain structures. Each CNV is associated with differences in cognitive, neurodevelopmental and neuropsychiatric traits, with characteristic patterns of brain structural abnormalities. Evidence of gene-dosage effects on distinct brain regions also emerged, providing further insight into genotype-phenotype relationships. Taken together, these results offer a more comprehensive picture of molecular mechanisms involved in typical and atypical brain development. This "genotype-first" approach also contributes to our understanding of the etiopathogenesis of brain disorders. Finally, we outline future directions to better understand effects of CNVs on brain structure and behavior.

Adaptation dynamics between copy-number and point mutations

Together, copy-number and point mutations form the basis for most evolutionary novelty, through the process of gene duplication and divergence. While a plethora of genomic data reveals the long-term fate of diverging coding sequences and their cis-regulatory elements, little is known about the early dynamics around the duplication event itself. In microorganisms, selection for increased gene expression often drives the expansion of gene copy-number mutations, which serves as a crude adaptation, prior to divergence through refining point mutations. Using a simple synthetic genetic reporter system that can distinguish between copy-number and point mutations, we study their early and transient adaptive dynamics in real time in Escherichia coli. We find two qualitatively different routes of adaptation, depending on the level of functional improvement needed. In conditions of high gene expression demand, the two mutation types occur as a combination. However, under low gene expression demand, copy-number and point mutations are mutually exclusive; here, owing to their higher frequency, adaptation is dominated by copy-number mutations, in a process we term amplification hindrance. Ultimately, due to high reversal rates and pleiotropic cost, copy-number mutations may not only serve as a crude and transient adaptation, but also constrain sequence divergence over evolutionary time scales.

Radiomic Analysis of Contrast-Enhanced MRI Predicts DNA Copy-Number Subtype and Outcome in Lower-Grade Gliomas

RATIONALE AND OBJECTIVES

DNA copy-number (CN)2-subtype impairs outcomes in patients with lower-grade gliomas (LGG). We aimed to determine the value of preoperative nomograms integrating radiomic and radiographic (RR) features in predicting DNA copy-number subtype.

METHODS

Data of 153 consecutive patients were retrospectively analyzed. A total of 1167 radiomics features were extracted from contrast-enhanced MR images. LASSO logistic regression was performed to choose the key features and construct a radiomics signature. Three CN-related RR model were built with multivariate logistic regression.

RESULTS

CN2-subtype was associated with shortest median PFS(p <0.001) and OS (p <0.001). The radiomics nomogram, which incorporated the signature (AUC:0.891, OR: 2.345; p = 0.001), extranodular growth (OR: 14.413; p <0.001) and width (OR: 0.194; p = 0.027), distinguished CN2-subtype with an AUC of 0.924(95%CI: 0.869-0.979).The radiomics nomogram, which incorporated the signature (AUC:0.730, OR: 2.408; p = 0.001), hemorrhage (OR: 0.100; p <0.001), poorly-defined margin (OR:4.433; p = 0.001) and volume>=60cm3 (OR: 4.195; p = 0.002) were associated with CN1-subtype (AUC:0.829,95%CI:0.765-0.892).The radiomics nomogram, which incorporated the signature (AUC:0.660, OR: 2.518; p = 0.003), necrosis/cystic(OR:6.975; p = 0.008), hemorrhage (OR:3.723; p = 0.024), poorly-defined margin (OR:0.124; p <0.001) and frontal lobe tumors (OR: 4.870; p <0.001) were associated with CN3-subtype (AUC: 0.837,95%CI: 0.767-0.909).All three RR models showed good discrimination and calibration. Decision curve analysis indicated that all RR models were clinically useful. The average accuracy of the ten-fold cross validation was 92.8% for CN2-subtype, 72.6% for CN1-subtype and 79.0% for CN3-subtype.

CONCLUSION

The shortest PFS and OS was observed in LGG patients with CN2-subtype. The RR models, integrating radiomic and radiographic features, demonstrates good performance for predicting DNA copy-number subtype and clinical outcomes.

Obesity in Children with Leptin Receptor Gene Polymorphisms

INTRODUCTION

The study of single nucleotide polymorphisms (SNPs) of the leptin receptor gene (LEPR) based on next generation genomic sequencing (NGS) data is becoming an increasingly important aspect of diagnosis, treatment and prevention of both metabolically healthy (MHO) and metabolically unhealthy obesity (MUO) phenotypes.

MATERIAL AND METHODS

35 obese children 6-18 years old were examined by the NGS method with bioinformatic analysis. The main group (n = 18) was formed by children with MUO, according to the recommendations of the expert group of the National Heart, Lung, and Blood Institute. The control group (n = 17) was represented by children with MHO. Statistical methods were used: analysis of variance, Wald's sequential analysis, Spearman's correlation analysis, analysis of nominal data and multiple discriminant analysis.

RESULTS

10 types of non-synonymous SNPs (rs3790435, rs1137100, rs2186248, rs70940803, rs79639154, rs1359482195, rs1137101, rs1805094, rs13306520, rs13306522) of the LEPR gene in obese children have been identified. Multiple discriminant analysis demonstrated that the following LEPR SNPs are of greatest importance in the development of MUO: rs3790435, rs13306522, rs13306520. Analysis of nominal data revealed significant differences in the groups for Copy number variation (CNV) rs3790435 of the LEPR gene. Wald's analysis allowed us to identify 6 important predictors of MUO (І ≥ 0.5): 2 CNV rs3790435 (Relative Risk, RR = 2, Prognostic coefficient, PC = +2.76); male gender of the child (RR = 1.3, PC = +1.35); rs3790435 (RR = 1.9, PC = +2.76); hyperleptinemia more than 40.56 ng/ml (RR = 2, PC = +3); CNV rs1359482195 ≥ 3 (RR = 1.9, PC = +5.8); SNP of the LEPR gene ≥4 (RR = 3.8, PC = +5.8).

CONCLUSION

Children with the genotype rs3790435 gene LEPR had signs of MUO more often.

Analysis of Genomic Copy Number Variation in Miscarriages During Early and Middle Pregnancy

The purpose of this study was to explore the copy number variations (CNVs) associated with miscarriage during early and middle pregnancy and provide useful genetic guidance for pregnancy and prenatal diagnosis. A total of 505 fetal specimens were collected and CNV sequencing (CNV-seq) analysis was performed to determine the types and clinical significance of CNVs, and relevant medical records were collected. The chromosomal abnormality rate was 54.3% (274/505), among which the numerical chromosomal abnormality rate was 40.0% (202/505) and structural chromosomal abnormality rate was 14.3% (72/505). Chromosomal monosomy mainly occurred on sex chromosomes, and chromosomal trisomy mainly occurred on chromosomes 16, 22, 21, 15, 13, and 9. The incidence of numerical chromosomal abnormalities in ≥35 year-old age pregnant women was significantly higher than <35 year-old age group. The highest incidence of pathogenic CNV (pCNV) was found in fetuses at ≤6 weeks of pregnancy (5.26%), and the incidence of variants of unknown significance (VOUS) CNVs decreased gradually with the increase of gestational age. The rate of chromosomal abnormalities of fetuses in early pregnancy (59.5%) was higher than that of fetuses in middle pregnancy (27.2%) (p < 0.001). There were 168 genes in VOUS + pCNV regions. 41 functions and 12 pathways (p < 0.05) were enriched of these genes by Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Some meaningful genetic etiology information such as genes and pathways has been obtained, it may provide useful genetic guidance for pregnancy and prenatal diagnosis.

A Comprehensive Multiomics Analysis Identified Ubiquilin 4 as a Promising Prognostic Biomarker of Immune-Related Therapy in Pan-Cancer

Recently, it was reported that ubiquilin 4 (UBQLN4) alteration was associated with genomic instability in some cancers. However, whether UBQLN4 is a valuable biomarker for the prognosis of immunotherapy in pan-cancer was not identified. We evaluated the biologic and oncologic significance of UBQLN4 in pan-cancer at multiomics level, such as expression, mutation, copy number variation (CNV), methylation, and N⁶-methyladenosine (m⁶A) methylation. These omics data were obtained from several public databases, including Oncomine, The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), the Genotype-Tissue Expression (GTEx), the Human Protein Atlas (HPA), Gene Set Cancer Analysis (GSCA), m⁶A-Atlas, CancerSEA, and RNAactDrug. We found that UBQLN4 mRNA and protein were overexpressed in most cancer types, and the expression, mutation, CNV, and methylation of UBQLN4 were associated with the prognosis of some cancers. Mechanistically, UBQLN4 was involved in angiogenesis, DNA damage, apoptosis, and the pathway of PI3K/AKT and TSC/mTOR. Moreover, UBQLN4 mRNA was significantly correlated with immune checkpoints, tumor mutational burden (TMB), microsatellite instability (MSI), and mismatch repair (MMR). And, the correlation among UBQLN4 mRNA, CNV, and methylation and immune microenvironment was also identified. Furthermore, UBQLN4 was associated with the sensitivity of chemotherapy and targeted drugs at multiomics level. In conclusion, UBQLN4 was a promising prognostic biomarker of immune-related therapy in pan-cancer.

Bridging Tumorigenesis and Therapy Resistance With a Non-Darwinian and Non-Lamarckian Mechanism of Adaptive Evolution

Although neo-Darwinian (and less often Lamarckian) dynamics are regularly invoked to interpret cancer's multifarious molecular profiles, they shine little light on how tumorigenesis unfolds and often fail to fully capture the frequency and breadth of resistance mechanisms. This uncertainty frames one of the most problematic gaps between science and practice in modern times. Here, we offer a theory of adaptive cancer evolution, which builds on a molecular mechanism that lies outside neo-Darwinian and Lamarckian schemes. This mechanism coherently integrates non-genetic and genetic changes, ecological and evolutionary time scales, and shifts the spotlight away from positive selection towards purifying selection, genetic drift, and the creative-disruptive power of environmental change. The surprisingly simple use-it or lose-it rationale of the proposed theory can help predict molecular dynamics during tumorigenesis. It also provides simple rules of thumb that should help improve therapeutic approaches in cancer.

Copy number variation characterization and possible candidate genes in miscarriage and stillbirth by next-generation sequencing analysis

Background

The present study aimed to explore the etiological relationship between miscarriage and stillbirth and copy number variations (CNVs), as well as provide useful genetic guidance for high‐risk pregnancy.

Methods

In total, 659 fetal samples were recruited and subjected to DNA extraction and CNV sequencing (CNV‐seq), relevant medical records were collected.

Results

There were 322 cases (48.86%) with chromosomal abnormalities, including 230 with numerical abnormalities and 92 with structural abnormalities. Chromosomal monosomy variations mainly occurred on sex chromosomes and trisomy variations mainly occurred on chromosomes 16, 22, 21, 18, 13 and 15. In total, 41 pathogenic CNVs (23 microdeletions and 18 microduplications) were detected in 27 fetal tissues. The rates of numerical chromosomal abnormalities were 29.30% (109/372), 32.39% (57/176) and 57.66% (64/111) in < 30‐year‐old, 30–34‐year‐old and ≥ 35‐year‐old age pregnant women, respectively, and increased with an increasing age (p < 0.001). There was statistically significant difference (χ² = 7.595, p = 0.022) in the rates of structural chromosomal abnormalities in these groups (13.71%, 18.75% and 7.21%, respectively). The rates of numerical chromosomal abnormalities were 45.44% (219/482), 7.80% (11/141) and 0% (0/36) in the ≤ 13 gestational weeks, 14–27 weeks and ≥ 28 weeks groups, respectively, and decreased with respect to the increasing gestational age of the fetuses (p < 0.001).

Conclusions

The present study has obtained useful and accurate genetic etiology information that will provide useful genetic guidance for high‐risk pregnancies.

HNRNPH1-stabilized LINC00662 promotes ovarian cancer progression by activating the GRP78/p38 pathway

Numerous studies suggest an important role for copy number alterations (CNAs) in cancer progression. However, CNAs of long intergenic noncoding RNAs (lincRNAs) in ovarian cancer (OC) and their potential functions have not been fully investigated. Here, based on analysis of The Cancer Genome Atlas (TCGA) database, we identified in this study an oncogenic lincRNA termed LINC00662 that exhibited a significant correlation between its CNA and its increased expression. LINC00662 overexpression is highly associated with malignant features in OC patients and is a prognostic indicator. LINC00662 significantly promotes OC cell proliferation and metastasis in vitro and in vivo. Mechanistically, LINC00662 is stabilized by heterogeneous nuclear ribonucleoprotein H1 (HNRNPH1). Moreover, LINC00662 exerts oncogenic effects by interacting with glucose-regulated protein 78 (GRP78) and preventing its ubiquitination in OC cells, leading to activation of the oncogenic p38 MAPK signaling pathway. Taken together, our results define an oncogenic role for LINC00662 in OC progression mediated via GRP78/p38 signaling, with potential implications regarding therapeutic targets for OC.

Copy number alteration profiling facilitates differential diagnosis between ossifying fibroma and fibrous dysplasia of the jaws

Ossifying fibroma (OF) and fibrous dysplasia (FD) are two fibro-osseous lesions with overlapping clinicopathological features, making diagnosis challenging. In this study, we applied a whole-genome shallow sequencing approach to facilitate differential diagnosis via precise profiling of copy number alterations (CNAs) using minute amounts of DNA extracted from morphologically correlated microdissected tissue samples. Freshly frozen tissue specimens from OF (n = 29) and FD (n = 28) patients were obtained for analysis. Lesion fibrous tissues and surrounding normal tissues were obtained by laser capture microdissection (LCM), with ~30–50 cells (5 000–10 000 µm²) per sample. We found that the rate of recurrent CNAs in OF cases was much higher (44.8%, 13 of 29) than that in FD cases (3.6%, 1 of 28). Sixty-nine percent (9 of 13) of the CNA-containing OF cases involved segmental amplifications and deletions on Chrs 7 and 12. We also identified eight CNA-associated genes (HILPDA, CALD1, C1GALT1, MICALL2, PHF14, AIMP2, MDM2, and CDK4) with amplified expression, which was consistent with the copy number changes. We further confirmed a jaw lesion with a previous uncertain diagnosis due to its ambiguous morphological features and the absence of GNAS mutation as OF based on the typical Chr 12 amplification pattern in its CNA profile. Moreover, analysis of a set of longitudinal samples collected from an individual with a cellular lesion in suspicion of OF at the first surgery, recurrence and the latest malignant transformation revealed identical CNA patterns at the three time points, suggesting that copy number profiling can be used as an important tool to identify borderline lesions or lesions with malignant potential. Overall, CNA profiling of fibro-osseous lesions can greatly improve differential diagnosis between OF and FD and help predict disease progression.

A Novel Computational Framework to Predict Disease-Related Copy Number Variations by Integrating Multiple Data Sources

Copy number variation (CNV) may contribute to the development of complex diseases. However, due to the complex mechanism of path association and the lack of sufficient samples, understanding the relationship between CNV and cancer remains a major challenge. The unprecedented abundance of CNV, gene, and disease label data provides us with an opportunity to design a new machine learning framework to predict potential disease-related CNVs. In this paper, we developed a novel machine learning approach, namely, IHI-BMLLR (Integrating Heterogeneous Information sources with Biweight Mid-correlation and L1-regularized Logistic Regression under stability selection), to predict the CNV-disease path associations by using a data set containing CNV, disease state labels, and gene data. CNVs, genes, and diseases are connected through edges and then constitute a biological association network. To construct a biological network, we first used a self-adaptive biweight mid-correlation (BM) formula to calculate correlation coefficients between CNVs and genes. Then, we used logistic regression with L1 penalty (LLR) function to detect genes related to disease. We added stability selection strategy, which can effectively reduce false positives, when using self-adaptive BM and LLR. Finally, a weighted path search algorithm was applied to find top D path associations and important CNVs. The experimental results on both simulation and prostate cancer data show that IHI-BMLLR is significantly better than two state-of-the-art CNV detection methods (i.e., CCRET and DPtest) under false-positive control. Furthermore, we applied IHI-BMLLR to prostate cancer data and found significant path associations. Three new cancer-related genes were discovered in the paths, and these genes need to be verified by biological research in the future.

Characterization of TERT and BRAF copy number variation in papillary thyroid carcinoma: An analysis of the cancer genome atlas study

Alterations in the genome, including mutations and copy number variation (CNV), can drive cancer progression. The Cancer Genome Atlas (TCGA) project studying papillary thyroid cancer (PTC) identified a number of recurrent arm-level copy number amplifications, some spanning genes that are also commonly mutated in thyroid cancer. Herein, we focus on the role of TERT and BRAF CNV in PTC, including its relation to mutation status, gene expression, and clinicopathological characteristics. Utilizing TCGA CNV data, we identified focal amplifications and deletions involving the TERT and BRAF loci. TERT amplifications are more frequent in later stage thyroid tumors; in contrast, BRAF amplifications are not associated with stage. Furthermore, TERT amplifications are more frequently found in tumors also harboring TERT mutations, the combination further increasing TERT expression. Conversely, BRAF amplifications are more frequently found in BRAF wildtype tumors, and are more common in the follicular subtype of PTC as well as classic PTCs associated with a high follicular component and a RAS-like expression profile (assessed by the BRAF/RAS score). This is the first study to examine the TCGA thyroid dataset for gene-level CNV of TERT and BRAF, and their relationship with mutation status, tumor type and tumor stage. Assessing the differences in patterns of TERT and BRAF amplifications in the context of the mutation status of these genes may provide insight into the differing roles CNV can play depending on tumor type, and may lead to a better understanding of cancer drivers in thyroid cancer.

Genome-wide copy number variations in a large cohort of bantu African children

BACKGROUND

Copy number variations (CNVs) account for a substantial proportion of inter-individual genomic variation. However, a majority of genomic variation studies have focused on single-nucleotide variations (SNVs), with limited genome-wide analysis of CNVs in large cohorts, especially in populations that are under-represented in genetic studies including people of African descent.

METHODS

We carried out a genome-wide copy number analysis in > 3400 healthy Bantu Africans from Tanzania. Signal intensity data from high density (> 2.5 million probes) genotyping arrays were used for CNV calling with three algorithms including PennCNV, DNAcopy and VanillaICE. Stringent quality metrics and filtering criteria were applied to obtain high confidence CNVs.

RESULTS

We identified over 400,000 CNVs larger than 1 kilobase (kb), for an average of 120 CNVs (SE = 2.57) per individual. We detected 866 large CNVs (≥ 300 kb), some of which overlapped genomic regions previously associated with multiple congenital anomaly syndromes, including Prader-Willi/Angelman syndrome (Type1) and 22q11.2 deletion syndrome. Furthermore, several of the common CNVs seen in our cohort (≥ 5%) overlap genes previously associated with developmental disorders.

CONCLUSIONS

These findings may help refine the phenotypic outcomes and penetrance of variations affecting genes and genomic regions previously implicated in diseases. Our study provides one of the largest datasets of CNVs from individuals of African ancestry, enabling improved clinical evaluation and disease association of CNVs observed in research and clinical studies in African populations.

DNA damage and growth hormone hypersecretion in pituitary somatotroph adenomas

Drivers of sporadic benign pituitary adenoma growth are largely unknown. Whole-exome sequencing of 159 prospectively resected pituitary adenomas showed that somatic copy number alteration (SCNA) rather than mutation is a hallmark of hormone-secreting adenomas and that SCNAs correlate with adenoma phenotype. Using single-gene SCNA pathway analysis, we observed that both cAMP and Fanconi anemia DNA damage repair pathways were affected by SCNAs in growth hormone-secreting (GH-secreting) somatotroph adenomas. As somatotroph differentiation and GH secretion are dependent on cAMP activation and we previously showed DNA damage, aneuploidy, and senescence in somatotroph adenomas, we studied links between cAMP signaling and DNA damage. Stimulation of cAMP in C57BL/6 mouse primary pituitary cultures using forskolin or a long-acting GH-releasing hormone (GHRH) analog increased GH production and DNA damage measured by H2AX phosphorylation and a comet assay. Octreotide, a somatostatin receptor ligand that targets somatotroph adenoma GH secretion in patients with acromegaly, inhibited cAMP and GH and reversed DNA damage induction. In vivo long-acting GHRH treatment also induced pituitary DNA damage in mice. We conclude that cAMP, which induces somatotroph proliferation and GH secretion, may concomitantly induce DNA damage, potentially linking hormone hypersecretion to SCNA and genome instability. These results elucidating somatotroph adenoma pathophysiology identify pathways for targeted treatment.