Copy number variants, diseases and gene expression

regNet: an R package for network-based propagation of gene expression alterations

SUMMARY

Gene expression alterations and potentially underlying gene copy number mutations can be measured routinely in the wet lab, but it is still extremely challenging to quantify impacts of altered genes on clinically relevant characteristics to predict putative driver genes. We developed the R package regNet that utilizes gene expression and copy number data to learn regulatory networks for the quantification of potential impacts of individual gene expression alterations on user-defined target genes via network propagation. We demonstrate the value of regNet by identifying putative major regulators that distinguish pilocytic from diffuse astrocytomas and by predicting putative impacts of glioblastoma-specific gene copy number alterations on cell cycle pathway genes and patient survival.

AVAILABILITY AND IMPLEMENTATION

regNet is available for download at https://github.com/seifemi/regNet under GNU GPL-3.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genomic landscape of copy number variation and copy neutral loss of heterozygosity events in equine sarcoids reveals increased instability of the sarcoid genome

Although they are the most common neoplasms in equids, sarcoids are not fully characterized at the molecular level. Therefore, the objective of this study was to characterize the landscape of structural rearrangements, such as copy number variation (CNV) and copy neutral loss of heterozygosity (cnLOH), in the genomes of sarcoid tumor cells. This information will not only broaden our understanding of the characteristics of this genome but will also improve the general knowledge of this tumor and the mechanisms involved in its generation. To this end, Equine SNP64K Illumina microarrays were applied along with bioinformatics tools dedicated for signal intensity analysis. The analysis revealed increased instability of the genome of sarcoid cells compared with unaltered skin tissue samples, which was manifested by the prevalence of CNV and cnLOH events. Many of the identified CNVs overlapped with the other research results, but the simultaneously observed variability in the number and sizes of detected aberrations indicated a need for further studies and the development of more reliable bioinformatics algorithms. The functional analysis of genes co-localized with the identified aberrations revealed that these genes are engaged in vital cellular processes. In addition, a number of these genes directly contribute to neoplastic transformation. Furthermore, large numbers of cnLOH events identified in the sarcoids suggested that they may play no less significant roles than CNVs in the carcinogenesis of this tumor. Thus, our results indicate the importance of cnLOH and CNV in equine sarcoid oncogenesis and present a direction of future research.

Impact of Zygosity on Bimodal Phenotype Distributions

Allele number, or zygosity, is a clear determinant of gene expression in diploid cells. However, the relationship between the number of copies of a gene and its expression can be hard to anticipate, especially when the gene in question is embedded in a regulatory circuit that contains feedback. Here, we study this question making use of the natural genetic variability of human populations, which allows us to compare the expression profiles of a receptor protein in natural killer cells among donors infected with human cytomegalovirus with one or two copies of the allele. Crucially, the distribution of gene expression in many of the donors is bimodal, which indicates the presence of a positive feedback loop somewhere in the regulatory environment of the gene. Three separate gene-circuit models differing in the location of the positive feedback loop with respect to the gene can all reproduce the homozygous data. However, when the resulting fitted models are applied to the hemizygous donors, one model (the one with the positive feedback located at the level of gene transcription) is superior in describing the experimentally observed gene-expression profile. In that way, our work shows that zygosity can help us relate the structure and function of gene regulatory networks.

An enhancer haplotype may influence BCL11A expression levels and the response to hydroxyurea in β-thalassemia patients

AIM

To identify the BCL11A intron-2 enhancer linkage disequilibrium (LD) block, harboring two previously identified SNPs, associating with the hydroxyurea response in β-thalassemia patients and the functional significance of this region.

MATERIALS & METHODS

Several neighboring SNPs were genotyped in our cohort. The associating LD block was identified, and its function studied in K562 erythroid cells via CRISPR/Cas9 genome editing.

RESULTS

A haplotype harboring three tag SNPs correlated significantly with the HU-response and BCL11A transcript levels in the patients' reticulocytes. Two deletions encompassing this LD block significantly reduced BCL11A transcript levels in K562 cells.

CONCLUSION

Our data suggest an essential role for this LD block in BCL11A expression levels and the response to hydroxyurea in β-thalassemia patients.

Aneuploidy Causes Non-genetic Individuality

Phenotypic variability is a hallmark of diseases involving chromosome gains and losses, such as Down syndrome and cancer. Allelic variances have been thought to be the sole cause of this heterogeneity. Here, we systematically examine the consequences of gaining and losing single or multiple chromosomes to show that the aneuploid state causes non-genetic phenotypic variability. Yeast cell populations harboring the same defined aneuploidy exhibit heterogeneity in cell-cycle progression and response to environmental perturbations. Variability increases with degree of aneuploidy and is partly due to gene copy number imbalances, suggesting that subtle changes in gene expression impact the robustness of biological networks and cause alternate behaviors when they occur across many genes. As inbred trisomic mice also exhibit variable phenotypes, we further propose that non-genetic individuality is a universal characteristic of the aneuploid state that may contribute to variability in presentation and treatment responses of diseases caused by aneuploidy.

A novel pathway-based distance score enhances assessment of disease heterogeneity in gene expression

BACKGROUND

Distance based unsupervised clustering of gene expression data is commonly used to identify heterogeneity in biologic samples. However, high noise levels in gene expression data and relatively high correlation between genes are often encountered, so traditional distances such as Euclidean distance may not be effective at discriminating the biological differences between samples. An alternative method to examine disease phenotypes is to use pre-defined biological pathways. These pathways have been shown to be perturbed in different ways in different subjects who have similar clinical features. We hypothesize that differences in the expressions of genes in a given pathway are more predictive of differences in biological differences compared to standard approaches and if integrated into clustering analysis will enhance the robustness and accuracy of the clustering method. To examine this hypothesis, we developed a novel computational method to assess the biological differences between samples using gene expression data by assuming that ontologically defined biological pathways in biologically similar samples have similar behavior.

RESULTS

Pre-defined biological pathways were downloaded and genes in each pathway were used to cluster samples using the Gaussian mixture model. The clustering results across different pathways were then summarized to calculate the pathway-based distance score between samples. This method was applied to both simulated and real data sets and compared to the traditional Euclidean distance and another pathway-based clustering method, Pathifier. The results show that the pathway-based distance score performs significantly better than the Euclidean distance, especially when the heterogeneity is low and genes in the same pathways are correlated. Compared to Pathifier, we demonstrated that our approach achieves higher accuracy and robustness for small pathways. When the pathway size is large, by downsampling the pathways into smaller pathways, our approach was able to achieve comparable performance.

CONCLUSIONS

We have developed a novel distance score that represents the biological differences between samples using gene expression data and pre-defined biological pathway information. Application of this distance score results in more accurate, robust, and biologically meaningful clustering results in both simulated data and real data when compared to traditional methods. It also has comparable or better performance compared to Pathifier.

Genetics of β-Amyloid Precursor Protein in Alzheimer's Disease

Alzheimer's disease (AD) is characterized neuropathologically by neuronal cell loss, extracellular neuritic plaques composed of β-amyloid (Aβ), and intracellular neurofibrillary tangles composed of hyperphosphorylated tau protein. Aβ is generated by proteolytic processing of the β-amyloid precursor protein (APP). Most individuals with Down syndrome (DS) have three copies of APP, leading to elevated APP expression, increased Aβ deposition, and characteristic AD neuropathology. Sequencing of APP in familial early-onset AD identified missense mutations that cause AD, while a recently discovered coding variant, APP A673T, reduces the risk for AD. Cellular and animal studies show that risk-associated mutations increase total Aβ levels, Aβ42 levels, or Aβ fibrillogenesis, while protective alleles reduce Aβ levels. Together, these studies provide compelling evidence for the Aβ hypothesis and suggest that therapeutics that reduces Aβ levels or Aβ fibrillogenesis should lower the risk for or prevent AD.

Characterization of Copy Number Variation's Potential Role in Marek's Disease

Marek’s Disease (MD) is a highly contagious pathogenic and oncogenic disease primarily affecting chickens. Chicken Lines 63 and 72, as well as their recombinant congenic strains (RCS) with varied susceptibility to MD, are ideal models to study the complex mechanisms of genetic resistance to MD. In this study, we investigated copy number variation (CNV) in these inbred chicken lines using the Affymetrix Axiom HD 600 K SNP genotyping array. We detected 393 CNV segments across all ten chicken lines, of which 12 CNVs were specifically identified in Line 72. We then assessed genetic structure based on CNV and observed markedly different patterns. Finally, we validated two deletion events in Line 72 and correlated them with genes expression using qPCR and RNA-seq, respectively. Our combined results indicated that these two CNV deletions were likely to contribute to MD susceptibility.

Extensive Copy Number Variation in Fermentation-Related Genes Among Saccharomyces cerevisiae Wine Strains

Due to the importance of Saccharomyces cerevisiae in wine-making, the genomic variation of wine yeast strains has been extensively studied. One of the major insights stemming from these studies is that wine yeast strains harbor low levels of genetic diversity in the form of single nucleotide polymorphisms (SNPs). Genomic structural variants, such as copy number (CN) variants, are another major type of variation segregating in natural populations. To test whether genetic diversity in CN variation is also low across wine yeast strains, we examined genome-wide levels of CN variation in 132 whole-genome sequences of S. cerevisiae wine strains. We found an average of 97.8 CN variable regions (CNVRs) affecting ∼4% of the genome per strain. Using two different measures of CN diversity, we found that gene families involved in fermentation-related processes such as copper resistance (CUP), flocculation (FLO), and glucose metabolism (HXT), as well as the SNO gene family whose members are expressed before or during the diauxic shift, showed substantial CN diversity across the 132 strains examined. Importantly, these same gene families have been shown, through comparative transcriptomic and functional assays, to be associated with adaptation to the wine fermentation environment. Our results suggest that CN variation is a substantial contributor to the genomic diversity of wine yeast strains, and identify several candidate loci whose levels of CN variation may affect the adaptation and performance of wine yeast strains during fermentation.

CNV discovery for milk composition traits in dairy cattle using whole genome resequencing

Background

Copy number variations (CNVs) are important and widely distributed in the genome. CNV detection opens a new avenue for exploring genes associated with complex traits in humans, animals and plants. Herein, we present a genome-wide assessment of CNVs that are potentially associated with milk composition traits in dairy cattle.

Results

In this study, CNVs were detected based on whole genome re-sequencing data of eight Holstein bulls from four half- and/or full-sib families, with extremely high and low estimated breeding values (EBVs) of milk protein percentage and fat percentage. The range of coverage depth per individual was 8.2–11.9×. Using CNVnator, we identified a total of 14,821 CNVs, including 5025 duplications and 9796 deletions. Among them, 487 differential CNV regions (CNVRs) comprising ~8.23 Mb of the cattle genome were observed between the high and low groups. Annotation of these differential CNVRs were performed based on the cattle genome reference assembly (UMD3.1) and totally 235 functional genes were found within the CNVRs. By Gene Ontology and KEGG pathway analyses, we found that genes were significantly enriched for specific biological functions related to protein and lipid metabolism, insulin/IGF pathway-protein kinase B signaling cascade, prolactin signaling pathway and AMPK signaling pathways. These genes included INS, IGF2, FOXO3, TH, SCD5, GALNT18, GALNT16, ART3, SNCA and WNT7A, implying their potential association with milk protein and fat traits. In addition, 95 CNVRs were overlapped with 75 known QTLs that are associated with milk protein and fat traits of dairy cattle (Cattle QTLdb).

Conclusions

In conclusion, based on NGS of 8 Holstein bulls with extremely high and low EBVs for milk PP and FP, we identified a total of 14,821 CNVs, 487 differential CNVRs between groups, and 10 genes, which were suggested as promising candidate genes for milk protein and fat traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3636-3) contains supplementary material, which is available to authorized users.

Prediction of years of life after diagnosis of breast cancer using omics and omic-by-treatment interactions

Breast cancer (BC) is the second most common type of cancer and a major cause of death for women. Commonly, BC patients are assigned to risk groups based on the combination of prognostic and prediction factors (eg, patient age, tumor size, tumor grade, hormone receptor status, etc). Although this approach is able to identify risk groups with different prognosis, patients are highly heterogeneous in their response to treatments. To improve the prediction of BC patients, we extended clinical models (including prognostic and prediction factors with whole-omic data) to integrate omics profiles for gene expression and copy number variants (CNVs). We describe a modeling framework that is able to incorporate clinical risk factors, high-dimensional omics profiles, and interactions between omics and non-omic factors (eg, treatment). We used the proposed modeling framework and data from METABRIC (Molecular Taxonomy of Breast Cancer Consortium) to assess the impact on the accuracy of BC patient survival predictions when omics and omic-by-treatment interactions are being considered. Our analysis shows that omics and omic-by-treatment interactions explain a sizable fraction of the variance on survival time that is not explained by commonly used clinical covariates. The sizable interaction effects observed, together with the increase in prediction accuracy, suggest that whole-omic profiles could be used to improve prognosis prediction among BC patients.

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.

Copy-number variation of housekeeping gene rpl13a in rat strains selected for nervous system excitability

We evaluated copy number variation (CNV) for four genes in rat strains differing in nervous system excitability. rpl13a copy number is significantly reduced in hippocampus and bone marrow in rats with a high excitability threshold and stress. The observed phenomenon may be associated with a role for rpl13a in lipid metabolism.

iGC-an integrated analysis package of gene expression and copy number alteration

BACKGROUND

With the advancement in high-throughput technologies, researchers can simultaneously investigate gene expression and copy number alteration (CNA) data from individual patients at a lower cost. Traditional analysis methods analyze each type of data individually and integrate their results using Venn diagrams. Challenges arise, however, when the results are irreproducible and inconsistent across multiple platforms. To address these issues, one possible approach is to concurrently analyze both gene expression profiling and CNAs in the same individual.

RESULTS

We have developed an open-source R/Bioconductor package (iGC). Multiple input formats are supported and users can define their own criteria for identifying differentially expressed genes driven by CNAs. The analysis of two real microarray datasets demonstrated that the CNA-driven genes identified by the iGC package showed significantly higher Pearson correlation coefficients with their gene expression levels and copy numbers than those genes located in a genomic region with CNA. Compared with the Venn diagram approach, the iGC package showed better performance.

CONCLUSION

The iGC package is effective and useful for identifying CNA-driven genes. By simultaneously considering both comparative genomic and transcriptomic data, it can provide better understanding of biological and medical questions. The iGC package's source code and manual are freely available at https://www.bioconductor.org/packages/release/bioc/html/iGC.html .

Copy Number Variation Analysis by Droplet Digital PCR

The health impact of many copy number variants in our genome remains still largely to be discovered. Detecting and genotyping this often complex variation presents a technical challenge. Here we describe a 96-well format droplet digital PCR (ddPCR) protocol for genotyping a common copy variant in the human haptoglobin gene. ddPCR allows for high-throughput and accurate quantitation of gene copy numbers.

Between-species differences in gene copy number are enriched among functions critical for adaptive evolution in Arabidopsis halleri

Background

Gene copy number divergence between species is a form of genetic polymorphism that contributes significantly to both genome size and phenotypic variation. In plants, copy number expansions of single genes were implicated in cultivar- or species-specific tolerance of high levels of soil boron, aluminium or calamine-type heavy metals, respectively. Arabidopsis halleri is a zinc- and cadmium-hyperaccumulating extremophile species capable of growing on heavy-metal contaminated, toxic soils. In contrast, its non-accumulating sister species A. lyrata and the closely related reference model species A. thaliana exhibit merely basal metal tolerance.

Results

For a genome-wide assessment of the role of copy number divergence (CND) in lineage-specific environmental adaptation, we conducted cross-species array comparative genome hybridizations of three plant species and developed a global signal scaling procedure to adjust for sequence divergence. In A. halleri, transition metal homeostasis functions are enriched twofold among the genes detected as copy number expanded. Moreover, biotic stress functions including mostly disease Resistance (R) gene-related genes are enriched twofold among genes detected as copy number reduced, when compared to the abundance of these functions among all genes.

Conclusions

Our results provide genome-wide support for a link between evolutionary adaptation and CND in A. halleri as shown previously for Heavy metal ATPase4. Moreover our results support the hypothesis that elemental defences, which result from the hyperaccumulation of toxic metals, allow the reduction of classical defences against biotic stress as a trade-off.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3319-5) contains supplementary material, which is available to authorized users.

MiRNAs and miRNA Polymorphisms Modify Drug Response

Differences in expression of drug response-related genes contribute to inter-individual variation in drugs’ biological effects. MicroRNAs (miRNAs) are small noncoding RNAs emerging as new players in epigenetic regulation of gene expression at post-transcriptional level. MiRNAs regulate the expression of genes involved in drug metabolism, drug transportation, drug targets and downstream signal molecules directly or indirectly. MiRNA polymorphisms, the genetic variations affecting miRNA expression and/or miRNA-mRNA interaction, provide a new insight into the understanding of inter-individual difference in drug response. Here, we provide an overview of the recent progress in miRNAs mediated regulation of biotransformation enzymes, drug transporters, and nuclear receptors. We also describe the implications of miRNA polymorphisms in cancer chemotherapy response.

Identification of gene-drug interactions that impact patient survival in TCGA

BACKGROUND

With the advent of large scale biological data collection for various diseases, data analysis pipelines and workflows need to be established to build frameworks for integrative analysis. Here the authors present a pipeline for identifying disease specific gene-drug interactions using CNV (Copy Number Variation) and clinical data from the TCGA (The Cancer Genome Atlas) project. Two cancer types were selected for analysis, LGG (Brain lower grade glioma) and GBM (Glioblastoma multiforme), due to the possible progression from LGG to GBM in some cases. The copy number and clinical data were then used to preform survival analysis on a gene by gene basis on sub-populations of patients exposed to a given drug.

RESULTS

Several gene-drug interactions are identified, where the copy number of a gene is associated to survival of a patient exposed to a certain drug. Both Irinotecan/HAS2 (Hyaluronan synthase 2) and Bevacizumab/PGAM1 (Phosphoglycerate mutase 1) are interactions found in this study with independent confirmation. Independent work in colon, breast cancer and leukemia (Györffy, Breast Cancer Res Treat 123:725-731, 2010; Mueller, Mol Cancer Ther 11:3024-3032, 2010; Hitosugi, Cancer Cell 13:585-600, 2012) showed these two interactions can lead to increased survival.

CONCLUSION

While the pipeline produced several possible interactions where increased survival is linked to normal or increased copy number of a given gene for patients treated with a given drug, no instance of low copy number or full deletion was linked to increased survival. The development of this pipeline shows a promising utility to identify possible beneficial gene-drug interactions that could improve patient survival and may illustrate some of the problems inherent in this kind of analysis on these data.

Major influence of repetitive elements on disease-associated copy number variants (CNVs)

Copy number variants (CNVs) are important contributors to the human pathogenic genetic diversity as demonstrated by a number of cases reported in the literature. The high homology between repetitive elements may guide genomic stability which will give rise to CNVs either by non-allelic homologous recombination (NAHR) or non-homologous end joining (NHEJ). Here, we present a short guide based on previously documented cases of disease-associated CNVs in order to provide a general view on the impact of repeated elements on the stability of the genomic sequence and consequently in the origin of the human pathogenic variome.

Copy number variation contributes to cryptic genetic variation in outbreak lineages of Cryptococcus gattii from the North American Pacific Northwest

Background

Copy number variants (CNVs) are a class of structural variants (SVs) and are defined as fragments of DNA that are present at variable copy number in comparison with a reference genome. Recent advances in bioinformatics methodologies and sequencing technologies have enabled the high-resolution quantification of genome-wide CNVs. In pathogenic fungi SVs have been shown to alter gene expression, influence host specificity, and drive fungicide resistance, but little attention has focused specifically on CNVs. Using publicly available sequencing data, we identified 90 isolates across 212 Cryptococcus gattii genomes that belong to the VGII subgroups responsible for the recent deadly outbreaks in the North American Pacific Northwest. We generated CNV profiles for each sample to investigate the prevalence and function of CNV in C. gattii.

Results

We identified eight genetic clusters among publicly available Illumina whole genome sequence data from 212 C. gattii isolates through population structure analysis. Three clusters represent the VGIIa, VGIIb, and VGIIc subgroups from the North American Pacific Northwest. CNV was bioinformatically predicted and affected ~300–400 Kilobases (Kb) of the C. gattii VGII subgroup genomes. Sixty-seven loci, encompassing 58 genes, showed highly divergent patterns of copy number variation between VGII subgroups. Analysis of PFam domains within divergent CN variable genes revealed enrichment of protein domains associated with transport, cell wall organization and external encapsulating structure.

Conclusions

CNVs may contribute to pathological and phenotypic differences observed between the C. gattii VGIIa, VGIIb, and VGIIc subpopulations. Genes overlapping with population differentiated CNVs were enriched for several virulence related functional terms. These results uncover novel candidate genes to examine the genetic and functional underpinnings of C. gattii pathogenicity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3044-0) contains supplementary material, which is available to authorized users.

Concordance of copy number loss and down-regulation of tumor suppressor genes: a pan-cancer study

Background

Tumor suppressor genes (TSGs) encode the guardian molecules to control cell growth. The genomic alteration of TSGs may cause tumorigenesis and promote cancer progression. So far, investigators have mainly studied the functional effects of somatic single nucleotide variants in TSGs. Copy number variation (CNV) is another important form of genetic variation, and is often involved in cancer biology and drug treatment, but studies of CNV in TSGs are less represented in literature. In addition, there is a lack of a combinatory analysis of gene expression and CNV in this important gene set. Such a study may provide more insights into the relationship between gene dosage and tumorigenesis. To meet this demand, we performed a systematic analysis of CNVs and gene expression in TSGs to provide a systematic view of CNV and gene expression change in TSGs in pan-cancer.

Results

We identified 1170 TSGs with copy number gain or loss in 5846 tumor samples. Among them, 207 TSGs tended to have copy number loss (CNL), from which fifteen CNL hotspot regions were identified. The functional enrichment analysis revealed that the 207 TSGs were enriched in cancer-related pathways such as P53 signaling pathway and the P53 interactome. We further performed integrative analyses of CNV with gene expression using the data from the matched tumor samples. We found 81 TSGs with concordant CNL events and decreased gene expression in the tumor samples we examined. Remarkably, seven TSGs displayed concordant CNL and gene down-regulation in at least 50 tumor samples: MTAP (212 samples), PTEN (139), MCPH1 (85), FBXO25 (67), SMAD4 (64), TRIM35 (57), and RB1 (54). Specifically to MTAP, this concordance was found in 14 cancer types, an observation that is not much reported in literature yet. Further network-based analysis revealed that these TSGs with concordant CNL and gene down-regulation were highly connected.

Conclusions

This study provides a draft landscape of CNV in pan-cancer. Our findings of systematic concordance between CNL and down-regulation of gene expression may help better understand the TSG biology in tumorigenesis and cancer progression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2904-y) contains supplementary material, which is available to authorized users.

Associations of genetic risk scores based on adult adiposity pathways with childhood growth and adiposity measures

BACKGROUND

Results from genome-wide association studies (GWAS) identified many loci and biological pathways that influence adult body mass index (BMI). We aimed to identify if biological pathways related to adult BMI also affect infant growth and childhood adiposity measures.

METHODS

We used data from a population-based prospective cohort study among 3,975 children with a mean age of 6 years. Genetic risk scores were constructed based on the 97 SNPs associated with adult BMI previously identified with GWAS and on 28 BMI related biological pathways based on subsets of these 97 SNPs. Outcomes were infant peak weight velocity, BMI at adiposity peak and age at adiposity peak, and childhood BMI, total fat mass percentage, android/gynoid fat ratio, and preperitoneal fat area. Analyses were performed using linear regression models.

RESULTS

A higher overall adult BMI risk score was associated with infant BMI at adiposity peak and childhood BMI, total fat mass, android/gynoid fat ratio, and preperitoneal fat area (all p-values < 0.05). Analyses focused on specific biological pathways showed that the membrane proteins genetic risk score was associated with infant peak weight velocity, and the genetic risk scores related to neuronal developmental processes, hypothalamic processes, cyclicAMP, WNT-signaling, membrane proteins, monogenic obesity and/or energy homeostasis, glucose homeostasis, cell cycle, and muscle biology pathways were associated with childhood adiposity measures (all p-values <0.05). None of the pathways were associated with childhood preperitoneal fat area.

CONCLUSIONS

A genetic risk score based on 97 SNPs related to adult BMI was associated with peak weight velocity during infancy and general and abdominal fat measurements at the age of 6 years. Risk scores based on genetic variants linked to specific biological pathways, including central nervous system and hypothalamic processes, influence body fat development from early life onwards.

Phylogenetic-Derived Insights into the Evolution of Sialylation in Eukaryotes: Comprehensive Analysis of Vertebrate β-Galactoside α2,3/6-Sialyltransferases (ST3Gal and ST6Gal)

Cell surface of eukaryotic cells is covered with a wide variety of sialylated molecules involved in diverse biological processes and taking part in cell–cell interactions. Although the physiological relevance of these sialylated glycoconjugates in vertebrates begins to be deciphered, the origin and evolution of the genetic machinery implicated in their biosynthetic pathway are poorly understood. Among the variety of actors involved in the sialylation machinery, sialyltransferases are key enzymes for the biosynthesis of sialylated molecules. This review focus on β-galactoside α2,3/6-sialyltransferases belonging to the ST3Gal and ST6Gal families. We propose here an outline of the evolutionary history of these two major ST families. Comparative genomics, molecular phylogeny and structural bioinformatics provided insights into the functional innovations in sialic acid metabolism and enabled to explore how ST-gene function evolved in vertebrates.

A Highly Polymorphic Copy Number Variant in the NSF Gene is Associated with Cocaine Dependence

Cocaine dependence is a complex psychiatric disorder involving both genetic and environmental factors. Several neurotransmitter systems mediate cocaine’s effects, dependence and relapse, being the components of the neurotransmitter release machinery good candidates for the disorder. Previously, we identified a risk haplotype for cocaine dependence in the NSF gene, encoding the protein N-Ethylmaleimide-Sensitive Factor essential for synaptic vesicle turnover. Here we examined the possible contribution to cocaine dependence of a large copy number variant (CNV) that encompasses part of the NSF gene. We performed a case-control association study in a discovery sample (359 cases and 356 controls) and identified an association between cocaine dependence and the CNV (P = 0.013), that was confirmed in the replication sample (508 cases and 569 controls, P = 7.1e-03) and in a pooled analysis (P = 1.8e-04), with an over-representation of low number of copies in cases. Subsequently, we studied the functional impact of the CNV on gene expression and found that the levels of two NSF transcripts were significantly increased in peripheral blood mononuclear cells (PBMC) along with the number of copies of the CNV. These results, together with a previous study from our group, support the role of NSF in the susceptibility to cocaine dependence.

Identification of copy number variations in three Chinese horse breeds using 70K single nucleotide polymorphism BeadChip array

Copy number variation (CNV), an essential form of genetic variation, has been increasingly recognized as one promising genetic marker in the analysis of animal genomes. Here, we used the Equine 70K single nucleotide polymorphism genotyping array for the genome-wide detection of CNVs in 96 horses from three diverse Chinese breeds: Debao pony (DB), Mongolian horse (MG) and Yili horse (YL). A total of 287 CNVs were determined and merged into 122 CNV regions (CNVRs) ranging from 199 bp to 2344 kb in size and distributed in a heterogeneous manner on chromosomes. These CNVRs were integrated with seven existing reports to generate a composite genome-wide dataset of 1558 equine CNVRs, revealing 69 (56.6%) novel CNVRs. The majority (69.7%) of the 122 CNVRs overlapped with 438 genes, whereas 30.3% were located in intergenic regions. Most of these genes were associated with common CNVRs, which were shared by divergent horse breeds. As many as 60, 42 and 91 genes overlapping with the breed-specific ss were identified in DB, MG and YL respectively. Among these genes, FGF11, SPEM1, PPARG, CIDEB, HIVEP1 and GALR may have potential relevance to breed-specific traits. These findings provide valuable information for understanding the equine genome and facilitating association studies of economically important traits with equine CNVRs in the future.

Further evidence for deletions in 7p14.1 contributing to nonsyndromic cleft lip with or without cleft palate

BACKGROUND

Nonsyndromic cleft with or without cleft palate (nsCL/P) is a common birth defect. Although genome-wide association studies (GWAS) have identified numerous risk variants, a considerable fraction of the genetic heritability remains unknown. The aim of the present study was to replicate a previous finding that de novo deletions in a 62 kb region of chromosome 7p14 are a risk factor for nsCL/P, using an independent cohort.

METHODS

Data from a published case-control GWAS cohort of 399 patients and 1318 controls were used. Copy number variant (CNV) detection in the 62 kb candidate region of 7p14 was performed using QuantiSNP. Putative CNVs in probands were verified and validated by quantitative polymerase chain reaction. Segregation analyses were performed in family members for whom DNA was available.

RESULTS

Within the 62 kb candidate region, a deletion of 7.4 kb showed association with nsCL/P (13/387 cases, 20/1300 controls, plowest = 0.024, odds ratio = 2.22). In all families with a sporadic case (n = 3), the deletion occurred de novo. In multiplex families, both incomplete segregation and incomplete penetrance were observed.

CONCLUSION

The present data support the hypothesis that deletions at 7p14 are a common risk factor for nsCL/P. Genome-wide CNV analyses in nsCL/P cohorts are warranted to explore the functional relevance of these deletions and their contribution to nsCL/P, and to determine exact breakpoints. Birth Defects Research (Part A) 106:767-772, 2016. © 2016 Wiley Periodicals, Inc.

Cytogenetic abnormalities and genomic copy number variations in EPO (7q22) and SEC-61(7p11) genes in primary myelodysplastic syndromes

Myelodysplastic syndromes (MDSs) are heterogeneous clonal haematopoeitic stem cell disorders characterized by ineffective haematopoeisis, cytopenias and risk of progression to AML. We studied 150 MDS patients for cytogenetic aberrations and 60 patients with normal karyotype and 40 patients harboring cytogenetic abnormalities for copy number variations (CNVs). Cytogenetic abnormalities were detected in 46% of patients with a majority of patients harboring abnormalities of chromosome 7 and del (20q) at frequencies of 16% and 12% respectively. We explored the potential of quantitative multiplex PCR assay of short fluorescent fragments (QMPSF) to identify CNVs and correlated the findings with cytogenetic data and disease prognosis. CNVs (n=31) were detected in 28.3% of karyotypically normal and 23% patients with abnormal karyotype. Genetic losses or deletions (n=26) were more frequent than duplications (n=5). EPO (7q22) and SEC-61(7p11) emerged as new candidate genes susceptible to genetic losses with 57.7% deletions identified in regions on chromosome 7. The CNVs correlated with International Prognostic Scoring System (IPSS) intermediate disease risk group. Our integrative cytogenetic and copy number variation study suggests that abnormalities of chromosome 7 are predominant in Indian population and that they may play a secondary role in disease progression and should be evaluated further for asserting their clinical significance and influence on disease prognosis.

Comparative analysis of copy number variations in ulcerative colitis associated and sporadic colorectal neoplasia

Background

The incidence of and mortality from colorectal cancers (CRC) can be reduced by early detection. Currently there is a lack of established markers to detect early neoplastic changes. We aimed to identify the copy number variations (CNVs) and the associated genes which could be potential markers for the detection of neoplasia in both ulcerative colitis-associated neoplasia (UC-CRN) and sporadic colorectal neoplasia (S-CRN).

Methods

We employed array comparative genome hybridization (aCGH) to identify CNVs in tissue samples of UC nonprogressor, progressor and sporadic CRC. Select genes within these CNV regions as a panel of markers were validated using quantitative real time PCR (qRT-PCR) method along with the microsatellite instability (MSI) in an independent cohort of samples. Immunohistochemistry (IHC) analysis was also performed.

Results

Integrated analysis showed 10 overlapping CNV regions between UC-Progressor and S-CRN, with the 8q and 12p regions showing greater overlap. The qRT-PCR based panel of MYC, MYCN, CCND1, CCND2, EGFR and FNDC3A was successful in detecting neoplasia with an overall accuracy of 54 % in S-CRN compared to that of 29 % in UC neoplastic samples. IHC study showed that p53 and CCND1 were significantly overexpressed with an increasing frequency from pre-neoplastic to neoplastic stages. EGFR and AMACR were expressed only in the neoplastic conditions.

Conclusion

CNVs that are common and unique to both UC-associated and sporadic colorectal neoplasm could be the key players driving carcinogenesis. Comparative analysis of CNVs provides testable driver aberrations but needs further evaluation in larger cohorts of samples. These markers may help in developing more effective neoplasia-detection strategies during screening and surveillance programs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-016-2303-4) contains supplementary material, which is available to authorized users.

Towards plant pangenomics

As an increasing number of genome sequences become available for a wide range of species, there is a growing understanding that the genome of a single individual is insufficient to represent the gene diversity within a whole species. Many studies examine the sequence diversity within genes, and this allelic variation is an important source of phenotypic variation which can be selected for by man or nature. However, the significant gene presence/absence variation that has been observed within species and the impact of this variation on traits is only now being studied in detail. The sum of the genes for a species is termed the pangenome, and the determination and characterization of the pangenome is a requirement to understand variation within a species. In this review, we explore the current progress in pangenomics as well as methods and approaches for the characterization of pangenomes for a wide range of plant species.

Noise propagation with interlinked feed-forward pathways

Functionally similar pathways are often seen in biological systems, forming feed-forward controls. The robustness in network motifs such as feed-forward loops (FFLs) has been reported previously. In this work, we studied noise propagation in a development network that has multiple interlinked FFLs. A FFL has the potential of asymmetric noise-filtering (i.e., it works at either the “ON” or the “OFF” state in the target gene). With multiple, interlinked FFLs, we show that the propagated noises are largely filtered regardless of the states in the input genes. The noise-filtering property of an interlinked FFL can be largely derived from that of the individual FFLs, and with interlinked FFLs, it is possible to filter noises in both “ON” and “OFF” states in the output. We demonstrated the noise filtering effect in the developmental regulatory network of Caenorhabditis elegans that controls the timing of distal tip cell (DTC) migration. The roles of positive feedback loops involving blmp-1 and the degradation regulation of DRE-1 also studied. Our analyses allow for better inference from network structures to noise-filtering properties, and provide insights into the mechanisms behind the precise DTC migration controls in space and time.

Drug Repurposing Is a New Opportunity for Developing Drugs against Neuropsychiatric Disorders

Better the drugs you know than the drugs you do not know. Drug repurposing is a promising, fast, and cost effective method that can overcome traditional de novo drug discovery and development challenges of targeting neuropsychiatric and other disorders. Drug discovery and development targeting neuropsychiatric disorders are complicated because of the limitations in understanding pathophysiological phenomena. In addition, traditional de novo drug discovery and development are risky, expensive, and time-consuming processes. One alternative approach, drug repurposing, has emerged taking advantage of off-target effects of the existing drugs. In order to identify new opportunities for the existing drugs, it is essential for us to understand the mechanisms of action of drugs, both biologically and pharmacologically. By doing this, drug repurposing would be a more effective method to develop drugs against neuropsychiatric and other disorders. Here, we review the difficulties in drug discovery and development in neuropsychiatric disorders and the extent and perspectives of drug repurposing.

miR-23a, miR-146a and miR-301a confer predisposition to Vogt-Koyanagi-Harada syndrome but not to Behcet's disease

Ninety-eight miRNAs are involved in the immune response. However, the genetic roles of these miRNAs remain unclear in Behcet's disease (BD) and Vogt-Koyanagi-Harada (VKH) syndrome. This study aimed to explore the association and functional roles of copy number variants (CNV) in several miRNAs with BD and VKH syndrome. Genotyping of CNVs was examined by TaqMan PCR. The expression of miR-23a, transfection efficiency and cytokine production were measured by real-time PCR, flow cytometry or ELISA. First, replication and combined studies for miR-23a, miR-146a and miR-301a demonstrated a similar association with VKH syndrome (Combined: P = 5.53 × 10(-8); P = 8.43 × 10(-31); P = 9.23 × 10(-8), respectively). No association of CNVs of the above mentioned miRNAs was observed in BD patients. mRNA expression of miR-23a showed a positive association with its copy numbers. Additionally, individuals with high copy number of miR-23a show an increased production of interleukin-6 (IL-6), but not IL-8 and monocyte chemoattractant protein-1 (MCP-1) by stimulated PBMCs. miR-23a transfected ARPE-19 cells modulated the production of IL-6 and IL-8, but not MCP-1. Our results suggest that CNVs of miR-146a, miR-23a and miR-301a confer susceptibility to VKH syndrome, but not to BD. The contribution of miR-23a to VKH syndrome may be mediated by increasing the production of IL-6.

Clinical and Genetic Heterogeneity of the 15q13.3 Microdeletion Syndrome

The 15q13.3 microdeletion is a recurrent CNV, presumably mediated by NAHR between segmental duplications in chromosome 15. The 15q13.3 deletion and duplication are associated with a wide range of clinical manifestations, such as intellectual deficits, seizures, autism, language and developmental delay, neuropsychiatric impairments, and behavioral problems illustrating incomplete penetrance and expressivity. This study comprises an evaluation of 106 symptomatic patients carrying the heterozygous deletion, as well as of 21 patients carrying the duplication, who have been described in previous studies. The analysis shows considerable heterogeneity for the manifestation of different key symptoms and familiar occurrence. Furthermore, 8 new patients are introduced. Convoluted familiar connections give new insights into the complexity of symptomatic manifestation. In previous studies, different opinions have been expressed as to the nature and precise location of the deletion breakpoints. Here, we show that not CHRNA7 and CHRFAM7A, but rather FAM7A or GOLGA8, serve as breakpoint regions concerning our patients. The deletion is described as heterogeneous in size. However, we assume that not only different breakpoints but also the imprecision of aCGH analysis on chromosome 15 due to segmental duplications accounts for the variability in size.

Structural Variant Detection by Large-scale Sequencing Reveals New Evolutionary Evidence on Breed Divergence between Chinese and European Pigs

In this study, we performed a genome-wide SV detection among the genomes of thirteen pigs from diverse Chinese and European originated breeds by next genetation sequencing, and constrcuted a single-nucleotide resolution map involving 56,930 putative SVs. We firstly identified a SV hotspot spanning 35 Mb region on the X chromosome specifically in the genomes of Chinese originated individuals. Further scrutinizing this region by large-scale sequencing data of extra 111 individuals, we obtained the confirmatory evidence on our initial finding. Moreover, thirty five SV-related genes within the hotspot region, being of importance for reproduction ability, rendered significant different evolution rates between Chinese and European originated breeds. The SV hotspot identified herein offers a novel evidence for assessing phylogenetic relationships, as well as likely explains the genetic difference of corresponding phenotypes and features, among Chinese and European pig breeds. Furthermore, we employed various SVs to infer genetic structure of individuls surveyed. We found SVs can clearly detect the difference of genetic background among individuals. This clues us that genome-wide SVs can capture majority of geneic variation and be applied into cladistic analyses. Characterizing whole genome SVs demonstrated that SVs are significantly enriched/depleted with various genomic features.

Genotyping Single Nucleotide Polymorphisms and Copy Number Variability of the FCGRs Expressed on NK Cells

Natural killer (NK) cells are one of the main effector immune cells involved in antibody-dependent cell-mediated cytotoxicity (ADCC). Upon recognition of cell-bound IgG antibodies, which occurs through Fc gamma receptors (FCGRs) expressed on the cell surface of NK cells, NK cells become activated and lyse target tumor or infected cells. The FCGRs, FCGR3A and FCGR2C, expressed on the surface of NK cells have single nucleotide polymorphisms (SNPs) that result in differential activity of NK cells. In addition to SNP genetic variation within each of these genes, the FCGRs are subject to copy number variation (CNV), which leads to variable protein expression levels on the cell surface. Studies have found that FCGR genotype for FCGR3A and FCGR2C is associated with variation in the response to immunotherapy.Due to high sequence homology within FCGR3 and FCGR2 families, there are difficulties associated with genotyping these specific receptors related to cross-amplification of non-targeted FCGRs. To improve specificity for both FCGR3A and FCGR2C, Rnase-H (RH) primers were designed to amplify specifically FCGR3A (while not co-amplifying FCGR3B) and FCGR2C (while not co-amplifying FCGR2B). In addition, fluorescently labeled locked nucleic acid (LNA) probes provide additional precision for determination of the SNPs within both FCGR3A and FCGR2C. For CNV determination, separate fluorescently labeled probes for FCGR3A, and for FCGR2C, can be used with the same RH primers for each gene. These probes can be combined in the same well with control primers/probe for a known diploid gene and used to calculate the copy number of both FCGR3A and FCGR2C. Here we provide new detailed methodology that allows for the specific amplification of these FCGRs in a single PCR reaction, allowing for genotyping of both the SNPs and CNVs using real-time PCR.

Changes in correlation between promoter methylation and gene expression in cancer

BACKGROUND

Methylation of high-density CpG regions known as CpG Islands (CGIs) has been widely described as a mechanism associated with gene expression regulation. Aberrant promoter methylation is considered a hallmark of cancer involved in silencing of tumor suppressor genes and activation of oncogenes. However, recent studies have also challenged the simple model of gene expression control by promoter methylation in cancer, and the precise mechanism of and role played by changes in DNA methylation in carcinogenesis remains elusive.

RESULTS

Using a large dataset of 672 matched cancerous and healthy methylomes, gene expression, and copy number profiles accross 3 types of tissues from The Cancer Genome Atlas (TCGA), we perform a detailed meta-analysis to clarify the interplay between promoter methylation and gene expression in normal and cancer samples. On the one hand, we recover the existence of a CpG island methylator phenotype (CIMP) with prognostic value in a subset of breast, colon and lung cancer samples, where a common subset of promoter CGIs hypomethylated in normal samples become hypermethylated. However, this hypermethylation is not accompanied by a decrease in expression of the corresponding genes, which are already lowly expressed in the normal genes. On the other hand, we identify tissue-specific sets of genes, different between normal and cancer samples, whose inter-individual variation in expression is significantly correlated with the variation in methylation of the 3' flanking regions of the promoter CGIs. These subsets of genes are not the same in the different tissues, nor between normal and cancerous samples, but transcription factors are over-represented in all subsets.

CONCLUSION

Our results suggest that epigenetic reprogramming in cancer does not contribute to cancer development via direct inhibition of gene expression through promoter hypermethylation. It may instead modify how the expression of a few specific genes, particularly transcription factors, are associated with DNA methylation variations in a tissue-dependent manner.

Rare genomic rearrangement in a boy with Williams-Beuren syndrome associated to XYY syndrome and intriguing behavior

Williams-Beuren syndrome (WBS) is caused by a hemizygous contiguous gene microdeletion of 1.55-1.84 Mb at 7q11.23 region. Approximately, 28 genes have been shown to contribute to classical phenotype of SWB with presence of dysmorphic facial features, supravalvular aortic stenosis (SVAS), intellectual disability, and overfriendliness. With the use of Microarray-based comparative genomic hybridization and other molecular cytogenetic techniques, is possible define with more accuracy partial or atypical deletion and refine the genotype-phenotype correlation. Here, we report on a rare genomic structural rearrangement in a boy with atypical deletion in 7q11.23 and XYY syndrome with characteristic clinical signs, but not sufficient for the diagnosis of WBS. Cytogenetic analysis of G-banding showed a karyotype 47,XYY. Analysis of DNA with the technique of MLPA (Multiplex Ligation-dependent Probe Amplification) using kits a combination of kits (P064, P036, P070, and P029) identified an atypical deletion on 7q11.23. In addition, high resolution SNP Oligonucleotide Microarray Analysis (SNP-array) confirmed the alterations found by MLPA and revealed others pathogenic CNVs, in the chromosomes 7 and X. The present report demonstrates an association not yet described in literature, between Williams-Beuren syndrome and 47,XYY. The identification of atypical deletion in 7q11.23 concomitant to additional pathogenic CNVs in others genomic regions allows a better comprehension of clinical consequences of atypical genomic rearrangements.

Integrated transcriptional profiling and genomic analyses reveal RPN2 and HMGB1 as promising biomarkers in colorectal cancer

Colorectal cancer (CRC) is a heterogeneous disease that is associated with a gradual accumulation of genetic and epigenetic alterations. Among all CRC stages, stage II tumors are highly heterogeneous with a high relapse rate in about 20–25 % of stage II CRC patients following surgery. Thus, a comprehensive analysis of gene signatures to identify aggressive and metastatic phenotypes in stage II CRC is desired for a more accurate disease classification and outcome prediction. By utilizing a Cancer Array, containing 440 oncogenes and tumor suppressors to profile mRNA expression, we identified a larger number of differentially expressed genes in poorly differentiated stage II colorectal adenocarcinoma tissues, compared to their matched normal tissues. Ontology and Ingenuity Pathway Analysis (IPA) indicated that these genes are involved in functional mechanisms associated with several transcription factors. Genomic alterations of these genes were also investigated through The Cancer Genome Atlas (TCGA) database, utilizing 195 published CRC specimens. The percentage of genomic alterations in these genes was ranked based on their mRNA expression, copy number variations and mutations. This data was further combined with published microarray studies from a large set of CRC tumors classified based on prognostic features. This led to the identification of eight candidate genes including RPN2, HMGB1, AARS, IGFBP3, STAT1, HYOU1, NQO1 and PEA15 that were associated with the progressive phenotype. In particular, RPN2 and HMGB1 displayed a higher genomic alteration frequency in CRC, compared to eight other major solid cancers. Immunohistochemistry was performed on additional 78 stage I–IV CRC samples, where RPN2 protein immunostaining exhibited a significant association with stage III/IV tumors, distant metastasis, and poor differentiation, indicating that RPN2 expression is associated with poor prognosis. Further, our study revealed significant transcriptional regulatory mechanisms, networks and gene signatures, underlying CRC malignant progression and phenotype warranting future clinical investigations.

Interactive analysis and assessment of single-cell copy-number variations

We present Ginkgo (http://qb.cshl.edu/ginkgo), a user-friendly, open-source web platform for the analysis of single-cell copy-number variations (CNVs). Ginkgo automatically constructs copy-number profiles of cells from mapped reads and constructs phylogenetic trees of related cells. We validated Ginkgo by reproducing the results of five major studies. After comparing three commonly used single-cell amplification techniques, we concluded that degenerate oligonucleotide-primed PCR is the most consistent for CNV analysis.

Deciphering the associations between gene expression and copy number alteration using a sparse double Laplacian shrinkage approach

MOTIVATION

Both gene expression levels (GEs) and copy number alterations (CNAs) have important biological implications. GEs are partly regulated by CNAs, and much effort has been devoted to understanding their relations. The regulation analysis is challenging with one gene expression possibly regulated by multiple CNAs and one CNA potentially regulating the expressions of multiple genes. The correlations among GEs and among CNAs make the analysis even more complicated. The existing methods have limitations and cannot comprehensively describe the regulation.

RESULTS

A sparse double Laplacian shrinkage method is developed. It jointly models the effects of multiple CNAs on multiple GEs. Penalization is adopted to achieve sparsity and identify the regulation relationships. Network adjacency is computed to describe the interconnections among GEs and among CNAs. Two Laplacian shrinkage penalties are imposed to accommodate the network adjacency measures. Simulation shows that the proposed method outperforms the competing alternatives with more accurate marker identification. The Cancer Genome Atlas data are analysed to further demonstrate advantages of the proposed method.

AVAILABILITY AND IMPLEMENTATION

R code is available at http://works.bepress.com/shuangge/49/.

Copy number variants associated with epilepsy from gene expression microarrays

We aimed to identify novel copy number variations (CNV) that might contribute to the pathogenesis of epilepsy. Epilepsy is a common brain disorder characterized by recurring seizures and various serious comorbidities, including respiratory, cardiovascular, and neurologic dysfunction. CNV have recently been considered as important risk factors for epilepsy. With public gene expression data from brain tissue of 23 epilepsy patients and 23 healthy controls, we detected CNV using the R language package CAFÉ. Real-time quantitative polymerase chain reaction validation was performed in a further nine patients and 10 controls. Functional analyses of the genes in the validated CNV were also carried out, using Ingenuity pathway analysis. Three copy number abnormalities (19q13.33, 19q13.11 and 4q35.1) were detected with the gene expression data. The duplication in 19q13.33 (approximately 1.22 million bases) was further validated in three additional epilepsy patients, and the deletion in 19q13.11 (approximately 855 kilobases) was further validated in another two epilepsy patients. The functional analyses of the genes in these two CNV suggested that they may be involved in the pathogenesis of epilepsy. The CNV that we detected may be common genetic etiological factors of epilepsy, and there is potential for the identification of a novel biomarker for treatment from these CNV regions.

Mechanisms of Origin, Phenotypic Effects and Diagnostic Implications of Complex Chromosome Rearrangements

Complex chromosome rearrangements (CCRs) are currently defined as structural genome variations that involve more than 2 chromosome breaks and result in exchanges of chromosomal segments. They are thought to be extremely rare, but their detection rate is rising because of improvements in molecular cytogenetic technology. Their population frequency is also underestimated, since many CCRs may not elicit a phenotypic effect. CCRs may be the result of fork stalling and template switching, microhomology-mediated break-induced repair, breakage-fusion-bridge cycles, or chromothripsis. Patients with chromosomal instability syndromes show elevated rates of CCRs due to impaired DNA double-strand break responses during meiosis. Therefore, the putative functions of the proteins encoded by ATM, BLM, WRN, ATR, MRE11, NBS1, and RAD51 in preventing CCRs are discussed. CCRs may exert a pathogenic effect by either (1) gene dosage-dependent mechanisms, e.g. haploinsufficiency, (2) mechanisms based on disruption of the genomic architecture, such that genes, parts of genes or regulatory elements are truncated, fused or relocated and thus their interactions disturbed - these mechanisms will predominantly affect gene expression - or (3) mixed mutation mechanisms in which a CCR on one chromosome is combined with a different type of mutation on the other chromosome. Such inferred mechanisms of pathogenicity need corroboration by mRNA sequencing. Also, future studies with in vitro models, such as inducible pluripotent stem cells from patients with CCRs, and transgenic model organisms should substantiate current inferences regarding putative pathogenic effects of CCRs. The ramifications of the growing body of information on CCRs for clinical and experimental genetics and future treatment modalities are briefly illustrated with 2 cases, one of which suggests KDM4C (JMJD2C) as a novel candidate gene for mental retardation.

A genome-wide scan for copy number variations using high-density single nucleotide polymorphism array in Simmental cattle

Copy number variations (CNVs) have recently been identified as promising sources of genetic variation, complementary to single nucleotide polymorphisms (SNPs). As a result, detection of CNVs has attracted a great deal of attention. In this study, we performed genome-wide CNV detection using Illumina Bovine HD BeadChip (770k) data on 792 Simmental cattle. A total of 263 CNV regions (CNVRs) were identified, which included 137 losses, 102 gains and 24 regions classified as both loss and gain, covering 35.48 Mb (1.41%) of the bovine genome. The length of these CNVRs ranged from 10.18 kb to 1.76 Mb, with an average length of 134.78 kb and a median length of 61.95 kb. In 136 of these regions, a total of 313 genes were identified related to biological functions such as transmembrane activity and olfactory transduction activity. To validate the results, we performed quantitative PCR to detect nine randomly selected CNVRs and successfully confirmed seven (77.6%) of them. Our results present a map of cattle CNVs derived from high-density SNP data, which expands the current CNV map of the cattle genome and provides useful information for investigation of genomic structural variation in cattle.

Genome-Wide Copy Number Variations Using SNP Genotyping in a Mixed Breed Swine Population

Copy number variations (CNVs) are increasingly understood to affect phenotypic variation. This study uses SNP genotyping of trios of mixed breed swine to add to the catalog of known genotypic variation in an important agricultural animal. PorcineSNP60 BeadChip genotypes were collected from 1802 pigs that combined to form 1621 trios. These trios were from the crosses of 50 boars with 525 sows producing 1621 piglets. The pigs were part of a population that was a mix of ¼ Duroc, ½ Landrace and ¼ Yorkshire breeds. Merging the overlapping CNVs that were observed in two or more individuals to form CNV regions (CNVRs) yielded 502 CNVRs across the autosomes. The CNVRs intersected genes, as defined by RefSeq, 84% of the time – 420 out of 502. The results of this study are compared and contrasted to other swine studies using similar and different methods of detecting CNVR. While progress is being made in this field, more work needs to be done to improve consistency and confidence in CNVR results.