The impact of human copy number variation on gene expression

Genomic structural variation contributes to evolved changes in gene expression in high-altitude Tibetan sheep

Tibetan sheep were introduced to the Qinghai Tibet plateau roughly 3,000 B.P., making this species a good model for investigating genetic mechanisms of high-altitude adaptation over a relatively short timescale. Here, we characterize genomic structural variants (SVs) that distinguish Tibetan sheep from closely related, low-altitude Hu sheep, and we examine associated changes in tissue-specific gene expression. We document differentiation between the two sheep breeds in frequencies of SVs associated with genes involved in cardiac function and circulation. In Tibetan sheep, we identified high-frequency SVs in a total of 462 genes, including EPAS1, PAPSS2, and PTPRD. Single-cell RNA-Seq data and luciferase reporter assays revealed that the SVs had cis-acting effects on the expression levels of these three genes in specific tissues and cell types. In Tibetan sheep, we identified a high-frequency chromosomal inversion that exhibited modified chromatin architectures relative to the noninverted allele that predominates in Hu sheep. The inversion harbors several genes with altered expression patterns related to heart protection, brown adipocyte proliferation, angiogenesis, and DNA repair. These findings indicate that SVs represent an important source of genetic variation in gene expression and may have contributed to high-altitude adaptation in Tibetan sheep.

Understanding genetic variability: exploring large-scale copy number variants through non-invasive prenatal testing in European populations

Large-scale copy number variants (CNVs) are structural alterations in the genome that involve the duplication or deletion of DNA segments, contributing to genetic diversity and playing a crucial role in the evolution and development of various diseases and disorders, as they can lead to the dosage imbalance of one or more genes. Massively parallel sequencing (MPS) has revolutionized the field of genetic analysis and contributed significantly to routine clinical diagnosis and screening. It offers a precise method for detecting CNVs with exceptional accuracy. In this context, a non-invasive prenatal test (NIPT) based on the sequencing of cell-free DNA (cfDNA) from pregnant women's plasma using a low-coverage whole genome MPS (WGS) approach represents a valuable source for population studies. Here, we analyzed genomic data of 12,732 pregnant women from the Slovak (9,230), Czech (1,583), and Hungarian (1,919) populations. We identified 5,062 CNVs ranging from 200 kbp and described their basic characteristics and differences between the subject populations. Our results suggest that re-analysis of sequencing data from routine WGS assays has the potential to obtain large-scale CNV population frequencies, which are not well known and may provide valuable information to support the classification and interpretation of this type of genetic variation. Furthermore, this could contribute to expanding knowledge about the central European genome without investing in additional laboratory work, as NIPTs are a relatively widely used screening method.

Investigating structural variant, indel and single nucleotide polymorphism differentiation between locally adapted Atlantic salmon populations

Genomic structural variants (SVs) are now recognized as an integral component of intraspecific polymorphism and are known to contribute to evolutionary processes in various organisms. However, they are inherently difficult to detect and genotype from readily available short-read sequencing data, and therefore remain poorly documented in wild populations. Salmonid species displaying strong interpopulation variability in both life history traits and habitat characteristics, such as Atlantic salmon (Salmo salar), offer a prime context for studying adaptive polymorphism, but the contribution of SVs to fine-scale local adaptation has yet to be explored. Here, we performed a comparative analysis of SVs, single nucleotide polymorphisms (SNPs) and small indels (<50 bp) segregating in the Romaine and Puyjalon salmon, two putatively locally adapted populations inhabiting neighboring rivers (Québec, Canada) and showing pronounced variation in life history traits, namely growth, fecundity, and age at maturity and smoltification. We first catalogued polymorphism using a hybrid SV characterization approach pairing both short- (16X) and long-read sequencing (20X) for variant discovery with graph-based genotyping of SVs across 60 salmon genomes, along with characterization of SNPs and small indels from short reads. We thus identified 115,907 SVs, 8,777,832 SNPs and 1,089,321 short indels, with SVs covering 4.8 times more base pairs than SNPs. All three variant types revealed a highly congruent population structure and similar patterns of F ST and density variation along the genome. Finally, we performed outlier detection and redundancy analysis (RDA) to identify variants of interest in the putative local adaptation of Romaine and Puyjalon salmon. Genes located near these variants were enriched for biological processes related to nervous system function, suggesting that observed variation in traits such as age at smoltification could arise from differences in neural development. This study therefore demonstrates the feasibility of large-scale SV characterization and highlights its relevance for salmonid population genomics.

Identify Candidate Genes Associated with the Weight and Egg Quality Traits in Wenshui Green Shell-Laying Chickens by the Copy Number Variation-Based Genome-Wide Association Study

Copy number variation (CNV), as an essential source of genetic variation, can have an impact on gene expression, genetic diversity, disease susceptibility, and species evolution in animals. To better understand the weight and egg quality traits of chickens, this paper aimed to detect CNVs in Wenshui green shell-laying chickens and conduct a copy number variation regions (CNVRs)-based genome-wide association study (GWAS) to identify variants and candidate genes associated with their weight and egg quality traits to support related breeding efforts. In our paper, we identified 11,035 CNVRs in Wenshui green shell-laying chickens, which collectively spanned a length of 13.1 Mb, representing approximately 1.4% of its autosomal genome. Out of these CNVRs, there were 10,446 loss types, 491 gain types, and 98 mixed types. Notably, two CNVRs showed significant correlations with egg quality, while four CNVRs exhibited significant associations with body weight. These significant CNVRs are located on chromosome 4. Further analysis identified potential candidate genes that influence weight and egg quality traits, including FAM184B, MED28, LAP3, ATOH8, ST3GAL5, LDB2, and SORCS2. In this paper, the CNV map of the Wenshui green shell-laying chicken genome was constructed for the first time through population genotyping. Additionally, CNVRs can be employed as molecular markers to genetically improve chickens' weight and egg quality traits.

Long-insert sequence capture detects high copy numbers in a defence-related beta-glucosidase gene βglu-1 with large variations in white spruce but not Norway spruce

Conifers are long-lived and slow-evolving, thus requiring effective defences against their fast-evolving insect natural enemies. The copy number variation (CNV) of two key acetophenone biosynthesis genes Ugt5/Ugt5b and βglu-1 may provide a plausible mechanism underlying the constitutively variable defence in white spruce (Picea glauca) against its primary defoliator, spruce budworm. This study develops a long-insert sequence capture probe set (Picea_hung_p1.0) for quantifying copy number of βglu-1-like, Ugt5-like genes and single-copy genes on 38 Norway spruce (Picea abies) and 40 P. glauca individuals from eight and nine provenances across Europe and North America respectively. We developed local assemblies (Piabi_c1.0 and Pigla_c.1.0), full-length transcriptomes (PIAB_v1 and PIGL_v1), and gene models to characterise the diversity of βglu-1 and Ugt5 genes. We observed very large copy numbers of βglu-1, with up to 381 copies in a single P. glauca individual. We observed among-provenance CNV of βglu-1 in P. glauca but not P. abies. Ugt5b was predominantly single-copy in both species. This study generates critical hypotheses for testing the emergence and mechanism of extreme CNV, the dosage effect on phenotype, and the varying copy number of genes with the same pathway. We demonstrate new approaches to overcome experimental challenges in genomic research in conifer defences.

Palmitoylation landscapes across human cancers reveal a role of palmitoylation in tumorigenesis

BACKGROUND

Protein palmitoylation, which is catalyzed by palmitoyl-transferase and de-palmitoyl-transferase, plays a crucial role in various biological processes. However, the landscape and dynamics of protein palmitoylation in human cancers are not well understood.

METHODS

We utilized 23 palmitoyl-acyltransferases and seven de-palmitoyl-acyltransferases as palmitoylation-related genes for protein palmitoylation analysis. Multiple publicly available datasets were employed to conduct pan-cancer analysis, examining the transcriptome, genomic alterations, clinical outcomes, and correlation with c-Myc (Myc) for palmitoylation-related genes. Real-time quantitative PCR and immunoblotting were performed to assess the expression of palmitoylation-related genes and global protein palmitoylation levels in cancer cells treated with Myc depletion or small molecule inhibitors. Protein docking and drug sensitivity analyses were employed to predict small molecules that target palmitoylation-related genes.

RESULTS

We identified associations between palmitoylation and cancer subtype, stage, and patient survival. We discovered that abnormal DNA methylation and oncogenic Myc-driven transcriptional regulation synergistically contribute to the dysregulation of palmitoylation-related genes. This dysregulation of palmitoylation was closely correlated with immune infiltration in the tumor microenvironment and the response to immunotherapy. Importantly, dysregulated palmitoylation was found to modulate canonical cancer-related pathways, thus influencing tumorigenesis. To support our findings, we performed a proof-of-concept experiment showing that depletion of Myc led to reduced expression of most palmitoylation-related genes, resulting in decreased global protein palmitoylation levels. Through mass spectrometry and enrichment analyses, we also identified palmitoyl-acyltransferases ZDHHC7 and ZDHHC23 as significant contributors to mTOR signaling, DNA repair, and immune pathways, highlighting their potential roles in tumorigenesis. Additionally, our study explored the potential of three small molecular (BI-2531, etoposide, and piperlongumine) to modulate palmitoylation by targeting the expression or activity of palmitoylation-related genes or enzymes.

CONCLUSIONS

Overall, our findings underscore the critical role of dysregulated palmitoylation in tumorigenesis and the response to immunotherapy, mediated through classical cancer-related pathways and immune cell infiltration. Additionally, we propose that the aforementioned three small molecule hold promise as potential therapeutics for modulating palmitoylation, thereby offering novel avenues for cancer therapy.

Targeted large fragment deletion in plants using paired crRNAs with type I CRISPR system

The CRISPR-Cas systems have been widely used as genome editing tools, with type II and V systems typically introducing small indels, and type I system mediating long-range deletions. However, the precision of type I systems for large fragment deletion is still remained to be optimized. Here, we developed a compact Cascade-Cas3 Dvu I-C system with Cas11c for plant genome editing. The Dvu I-C system was efficient to introduce controllable large fragment deletion up to at least 20 kb using paired crRNAs. The paired-crRNAs design also improved the controllability of deletions for the type I-E system. Dvu I-C system was sensitive to spacer length and mismatch, which was benefit for target specificity. In addition, we showed that the Dvu I-C system was efficient for generating stable transgenic lines in maize and rice with the editing efficiency up to 86.67%. Overall, Dvu I-C system we developed here is powerful for achieving controllable large fragment deletions.

Copy Number Variation Regions Differing in Segregation Patterns Span Different Sets of Genes

Copy number variations regions (CNVRs) can be classified either as segregating, when found in both parents, and offspring, or non-segregating. A total of 65 segregating and 31 non-segregating CNVRs identified in at least 10 individuals within a dense pedigree of the Gochu Asturcelta pig breed was subjected to enrichment and functional annotation analyses to ascertain their functional independence and importance. Enrichment analyses allowed us to annotate 1018 and 351 candidate genes within the bounds of the segregating and non-segregating CNVRs, respectively. The information retrieved suggested that the candidate genes spanned by segregating and non-segregating CNVRs were functionally independent. Functional annotation analyses allowed us to identify nine different significantly enriched functional annotation clusters (ACs) in segregating CNVR candidate genes mainly involved in immunity and regulation of the cell cycle. Up to five significantly enriched ACs, mainly involved in reproduction and meat quality, were identified in non-segregating CNVRs. The current analysis fits with previous reports suggesting that segregating CNVRs would explain performance at the population level, whereas non-segregating CNVRs could explain between-individuals differences in performance.

Association of PKLR gene copy number, expression levels and enzyme activity with 2,3,7,8-TCDD exposure in individuals exposed to Agent Orange/Dioxin in Vietnam

2,3,7,8-tetrachlorodibenzo-p-dioxin (2,3,7,8-TCDD) is the most toxic congener of dioxin and has serious long-term effects on the environment and human health. Pyruvate Kinase L/R (PKLR) gene expression levels and gene variants are associated with pyruvate kinase enzyme deficiency, which has been identified as the cause of several diseases linked to dioxin exposure. In this study, we estimated PKLR gene copy number and gene expression levels using real-time quantitative PCR (RT-qPCR) assays, genotyped PKLR SNP rs3020781 by Sanger sequencing, and quantified plasma pyruvate kinase enzyme activity in 100 individuals exposed to Agent Orange/Dioxin near Bien Hoa and Da Nang airfields in Vietnam and 100 healthy controls. The means of PKLR copy numbers and PKLR gene expression levels were significantly higher, while pyruvate kinase enzyme activity was significantly decreased in Agent Orange/Dioxin-exposed individuals compared to healthy controls (P < 0.0001). Positive correlations of PKLR gene copy number and gene expression with 2,3,7,8-TCDD concentrations were observed (r = 0.2, P = 0.045 and r = 0.54, P < 0.0001, respectively). In contrast, pyruvate kinase enzyme activity was inversely correlated with 2,3,7,8-TCDD concentrations (r = -0.52, P < 0.0001). PKLR gene copy number and gene expression levels were also inversely correlated with pyruvate kinase enzyme activity. Additionally, PKLR SNP rs3020781 was found to be associated with 2,3,7,8-TCDD concentrations and PKLR gene expression. In conclusion, PKLR copy number, gene expression levels, and pyruvate kinase enzyme activity are associated with 2,3,7,8-TCDD exposure in individuals living in Agent Orange/Dioxin-contaminated areas.

A comprehensive analysis of copy number variations in diverse apple populations

BACKGROUND

As an important source of genetic variation, copy number variation (CNV) can alter the dosage of DNA segments, which in turn may affect gene expression level and phenotype. However, our knowledge of CNV in apple is still limited. Here, we obtained high-confidence CNVs and investigated their functional impact based on genome resequencing data of two apple populations, cultivars and wild relatives.

RESULTS

In this study, we identified 914,610 CNVs comprising 14,839 CNV regions (CNVRs) from 346 apple accessions, including 289 cultivars and 57 wild relatives. CNVRs summed to 71.19 Mb, accounting for 10.03% of the apple genome. Under the low linkage disequilibrium (LD) with nearby SNPs, they could also accurately reflect the population structure of apple independent of SNPs. Furthermore, A total of 3,621 genes were covered by CNVRs and functionally involved in biological processes such as defense response, reproduction and metabolic processes. In addition, the population differentiation index ([Formula: see text]) analysis between cultivars and wild relatives revealed 127 CN-differentiated genes, which may contribute to trait differences in these two populations.

CONCLUSIONS

This study was based on identification of CNVs from 346 diverse apple accessions, which to our knowledge was the largest dataset for CNV analysis in apple. Our work presented the first comprehensive CNV map and provided valuable resources for understanding genomic variations in apple.

A crowdsourcing database for the copy-number variation of the Spanish population

BACKGROUND

Despite being a very common type of genetic variation, the distribution of copy-number variations (CNVs) in the population is still poorly understood. The knowledge of the genetic variability, especially at the level of the local population, is a critical factor for distinguishing pathogenic from non-pathogenic variation in the discovery of new disease variants.

RESULTS

Here, we present the SPAnish Copy Number Alterations Collaborative Server (SPACNACS), which currently contains copy number variation profiles obtained from more than 400 genomes and exomes of unrelated Spanish individuals. By means of a collaborative crowdsourcing effort whole genome and whole exome sequencing data, produced by local genomic projects and for other purposes, is continuously collected. Once checked both, the Spanish ancestry and the lack of kinship with other individuals in the SPACNACS, the CNVs are inferred for these sequences and they are used to populate the database. A web interface allows querying the database with different filters that include ICD10 upper categories. This allows discarding samples from the disease under study and obtaining pseudo-control CNV profiles from the local population. We also show here additional studies on the local impact of CNVs in some phenotypes and on pharmacogenomic variants. SPACNACS can be accessed at: http://csvs.clinbioinfosspa.es/spacnacs/ .

CONCLUSION

SPACNACS facilitates disease gene discovery by providing detailed information of the local variability of the population and exemplifies how to reuse genomic data produced for other purposes to build a local reference database.

Pyroptosis relates to tumor microenvironment remodeling and prognosis: A pan-cancer perspective

Background and aim

Pyroptosis is an inflammatory form of programmed cell death implicated in inflammation and disease. Moreover, inducing pyroptosis has been appreciated as anti-cancer therapy for its ability to unleash anti-cancer immune responses.

Methods

Utilizing the data available in The Cancer Genome Atlas (TCGA), pyroptosis-related genes' (PRGs) expression, genomic aberrations, and clinical significance were systematically analyzed in pan-cancer. A GSVA score was obtained to rate pyroptosis level and divide the cancers into pyroptosis-low and pyroptosis-high groups. Immunohistochemistry (IHC) was used to evaluate the differential expression of major PRGs (GSDMC, GSDMD, GSDME, NLRP3, NLRC4, IL1B) in selected tumor types (COAD, HNSC, KIRC, LIHC, LUAD, LUSC). Selection of tumors for immunohistochemistry (IHC) was based on their expression pattern in TCGA cancers, clinical relevance, tumor epidemiology, and sample availability.

Results

Differential expression of PRGs was evident in various cancers and associated with prognosis which was driven by genomic variations and epigenetic abnormalities, such as single nucleotide variations (SNVs), copy number variation (CNV) and DNA methylation level. For example, methylation of PRGs in lower grade glioma (LGG), uveal melanoma (UVM) and kidney renal clear cell carcinoma (KIRC) were predictive of improved survival as upregulation of PRGs was risky in these cancers. Pyroptosis level significantly differentiated tumor from normal samples in 15 types of cancers, exhibited a progressive trend with cancer stage, observed variation among cancer subtypes, and showed a significant association with cancer prognosis. Higher pyroptosis level was associated with worst prognosis in majority of the cancers in terms of OS (KIRC, LGG, and UVM), PFS (GBM, KIRC, LGG, PRAD, THCA, and THYM) and DSS (KIRC and LGG) as estimated by Kaplan-Meier survival curves. Moreover, Pyroptosis level was strongly indicative of a hot tumor immune microenvironment with high presence of CD8+ T cell and other T cell subtypes. Several oncogenic pathways, such as P53 pathway, DNA repair, KRAS signaling, epithelial-mesenchymal transition (EMT), IL6 JAK STAT3 signaling, IL2 STAT5 signaling, PI3K AKT MTOR signaling and angiogenesis, were enriched in pyroptosis-hi subgroups across cancers.

Conclusions

Genetic alterations in PRGs greatly influence the pyroptosis level and cancer prognosis. A relatively hot tumor immune microenvironment was associated with pyroptosis irrespective of the cancer prognosis. Overall, our study reveals the critical role of pyroptosis in cancer and highlights pyroptosis-based therapeutic vulnerabilities.

Coexisting Molecular Alterations Increase the Risk of Malignancy in Thyroid Nodules with Copy Number Alterations

Molecular mutations and alterations play a role in thyroid tumorigenesis. Different alterations are associated with different clinical and pathological characteristics. Copy number alterations (CNAs) are known to be present in some thyroid tumors; however, their idiosyncratic clinicopathological implications are not yet well elucidated. A retrospective chart review was performed to identify patients with CNAs on pre-operative molecular testing results who subsequently underwent surgical treatment between January 2016 and April 2022 at McGill University teaching hospitals. Of the 316 patients with thyroid nodules who opted for molecular testing with ThyroSeqV3 followed by surgery, 67 (21.2%) nodules were positive for CNAs, including 23 Bethesda III, 31 Bethesda IV, 12 Bethesda V and 1 Bethesda VI nodules. On surgical pathology, 29.9% were benign and 70.1% were malignant or non-invasive follicular thyroid neoplasm with papillary-like nuclear features (NIFTP). Among those that were malignant/NIFTP, 17.02% were considered to be aggressive cancers. The presence of other molecular alterations was found to be an independent predictor of malignancy in multivariate analysis (OR = 5.087, 95% C.I. = 1.12-23.04, p = 0.035). No unique factor was correlated with aggressiveness; however, CNA-positive thyroid nodules that were associated with high-risk mutations such as BRAF V600E, TP53, NTRK1/3 fusion, or PTEN mutation with high allele frequency (AF) ended up being aggressive cancers. Most of the CNA-positive thyroid nodules resulted in follicular patterned tumors in 41 (65.2%) cases and oncocytic tumors in 20 (29.9%) cases. This study demonstrates that 70.1% of surgically resected thyroid nodules with CNAs were malignant/NIFTP. Most CNA-positive thyroid nodules were either oncocytic patterned tumors or follicular patterned tumors. Furthermore, CNA-positive thyroid nodules were more likely to be malignant if they were associated with other molecular alterations or mutations.

Copy Number Variation of Circulating Tumor DNA (ctDNA) Detected Using NIPT in Neoadjuvant Chemotherapy-Treated Ovarian Cancer Patients

Analysis of circulating tumor DNA (ctDNA) can be used to characterize and monitor cancers. Recently, non-invasive prenatal testing (NIPT) as a new next-generation sequencing (NGS)-based approach has been applied for detecting ctDNA. This study aimed to investigate the copy number variations (CNVs) utilizing the non-invasive prenatal testing in plasma ctDNA from ovarian cancer (OC) patients who were treated with neoadjuvant chemotherapy (NAC). The plasma samples of six patients, including stages II–IV, were collected during the pre- and post-NAC treatment that were divided into NAC-sensitive and NAC-resistant groups during the follow-up time. CNV analysis was performed using the NIPT via two methods “an open-source algorithm WISECONDORX and NextGENe software.” Results of these methods were compared in pre- and post-NAC of OC patients. Finally, bioinformatics tools were used for data mining from The Cancer Genome Atlas (TCGA) to investigate CNVs in OC patients. WISECONDORX analysis indicated fewer CNV changes on chromosomes before treatment in the NAC-sensitive rather than NAC-resistant patients. NextGENe data indicated that CNVs are not only observed in the coding genes but also in non-coding genes. CNVs in six genes were identified, including HSF1, TMEM249, MROH1, GSTT2B, ABR, and NOMO2, only in NAC-resistant patients. The comparison of these six genes in NAC-resistant patients with The Cancer Genome Atlas data illustrated that the total alteration frequency is amplification, and the highest incidence of the CNVs (≥35% based on TCGA data) is found in MROH1, TMEM249, and HSF1 genes on the chromosome (Chr) 8. Based on TCGA data, survival analysis showed a significant reduction in the overall survival among chemotherapy-resistant patients as well as a high expression level of these three genes compared to that of sensitive samples (all, p < 0.0001). The continued Chr8 study using WISECONDORX revealed CNV modifications in NAC-resistant patients prior to NAC therapy, but no CNV changes were observed in NAC-sensitive individuals. Our findings showed that low coverage whole-genome sequencing analysis used for NIPT could identify CNVs in ctDNA of OC patients before and after chemotherapy. These CNVs are different in NAC-sensitive and -resistant patients highlighting the potential application of this approach in cancer patient management.

Response Best-subset Selector for Multivariate Regression with High-dimensional Response Variables

This article investigates the statistical problem of response-variable selection with high-dimensional response variables and a diverging number of predictor variables with respect to the sample size in the framework of multivariate linear regression. A response best-subset selection model is proposed by introducing a 0–1 selection indictor for each response variable, then a response best-subset selector is developed by introducing a separation parameter and a novel penalized least-squares function. The developed procedure can perform response-variable selection and regression-coefficient estimation simultaneously, and the proposed response best-subset selector has model consistency under mild conditions for both fixed and diverging numbers of predictor variables. Also, consistency and asymptotic normality of regression-coefficient estimators are presented for cases with a fixed dimension, and it is discovered that the Bonferroni test is a special response best-subset selector. Finite-sample simulations show that the response best-subset selector has strong advantages over existing competitors in terms of the Matthews correlation coefficient, a criterion aimed at balancing accuracies for both true and false response variables. An analysis of actual data demonstrates the effectiveness of the response best-subset selector in an application involving the identification of dosage-sensitive genes.

Overexpression of ERCC6L correlates with poor prognosis and confers malignant phenotypes of lung adenocarcinoma

Excision repair cross‑complementation group 6 like (ERCC6L) has been reported to be upregulated in a variety of malignant tumors and plays a critical oncogenic role. However, the role and molecular mechanism of ERCC6L in lung adenocarcinoma (LUAD) remain unclear, and were therefore investigated in the present study. Clinical data of patients with LUAD were obtained and bioinformatics analysis was performed to investigate the expression characteristics, prognostic value, and biological function of ERCC6L. In addition, cell function experiments were performed to detect the effect of ERCC6L silencing on the biological behavior of LUAD cells. The results revealed that ERCC6L expression was significantly higher in LUAD tissues vs. normal lung tissues and closely associated with nodal invasion, advanced clinical stage and survival in LUAD. Overexpression of ERCC6L was an independent prognostic biomarker of overall survival, progression‑free interval, and disease‑specific survival in patients with LUAD. DNA amplification and low methylation levels of ERCC6L suggested regulation at both the genetic and epigenetic levels. The most significant positive genes co‑expressed with ERCC6L were mainly enriched in the cell cycle signaling pathway. The major functions of ERCC6L in LUAD cells were positively correlated with the cell cycle, DNA damage, DNA repair, proliferation, invasion and epithelial‑mesenchymal transition (EMT). Knockdown of ERCC6L inhibited the proliferative, migratory and invasive abilities of A549 and PC9 cells. It also promoted cell apoptosis, and led to cell cycle arrest in the S phase. ERCC6L may regulate the EMT process through the Wnt/β‑catenin and Wnt/Notch 3 signaling pathways, thus regulating the tumorigenesis and progression of LUAD. The overexpression of ERCC6L may be a biological indicator for the diagnosis and prognosis of LUAD. ERCC6L may be a novel molecular target for the treatment of lung cancer.

Intra-specific copy number variation of MHC class II genes in the Siamese fighting fish

Duplicates of genes for major histocompatibility complex (MHC) molecules can be subjected to selection independently and vary markedly in their evolutionary rates, sequence polymorphism, and functional roles. Therefore, without a thorough understanding of their copy number variation (CNV) in the genome, the MHC-dependent fitness consequences within a species could be misinterpreted. Studying the intra-specific CNV of this highly polymorphic gene, however, has long been hindered by the difficulties in assigning alleles to loci and the lack of high-quality genomic data. Here, using the high-quality genome of the Siamese fighting fish (Betta splendens), a model for mate choice studies, and the whole-genome sequencing (WGS) data of 17 Betta species, we achieved locus-specific amplification of their three classical MHC class II genes - DAB1, DAB2, and DAB3. By performing quantitative PCR and depth-of-coverage analysis using the WGS data, we revealed intra-specific CNV at the DAB3 locus. We identified individuals that had two allelic copies (i.e., heterozygous or homozygous) or one allele (i.e., hemizygous) and individuals without this gene. The CNV was due to the deletion of a 20-kb-long genomic region harboring both the DAA3 and DAB3 genes. We further showed that the three DAB genes were under different modes of selection, which also applies to their corresponding DAA genes that share similar pattern of polymorphism. Our study demonstrates a combined approach to study CNV within a species, which is crucial for the understanding of multigene family evolution and the fitness consequences of CNV.

Copy Number Analyses Identified a Novel Gene: APOBEC3A Related to Lipid Metabolism in the Pathogenesis of Preeclampsia

Background

Preeclampsia is a heterogeneous and complex disease with its pathogenesis mechanism not fully elucidated. A certain subset of patients with preeclampsia exhibit disturbances in lipid metabolism before clinical symptoms. Moreover, there is a tendency for preeclampsia to run in families. Whether genetic factors play a role in abnormal lipid metabolism during the incidence of preeclampsia has not been well investigated.

Methods

Preeclampsia patients (n = 110) and healthy age- and gravidity-matched pregnant women (n = 110) were enrolled in this study. Peripheral blood specimens were used for genomic analysis (n = 10/group) or laboratory validation (n = 100/group). We retrospectively obtained the baseline clinical characteristics of 68 preeclampsia patients and 107 controls in early pregnancy (12–14 gestational weeks). Correlation analyses between differential genes and baseline lipid profiles were performed to identify candidate genes. In vitro and in vivo gain-of-function models were constructed with lentivirus and adeno-associated virus systems, respectively, to investigate the role of candidate genes in regulating lipid metabolism and the development of preeclampsia.

Results

We observed that preeclampsia patients exhibited significantly elevated plasma TC (P = 0.037) and TG (P < 0.001) levels and increased body mass index (P = 0.006) before the disease onset. Within the region of 27 differential copy number variations, six genes potentially connected with lipid metabolism were identified. The aberrant copies of APOBEC3A, APOBEC3A_B, BTNL3, and LMF1 between preeclampsia patients and controls were verified by quantitative polymerase chain reaction. Especially, APOBEC3A showed a significant positive correlation with TC (P < 0.001) and LDL (P = 0.048) in early pregnancy. Then, our in vitro data revealed that overexpression of APOBEC3A disrupted lipid metabolism in HepG2 cells and affected both cholesterol and fatty acid metabolisms. Finally, in vivo study in a hepatic-specific overexpressed APOBEC3A mouse model revealed abnormal parameters related to lipid metabolism. Pregnant mice of the same model at the end of pregnancy showed changes related to preeclampsia-like symptoms, such as increases in sFlt-1 levels and sFlt-1/PLGF ratios in the placenta and decreases in fetal weight.

Conclusion

Our findings established a new link between genetics and lipid metabolism in the pathogenesis of preeclampsia and could contribute to a better understanding of the molecular mechanisms of preeclampsia.

Modified eQTL and Somatic DNA Segment Alterations in Esophageal Squamous Cell Carcinoma for Genes Related to Immunity, DNA Repair, and Inflammation

We integrated ESCC expression and GWAS genotyping, to investigate eQTL and somatic DNA segment alterations, including somatic copy number alteration, allelic imbalance (AI), and loss of heterozygosity (LOH) in ESCC. First, in eQTL analysis, we used a classical approach based on genotype data from GWAS and expression signals in normal tissue samples, and then used a modified approach based on fold change in the tumor vs. normal samples. We focused on the genes in three pathways: inflammation, DNA repair, and immunity. Among the significant (p < 0.05) SNP-probe pairs from classical and modified eQTL analyses, 24 genes were shared by the two approaches, including 18 genes that showed the same numbers of SNPs and probes and 6 genes that had the different numbers of SNPs and probes. For these 18 genes, we found 28 SNP−probe pairs were correlated in opposite directions in the two approaches, indicating an intriguing difference between the classical and modified eQTL approaches. Second, we analyzed the somatic DNA segment alterations. Across the 24 genes, abnormal gene expression on mRNA arrays was seen in 19−95% of cases and 26−78% showed somatic DNA segment alterations on Affymetrix GeneChip Human Mapping Arrays. The results suggested that this strategy could identify gene expression and somatic DNA segment alterations for biological markers (genes) by combining classical and modified eQTLs and somatic DNA evaluation on SNP arrays. Thus, this study approach may allow us to understand functionality indicative of potentially relevant biomarkers in ESCC.

Gene Variants Involved in Nonsense-Mediated mRNA Decay Suggest a Role in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) is a heterogeneous neurodevelopmental condition with unclear etiology. Many genes have been associated with ASD risk, but the underlying mechanisms are still poorly understood. An important post-transcriptional regulatory mechanism that plays an essential role during neurodevelopment, the Nonsense-Mediated mRNA Decay (NMD) pathway, may contribute to ASD risk. In this study, we gathered a list of 46 NMD factors and regulators and investigated the role of genetic variants in these genes in ASD. By conducting a comprehensive search for Single Nucleotide Variants (SNVs) in NMD genes using Whole Exome Sequencing data from 1828 ASD patients, we identified 270 SNVs predicted to be damaging in 28.7% of the population. We also analyzed Copy Number Variants (CNVs) from two cohorts of ASD patients (N = 3570) and discovered 38 CNVs in 1% of cases. Importantly, we discovered 136 genetic variants (125 SNVs and 11 CNVs) in 258 ASD patients that were located within protein domains required for NMD. These gene variants are classified as damaging using in silico prediction tools, and therefore may interfere with proper NMD function in ASD. The discovery of NMD genes as candidates for ASD in large patient genomic datasets provides evidence supporting the involvement of the NMD pathway in ASD pathophysiology.

What genes are differentially expressed in individuals with schizophrenia? A systematic review

Schizophrenia is a severe, complex mental disorder characterized by a combination of positive symptoms, negative symptoms, and impaired cognitive function. Schizophrenia is highly heritable (~80%) with multifactorial etiology and complex polygenic genetic architecture. Despite the large number of genetic variants associated with schizophrenia, few causal variants have been established. Gaining insight into the mechanistic influences of these genetic variants may facilitate our ability to apply these findings to prevention and treatment. Though there have been more than 300 studies of gene expression in schizophrenia over the past 15 years, none of the studies have yielded consistent evidence for specific genes that contribute to schizophrenia risk. The aim of this work is to conduct a systematic review and synthesis of case-control studies of genome-wide gene expression in schizophrenia. Comprehensive literature searches were completed in PubMed, EmBase, and Web of Science, and after a systematic review of the studies, data were extracted from those that met the following inclusion criteria: human case-control studies comparing the genome-wide transcriptome of individuals diagnosed with schizophrenia to healthy controls published between January 1, 2000 and June 30, 2020 in the English language. Genes differentially expressed in cases were extracted from these studies, and overlapping genes were compared to previous research findings from the genome-wide association, structural variation, and tissue-expression studies. The transcriptome-wide analysis identified different genes than those previously reported in genome-wide association, exome sequencing, and structural variation studies of schizophrenia. Only one gene, GBP2, was replicated in five studies. Previous work has shown that this gene may play a role in immune function in the etiology of schizophrenia, which in turn could have implications for risk profiling, prevention, and treatment. This review highlights the methodological inconsistencies that impede valid meta-analyses and synthesis across studies. Standardization of the use of covariates, gene nomenclature, and methods for reporting results could enhance our understanding of the potential mechanisms through which genes exert their influence on the etiology of schizophrenia. Although these results are promising, collaborative efforts with harmonization of methodology will facilitate the identification of the role of genes underlying schizophrenia.

Genome-wide analysis of deletions in maize population reveals abundant genetic diversity and functional impact

Two read depth methods were jointly used in next-generation sequencing data to identify deletions in maize population. GWAS by deletions were analyzed for gene expression pattern and classical traits, respectively. Many studies have confirmed that structural variation (SV) is pervasive throughout the maize genome. Deletion is one type of SV that may impact gene expression and cause phenotypic changes in quantitative traits. In this study, two read count approaches were used to analyze the deletions in the whole-genome sequencing data of 270 maize inbred lines. A total of 19,754 deletion windows overlapped 12,751 genes, which were unevenly distributed across the genome. The deletions explained population structure well and correlated with genomic features. The deletion proportion of genes was determined to be negatively correlated with its expression. The detection of gene expression quantitative trait loci (eQTL) indicated that local eQTL were fewer but had larger effects than distant ones. The common associated genes were related to basic metabolic processes, whereas unique associated genes with eQTL played a role in the stress or stimulus responses in multiple tissues. Compared with the eQTL detected by SNPs derived from the same sequencing data, 89.4% of the associated genes could be detected by both markers. The effect of top eQTL detected by SNPs was usually larger than that detected by deletions for the same gene. A genome-wide association study (GWAS) on flowering time and plant height illustrated that only a few loci could be consistently captured by SNPs, suggesting that combining deletion and SNP for GWAS was an excellent strategy to dissect trait architecture. Our findings will provide insights into characteristic and biological function of genome-wide deletions in maize.

CNVs with adaptive potential in Rangifer tarandus: genome architecture and new annotated assembly

Rangifer tarandus has experienced recent drastic population size reductions throughout its circumpolar distribution and preserving the species implies genetic diversity conservation. To facilitate genomic studies of the species populations, we improved the genome assembly by combining long read and linked read and obtained a new highly accurate and contiguous genome assembly made of 13,994 scaffolds (L90 = 131 scaffolds). Using de novo transcriptome assembly of RNA-sequencing reads and similarity with annotated human gene sequences, 17,394 robust gene models were identified. As copy number variations (CNVs) likely play a role in adaptation, we additionally investigated these variations among 20 genomes representing three caribou ecotypes (migratory, boreal and mountain). A total of 1,698 large CNVs (length > 1 kb) showing a genome distribution including hotspots were identified. 43 large CNVs were particularly distinctive of the migratory and sedentary ecotypes and included genes annotated for functions likely related to the expected adaptations. This work includes the first publicly available annotation of the caribou genome and the first assembly allowing genome architecture analyses, including the likely adaptive CNVs reported here.

Integrated multi-omics analysis of the clinical relevance and potential regulatory mechanisms of splicing factors in hepatocellular carcinoma

Splicing factors (SFs) have been increasingly documented to perturb the genome of cancers. However, little is known about the alterations of SFs in hepatocellular carcinoma (HCC). This study comprehensively delineated the genomic and epigenomic characteristics of 404 SFs in HCC based on the multi-omics data from the Cancer Genome Atlas database. The analysis revealed several clinically relevant SFs that could be effective biomarkers for monitoring the onset and prognosis of HCC (such as, HSPB1, DDX39A, and NELFE, which were the three most significant clinically relevant SFs). Functional enrichment analysis of these indicators showed the enrichment of pathways related to splicing and mRNA processes. Furthermore, the study found that SF copy number variation is common in HCC and could be a typical characteristic of hepato-carcinogenesis; the complex expression regulation of SFs was significantly affected by copy number variant and methylation. Several SFs with significant mutation patterns were identified (such as, RNF213, SF3B1, SPEN, NOVA1, and EEF1A1), and the potential regulatory network of SFs was constructed to identify their potential mechanisms for regulating clinically relevant alternative splicing events. Therefore, this study established a foundation to uncover the broad molecular spectrum of SFs for future functional and therapeutic studies of HCC.

Protein-gene Expression Nexus: Comprehensive characterization of human cancer cell lines with proteogenomic analysis

Researchers have gained new therapeutic insights using multi-omics platform approaches to study DNA, RNA, and proteins of comprehensively characterized human cancer cell lines. To improve our understanding of the molecular features associated with oncogenic modulation in cancer, we proposed a proteogenomic database for human cancer cell lines, called Protein-gene Expression Nexus (PEN). We have expanded the characterization of cancer cell lines to include genetic, mRNA, and protein data of 145 cancer cell lines from various public studies. PEN contains proteomic and phosphoproteomic data on 4,129,728 peptides, 13,862 proteins, 7,138 phosphorylation site-associated genomic variations, 117 studies, and 12 cancer. We analyzed functional characterizations along with the integrated datasets, such as cis/trans association for copy number alteration (CNA), single amino acid variation for coding genes, post-translation modification site variation for Single Amino Acid Variation, and novel peptide expression for noncoding regions and fusion genes. PEN provides a user-friendly interface for searching, browsing, and downloading data and also supports the visualization of genome-wide association between CNA and expression, novel peptide landscape, mRNA-protein abundance, and functional annotation. Together, this dataset and PEN data portal provide a resource to accelerate cancer research using model cancer cell lines. PEN is freely accessible at http://combio.snu.ac.kr/pen.

Screening for Copy Number Variations of the 15q13.3 Hotspot in CHRNA7 Gene and Expression in Patients with Migraines

Background: a migraine is a neurological disease. Copy number variation (CNV) is a phenomenon in which parts of the genome are repeated. We investigated the effects of the CNV and gene expression at the location 15q13.3 in the Cholinergic Receptor Nicotinic Alpha 7 Subunit (CHRNA7) gene, which we believe to be effective in the migraine clinic. Methods: we evaluated changes in CHRNA7 gene expression levels and CNV of 15q13.3 in patients with migraine (n = 102, with aura, n = 43; without aura, n = 59) according to healthy controls (n = 120) by q-PCR. The data obtained were analyzed against the reference telomerase reverse transcriptase (TERT) gene with the double copy number by standard curve analysis. Copy numbers were graded as a normal copy (2), gain (2>), and loss (<2). Results: we analyzed using the 2^−ΔΔCT calculation method. The CHRNA7 gene was significantly downregulated in patients (p < 0.05). The analysis of CNV in the CHRNA7 gene was statistically significant in the patient group, according to healthy controls (p < 0.05). A decreased copy number indicates a dosage loss. However, no significant difference was observed among gain, normal, and loss copy numbers and expression values in patients (p > 0.05). The change in CNV was not associated with the downregulation of the CHRNA7 gene. Conclusion: Downregulation of the CHRNA7 gene may contribute to the formation of migraine by inactivation of the alpha-7 nicotinic receptor (α7nAChR). The association of CNV gains and losses with migraines will lead to better understanding of the molecular mechanisms and pathogenesis, to better define the disease, to be used as a treatment target.

Copy number variation and expression of exportin-4 associates with severity of fibrosis in metabolic associated fatty liver disease

Background

Liver fibrosis risk is a heritable trait, the outcome of which is the net deposition of extracellular matrix by hepatic stellate cell-derived myofibroblasts. Whereas nucleotide sequence variations have been extensively studied in liver fibrosis, the role of copy number variations (CNV) in which genes exist in abnormal numbers of copies (mostly due to duplication or deletion) has had limited exploration.

Methods

The impact of the XPO4 CNV on histological liver damage was examined in a cohort comprised 646 Caucasian patients with biopsy-proven MAFLD and 170 healthy controls. XPO4 expression was modulated and function was examined in human and animal models.

Findings

Here we demonstrate in a cohort of 816 subjects, 646 with biopsy-proven metabolic associated liver disease (MAFLD) and 170 controls, that duplication in the exportin 4 (XPO4) CNV is associated with the severity of liver fibrosis. Functionally, this occurs via reduced expression of hepatic XPO4 that maintains sustained activation of SMAD3/SMAD4 and promotes TGF-β1-mediated HSC activation and fibrosis. This effect was mediated through termination of nuclear SMAD3 signalling. XPO4 demonstrated preferential binding to SMAD3 compared to other SMADs and led to reduced SMAD3-mediated responses as shown by attenuation of TGFβ1 induced SMAD transcriptional activity, reductions in the recruitment of SMAD3 to target gene promoters following TGF-β1, as well as attenuation of SMAD3 phosphorylation and disturbed SMAD3/SMAD4 complex formation.

Interpretation

We conclude that a CNV in XPO4 is a critical mediator of fibrosis severity and can be exploited as a therapeutic target for liver fibrosis.

Funding

ME and JG are supported by the Robert W. Storr Bequest to the Sydney Medical Foundation, University of Sydney; a National Health and Medical Research Council of Australia (NHMRC) Program Grant (APP1053206) and Project and ideas grants (APP2001692, APP1107178 and APP1108422). AB is supported by an Australian Government Research Training Program (RTP) scholarship. EB is supported by Horizon 2020 under grant 634413 for the project EPoS.

Quantitative neurogenetics: applications in understanding disease

Neurodevelopmental and neurodegenerative disorders (NNDs) are a group of conditions with a broad range of core and co-morbidities, associated with dysfunction of the central nervous system. Improvements in high throughput sequencing have led to the detection of putative risk genetic loci for NNDs, however, quantitative neurogenetic approaches need to be further developed in order to establish causality and underlying molecular genetic mechanisms of pathogenesis. Here, we discuss an approach for prioritizing the contribution of genetic risk loci to complex-NND pathogenesis by estimating the possible impacts of these loci on gene regulation. Furthermore, we highlight the use of a tissue-specificity gene expression index and the application of artificial intelligence (AI) to improve the interpretation of the role of genetic risk elements in NND pathogenesis. Given that NND symptoms are associated with brain dysfunction, risk loci with direct, causative actions would comprise genes with essential functions in neural cells that are highly expressed in the brain. Indeed, NND risk genes implicated in brain dysfunction are disproportionately enriched in the brain compared with other tissues, which we refer to as brain-specific expressed genes. In addition, the tissue-specificity gene expression index can be used as a handle to identify non-brain contexts that are involved in NND pathogenesis. Lastly, we discuss how using an AI approach provides the opportunity to integrate the biological impacts of risk loci to identify those putative combinations of causative relationships through which genetic factors contribute to NND pathogenesis.

Relationship between TSHR, BRAF and PIK3CA gene copy number variations and thyroid nodules

PURPOSE

This study aims to investigate the relationship between the TSHR, BRAF, and PIK3CA gene copy number variations (CNVs) and thyroid nodules by analyzing gene CNVs, and to explore the interaction between iodine status and the above genes CNVs in the occurrence of thyroid nodules.

METHODS

Three hundred and ninety-five subjects were selected from 3 regions with different iodine status in Shanxi Province of China, including 192 patients with thyroid nodules and 203 healthy controls. The basic information about subjects had been obtained through a questionnaire. B ultrasound was utilized to check thyroid nodules. Blood and urine samples were harvested to detect the thyroid function and urinary iodine concentration. Real-time quantitative polymerase chains reaction (RT-PCR) served to detect CNVs in DNA from human blood.

RESULTS

There was an association between TSHR gene CNV and thyroid nodules (χ² = 8.403, P = 0.004). The prevalence of BRAF and PIK3CA gene CNVs was not statistically significant between the case group and the control group. Differences in the TSHR gene CNV rates for cases of the 3 areas were statistically significant (χ² = 10.072, P = 0.007). No statistical difference in the prevalence rates of the 3 genes CNVs between diverse characteristics of thyroid nodules was observed. UIC > 300 μg/L (OR = 1.74, 95% CI: 1.02-2.96, P = 0.041) and TSHR gene CNV (OR = 3.53, 95% CI: 1.40-8.92, P = 0.008) were risk factors for thyroid nodules. There was no significant interaction between the UIC and the examined genes CNVs.

CONCLUSIONS

TSHR gene CNV and high urinary iodine levels can increase the risk of thyroid nodules. But the interactions between the 3 above genes CNVs and iodine nutrition were not found in the occurrence of thyroid nodules.

From Stem Cells to Populations-Using hiPSC, Next-Generation Sequencing, and GWAS to Explore the Genetic and Molecular Mechanisms of Congenital Heart Defects

Congenital heart defects (CHD) are developmental malformations affecting the heart and the great vessels. Early heart development requires temporally regulated crosstalk between multiple cell types, signaling pathways, and mechanical forces of early blood flow. While both genetic and environmental factors have been recognized to be involved, identifying causal genes in non-syndromic CHD has been difficult. While variants following Mendelian inheritance have been identified by linkage analysis in a few families with multiple affected members, the inheritance pattern in most familial cases is complex, with reduced penetrance and variable expressivity. Furthermore, most non-syndromic CHD are sporadic. Improved sequencing technologies and large biobank collections have enabled genome-wide association studies (GWAS) in non-syndromic CHD. The ability to generate human to create human induced pluripotent stem cells (hiPSC) and further differentiate them to organotypic cells enables further exploration of genotype–phenotype correlations in patient-derived cells. Here we review how these technologies can be used in unraveling the genetics and molecular mechanisms of heart development.

Genomics of Adaptations in Ungulates

Ungulates are a group of hoofed animals that have long interacted with humans as essential sources of food, labor, clothing, and transportation. These consist of domesticated, feral, and wild species raised in a wide range of habitats and biomes. Given the diverse and extreme environments inhabited by ungulates, unique adaptive traits are fundamental for fitness. The documentation of genes that underlie their genomic signatures of selection is crucial in this regard. The increasing availability of advanced sequencing technologies has seen the rapid growth of ungulate genomic resources, which offers an exceptional opportunity to understand their adaptive evolution. Here, we summarize the current knowledge on evolutionary genetic signatures underlying the adaptations of ungulates to different habitats.

The Role of Structural Variation in Adaptation and Evolution of Yeast and Other Fungi

Mutations in DNA can be limited to one or a few nucleotides, or encompass larger deletions, insertions, duplications, inversions and translocations that span long stretches of DNA or even full chromosomes. These so-called structural variations (SVs) can alter the gene copy number, modify open reading frames, change regulatory sequences or chromatin structure and thus result in major phenotypic changes. As some of the best-known examples of SV are linked to severe genetic disorders, this type of mutation has traditionally been regarded as negative and of little importance for adaptive evolution. However, the advent of genomic technologies uncovered the ubiquity of SVs even in healthy organisms. Moreover, experimental evolution studies suggest that SV is an important driver of evolution and adaptation to new environments. Here, we provide an overview of the causes and consequences of SV and their role in adaptation, with specific emphasis on fungi since these have proven to be excellent models to study SV.

The contribution of CNVs to the most common aging-related neurodegenerative diseases

Alzheimer and Parkinson's diseases are neurodegenerative aging-related pathological conditions, mainly caused by the interplay of genetic and non-genetic factors and whose incidence rate is going to drastically increase given the growing life expectancy. To address these complex multifactorial traits, a systems biology strategy is needed to highlight genotype-phenotype correlations as well as overlapping gene signatures. Copy number variants (CNVs) are structural chromosomal imbalances that can have pathogenic nature causing or contributing to the disease onset or progression. Moreover, neurons affected by CNVs have been found to decline in number depending on age in healthy controls and may be selectively vulnerable to aging-related cell-death. In this review, we aim to update the reader on the role of these variations in the pathogenesis of Alzheimer and Parkinson diseases. To widen the comprehension of pathogenic mechanisms underlying them, we discuss variations detected from blood or brain specimens, as well as overlapped signatures between the two pathologies.

Extensive genome-wide duplications in the eastern oyster (Crassostrea virginica)

Genomic structural variation is an important source of genetic and phenotypic diversity, playing a critical role in evolution. The recent availability of a high-quality reference genome for the eastern oyster, Crassostrea virginica, and whole-genome sequence data of samples from across the species range in the USA, provides an opportunity to explore structural variation across the genome of this species. Our analysis shows significantly greater individual-level duplications of regions across the genome than that of most model vertebrate species. Duplications are widespread across all ten chromosomes with variation in frequency per chromosome. The eastern oyster shows a large interindividual variation in duplications as well as particular chromosomal regions with a higher density of duplications. A high percentage of duplications seen in C. virginica lie completely within genes and exons, suggesting the potential for impacts on gene function. These results support the hypothesis that structural changes may play a significant role in standing genetic variation in C. virginica, and potentially have a role in their adaptive and evolutionary success. Altogether, these results suggest that copy number variation plays an important role in the genomic variation of C. virginica. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Genomic and transcriptomic characterization of the human glioblastoma cell line AHOL1

Cancer cell lines are widely used as in vitro models of tumorigenesis, facilitating fundamental discoveries in cancer biology and translational medicine. Currently, there are few options for glioblastoma (GBM) treatment and limited in vitro models with accurate genomic and transcriptomic characterization. Here, a detailed characterization of a new GBM cell line, namely AHOL1, was conducted in order to fully characterize its molecular composition based on its karyotype, copy number alteration (CNA), and transcriptome profiling, followed by the validation of key elements associated with GBM tumorigenesis. Large numbers of CNAs and differentially expressed genes (DEGs) were identified. CNAs were distributed throughout the genome, including gains at Xq11.1-q28, Xp22.33-p11.1, Xq21.1-q21.33, 4p15.1-p14, 8q23.2-q23.3 and losses at Yq11.21-q12, Yp11.31-p11.2, and 15q11.1-q11.2 positions. Nine druggable genes were identified, including HCRTR2, ETV1, PTPRD, PRKX, STS, RPS6KA6, ZFY, USP9Y, and KDM5D. By integrating DEGs and CNAs, we identified 57 overlapping genes enriched in fourteen pathways. Altered expression of several cancer-related candidates found in the DEGs-CNA dataset was confirmed by RT-qPCR. Taken together, this first comprehensive genomic and transcriptomic landscape of AHOL1 provides unique resources for further studies and identifies several druggable targets that may be useful for therapeutics and biologic and molecular investigation of GBM.

Cell Adhesion Molecules Involved in Neurodevelopmental Pathways Implicated in 3p-Deletion Syndrome and Autism Spectrum Disorder

Autism spectrum disorder (ASD) is characterized by impaired social interaction, language delay and repetitive or restrictive behaviors. With increasing prevalence, ASD is currently estimated to affect 0.5–2.0% of the global population. However, its etiology remains unclear due to high genetic and phenotypic heterogeneity. Copy number variations (CNVs) are implicated in several forms of syndromic ASD and have been demonstrated to contribute toward ASD development by altering gene dosage and expression. Increasing evidence points toward the p-arm of chromosome 3 (chromosome 3p) as an ASD risk locus. Deletions occurring at chromosome 3p result in 3p-deletion syndrome (Del3p), a rare genetic disorder characterized by developmental delay, intellectual disability, facial dysmorphisms and often, ASD or ASD-associated behaviors. Therefore, we hypothesize that overlapping molecular mechanisms underlie the pathogenesis of Del3p and ASD. To investigate which genes encoded in chromosome 3p could contribute toward Del3p and ASD, we performed a comprehensive literature review and collated reports investigating the phenotypes of individuals with chromosome 3p CNVs. We observe that high frequencies of CNVs occur in the 3p26.3 region, the terminal cytoband of chromosome 3p. This suggests that CNVs disrupting genes encoded within the 3p26.3 region are likely to contribute toward the neurodevelopmental phenotypes observed in individuals affected by Del3p. The 3p26.3 region contains three consecutive genes encoding closely related neuronal immunoglobulin cell adhesion molecules (IgCAMs): Close Homolog of L1 (CHL1), Contactin-6 (CNTN6), and Contactin-4 (CNTN4). CNVs disrupting these neuronal IgCAMs may contribute toward ASD phenotypes as they have been associated with key roles in neurodevelopment. CHL1, CNTN6, and CNTN4 have been observed to promote neurogenesis and neuronal survival, and regulate neuritogenesis and synaptic function. Furthermore, there is evidence that these neuronal IgCAMs possess overlapping interactomes and participate in common signaling pathways regulating axon guidance. Notably, mouse models deficient for these neuronal IgCAMs do not display strong deficits in axonal migration or behavioral phenotypes, which is in contrast to the pronounced defects in neuritogenesis and axon guidance observed in vitro. This suggests that when CHL1, CNTN6, or CNTN4 function is disrupted by CNVs, other neuronal IgCAMs may suppress behavioral phenotypes by compensating for the loss of function.

Genetic alterations in the 3q26.31-32 locus confer an aggressive prostate cancer phenotype

Large-scale genetic aberrations that underpin prostate cancer development and progression, such as copy-number alterations (CNAs), have been described but the consequences of specific changes in many identified loci is limited. Germline SNPs in the 3q26.31 locus are associated with aggressive prostate cancer, and is the location of NAALADL2, a gene overexpressed in aggressive disease. The closest gene to NAALADL2 is TBL1XR1, which is implicated in tumour development and progression. Using publicly-available cancer genomic data we report that NAALADL2 and TBL1XR1 gains/amplifications are more prevalent in aggressive sub-types of prostate cancer when compared to primary cohorts. In primary disease, gains/amplifications occurred in 15.99% (95% CI: 13.02–18.95) and 14.96% (95% CI: 12.08–17.84%) for NAALADL2 and TBL1XR1 respectively, increasing in frequency in higher Gleason grade and stage tumours. Gains/amplifications result in transcriptional changes and the development of a pro-proliferative and aggressive phenotype. These results support a pivotal role for copy-number gains in this genetic region.

Benjamin Simpson et al. use publicly available cancer genomic data to investigate copy number changes at the 3q26.31–32 locus, which has been associated with aggressive prostate cancer based on single-nucleotide polymorphisms. They find that gains of NAALADL2 and TBL1XR1 in this locus are associated with more aggressive subtypes of prostate cancer and the transcription of pro-proliferative signalling processes.

Systematic identification of genetic systems associated with phenotypes in patients with rare genomic copy number variations

Copy number variation (CNV) related disorders tend to show complex phenotypic profiles that do not match known diseases. This makes it difficult to ascertain their underlying molecular basis. A potential solution is to compare the affected genomic regions for multiple patients that share a pathological phenotype, looking for commonalities. Here, we present a novel approach to associate phenotypes with functional systems, in terms of GO categories and KEGG and Reactome pathways, based on patient data. The approach uses genomic and phenomic data from the same patients, finding shared genomic regions between patients with similar phenotypes. These regions are mapped to genes to find associated functional systems. We applied the approach to analyse patients in the DECIPHER database with de novo CNVs, finding functional systems associated with most phenotypes, often due to mutations affecting related genes in the same genomic region. Manual inspection of the ten top-scoring phenotypes found multiple FunSys connections supported by the previous studies for seven of them. The workflow also produces reports focussed on the genes and FunSys connected to the different phenotypes, alongside patient-specific reports, which give details of the associated genes and FunSys for each individual in the cohort. These can be run in "confidential" mode, preserving patient confidentiality. The workflow presented here can be used to associate phenotypes with functional systems using data at the level of a whole cohort of patients, identifying important connections that could not be found when considering them individually. The full workflow is available for download, enabling it to be run on any patient cohort for which phenotypic and CNV data are available.

Copy number variation of IL17RA gene and its association with the ankylosing spondylitis risk in Iranian patients: a case-control study

Background

Ankylosing spondylitis (AS) is considered as a subtype of spondyloarthritis (SpA) that mainly leads to fatigue, stiffness, spinal ankylosis, and impaired physical functions with reduced quality of life. Interleukin (IL)-17A provokes additional inflammatory mediators and recruits immune cells to the inflamed site. IL17 expression increased in various inflammatory disorders including psoriasis, rheumatoid arthritis, multiple sclerosis, crohn’s disease, and ankylosing spondylitis. The current study aimed to evaluate the association of IL17RA copy number changes with the susceptibility to AS and their correlation to IL17RA expression in Iranian population.

Methods

IL17RA copy number genotyping assessments were carried out in 455 AS patients and 450 healthy controls, using custom TaqMan CNV assays. TaqMan primers and probe were located in Chr.22:17109553 based on pre-designed IL17RA Copy Number Assay ID, Hs02339506_cn. mRNA expression of IL17RA was also measured by SYBR Green real-time polymerase chain reaction (PCR).

Results

A IL17RA copy number loss (< 2) was associated with AS compared to 2 copies as reference (OR:2.18, 95% CI: (1.38–3.44), P-value < 0.001) and increased the risk of AS. IL17RA mRNA expression showed a significant increase in peripheral blood mononuclear cells (PBMCs) of all AS individuals than controls. The mRNA expression level of 2 copies was significantly higher in AS patients.

Conclusions

Our findings revealed that a low copy number of IL17RA might confer a susceptibility risk to AS. However, it is probably not directly involved in the regulation of IL17RA mRNA expression. Epigenetic mechanisms like DNA methylation, post-transcriptional, and -translational modifications that regulate the expression of the genes may contribute in upregulation of IL17RA mRNA expression in the loss of gene copy number condition.

Functional annotation of rare structural variation in the human brain

Structural variants (SVs) contribute to many disorders, yet, functionally annotating them remains a major challenge. Here, we integrate SVs with RNA-sequencing from human post-mortem brains to quantify their dosage and regulatory effects. We show that genic and regulatory SVs exist at significantly lower frequencies than intergenic SVs. Functional impact of copy number variants (CNVs) stems from both the proportion of genic and regulatory content altered and loss-of-function intolerance of the gene. We train a linear model to predict expression effects of rare CNVs and use it to annotate regulatory disruption of CNVs from 14,891 independent genome-sequenced individuals. Pathogenic deletions implicated in neurodevelopmental disorders show significantly more extreme regulatory disruption scores and if rank ordered would be prioritized higher than using frequency or length alone. This work shows the deleteriousness of regulatory SVs, particularly those altering CTCF sites and provides a simple approach for functionally annotating the regulatory consequences of CNVs.

De novo duplication on Chromosome 19 observed in nuclear family displaying neurodevelopmental disorders

Pleiotropy and variable expressivity have been cited to explain the seemingly distinct neurodevelopmental disorders due to a common genetic etiology within the same family. Here we present a family with a de novo 1-Mb duplication involving 18 genes on Chromosome 19. Within the family there are multiple cases of neurodevelopmental disorders including autism spectrum disorder, attention deficit/hyperactivity disorder, intellectual disability, and psychiatric disease in individuals carrying this copy-number variant (CNV). Quantitative polymerase chain reaction (PCR) confirmed the CNV was de novo in the mother and inherited by both sons. Whole-exome sequencing did not uncover further genetic risk factors segregating within the family. Transcriptome analysis of peripheral blood demonstrated a ∼1.5-fold increase in RNA transcript abundance in 12 of the 15 detected genes within the CNV region for individuals carrying the CNV compared with their noncarrier relatives. Examination of transcript abundance across the rest of the transcriptome identified 407 differentially expressed genes (P-value < 0.05; adjusted P-value < 0.1) mapping to immune response, response to endoplasmic reticulum stress, and regulation of epithelial cell proliferation pathways. 16S microbiome profiling demonstrated compositional difference in the gut bacteria between the half-brothers. These results raise the possibility that the observed CNV may contribute to the varied phenotypic characteristics in family members through alterations in gene expression and/or dysbiosis of the gut microbiome. More broadly, there is growing evidence that different neurodevelopmental and psychiatric disorders can share the same genetic variant, which lays a framework for later neurodevelopmental and psychiatric manifestations.

Copy number variation in human genomes from three major ethno-linguistic groups in Africa

Background

Copy number variation is an important class of genomic variation that has been reported in 75% of the human genome. However, it is underreported in African populations. Copy number variants (CNVs) could have important impacts on disease susceptibility and environmental adaptation. To describe CNVs and their possible impacts in Africans, we sequenced genomes of 232 individuals from three major African ethno-linguistic groups: (1) Niger Congo A from Guinea and Côte d’Ivoire, (2) Niger Congo B from Uganda and the Democratic Republic of Congo and (3) Nilo-Saharans from Uganda. We used GenomeSTRiP and cn.MOPS to identify copy number variant regions (CNVRs).

Results

We detected 7608 CNVRs, of which 2172 were only deletions, 2384 were only insertions and 3052 had both. We detected 224 previously un-described CNVRs. The majority of novel CNVRs were present at low frequency and were not shared between populations. We tested for evidence of selection associated with CNVs and also for population structure. Signatures of selection identified previously, using SNPs from the same populations, were overrepresented in CNVRs. When CNVs were tagged with SNP haplotypes to identify SNPs that could predict the presence of CNVs, we identified haplotypes tagging 3096 CNVRs, 372 CNVRs had SNPs with evidence of selection (iHS > 3) and 222 CNVRs had both. This was more than expected (p < 0.0001) and included loci where CNVs have previously been associated with HIV, Rhesus D and preeclampsia. When integrated with 1000 Genomes CNV data, we replicated their observation of population stratification by continent but no clustering by populations within Africa, despite inclusion of Nilo-Saharans and Niger-Congo populations within our dataset.

Conclusions

Novel CNVRs in the current study increase representation of African diversity in the database of genomic variants. Over-representation of CNVRs in SNP signatures of selection and an excess of SNPs that both tag CNVs and are subject to selection show that CNVs may be the actual targets of selection at some loci. However, unlike SNPs, CNVs alone do not resolve African ethno-linguistic groups. Tag haplotypes for CNVs identified may be useful in predicting African CNVs in future studies where only SNP data is available.

A Roadmap for Understanding the Evolutionary Significance of Structural Genomic Variation

Structural genomic variants (SVs) are ubiquitous and play a major role in adaptation and speciation. Yet, comparative and population genomics have focused predominantly on gene duplications and large-effect inversions. The lack of a common framework for studying all SVs is hampering progress towards a more systematic assessment of their evolutionary significance. Here we (i) review how different types of SVs affect ecological and evolutionary processes; (ii) suggest unifying definitions and recommendations for future studies; and (iii) provide a roadmap for the integration of SVs in ecoevolutionary studies. In doing so, we lay the foundation for population genomics, theoretical, and experimental approaches to understand how the full spectrum of SVs impacts ecological and evolutionary processes.

Multi-Omics Analysis Reveals the Pan-Cancer Landscape of Bone Morphogenetic Proteins

Background

Bone morphogenetic proteins (BMPs) are widely involved in cancer development. However, a wealth of conflicting data raises the question of whether BMPs serve as oncogenes or as cancer suppressors.

Material/Methods

By integrating multi-omics data across cancers, we comprehensively analyzed the genomic and pharmacogenomic landscape of BMP genes across cancers.

Results

Surprisingly, our data indicate that BMPs are globally downregulated in cancers. Further genetics and epigenetics analyses show that this abnormal expression is driven by copy number variations, especially heterozygous amplification. We next assessed the BMP-associated pathways and demonstrated that they suppress cell cycle and estrogen hormone pathways. Bone morphogenetic protein interacts with 58 compounds, and their dysfunction can induce drug sensitivity.

Conclusions

Our results define the landscape of the BMP family at a systems level and open potential therapeutic opportunities for cancer patients.

Integrative analysis of genomic amplification-dependent expression and loss-of-function screen identifies ASAP1 as a driver gene in triple-negative breast cancer progression

The genetically heterogeneous triple-negative breast cancer (TNBC) continues to be an intractable disease, due to lack of effective targeted therapies. Gene amplification is a major event in tumorigenesis. Genes with amplification-dependent expression are being explored as therapeutic targets for cancer treatment. In this study, we have applied Analytical Multi-scale Identification of Recurring Events analysis and transcript quantification in the TNBC genome across 222 TNBC tumors and identified 138 candidate genes with positive correlation in copy number gain (CNG) and gene expression. siRNA-based loss-of-function screen of the candidate genes has validated EGFR, MYC, ASAP1, IRF2BP2, and CCT5 genes as drivers promoting proliferation in different TNBC cells. MYC, ASAP1, IRF2BP2, and CCT5 display frequent CNG and concurrent expression over 2173 breast cancer tumors (cBioPortal dataset). More frequently are MYC and ASAP1 amplified in TNBC tumors (>30%, n = 320). In particular, high expression of ASAP1, the ADP-ribosylation factor GTPase-activating protein, is significantly related to poor metastatic relapse-free survival of TNBC patients (n = 257, bc-GenExMiner). Furthermore, we have revealed that silencing of ASAP1 modulates numerous cytokine and apoptosis signaling components, such as IL1B, TRAF1, AIFM2, and MAP3K11 that are clinically relevant to survival outcomes of TNBC patients. ASAP1 has been reported to promote invasion and metastasis in various cancer cells. Our findings that ASAP1 is an amplification-dependent TNBC driver gene promoting TNBC cell proliferation, functioning upstream apoptosis components, and correlating to clinical outcomes of TNBC patients, support ASAP1 as a potential actionable target for TNBC treatment.

Resources for functional genomic studies of health and development in nonhuman primates

Primates display a wide range of phenotypic variation underlaid by complex genetically regulated mechanisms. The links among DNA sequence, gene function, and phenotype have been of interest from an evolutionary perspective, to understand functional genome evolution and its phenotypic consequences, and from a biomedical perspective to understand the shared and human-specific roots of health and disease. Progress in methods for characterizing genetic, transcriptomic, and DNA methylation (DNAm) variation is driving the rapid development of extensive omics resources, which are now increasingly available from humans as well as a growing number of nonhuman primates (NHPs). The fast growth of large-scale genomic data is driving the emergence of integrated tools and databases, thus facilitating studies of gene functionality across primates. This review describes NHP genomic resources that can aid in exploration of how genes shape primate phenotypes. It focuses on the gene expression trajectories across development in different tissues, the identification of functional genetic variation (including variants deleterious for protein function and regulatory variants modulating gene expression), and DNAm profiles as an emerging tool to understand the process of aging. These resources enable comparative functional genomics approaches to identify species-specific and primate-shared gene functionalities associated with health and development.