The impact of copy number variation on local gene expression in mouse hematopoietic stem and progenitor cells

Molecular Assessment of Methylglyoxal-Induced Toxicity and Therapeutic Approaches in Various Diseases: Exploring the Interplay with the Glyoxalase System

This comprehensive exploration delves into the intricate interplay of methylglyoxal (MG) and glyoxalase 1 (GLO I) in various physiological and pathological contexts. The linchpin of the narrative revolves around the role of these small molecules in age-related issues, diabetes, obesity, cardiovascular diseases, and neurodegenerative disorders. Methylglyoxal, a reactive dicarbonyl metabolite, takes center stage, becoming a principal player in the development of AGEs and contributing to cell and tissue dysfunction. The dual facets of GLO I-activation and inhibition-unfold as potential therapeutic avenues. Activators, spanning synthetic drugs like candesartan to natural compounds like polyphenols and isothiocyanates, aim to restore GLO I function. These molecular enhancers showcase promising outcomes in conditions such as diabetic retinopathy, kidney disease, and beyond. On the contrary, GLO I inhibitors emerge as crucial players in cancer treatment, offering new possibilities in diseases associated with inflammation and multidrug resistance. The symphony of small molecules, from GLO I activators to inhibitors, presents a nuanced understanding of MG regulation. From natural compounds to synthetic drugs, each element contributes to a molecular orchestra, promising novel interventions and personalized approaches in the pursuit of health and wellbeing. The abstract concludes with an emphasis on the necessity of rigorous clinical trials to validate these findings and acknowledges the importance of individual variability in the complex landscape of health.

Comprehensive Analysis of Copy Number Variations on Glycoside Hydrolase 45 Genes among Different Bursaphelenchus xylophilus Strains

Bursaphelenchus xylophilus is considered the most dangerous quarantine pest in China. It causes enormous economic and ecological losses in many countries from Asia and Europe. The glycoside hydrolase 45 gene family has been demonstrated in early studies to contribute to the cell wall degradation ability of B. xylophilus during its infection. However, the copy number variation (CNV) of the GH45 gene and its association with B. xylophilus pathogenicity were not fully elucidated. In this study, we found that the GH45 gene with two copies is the most predominant type among 259 B. xylophilus strains collected from China and Japan. Additionally, 18 strains are identified as GH45 genes with a single copy, and only two strains are verified to have three copies. Subsequent expression analysis and inoculation test suggest that the copy numbers of the GH45 gene are correlated with gene expression as well as the B. xylophilus pathogenicity. B. xylophilus strains with more copies of the GH45 gene usually exhibit more abundant expression and cause more severe wilt symptoms on pine trees. The aforementioned results indicated the potential regulatory effects of CNV in B. xylophilus and provided novel information to better understand the molecular pathogenesis of this devastating pest.

Virulence and Antibiotic Resistance Genes in Listeria monocytogenes Strains Isolated From Ready-to-Eat Foods in Chile

Listeria monocytogenes is causing listeriosis, a rare but severe foodborne infection. Listeriosis affects pregnant women, newborns, older adults, and immunocompromised individuals. Ready-to-eat (RTE) foods are the most common sources of transmission of the pathogen This study explored the virulence factors and antibiotic resistance in L. monocytogenes strains isolated from ready-to-eat (RTE) foods through in vitro and in silico testing by whole-genome sequencing (WGS). The overall positivity of L. monocytogenes in RTE food samples was 3.1% and 14 strains were isolated. L. monocytogenes ST8, ST2763, ST1, ST3, ST5, ST7, ST9, ST14, ST193, and ST451 sequence types were identified by average nucleotide identity, ribosomal multilocus sequence typing (rMLST), and core genome MLST. Seven isolates had serotype 1/2a, five 1/2b, one 4b, and one 1/2c. Three strains exhibited in vitro resistance to ampicillin and 100% of the strains carried the fosX, lin, norB, mprF, tetA, and tetC resistance genes. In addition, the arsBC, bcrBC, and clpL genes were detected, which conferred resistance to stress and disinfectants. All strains harbored hlyA, prfA, and inlA genes almost thirty-two the showed the bsh, clpCEP, hly, hpt, iap/cwhA, inlA, inlB, ipeA, lspA, mpl, plcA, pclB, oat, pdgA, and prfA genes. One isolate exhibited a type 11 premature stop codon (PMSC) in the inlA gene and another isolate a new mutation (deletion of A in position 819). The Inc18(rep25), Inc18(rep26), and N1011A plasmids and MGEs were found in nine isolates. Ten isolates showed CAS-Type II-B systems; in addition, Anti-CRISPR AcrIIA1 and AcrIIA3 phage-associated systems were detected in three genomes. These virulence and antibiotic resistance traits in the strains isolated in the RTE foods indicate a potential public health risk for consumers.

Distribution of Copy Number Variation in SYT11 Gene and Its Association with Growth Conformation Traits in Chinese Cattle

Simple Summary

It is known that many different breeds of cattle are widely distributed in China. However, due to a lengthy selection of draught direction, there are obvious shortcomings in Chinese cattle, such as less meat production, slow weight gain, poor meat quality, and a lack of specialized beef cattle breeds. Animal breeding heavily benefits from molecular technologies, among which molecular genetic markers were widely used to improve the economic traits of beef cattle. Because the copy number variation (CNV) involves a longer DNA sequence or even the entire functional gene, it may have a greater impact on the phenotype. Recent studies have indicated that CNVs are widespread in the Chinese cattle genome. By investigating the effects of CNVs on gene expression and cattle traits, we aim to find those genomic variations which could significantly affect cattle traits, and which could provide a basis for genetic selection and molecular breeding of local Chinese cattle.

Abstract

Currently, studies of the SYT11 gene mainly focus on neurological diseases such as schizophrenia and Parkinson’s disease. However, some studies have shown that the C2B domain of SYT11 can interact with RISC components and affect the gene regulation of miRNA, which is important for cell differentiation, proliferation, and apoptosis, and therefore has an impact on muscle growth and development in animals. The whole-genome resequencing data detected a CNV in the SYT11 gene, and this may affect cattle growth traits. In this study, CNV distribution of 672 individuals from four cattle breeds, Yunling, Pinan, Xianan, and Qinchuan, were detected by qPCR. The relationship between CNV, gene expression and growth traits was further investigated. The results showed that the proportion of multiple copy types was the largest in all cattle breeds, but there were some differences among different breeds. The normal type had higher gene expression than the abnormal copy type. The CNVs of the SYT11 gene were significantly correlated with body length, cannon circumference, chest depth, rump length, and forehead size of Yunling cattle, and was significantly correlated with the bodyweight of Xianan cattle, respectively. These data improve our understanding of the effects of CNV on cattle growth traits. Our results suggest that the CNV of SYT11 gene is a protentional molecular marker, which may be used to improve growth traits in Chinese cattle.

The Glyoxalase System in Age-Related Diseases: Nutritional Intervention as Anti-Ageing Strategy

The glyoxalase system is critical for the detoxification of advanced glycation end-products (AGEs). AGEs are toxic compounds resulting from the non-enzymatic modification of biomolecules by sugars or their metabolites through a process called glycation. AGEs have adverse effects on many tissues, playing a pathogenic role in the progression of molecular and cellular aging. Due to the age-related decline in different anti-AGE mechanisms, including detoxifying mechanisms and proteolytic capacities, glycated biomolecules are accumulated during normal aging in our body in a tissue-dependent manner. Viewed in this way, anti-AGE detoxifying systems are proposed as therapeutic targets to fight pathological dysfunction associated with AGE accumulation and cytotoxicity. Here, we summarize the current state of knowledge related to the protective mechanisms against glycative stress, with a special emphasis on the glyoxalase system as the primary mechanism for detoxifying the reactive intermediates of glycation. This review focuses on glyoxalase 1 (GLO1), the first enzyme of the glyoxalase system, and the rate-limiting enzyme of this catalytic process. Although GLO1 is ubiquitously expressed, protein levels and activities are regulated in a tissue-dependent manner. We provide a comparative analysis of GLO1 protein in different tissues. Our findings indicate a role for the glyoxalase system in homeostasis in the eye retina, a highly oxygenated tissue with rapid protein turnover. We also describe modulation of the glyoxalase system as a therapeutic target to delay the development of age-related diseases and summarize the literature that describes the current knowledge about nutritional compounds with properties to modulate the glyoxalase system.

Genomic Copy Number Variation Study of Nine Macaca Species Provides New Insights into Their Genetic Divergence, Adaptation, and Biomedical Application

Copy number variation (CNV) can promote phenotypic diversification and adaptive evolution. However, the genomic architecture of CNVs among Macaca species remains scarcely reported, and the roles of CNVs in adaptation and evolution of macaques have not been well addressed. Here, we identified and characterized 1,479 genome-wide hetero-specific CNVs across nine Macaca species with bioinformatic methods, along with 26 CNV-dense regions and dozens of lineage-specific CNVs. The genes intersecting CNVs were overrepresented in nutritional metabolism, xenobiotics/drug metabolism, and immune-related pathways. Population-level transcriptome data showed that nearly 46% of CNV genes were differentially expressed across populations and also mainly consisted of metabolic and immune-related genes, which implied the role of CNVs in environmental adaptation of Macaca. Several CNVs overlapping drug metabolism genes were verified with genomic quantitative polymerase chain reaction, suggesting that these macaques may have different drug metabolism features. The CNV-dense regions, including 15 first reported here, represent unstable genomic segments in macaques where biological innovation may evolve. Twelve gains and 40 losses specific to the Barbary macaque contain genes with essential roles in energy homeostasis and immunity defense, inferring the genetic basis of its unique distribution in North Africa. Our study not only elucidated the genetic diversity across Macaca species from the perspective of structural variation but also provided suggestive evidence for the role of CNVs in adaptation and genome evolution. Additionally, our findings provide new insights into the application of diverse macaques to drug study.

Structural variation at the CYP2C locus: Characterization of deletion and duplication alleles

The human CYP2C locus harbors the polymorphic CYP2C18, CYP2C19, CYP2C9, and CYP2C8 genes, and of these, CYP2C19 and CYP2C9 are directly involved in the metabolism of ~15% of all medications. All variant CYP2C19 and CYP2C9 star (*) allele haplotypes currently cataloged by the Pharmacogene Variation (PharmVar) Consortium are defined by sequence variants. To determine if structural variation also occurs at the CYP2C locus, the 10q23.33 region was interrogated across deidentified clinical chromosomal microarray (CMA) data from 20,642 patients tested at two academic medical centers. Fourteen copy number variants that affected the coding region of CYP2C genes were detected in the clinical CMA cohorts, which ranged in size from 39.2 to 1,043.3 kb. Selected deletions and duplications were confirmed by MLPA or ddPCR. Analysis of the clinical CMA and an additional 78,839 cases from the Database of Genomic Variants (DGV) and ClinGen (total n = 99,481) indicated that the carrier frequency of a CYP2C structural variant is ~1 in 1,000, with ~1 in 2,000 being a CYP2C19 full gene or partial-gene deletion carrier, designated by PharmVar as CYP2C19*36 and *37, respectively. Although these structural variants are rare in the general population, their detection will likely improve metabolizer phenotype prediction when interrogated for research and/or clinical testing.

Methylglyoxal, a Highly Reactive Dicarbonyl Compound, in Diabetes, Its Vascular Complications, and Other Age-Related Diseases

The formation and accumulation of methylglyoxal (MGO), a highly reactive dicarbonyl compound, has been implicated in the pathogenesis of type 2 diabetes, vascular complications of diabetes, and several other age-related chronic inflammatory diseases such as cardiovascular disease, cancer, and disorders of the central nervous system. MGO is mainly formed as a byproduct of glycolysis and, under physiological circumstances, detoxified by the glyoxalase system. MGO is the major precursor of nonenzymatic glycation of proteins and DNA, subsequently leading to the formation of advanced glycation end products (AGEs). MGO and MGO-derived AGEs can impact on organs and tissues affecting their functions and structure. In this review we summarize the formation of MGO, the detoxification of MGO by the glyoxalase system, and the biochemical pathways through which MGO is linked to the development of diabetes, vascular complications of diabetes, and other age-related diseases. Although interventions to treat MGO-associated complications are not yet available in the clinical setting, several strategies to lower MGO have been developed over the years. We will summarize several new directions to target MGO stress including glyoxalase inducers and MGO scavengers. Targeting MGO burden may provide new therapeutic applications to mitigate diseases in which MGO plays a crucial role.

Integrative analysis of rare copy number variants and gene expression data in alopecia areata implicates an aetiological role for autophagy

Alopecia areata (AA) is a highly prevalent autoimmune disease that attacks the hair follicle and leads to hair loss that can range from small patches to complete loss of scalp and body hair. Our previous linkage and genome-wide association studies (GWAS) generated strong evidence for aetiological contributions from inherited genetic variants at different population frequencies, including both rare mutations and common polymorphisms. Additionally, we conducted gene expression (GE) studies on scalp biopsies of 96 patients and controls to establish signatures of active disease. In this study, we performed an integrative analysis on these two datasets to test the hypothesis that rare CNVs in patients with AA could be leveraged to identify drivers of disease in our AA GE signatures. We analysed copy number variants (CNVs) in a case-control cohort of 673 patients with AA and 16 311 controls independent of the case-control cohort of 96 research participants used in our GE study. Using an integrative computational analysis, we identified 14 genes whose expression levels were altered by CNVs in a consistent direction of effect, corresponding to gene expression changes in lesional skin of patients. Four of these genes were affected by CNVs in three or more unrelated patients with AA, including ATG4B and SMARCA2, which are involved in autophagy and chromatin remodelling, respectively. Our findings identified new classes of genes with potential contributions to AA pathogenesis.

Comparative sequence alignment reveals River Buffalo genomic structural differences compared with cattle

This study sought to characterize differences in gene content, regulation and structure between taurine cattle and river buffalo (one subspecies of domestic water buffalo) using the extensively annotated UMD3.1 cattle reference genome as a basis for comparisons. We identified 127 deletion CNV regions in river buffalo representing 5 annotated cattle genes. We also characterized 583 merged mobile element insertion (MEI) events within the upstream regions of annotated cattle genes. Transcriptome analysis in various tissue types on river buffalo confirmed the absence of four cattle genes. Four genes which may be related to phenotypic differences in meat quality and color, had upstream MEI predictions and were found to have significantly elevated expression in river buffalo compared with cattle. Our comparative alignment approach and gene expression analyses suggested a functional role for many genomic structural variations, which may contribute to the unique phenotypes of river buffalo.

Widespread modulation of gene expression by copy number variation in skeletal muscle

Copy number variation (CNV) is a frequently observed deviation from the diploid state due to duplication or deletion of genomic regions. Although intensively analyzed for association with diseases and production traits, the specific mechanisms and extent by which such variations affect the phenotype are incompletely understood. We present an integrative study on CNV and genome-wide gene expression in Brazilian Bos indicus cattle. We analyzed CNVs inferred from SNP-chip data for effects on gene expression measured with RNA-seq in skeletal muscle samples of 183 steers. Local effects, where expression changes coincided with CNVs in the respective genes, were restricted to immune genes. Distal effects were attributable to several high-impact CNVs that modulated remote expression in an orchestrated and intertwined fashion. These CNVs were located in the vicinity of major skeletal muscle pathway regulators and associated genes were enriched for proteolysis, autophagy, and muscle structure development. From association analysis between CNVs and several meat quality and production traits, we found CNV-associated expression effects to also manifest at the phenotype level. Based on genome sequences of the population founders, we further demonstrate that CNVs with impact on expression and phenotype are passed on from one generation to another.

Copy number variations of five Y chromosome genes in donkeys

Abstract. In mammals, the Y chromosome plays a pivotal role in male sex determination and is essential for normal sperm production. A number of studies were conducted on Y chromosome genes of various species and identified single-copy and multi-copy genes. However, limited studies about donkey Y chromosome genes have been done. In this study, 263 male samples from 13 Chinese donkey breeds were collected to analyze the copy number variations (CNVs) of five Y chromosome genes using the quantitative PCR (qPCR) method. These five genes (cullin 4 B Y (CUL4BY), equus testis-specific transcript y1 (ETSTY1), equus testis-specific transcript y4 (ETSTY4), equus testis-specific transcript Y 5 (ETSTY5), and sex-determining region Y (SRY) were identified as multi-copy, whose median copy numbers (MCNs) were 5, 45, 2, and 2, and 13 with CNV ranges of 1–57, 1–227, 1–37, 1–86 and 1–152, respectively. The CNVs of these five genes were shared in different breeds. Compared to previous studies, the copy numbers of five genes showed some distinct consequences in this study. In particular, the well-known single-copy SRY gene showed CNVs in donkeys. Our results provided genetic variations of donkey Y chromosome genes.

Comparative genomic hybridization and transcriptome sequencing reveal that two genes, OsI_14279 (LOC_Os03g62620) and OsI_10794 (LOC_Os03g14950) regulate the mutation in the γ-rl rice mutant

We previously established the genetic locus of the rolled-leaf mutant, γ-rl, to chromosome 3. In this study, we performed a comparative genomic hybridization (CGH) analysis to identify the genes responsible for the γ-rl mutant phenotype. This was combined with RNA transcriptome sequencing (RNA-seq) to analyze differences in the mRNA expression in seeds 12 h after germination. Using the reference genome of the "indica type" rice from GenBank, we created a chip with 386,000 high density DNA probes designed to target chromosome 3. The genomic DNA from γ-rl and Qinghuazhan (the wild-type) was used for hybridization against the chip to compare signal differences. We uncovered 49 regions with significant differences in hybridization signals including deletions and insertions. RNA-seq analysis between γ-rl and QHZ identified 1060 differentially expressed genes, which potentially regulate numerous biological activities. Moreover, we identified 72 annotated genes in the 49 regions discovered in CGH. Among these, 44 genes showed differential expression in RNA-seq. qRT-PCR validation of the candidate genes confirmed that seven of the 44 genes showed a significant change in their expression levels. Among these, four genes [OsI_10125 (LOC_Os03g06654), OsI_14045 (LOC_Os03g62490), OsI_14279 (LOC_Os03g62620) and OsI_14326 (LOC_Os03g63250)] were down regulated and three genes [(OsI_10794 (LOC_Os03g14950), OsI_11412 (LOC_Os03g21250) and OsI_14152 (LOC_Os03g61360)] were up regulated with a fold change ≥2.0 and a P value ≤ 0.01. Finally, we constructed transgenic plants to study the in vivo functions of these genes. RNAi knock down of LOC_Os03g62620 resulted in rolled-leaf, lower seed-setting and decreased seed growth phenotypes. Transgenic plants with LOC_Os03g14950 over-expression showed dwarf plants with a shortened leaf phenotype. Our results, LOC_Os03g62620 and LOC_Os03g14950 as the essential genes responsible for creating the γ-rl mutant phenotypes suggested that these genes may play crucial roles in regulating rice leaf development and seed growth.

Glyoxalase 1 copy number variation in patients with well differentiated gastro-entero-pancreatic neuroendocrine tumours (GEP-NET)

Background

The glyoxalase-1 gene (GLO1) is a hotspot for copy-number variation (CNV) in human genomes. Increased GLO1 copy-number is associated with multidrug resistance in tumour chemotherapy, but prevalence of GLO1 CNV in gastro-entero-pancreatic neuroendocrine tumours (GEP-NET) is unknown.

Methods

GLO1 copy-number variation was measured in 39 patients with GEP-NET (midgut NET, n = 25; pancreatic NET, n = 14) after curative or debulking surgical treatment. Primary tumour tissue, surrounding healthy tissue and, where applicable, additional metastatic tumour tissue were analysed, using real time qPCR. Progression and survival following surgical treatment were monitored over 4.2 ± 0.5 years.

Results

In the pooled GEP-NET cohort, GLO1 copy-number in healthy tissue was 2.0 in all samples but significantly increased in primary tumour tissue in 43% of patients with pancreatic NET and in 72% of patients with midgut NET, mainly driven by significantly higher GLO1 copy-number in midgut NET. In tissue from additional metastases resection (18 midgut NET and one pancreatic NET), GLO1 copy number was also increased, compared with healthy tissue; but was not significantly different compared with primary tumour tissue. During mean 3 - 5 years follow-up, 8 patients died and 16 patients showed radiological progression. In midgut NET, a high GLO1 copy-number was associated with earlier progression. In NETs with increased GLO1 copy number, there was increased Glo1 protein expression compared to non-malignant tissue.

Conclusions

GLO1 copy-number was increased in a large percentage of patients with GEP-NET and correlated positively with increased Glo1 protein in tumour tissue. Analysis of GLO1 copy-number variation particularly in patients with midgut NET could be a novel prognostic marker for tumour progression.

Methylglyoxal-induced dicarbonyl stress in aging and disease: first steps towards glyoxalase 1-based treatments

Dicarbonyl stress is the abnormal accumulation of dicarbonyl metabolites leading to increased protein and DNA modification contributing to cell and tissue dysfunction in aging and disease. It is produced by increased formation and/or decreased metabolism of dicarbonyl metabolites. MG (methylglyoxal) is a dicarbonyl metabolite of relatively high flux of formation and precursor of the most quantitatively and functionally important spontaneous modifications of protein and DNA clinically. Major MG-derived adducts are arginine-derived hydroimidazolones of protein and deoxyguanosine-derived imidazopurinones of DNA. These are formed non-oxidatively. The glyoxalase system provides an efficient and essential basal and stress-response-inducible enzymatic defence against dicarbonyl stress by the reduced glutathione-dependent metabolism of methylglyoxal by glyoxalase 1. The GLO1 gene encoding glyoxalase 1 has low prevalence duplication and high prevalence amplification in some tumours. Dicarbonyl stress contributes to aging, disease and activity of cytotoxic chemotherapeutic agents. It is found at a low, moderate and severe level in obesity, diabetes and renal failure respectively, where it contributes to the development of metabolic and vascular complications. Increased glyoxalase 1 expression confers multidrug resistance to cancer chemotherapy and has relatively high prevalence in liver, lung and breast cancers. Studies of dicarbonyl stress are providing improved understanding of aging and disease and the basis for rational design of novel pharmaceuticals: glyoxalase 1 inducers for obesity, diabetes and cardiovascular disease and glyoxalase 1 inhibitors for multidrug-resistant tumours. The first clinical trial of a glyoxalase 1 inducer in overweight and obese subjects showed improved glycaemic control, insulin resistance and vascular function.

Latexin Inactivation Enhances Survival and Long-Term Engraftment of Hematopoietic Stem Cells and Expands the Entire Hematopoietic System in Mice

Natural genetic diversity offers an important yet largely untapped resource to decipher the molecular mechanisms regulating hematopoietic stem cell (HSC) function. Latexin (Lxn) is a negative stem cell regulatory gene identified on the basis of genetic diversity. By using an Lxn knockout mouse model, we found that Lxn inactivation in vivo led to the physiological expansion of the entire hematopoietic hierarchy. Loss of Lxn enhanced the competitive repopulation capacity and survival of HSCs in a cell-intrinsic manner. Gene profiling of Lxn-null HSCs showed altered expression of genes enriched in cell-matrix and cell-cell interactions. Thrombospondin 1 (Thbs1) was a potential downstream target with a dramatic downregulation in Lxn-null HSCs. Enforced expression of Thbs1 restored the Lxn inactivation-mediated HSC phenotypes. This study reveals that Lxn plays an important role in the maintenance of homeostatic hematopoiesis, and it may lead to development of safe and effective approaches to manipulate HSCs for clinical benefit.

A vast genomic deletion in the C56BL/6 genome affects different genes within the Ifi200 cluster on chromosome 1 and mediates obesity and insulin resistance

Background

Obesity, the excessive accumulation of body fat, is a highly heritable and genetically heterogeneous disorder. The complex, polygenic basis for the disease consisting of a network of different gene variants is still not completely known.

Results

In the current study we generated a BAC library of the obese-prone NZO strain to clarify the genomic alteration within the gene cluster Ifi200 on chr.1 including Ifi202b, an obesity gene that is in contrast to NZO not expressed in the lean B6 mouse. With the PacBio sequencing data of NZO BAC clones we identified a deletion spanning approximately 261.8 kb in the B6 reference genome. The deletion affects different members of the Ifi200 gene family which also includes the original first exon and 5′-regulatory parts of the Ifi202b gene and suggests to be the relevant cause of its expression deficiency in B6. In addition, the generation and characterization of congenic mice carrying the critical fragment on the B6 background demonstrate its crucial role for obesity and insulin resistance.

Conclusions

Our data reveal the reconstruction of a complex genomic region on mouse chr.1 resulting from deletions and duplications of Ifi200 genes and suggest to be relevant for the development of obesity. The results further demonstrate the complexity of the disease and highlight the importance for studying rare genetic variants as they can be causal for large effects.

Diversity and regulatory impact of copy number variation in the primate Macaca fascicularis

BACKGROUND

Copy number variations (CNVs) are a significant source of genetic diversity and commonly found in mammalian genomes. We have generated a genome-wide CNV map for Cynomolgus monkeys (Macaca fascicularis). This crab-eating macaque is the closest animal model to humans that is used in biomedical research.

RESULTS

We show that Cynomolgus monkey CNVs are in general much smaller in size than gene loci and are specific to the population of origin. Genome-wide expression data from five vitally important organs demonstrates that CNVs in close proximity to transcription start sites associate strongly with expression changes. Among these eQTL genes we find an overrepresentation of genes involved in metabolism, receptor activity, and transcription.

CONCLUSION

These results provide evidence that CNVs shape tissue transcriptomes in monkey populations, potentially offering an adaptive advantage. We suggest that this genetic diversity should be taken into account when using Cynomolgus macaques as models.

Methylglyoxal-Glyoxalase 1 Balance: The Root of Vascular Damage

The highly reactive dicarbonyl methylglyoxal (MGO) is mainly formed as byproduct of glycolysis. Therefore, high blood glucose levels determine increased MGO accumulation. Nonetheless, MGO levels are also increased as consequence of the ineffective action of its main detoxification pathway, the glyoxalase system, of which glyoxalase 1 (Glo1) is the rate-limiting enzyme. Indeed, a physiological decrease of Glo1 transcription and activity occurs not only in chronic hyperglycaemia but also with ageing, during which MGO accumulation occurs. MGO and its advanced glycated end products (AGEs) are associated with age-related diseases including diabetes, vascular dysfunction and neurodegeneration. Endothelial dysfunction is the first step in the initiation, progression and clinical outcome of vascular complications, such as retinopathy, nephropathy, impaired wound healing and macroangiopathy. Because of these considerations, studies have been centered on understanding the molecular basis of endothelial dysfunction in diabetes, unveiling a central role of MGO-Glo1 imbalance in the onset of vascular complications. This review focuses on the current understanding of MGO accumulation and Glo1 activity in diabetes, and their contribution on the impairment of endothelial function leading to diabetes-associated vascular damage.

Identifying anti-cancer drug response related genes using an integrative analysis of transcriptomic and genomic variations with cell line-based drug perturbations

Background

Clinical responses to anti-cancer therapies often only benefit a defined subset of patients. Predicting the best treatment strategy hinges on our ability to effectively translate genomic data into actionable information on drug responses.

Results

To achieve this goal, we compiled a comprehensive collection of baseline cancer genome data and drug response information derived from a large panel of cancer cell lines. This data set was applied to identify the signature genes relevant to drug sensitivity and their resistance by integrating CNVs and the gene expression of cell lines with in vitro drug responses. We presented an efficient in-silico pipeline for integrating heterogeneous cell line data sources with the simultaneous modeling of drug response values across all the drugs and cell lines. Potential signature genes correlated with drug response (sensitive or resistant) in different cancer types were identified. Using signature genes, our collaborative filtering-based drug response prediction model outperformed the 44 algorithms submitted to the DREAM competition on breast cancer cells. The functions of the identified drug response related signature genes were carefully analyzed at the pathway level and the synthetic lethality level. Furthermore, we validated these signature genes by applying them to the classification of the different subtypes of the TCGA tumor samples, and further uncovered their in vivo implications using clinical patient data.

Conclusions

Our work may have promise in translating genomic data into customized marker genes relevant to the response of specific drugs for a specific cancer type of individual patients.

Genome-Wide Deletion Screening with the Array CGH Method in Mouse Offspring Derived from Irradiated Spermatogonia Indicates that Mutagenic Responses are Highly Variable among Genes

Until the end of the 20th century, mouse germ cell data on induced mutation rates, which were collected using classical genetic methods at preselected specific loci, provided the principal basis for estimates of genetic risks from radiation in humans. The work reported on here is an extension of earlier efforts in this area using molecular methods. It focuses on validating the use of array comparative genomic hybridization (array CGH) methods for identifying radiation-induced copy number variants (CNVs) and specifically for DNA deletions. The emphasis on deletions stems from the view that it constitutes the predominant type of radiation-induced genetic damage, which is relevant for estimating genetic risks in humans. In the current study, deletion mutations were screened in the genomes of F1 mice born to unirradiated or 4 Gy irradiated sires at the spermatogonia stage (100 offspring each). The array CGH analysis was performed using a "2M array" with over 2 million probes with a mean interprobe distance of approximately 1 kb. The results provide evidence of five molecularly-confirmed paternally-derived deletions in the irradiated group (5/100) and one in the controls (1/100). These data support a calculation, which estimates that the mutation rate is 1 × 10^-2/Gy per genome for induced deletions; this is much lower than would be expected if one assumes that the specific locus rate of 1 × 10^-5/locus per Gy (at 34 loci) is applicable to other genes in the genome. The low observed rate of induced deletions suggests that the effective number of genes/genomic regions at which recoverable deletions could be induced would be only approximately 1,000. This estimate is far lower than expected from the size of the mouse genome (>20,000 genes). Such a discrepancy between observation and expectation can occur if the genome contains numerous genes that are far less sensitive to radiation-induced deletions, if many deletion-bearing offspring are not viable or if the current method is substandard for detecting small deletions.

Suppression of cytochrome P450 4B1: An early event in adrenocortical tumorigenesis

BACKGROUND

Adrenocortical carcinoma is a rare neoplasm with a poor prognosis. Conversely, adrenocortical adenomas are common and benign. Despite their shared histologic origin, little evidence exists to suggest that adrenocortical adenoma arises from adrenocortical carcinoma. Recent genetic analyses of adrenocortical carcinoma have shown recurrent gene copy deletion of CYP4B1, a cytochrome P450 isozyme. This study investigates a potential role for CYP4B1 in modulating adrenocortical tumorigenesis and/or conferring chemoresistance to adrenocortical carcinomas.

METHODS

Using TaqMan, real-time quantitative polymerase chain reaction techniques, we investigated CYP4B1 expression in normal adrenal cortex (n = 10), histologically confirmed adrenocortical adenomas (n = 10), and adrenocortical carcinomas (n = 10). Adrenocortical carcinoma cell lines were enforced to express CYP4B1, and effects on cell death and enhanced mitotane and cisplatin sensitivity were tested.

RESULTS

Gene expression analyses demonstrated suppression of CYP4B1 in 100% of both the adrenocortical adenomas (10/10) and adrenocortical carcinomas (10/10) tested. Average relative expression of CYP4B1 was decreased at 0.19 (0.01-0.50; P < .01) in adrenocortical adenomas and nearly absent in adrenocortical carcinomas (0.01; 0.00-0.05; P < .01). Protein expression correlated with mRNA expression. Ectopic expression of CYP4B1 promoted cytotoxicity and increased chemosensitivity in adrenocortical carcinoma cell lines.

CONCLUSION

CYP4B1 is silenced in both benign and malignant adrenocortical tumors and may contribute to tumorigenesis and chemoresistance. Sensitization of adrenocortical carcinoma cells engineered to overexpress CYP4B1 further supports this notion.

Reappraisal of putative glyoxalase 1-deficient mouse and dicarbonyl stress on embryonic stem cells in vitro

Glyoxalase 1 (Glo1) is a cytoplasmic enzyme with a cytoprotective function linked to metabolism of the cytotoxic side product of glycolysis, methylglyoxal (MG). It prevents dicarbonyl stress - the abnormal accumulation of reactive dicarbonyl metabolites, increasing protein and DNA damage. Increased Glo1 expression delays ageing and suppresses carcinogenesis, insulin resistance, cardiovascular disease and vascular complications of diabetes and renal failure. Surprisingly, gene trapping by the International Mouse Knockout Consortium (IMKC) to generate putative Glo1 knockout mice produced a mouse line with the phenotype characterised as normal and healthy. Here, we show that gene trapping mutation was successful, but the presence of Glo1 gene duplication, probably in the embryonic stem cells (ESCs) before gene trapping, maintained wild-type levels of Glo1 expression and activity and sustained the healthy phenotype. In further investigation of the consequences of dicarbonyl stress in ESCs, we found that prolonged exposure of mouse ESCs in culture to high concentrations of MG and/or hypoxia led to low-level increase in Glo1 copy number. In clinical translation, we found a high prevalence of low-level GLO1 copy number increase in renal failure where there is severe dicarbonyl stress. In conclusion, the IMKC Glo1 mutant mouse is not deficient in Glo1 expression through duplication of the Glo1 wild-type allele. Dicarbonyl stress and/or hypoxia induces low-level copy number alternation in ESCs. Similar processes may drive rare GLO1 duplication in health and disease.

Dicarbonyls and glyoxalase in disease mechanisms and clinical therapeutics

The reactive dicarbonyl metabolite methylglyoxal (MG) is the precursor of the major quantitative advanced glycation endproducts (AGEs) in physiological systems - arginine-derived hydroimidazolones and deoxyguanosine-derived imidazopurinones. The glyoxalase system in the cytoplasm of cells provides the primary defence against dicarbonyl glycation by catalysing the metabolism of MG and related reactive dicarbonyls. Dicarbonyl stress is the abnormal accumulation of dicarbonyl metabolites leading to increased protein and DNA modification contributing to cell and tissue dysfunction in ageing and disease. It is produced endogenously by increased formation and/or decreased metabolism of dicarbonyl metabolites. Dicarbonyl stress contributes to ageing, disease and activity of cytotoxic chemotherapeutic agents. It contributes to ageing through age-related decline in glyoxalase 1 (Glo-1) activity. Glo-1 has a dual role in cancer as a tumour suppressor protein prior to tumour development and mediator of multi-drug resistance in cancer treatment, implicating dicarbonyl glycation of DNA in carcinogenesis and dicarbonyl-driven cytotoxicity in mechanism of action of anticancer drugs. Glo-1 is a driver of cardiovascular disease, likely through dicarbonyl stress-driven dyslipidemia and vascular cell dysfunction. Dicarbonyl stress is also a contributing mediator of obesity and vascular complications of diabetes. There are also emerging roles in neurological disorders. Glo-1 responds to dicarbonyl stress to enhance cytoprotection at the transcriptional level through stress-responsive increase of Glo-1 expression. Small molecule Glo-1 inducers are in clinical development for improved metabolic, vascular and renal health and Glo-1 inhibitors in preclinical development for multidrug resistant cancer chemotherapy.

Genome-Wide Detection of CNVs and Their Association with Meat Tenderness in Nelore Cattle

Brazil is one of the largest beef producers and exporters in the world with the Nelore breed representing the vast majority of Brazilian cattle (Bos taurus indicus). Despite the great adaptability of the Nelore breed to tropical climate, meat tenderness (MT) remains to be improved. Several factors including genetic composition can influence MT. In this article, we report a genome-wide analysis of copy number variation (CNV) inferred from Illumina® High Density SNP-chip data for a Nelore population of 723 males. We detected >2,600 CNV regions (CNVRs) representing ≈6.5% of the genome. Comparing our results with previous studies revealed an overlap in ≈1400 CNVRs (>50%). A total of 1,155 CNVRs (43.6%) overlapped 2,750 genes. They were enriched for processes involving guanosine triphosphate (GTP), previously reported to influence skeletal muscle physiology and morphology. Nelore CNVRs also overlapped QTLs for MT reported in other breeds (8.9%, 236 CNVRs) and from a previous study with this population (4.1%, 109 CNVRs). Two CNVRs were also proximal to glutathione metabolism genes that were previously associated with MT. Genome-wide association study of CN state with estimated breeding values derived from meat shear force identified 6 regions, including a region on BTA3 that contains genes of the cAMP and cGMP pathway. Ten CNVRs that overlapped regions associated with MT were successfully validated by qPCR. Our results represent the first comprehensive CNV study in Bos taurus indicus cattle and identify regions in which copy number changes are potentially of importance for the MT phenotype.

Mitochondrial Dysfunction may explain symptom variation in Phelan-McDermid Syndrome

Phelan-McDermid Syndrome (PMS), which is defined by a deletion within 22q13, demonstrates significant phenotypic variation. Given that six mitochondrial genes are located within 22q13, including complex I and IV genes, we hypothesize that mitochondrial complex activity abnormalities may explain phenotypic variation in PMS symptoms. Complex I, II, II + III and IV activity was measured in 51 PMS participants. Caretakers completed questionnaires and provided genetic information through the PMS foundation registry. Complex activity was abnormal in 59% of PMS participants. Abnormalities were found in complex I and IV but not complex II + III and II activity, consistent with disruption of genes within the 22q13 region. However, complex activity abnormalities were not related to specific gene deletions suggesting a "neighboring effect" of regional deletions on adjacent gene expression. A specific combination of symptoms (autism spectrum disorder, developmental regression, failure-to-thrive, exercise intolerance/fatigue) was associated with complex activity abnormalities. 64% of 106 individuals in the PMS foundation registry who did not have complex activity measured also endorsed this pattern of symptoms. These data suggest that mitochondrial abnormalities, specifically abnormalities in complex I and IV activity, may explain some phenotypic variation in PMS individuals. These results point to novel pathophysiology mechanisms and treatment targets for PMS patients.

The Mouse Universal Genotyping Array: From Substrains to Subspecies

Genotyping microarrays are an important resource for genetic mapping, population genetics, and monitoring of the genetic integrity of laboratory stocks. We have developed the third generation of the Mouse Universal Genotyping Array (MUGA) series, GigaMUGA, a 143,259-probe Illumina Infinium II array for the house mouse (Mus musculus). The bulk of the content of GigaMUGA is optimized for genetic mapping in the Collaborative Cross and Diversity Outbred populations, and for substrain-level identification of laboratory mice. In addition to 141,090 single nucleotide polymorphism probes, GigaMUGA contains 2006 probes for copy number concentrated in structurally polymorphic regions of the mouse genome. The performance of the array is characterized in a set of 500 high-quality reference samples spanning laboratory inbred strains, recombinant inbred lines, outbred stocks, and wild-caught mice. GigaMUGA is highly informative across a wide range of genetically diverse samples, from laboratory substrains to other Mus species. In addition to describing the content and performance of the array, we provide detailed probe-level annotation and recommendations for quality control.

Genome Wide Distributions and Functional Characterization of Copy Number Variations between Chinese and Western Pigs

Copy number variations (CNVs) refer to large insertions, deletions and duplications in the genomic structure ranging from one thousand to several million bases in size. Since the development of next generation sequencing technology, several methods have been well built for detection of copy number variations with high credibility and accuracy. Evidence has shown that CNV occurring in gene region could lead to phenotypic changes due to the alteration in gene structure and dosage. However, it still remains unexplored whether CNVs underlie the phenotypic differences between Chinese and Western domestic pigs. Based on the read-depth methods, we investigated copy number variations using 49 individuals derived from both Chinese and Western pig breeds. A total of 3,131 copy number variation regions (CNVRs) were identified with an average size of 13.4 Kb in all individuals during domestication, harboring 1,363 genes. Among them, 129 and 147 CNVRs were Chinese and Western pig specific, respectively. Gene functional enrichments revealed that these CNVRs contribute to strong disease resistance and high prolificacy in Chinese domestic pigs, but strong muscle tissue development in Western domestic pigs. This finding is strongly consistent with the morphologic characteristics of Chinese and Western pigs, indicating that these group-specific CNVRs might have been preserved by artificial selection for the favored phenotypes during independent domestication of Chinese and Western pigs. In this study, we built high-resolution CNV maps in several domestic pig breeds and discovered the group specific CNVs by comparing Chinese and Western pigs, which could provide new insight into genomic variations during pigs’ independent domestication, and facilitate further functional studies of CNV-associated genes.

Genomic copy number variation in Mus musculus

BACKGROUND

Copy number variation is an important dimension of genetic diversity and has implications in development and disease. As an important model organism, the mouse is a prime candidate for copy number variant (CNV) characterization, but this has yet to be completed for a large sample size. Here we report CNV analysis of publicly available, high-density microarray data files for 351 mouse tail samples, including 290 mice that had not been characterized for CNVs previously.

RESULTS

We found 9634 putative autosomal CNVs across the samples affecting 6.87% of the mouse reference genome. We find significant differences in the degree of CNV uniqueness (single sample occurrence) and the nature of CNV-gene overlap between wild-caught mice and classical laboratory strains. CNV-gene overlap was associated with lipid metabolism, pheromone response and olfaction compared to immunity, carbohydrate metabolism and amino-acid metabolism for wild-caught mice and classical laboratory strains, respectively. Using two subspecies of wild-caught Mus musculus, we identified putative CNVs unique to those subspecies and show this diversity is better captured by wild-derived laboratory strains than by the classical laboratory strains. A total of 9 genic copy number variable regions (CNVRs) were selected for experimental confirmation by droplet digital PCR (ddPCR).

CONCLUSION

The analysis we present is a comprehensive, genome-wide analysis of CNVs in Mus musculus, which increases the number of known variants in the species and will accelerate the identification of novel variants in future studies.

Activity, regulation, copy number and function in the glyoxalase system

Molecular, catalytic and structural properties of glyoxalase pathway enzymes of many species are now known. Current research has focused on the regulation of activity and expression of Glo1 (glyoxalase I) and Glo2 (glyoxalase II) and their role in health and disease. Human GLO1 has MRE (metal-response element), IRE (insulin-response element), E2F4 (early gene 2 factor isoform 4), AP-2α (activating enhancer-binding protein 2α) and ARE (antioxidant response-element) regulatory elements and is a hotspot for copy number variation. The human Glo2 gene, HAGH (hydroxyacylglutathione hydrolase), has a regulatory p53-response element. Glo1 is linked to healthy aging, obesity, diabetes and diabetic complications, chronic renal disease, cardiovascular disease, other disorders and multidrug resistance in cancer chemotherapy. Mathematical modelling of the glyoxalase pathway predicts that pharmacological levels of increased Glo1 activity markedly decrease cellular methylglyoxal and related glycation, and pharmacological Glo1 inhibition markedly increases cellular methylglyoxal and related glycation. Glo1 inducers are in development to sustain healthy aging and for treatment of vascular complications of diabetes and other disorders, and cell-permeant Glo1 inhibitors are in development for treatment of multidrug-resistant tumours, malaria and potentially pathogenic bacteria and fungi.

Connecting Anxiety and Genomic Copy Number Variation: A Genome-Wide Analysis in CD-1 Mice

Genomic copy number variants (CNVs) have been implicated in multiple psychiatric disorders, but not much is known about their influence on anxiety disorders specifically. Using next-generation sequencing (NGS) and two additional array-based genotyping approaches, we detected CNVs in a mouse model consisting of two inbred mouse lines showing high (HAB) and low (LAB) anxiety-related behavior, respectively. An influence of CNVs on gene expression in the central (CeA) and basolateral (BLA) amygdala, paraventricular nucleus (PVN), and cingulate cortex (Cg) was shown by a two-proportion Z-test (p = 1.6 x 10^-31), with a positive correlation in the CeA (p = 0.0062), PVN (p = 0.0046) and Cg (p = 0.0114), indicating a contribution of CNVs to the genetic predisposition to trait anxiety in the specific context of HAB/LAB mice. In order to confirm anxiety-relevant CNVs and corresponding genes in a second mouse model, we further examined CD-1 outbred mice. We revealed the distribution of CNVs by genotyping 64 CD 1 individuals using a high-density genotyping array (Jackson Laboratory). 78 genes within those CNVs were identified to show nominally significant association (48 genes), or a statistical trend in their association (30 genes) with the time animals spent on the open arms of the elevated plus-maze (EPM). Fifteen of them were considered promising candidate genes of anxiety-related behavior as we could show a significant overlap (permutation test, p = 0.0051) with genes within HAB/LAB CNVs. Thus, here we provide what is to our knowledge the first extensive catalogue of CNVs in CD-1 mice and potential corresponding candidate genes linked to anxiety-related behavior in mice.

Adaptive potential of genomic structural variation in human and mammalian evolution

Because phenotypic innovations must be genetically heritable for biological evolution to proceed, it is natural to consider new mutation events as well as standing genetic variation as sources for their birth. Previous research has identified a number of single-nucleotide polymorphisms that underlie a subset of adaptive traits in organisms. However, another well-known class of variation, genomic structural variation, could have even greater potential to produce adaptive phenotypes, due to the variety of possible types of alterations (deletions, insertions, duplications, among others) at different genomic positions and with variable lengths. It is from these dramatic genomic alterations, and selection on their phenotypic consequences, that adaptations leading to biological diversification could be derived. In this review, using studies in humans and other mammals, we highlight examples of how phenotypic variation from structural variants might become adaptive in populations and potentially enable biological diversification. Phenotypic change arising from structural variants will be described according to their immediate effect on organismal metabolic processes, immunological response and physical features. Study of population dynamics of segregating structural variation can therefore provide a window into understanding current and historical biological diversification.

CNV instability associated with DNA replication dynamics: evidence for replicative mechanisms in CNV mutagenesis

Copy number variation (CNV) in the human genome is of vital importance to human health and evolution of our species. However, much of the molecular basis of CNV mutagenesis remains to be elucidated. Considering the DNA replication model of 'fork stalling and template switching' for CNV formation, we hypothesized that replication fork progression could be important for CNV mutagenesis. However, molecular assays of replication fork progression at the genome level are technically challenging. Instead, we conducted an estimation of DNA replication dynamics, as the statistic R, using the readily available data of replication timing. Small R-values can reflect 'slowed' replication, which could result from less fork initiation, reduced fork speed or fork barriers. We generated genome-wide profiles of R in the genomes of human, mouse and Drosophila. Intriguingly, the CNV breakpoints in all three genomes showed significantly biased distributions toward the genomic regions with small R-values, suggesting potential replication stress-induced CNV instability. Notably, among the human CNVs with distinct breakpoint junction characteristics, the homology-mediated and VNTR-mediated CNVs contribute the most to the correlation between CNV instability and the statistic R, consistent with the recent findings in the C. elegans and yeast genomes of repeat-induced DNA replication error and consequent CNV formation. The statistic R may reflect both replication stress and the effect of local genome architecture on fork progression. Our concordant observations suggest an important role for DNA replicative mechanisms in CNV mutagenesis and genome instability.

Identification of structural variation in mouse genomes

Structural variation is variation in structure of DNA regions affecting DNA sequence length and/or orientation. It generally includes deletions, insertions, copy-number gains, inversions, and transposable elements. Traditionally, the identification of structural variation in genomes has been challenging. However, with the recent advances in high-throughput DNA sequencing and paired-end mapping (PEM) methods, the ability to identify structural variation and their respective association to human diseases has improved considerably. In this review, we describe our current knowledge of structural variation in the mouse, one of the prime model systems for studying human diseases and mammalian biology. We further present the evolutionary implications of structural variation on transposable elements. We conclude with future directions on the study of structural variation in mouse genomes that will increase our understanding of molecular architecture and functional consequences of structural variation.

A genome-wide association study on copy-number variation identifies a 11q11 loss as a candidate susceptibility variant for colorectal cancer

Colorectal cancer (CRC) is a complex disease, and therefore its development is determined by the combination of both environmental factors and genetic variants. Although genome-wide association studies (GWAS) of SNP variation have conveniently identified 20 genetic variants so far, a significant proportion of the observed heritability is yet to be explained. Common copy-number variants (CNVs) are one of the most important genomic sources of variability, and hence a potential source to explain part of this missing genetic fraction. Therefore, we have performed a GWAS on CNVs to explore the relationship between common structural variation and CRC development. Phase 1 of the GWAS consisted of 881 cases and 667 controls from a Spanish cohort. Copy-number status was validated by quantitative PCR for each of those common CNVs potentially associated with CRC in phase I. Subsequently, SNPs were chosen as proxies for the validated CNVs for phase II replication (1,342 Spanish cases and 1,874 Spanish controls). Four common CNVs were found to be associated with CRC and were further replicated in Phase II. Finally, we found that SNP rs1944682, tagging a 11q11 CNV, was nominally associated with CRC susceptibility (p value = 0.039; OR = 1.122). This locus has been previously related to extreme obesity phenotypes, which could suggest a relationship between body weight and CRC susceptibility.

Increased genomic integrity of an improved protein-based mouse induced pluripotent stem cell method compared with current viral-induced strategies

It has recently been shown that genomic integrity (with respect to copy number variants [CNVs]) is compromised in human induced pluripotent stem cells (iPSCs) generated by viral-based ectopic expression of specific transcription factors (e.g., Oct4, Sox2, Klf4, and c-Myc). However, it is unclear how different methods for iPSC generation compare with one another with respect to CNV formation. Because array-based methods remain the gold standard for detecting unbalanced structural variants (i.e., CNVs), we have used this approach to comprehensively identify CNVs in iPSC as a proxy for determining whether our modified protein-based method minimizes genomic instability compared with retro- and lentiviral methods. In this study, we established an improved method for protein reprogramming by using partially purified reprogramming proteins, resulting in more efficient generation of iPSCs from C57/BL6J mouse hepatocytes than using protein extracts. We also developed a robust and unbiased 1 M custom array CGH platform to identify novel CNVs and previously described hot spots for CNV formation, allowing us to detect CNVs down to the size of 1.9 kb. The genomic integrity of these protein-based mouse iPSCs (p-miPSCs) was compared with miPSCs developed from viral-based strategies (i.e., retroviral: retro-miPSCs or lentiviral: lenti-miPSCs). We identified an increased CNV content in lenti-miPSCs and retro-miPSCs (29∼53 CNVs) compared with p-miPSCs (9∼10 CNVs), indicating that our improved protein-based reprogramming method maintains genomic integrity better than current viral reprogramming methods. Thus, our study, for the first time to our knowledge, demonstrates that reprogramming methods significantly influence the genomic integrity of resulting iPSCs.

Copy number variations among silkworms

Background

Copy number variations (CNVs), which are important source for genetic and phenotypic variation, have been shown to be associated with disease as well as important QTLs, especially in domesticated animals. However, little is known about the CNVs in silkworm.

Results

In this study, we have constructed the first CNVs map based on genome-wide analysis of CNVs in domesticated silkworm. Using next-generation sequencing as well as quantitative PCR (qPCR), we identified ~319 CNVs in total and almost half of them (~ 49%) were distributed on uncharacterized chromosome. The CNVs covered 10.8 Mb, which is about 2.3% of the entire silkworm genome. Furthermore, approximately 61% of CNVs directly overlapped with SDs in silkworm. The genes in CNVs are mainly related to reproduction, immunity, detoxification and signal recognition, which is consistent with the observations in mammals.

Conclusions

An initial CNVs map for silkworm has been described in this study. And this map provides new information for genetic variations in silkworm. Furthermore, the silkworm CNVs may play important roles in reproduction, immunity, detoxification and signal recognition. This study provided insight into the evolution of the silkworm genome and an invaluable resource for insect genomics research.

Copy number variation of glyoxalase I

The glyoxalase I gene GLO1 is a hotspot for copy number variation in the human and mouse genomes. The additional copies are often functional, giving rise to 2–4-fold increased glyoxalase I expression and activity. The prevalence of GLO1 copy number increase in the human population appears to be approximately 2% and may be linked to a risk of obesity, diabetes and aging. Increased GLO1 copy number has been found in human tumour cell lines and primary human tumours. The minimum common copy number increase region was approximately 1 Mb and it contained GLO1 and seven other genes. The increased copy number was generally functional, being associated with increased glyoxalase I protein and multidrug resistance in cancer chemotherapy. Glo1 duplication in the mouse genome is found within approximately 0.5 Mb of duplicated DNA. It was claimed to be linked to anxiety phenotypes, but other related discordant findings have doubted the association with glyoxalase I and further investigation is required.

Incorporating inter-relationships between different levels of genomic data into cancer clinical outcome prediction

In order to improve our understanding of cancer and develop multi-layered theoretical models for the underlying mechanism, it is essential to have enhanced understanding of the interactions between multiple levels of genomic data that contribute to tumor formation and progression. Although there exist recent approaches such as a graph-based framework that integrates multi-omics data including copy number alteration, methylation, gene expression, and miRNA data for cancer clinical outcome prediction, most of previous methods treat each genomic data as independent and the possible interplay between them is not explicitly incorporated to the model. However, cancer is dysregulated by multiple levels in the biological system through genomic, epigenomic, transcriptomic, and proteomic level. Thus, genomic features are likely to interact with other genomic features in the different genomic levels. In order to deepen our knowledge, it would be desirable to incorporate such inter-relationship information when integrating multi-omics data for cancer clinical outcome prediction. In this study, we propose a new graph-based framework that integrates not only multi-omics data but inter-relationship between them for better elucidating cancer clinical outcomes. In order to highlight the validity of the proposed framework, serous cystadenocarcinoma data from TCGA was adopted as a pilot task. The proposed model incorporating inter-relationship between different genomic features showed significantly improved performance compared to the model that does not consider inter-relationship when integrating multi-omics data. For the pair between miRNA and gene expression data, the model integrating miRNA, for example, gene expression, and inter-relationship between them with an AUC of 0.8476 (REI) outperformed the model combining miRNA and gene expression data with an AUC of 0.8404. Similar results were also obtained for other pairs between different levels of genomic data. Integration of different levels of data and inter-relationship between them can aid in extracting new biological knowledge by drawing an integrative conclusion from many pieces of information collected from diverse types of genomic data, eventually leading to more effective screening strategies and alternative therapies that may improve outcomes.

Relative impact of multi-layered genomic data on gene expression phenotypes in serous ovarian tumors

Background

The emerging multi-layers of genomic data have provided unprecedented opportunities for cancer research, especially for the association study between gene expressions and other types of genomic features. No previous approaches, however, provide an adequate statistical framework for or global analysis on the relative impact of different genomic feature layers to gene expression phenotypes.

Methods

We propose an integrative statistical framework based on a sparse regression to model the impact of multi-layered genomic features on gene expression traits. The proposed approach can be regarded as an integrative expression Quantitative Traits Loci approach in which not only the genetic variations of SNPs or copy number variations but also other features in both genomic and epigenomic levels are used to explain the expression of genes. To highlight the validity of the proposed approach, the TCGA ovarian cancer dataset was analysed as a pilot task.

Results

The analysis shows that our integrative approach has consistently superior power in predicting gene expression levels compared to that from each single data type-based analysis. Moreover, the proposed method has the advantage of producing a substantially reduced number of spurious associations. We provide an interesting characterization of genes in terms of its genomic association patterns. Important genomic features reported in previous ovarian cancer research are successfully identified as major hubs in the resulting association network between heterogeneous types of genomic features and genes.

Conclusions

In this paper, we model the gene expression phenotypes with respect to multiple different types of genomic data in an integrative framework. Our analysis reveals the global view on the relative contribution of different genomic feature types to gene expression phenotypes in ovarian cancer.

Structural variation-associated expression changes are paralleled by chromatin architecture modifications

Copy number variants (CNVs) influence the expression of genes that map not only within the rearrangement, but also to its flanks. To assess the possible mechanism(s) underlying this “neighboring effect”, we compared intrachromosomal interactions and histone modifications in cell lines of patients affected by genomic disorders and control individuals. Using chromosome conformation capture (4C-seq), we observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together. The newly identified interacting genes include AUTS2, mutations of which are associated with autism and intellectual disabilities. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We also pinpointed concomitant changes in histone modifications between samples.

We conclude that large genomic rearrangements can lead to chromatin conformation changes that extend far away from the structural variant, thereby possibly modulating expression globally and modifying the phenotype.

GEO Series accession number: GSE33784, GSE33867.

Identification of both copy number variation-type and constant-type core elements in a large segmental duplication region of the mouse genome

Background

Copy number variation (CNV), an important source of diversity in genomic structure, is frequently found in clusters called CNV regions (CNVRs). CNVRs are strongly associated with segmental duplications (SDs), but the composition of these complex repetitive structures remains unclear.

Results

We conducted self-comparative-plot analysis of all mouse chromosomes using the high-speed and large-scale-homology search algorithm SHEAP. For eight chromosomes, we identified various types of large SD as tartan-checked patterns within the self-comparative plots. A complex arrangement of diagonal split lines in the self-comparative-plots indicated the presence of large homologous repetitive sequences. We focused on one SD on chromosome 13 (SD13M), and developed SHEPHERD, a stepwise ab initio method, to extract longer repetitive elements and to characterize repetitive structures in this region. Analysis using SHEPHERD showed the existence of 60 core elements, which were expected to be the basic units that form SDs within the repetitive structure of SD13M. The demonstration that sequences homologous to the core elements (>70% homology) covered approximately 90% of the SD13M region indicated that our method can characterize the repetitive structure of SD13M effectively. Core elements were composed largely of fragmented repeats of a previously identified type, such as long interspersed nuclear elements (LINEs), together with partial genic regions. Comparative genome hybridization array analysis showed that whereas 42 core elements were components of CNVR that varied among mouse strains, 8 did not vary among strains (constant type), and the status of the others could not be determined. The CNV-type core elements contained significantly larger proportions of long terminal repeat (LTR) types of retrotransposon than the constant-type core elements, which had no CNV. The higher divergence rates observed in the CNV-type core elements than in the constant type indicate that the CNV-type core elements have a longer evolutionary history than constant-type core elements in SD13M.

Conclusions

Our methodology for the identification of repetitive core sequences simplifies characterization of the structures of large SDs and detailed analysis of CNV. The results of detailed structural and quantitative analyses in this study might help to elucidate the biological role of one of the SDs on chromosome 13.

Evolutionary dynamics of copy number variation in pig genomes in the context of adaptation and domestication

Background

Copy number variable regions (CNVRs) can result in drastic phenotypic differences and may therefore be subject to selection during domestication. Studying copy number variation in relation to domestication is highly relevant in pigs because of their very rich natural and domestication history that resulted in many different phenotypes. To investigate the evolutionary dynamic of CNVRs, we applied read depth method on next generation sequence data from 16 individuals, comprising wild boars and domestic pigs from Europe and Asia.

Results

We identified 3,118 CNVRs with an average size of 13 kilobases comprising a total of 39.2 megabases of the pig genome and 545 overlapping genes. Functional analyses revealed that CNVRs are enriched with genes related to sensory perception, neurological process and response to stimulus, suggesting their contribution to adaptation in the wild and behavioral changes during domestication. Variations of copy number (CN) of antimicrobial related genes suggest an ongoing process of evolution of these genes to combat food-borne pathogens. Likewise, some genes related to the omnivorous lifestyle of pigs, like genes involved in detoxification, were observed to be CN variable. A small portion of CNVRs was unique to domestic pigs and may have been selected during domestication. The majority of CNVRs, however, is shared between wild and domesticated individuals, indicating that domestication had minor effect on the overall diversity of CNVRs. Also, the excess of CNVRs in non-genic regions implies that a major part of these variations is likely to be (nearly) neutral. Comparison between different populations showed that larger populations have more CNVRs, highlighting that CNVRs are, like other genetic variation such as SNPs and microsatellites, reflecting demographic history rather than phenotypic diversity.

Conclusion

CNVRs in pigs are enriched for genes related to sensory perception, neurological process, and response to stimulus. The majority of CNVRs ascertained in domestic pigs are also variable in wild boars, suggesting that the domestication of the pig did not result in a change in CNVRs in domesticated pigs. The majority of variable regions were found to reflect demographic patterns rather than phenotypic.

Copy number variants in locally raised Chinese chicken genomes determined using array comparative genomic hybridization

Background

Copy number variants contribute to genetic variation in birds. Analyses of copy number variants in chicken breeds had focused primarily on those from commercial varieties with nothing known about the occurrence and diversity of copy number variants in locally raised Chinese chicken breeds. To address this deficiency, we characterized copy number variants in 11 chicken breeds and compared the variation among these breeds.

Results

We presented a detailed analysis of the copy number variants in locally raised Chinese chicken breeds identified using a customized comparative genomic hybridization array. We identified 833 copy number variants contained within 308 copy number variant regions. The median and mean sizes of the copy number variant regions were 14.6 kb and 35.1 kb, respectively. Of the copy number variant regions, 138 (45%) involved gain of DNA, 159 (52%) involved loss of DNA, and 11 (3%) involved both gain and loss of DNA. Principal component analysis and agglomerative hierarchical clustering revealed the close relatedness of the four locally raised chicken breeds, Shek-Ki, Langshan, Qingyuan partridge, and Wenchang. Biological process enrichment analysis of the copy number variant regions confirmed the greater variation among the four aforementioned varieties than among the seven other breeds studied.

Conclusion

Our description of the distribution of the copy number variants and comparison of the differences among the copy number variant regions of the 11 chicken breeds supplemented the information available concerning the copy number variants of other Chinese chicken breeds. In addition to its relevance for functional analysis, our results provided the first insight into how chicken breeds can be clustered on the basis of their genomic copy number variation.

Phenotypic impact of genomic structural variation: insights from and for human disease

Genomic structural variants have long been implicated in phenotypic diversity and human disease, but dissecting the mechanisms by which they exert their functional impact has proven elusive. Recently however, developments in high-throughput DNA sequencing and chromosomal engineering technology have facilitated the analysis of structural variants in human populations and model systems in unprecedented detail. In this Review, we describe how structural variants can affect molecular and cellular processes, leading to complex organismal phenotypes, including human disease. We further present advances in delineating disease-causing elements that are affected by structural variants, and we discuss future directions for research on the functional consequences of structural variants.

Glyoxalase 1 enzyme activity in erythrocytes and Ala111Glu polymorphism in type 1-diabetes patients

AIMS/HYPOTHESIS

The enzyme glyoxalase 1 (GLO1) can inactivate the glycoxidation product methylglyoxal that is thought to be an important contributor to the pathogenesis of vascular complications in diabetes. We aimed to study erythrocyte GLO1 activity and whether the Ala111Glu GLO1 gene polymorphism affected GLO1 activity.

METHODS

Fasting erythrocyte GLO1 activity was measured spectrophotometrically. The A111G gene polymorphism, assessed in DNA from leucocytes was analyzed in patients with type 1-diabetes and normal kidney function and compared with a control group.

RESULTS

Sixty-one patients with type 1-diabetes duration of 26.1 (10.7) years, mean (SD) with a HbA1c of 7.8 (0.9)%, 61.7 (9.9) mmol/mol and normal glomerular filtration rate were compared with 62 age- and sex-matched healthy controls. GLO1 activity was 0.206 (0.183-0.231) median (25-75% percentiles) U/mg Hb in the control group vs. 0.192 (0.165-0.224) in the diabetes group, (p = 0.149). In the diabetes group GLO1 correlated with HbA1c (r = 0.33, p < 0.01) and oxidized glutathione (GSSG) (r = - 0.34, p < 0.01) and in the control group with GSH (r = 0.37, p < 0.005) and fasting glucose (r = 0.26, p < 0.04). In a multiple regression analysis with GLO1 activity as the dependent variable, including the Ala111Glu polymorphism, the significant independent variables were log GSSG (ß - 0.318, p = 0.02) and HbA1c (ß 0.285, p = 0.041) in the diabetes group and log GSH, (ß 0.407, p = 0.004) in the control group.

CONCLUSIONS/INTERPRETATION

Erythrocyte glyoxalase 1 activity did not differ between patients with type 1-diabetes and controls. The Ala111Glu glyoxalase gene polymorphism did not have an effect on glyoxalase 1 activity in either group.