Copy number variants and genetic traits: closer to the resolution of phenotypic to genotypic variability

DeviCNV: detection and visualization of exon-level copy number variants in targeted next-generation sequencing data

BACKGROUND

Targeted next-generation sequencing (NGS) is increasingly being adopted in clinical laboratories for genomic diagnostic tests.

RESULTS

We developed a new computational method, DeviCNV, intended for the detection of exon-level copy number variants (CNVs) in targeted NGS data. DeviCNV builds linear regression models with bootstrapping for every probe to capture the relationship between read depth of an individual probe and the median of read depth values of all probes in the sample. From the regression models, it estimates the read depth ratio of the observed and predicted read depth with confidence interval for each probe which is applied to a circular binary segmentation (CBS) algorithm to obtain CNV candidates. Then, it assigns confidence scores to those candidates based on the reliability and strength of the CNV signals inferred from the read depth ratios of the probes within them. Finally, it also provides gene-centric plots with confidence levels of CNV candidates for visual inspection. We applied DeviCNV to targeted NGS data generated for newborn screening and demonstrated its ability to detect novel pathogenic CNVs from clinical samples.

CONCLUSIONS

We propose a new pragmatic method for detecting CNVs in targeted NGS data with an intuitive visualization and a systematic method to assign confidence scores for candidate CNVs. Since DeviCNV was developed for use in clinical diagnosis, sensitivity is increased by the detection of exon-level CNVs.

Copy number variations in Friesian horses and genetic risk factors for insect bite hypersensitivity

Background

Many common and relevant diseases affecting equine welfare have yet to be tested regarding structural variants such as copy number variations (CNVs). CNVs make up a substantial proportion of total genetic variability in populations of many species, resulting in more sequence differences between individuals than SNPs. Associations between CNVs and disease phenotypes have been established in several species, but equine CNV studies have been limited. Aim of this study was to identify CNVs and to perform a genome-wide association (GWA) study in Friesian horses to identify genomic loci associated with insect bite hypersensitivity (IBH), a common seasonal allergic dermatitis observed in many horse breeds worldwide.

Results

Genotypes were obtained using the Axiom® Equine Genotyping Array containing 670,796 SNPs. After quality control of genotypes, 15,041 CNVs and 5350 CNV regions (CNVRs) were identified in 222 Friesian horses. Coverage of the total genome by CNVRs was 11.2% with 49.2% of CNVRs containing genes. 58.0% of CNVRs were novel (i.e. so far only identified in Friesian horses). A SNP- and CNV-based GWA analysis was performed, where about half of the horses were affected by IBH. The SNP-based analysis showed a highly significant association between the MHC region on ECA20 and IBH in Friesian horses. Associations between the MHC region on ECA20 and IBH were also detected based on the CNV-based analysis. However, CNVs associated with IBH in Friesian horses were not often in close proximity to SNPs identified to be associated with IBH.

Conclusions

CNVs were identified in a large sample of the Friesian horse population, thereby contributing to our knowledge on CNVs in horses and facilitating our understanding of the equine genome and its phenotypic expression. A clear association was identified between the MHC region on ECA20 and IBH in Friesian horses based on both SNP- and CNV-based GWA studies. These results imply that MHC contributes to IBH sensitivity in Friesian horses. Although subsequent analyses are needed for verification, nucleotide differences, as well as more complex structural variations like CNVs, seem to contribute to IBH sensitivity. IBH should be considered as a common disease with a complex genomic architecture.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0657-0) contains supplementary material, which is available to authorized users.

Genome-wide significant risk factors on chromosome 19 and the APOE locus

The apolipoprotein E (APOE) gene on chromosome 19q13.32, was the first, and remains the strongest, genetic risk factor for Alzheimer’s disease (AD). Additional signals associated with AD have been located in chromosome 19, including ABCA7 (19p13.3) and CD33 (19q13.41). The ABCA7 gene has been replicated in most populations. However, the contribution to AD of other signals close to APOE gene remains controversial. Possible explanations for inconsistency between reports include long range linkage disequilibrium (LRLD). We analysed the contribution of ABCA7 and CD33 loci to AD risk and explore LRLD patterns across APOE region. To evaluate AD risk conferred by ABCA7 rs4147929:G>A and CD33 rs3865444:C>A, we used a large Spanish population (1796 AD cases, 2642 controls). The ABCA7 rs4147929:G>A SNP effect was nominally replicated in the Spanish cohort and reached genome-wide significance after meta-analysis (odds ratio (OR)=1.15, 95% confidence interval (95% CI)=1.12–1.19; P = 1.60 x 10^-19). CD33 rs3865444:C>A was not associated with AD in the dataset. The meta-analysis was also negative (OR=0.98, 95% CI=0.93–1.04; P=0.48). After exploring LRLD patterns between APOE and CD33 in several datasets, we found significant LD (D’ >0.20; P <0.030) between APOE-Ɛ2 and CD33 rs3865444C>A in two of five datasets, suggesting the presence of a non-universal long range interaction between these loci affecting to some populations. In conclusion, we provide here evidence of genetic association of the ABCA7 locus in the Spanish population and also propose a plausible explanation for the controversy on the contribution of CD33 to AD susceptibility.

Noise-robust assessment of SNP array based CNV calls through local noise estimation of log R ratios

Arrays based on single nucleotide polymorphisms (SNPs) have been successful for the large scale discovery of copy number variants (CNVs). However, current CNV calling algorithms still have limitations in detecting CNVs with high specificity and sensitivity, especially in case of small (<100 kb) CNVs. Therefore, this study presents a simple statistical analysis to evaluate CNV calls from SNP arrays in order to improve the noise-robustness of existing CNV calling algorithms. The proposed approach estimates local noise of log R ratios and returns the probability that a certain observation is different from this log R ratio noise level. This probability can be triggered at different thresholds to tailor specificity and/or sensitivity in a flexible way. Moreover, a comparison based on qPCR experiments showed that the proposed noise-robust CNV calls outperformed original ones for multiple threshold values.

Role of DNA copy number variation in dyslipidemias

PURPOSE OF REVIEW

DNA copy number variations (CNVs) are quantitative structural rearrangements that include deletions, duplications, and higher order amplifications. Because of technical limitations, the contribution of this common form of genetic variation to regulation of lipid metabolism and dyslipidemia has been underestimated.

RECENT FINDINGS

Recent literature involving CNVs and dyslipidemias has focused mainly on rare CNVs causing familial hypercholesterolemia, and a common CNV polymorphism as the major determinant of lipoprotein(a) plasma concentrations. Additionally, there is tantalizing evidence of largely uninvestigated but plausible presence of CNVs underlying other dyslipidemias. We also discuss the future role of improved technologies in facilitating more economic, routine CNV assessment in dyslipidemias.

SUMMARY

CNVs account for large proportion of human genetic variation and are already known to contribute to susceptibility of dyslipidemias, particularly in about 10% of familial hypercholesterolemia patients. Increasing availability of clinical next-generation sequencing and bioinformatics presents a cost-effective opportunity for novel CNV discoveries in dyslipidemias.

AMYCNE: Confident copy number assessment using whole genome sequencing data

Copy number variations (CNVs) within the human genome have been linked to a diversity of inherited diseases and phenotypic traits. The currently used methodology to measure copy numbers has limited resolution and/or precision, especially for regions with more than 4 copies. Whole genome sequencing (WGS) offers an alternative data source to allow for the detection and characterization of the copy number across different genomic regions in a single experiment. A plethora of tools have been developed to utilize WGS data for CNV detection. None of these tools are designed specifically to accurately estimate copy numbers of complex regions in a small cohort or clinical setting. Herein, we present AMYCNE (automatic modeling functionality for copy number estimation), a CNV analysis tool using WGS data. AMYCNE is multifunctional and performs copy number estimation of complex regions, annotation of VCF files, and CNV detection on individual samples. The performance of AMYCNE was evaluated using AMY1A ddPCR measurements from 86 unrelated individuals. In addition, we validated the accuracy of AMYCNE copy number predictions on two additional genes (FCGR3A and FCGR3B) using datasets available through the 1000 genomes consortium. Finally, we simulated levels of mosaic loss and gain of chromosome X and used this dataset for benchmarking AMYCNE. The results show a high concordance between AMYCNE and ddPCR, validating the use of AMYCNE to measure tandem AMY1 repeats with high accuracy. This opens up new possibilities for the use of WGS for accurate copy number determination of other complex regions in the genome in small cohorts or single individuals.

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.

Dosage-Dependent Expression Variation Suppressed on the Drosophila Male X Chromosome

DNA copy number variation is associated with many high phenotypic heterogeneity disorders. We systematically examined the impact of Drosophila melanogaster deletions on gene expression profiles to ask whether increased expression variability owing to reduced gene dose might underlie this phenotypic heterogeneity. Indeed, we found that one-dose genes have higher gene expression variability relative to two-dose genes. We then asked whether this increase in variability could be explained by intrinsic noise within cells due to stochastic biochemical events, or whether expression variability is due to extrinsic noise arising from more complex interactions. Our modeling showed that intrinsic gene expression noise averages at the organism level and thus cannot explain increased variation in one-dose gene expression. Interestingly, expression variability was related to the magnitude of expression compensation, suggesting that regulation, induced by gene dose reduction, is noisy. In a remarkable exception to this rule, the single X chromosome of males showed reduced expression variability, even compared with two-dose genes. Analysis of sex-transformed flies indicates that X expression variability is independent of the male differentiation program. Instead, we uncovered a correlation between occupancy of the chromatin-modifying protein encoded by males absent on the first (mof) and expression variability, linking noise suppression to the specialized X chromosome dosage compensation system. MOF occupancy on autosomes in both sexes also lowered transcriptional noise. Our results demonstrate that gene dose reduction can lead to heterogeneous responses, which are often noisy. This has implications for understanding gene network regulatory interactions and phenotypic heterogeneity. Additionally, chromatin modification appears to play a role in dampening transcriptional noise.

Copy Number Variation in SOX6 Contributes to Chicken Muscle Development

Copy number variations (CNVs), which cover many functional genes, are associated with complex diseases, phenotypic diversity and traits that are economically important to raising chickens. The sex-determining region Y-box 6 (Sox6) plays a key role in fast-twitch muscle fiber differentiation of zebrafish and mice, but it is still unknown whether SOX6 plays a role in chicken skeletal muscle development. We identified two copy number polymorphisms (CNPs) which were significantly related to different traits on the genome level in chickens by AccuCopy^® and CNVplex^® analyses. Notably, five white recessive rock (CN = 1, CN = 3) variant individuals and two Xinghua (CN = 3) variant individuals contain a CNP13 (chromosome5: 10,500,294-10,675,531) which overlaps with SOX6. There is a disordered region in SOX6 proteins 265-579 aa coded by a partial CNV overlapping region. A quantitative real-time polymerase chain reaction showed that the expression level of SOX6 mRNA was positively associated with CNV and highly expressed during the skeletal muscle cell differentiation in chickens. After the knockdown of the SOX6, the expression levels of IGFIR1, MYF6, SOX9, SHOX and CCND1 were significantly down-regulated. All of them directly linked to muscle development. These results suggest that the number of CNVs in the CNP13 is positively associated with the expression level of SOX6, which promotes the proliferation and differentiation of skeletal muscle cells by up-regulating the expression levels of the muscle-growth-related genes in chickens as in other animal species.

Genome-wide DNA copy number analysis in clonally expanded human ovarian cancer cells with distinct invasive/migratory capacities

Ovarian cancer has the worst prognosis of any gynecological malignancy, and generally presents with metastasis at advanced stages. Copy number variation (CNV) frequently contributes to the alteration of oncogenic drivers. In this study, we sought to identify genetic targets in heterogeneous clones from human ovarian cancers cells. We used array-based technology to systematically assess all the genes with CNVs in cell models clonally expanded from A2780 and SKOV3 ovarian cancer cell lines with distinct highly and minimally invasive/migratory capacities. We found that copy number alterations differed between matched highly and minimally invasive/migratory subclones, differentially affecting specific functional processes including immune response processes, DNA damage repair, cell cycle and cell proliferation. We also identified seven genes as strong candidates, including DDB1, ERCC1, ERCC2, PRPF19, BCAT1, CDKN1B and MARK4, by integrating the above data with gene expression and clinical outcome data. Thus, by determining the molecular signatures of heterogeneous invasive/migratory ovarian cancer cells, we identified genes that could be specifically targeted for the treatment and prognosis of advanced ovarian cancers.

Epigenetic variation between urban and rural populations of Darwin's finches

Background

The molecular basis of evolutionary change is assumed to be genetic variation. However, growing evidence suggests that epigenetic mechanisms, such as DNA methylation, may also be involved in rapid adaptation to new environments. An important first step in evaluating this hypothesis is to test for the presence of epigenetic variation between natural populations living under different environmental conditions.

Results

In the current study we explored variation between populations of Darwin’s finches, which comprise one of the best-studied examples of adaptive radiation. We tested for morphological, genetic, and epigenetic differences between adjacent “urban” and “rural” populations of each of two species of ground finches, Geospiza fortis and G. fuliginosa, on Santa Cruz Island in the Galápagos. Using data collected from more than 1000 birds, we found significant morphological differences between populations of G. fortis, but not G. fuliginosa. We did not find large size copy number variation (CNV) genetic differences between populations of either species. However, other genetic variants were not investigated. In contrast, we did find dramatic epigenetic differences between the urban and rural populations of both species, based on DNA methylation analysis. We explored genomic features and gene associations of the differentially DNA methylated regions (DMR), as well as their possible functional significance.

Conclusions

In summary, our study documents local population epigenetic variation within each of two species of Darwin’s finches.

Electronic supplementary material

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

Next generation sequencing: implications in personalized medicine and pharmacogenomics

A breakthrough in next generation sequencing (NGS) in the last decade provided an unprecedented opportunity to investigate genetic variations in humans and their roles in health and disease. NGS offers regional genomic sequencing such as whole exome sequencing of coding regions of all genes, as well as whole genome sequencing. RNA-seq offers sequencing of the entire transcriptome and ChIP-seq allows for sequencing the epigenetic architecture of the genome. Identifying genetic variations in individuals can be used to predict disease risk, with the potential to halt or retard disease progression. NGS can also be used to predict the response to or adverse effects of drugs or to calculate appropriate drug dosage. Such a personalized medicine also provides the possibility to treat diseases based on the genetic makeup of the patient. Here, we review the basics of NGS technologies and their application in human diseases to foster human healthcare and personalized medicine.

Next generation sequencing and array-based comparative genomic hybridization for molecular diagnosis of pediatric endocrine disorders

Next-generation sequencing (NGS) and array-based comparative genomic hybridization (array CGH) have enabled us to perform high-throughput mutation screening and genome-wide copy number analysis, respectively. These methods can be used for molecular diagnosis of pediatric endocrine disorders. NGS has determined the frequency and phenotypic variation of mutations in several disease-associated genes. Furthermore, whole exome analysis using NGS has successfully identified several novel causative genes for endocrine disorders. Array CGH is currently used as the standard procedure for molecular cytogenetic analysis. Array CGH can detect various submicroscopic genomic rearrangements involving exons or enhancers of disease-associated genes. This review introduces some examples of the use of NGS and array CGH for the molecular diagnosis of pediatric endocrine disorders.

Extensive homoeologous genome exchanges in allopolyploid crops revealed by mRNAseq-based visualization

Polyploidy, the possession of multiple sets of chromosomes, has been a predominant factor in the evolution and success of the angiosperms. Although artificially formed allopolyploids show a high rate of genome rearrangement, the genomes of cultivars and germplasm used for crop breeding were assumed stable and genome structural variation under the artificial selection process of commercial breeding has remained little studied. Here, we show, using a repurposed visualization method based on transcriptome sequence data, that genome structural rearrangement occurs frequently in varieties of three polyploid crops (oilseed rape, mustard rape and bread wheat), meaning that the extent of genome structural variation present in commercial crops is much higher than expected. Exchanges were found to occur most frequently where homoeologous chromosome segments are collinear to telomeres and in material produced as doubled haploids. The new insights into genome structural evolution enable us to reinterpret the results of recent studies and implicate homoeologous exchanges, not deletions, as being responsible for variation controlling important seed quality traits in rapeseed. Having begun to identify the extent of genome structural variation in polyploid crops, we can envisage new strategies for the global challenge of broadening crop genetic diversity and accelerating adaptation, such as the molecular identification and selection of genome deletions or duplications encompassing genes with trait-controlling dosage effects.

A genome-wide study of Hardy-Weinberg equilibrium with next generation sequence data

Statistical tests for Hardy–Weinberg equilibrium have been an important tool for detecting genotyping errors in the past, and remain important in the quality control of next generation sequence data. In this paper, we analyze complete chromosomes of the 1000 genomes project by using exact test procedures for autosomal and X-chromosomal variants. We find that the rate of disequilibrium largely exceeds what might be expected by chance alone for all chromosomes. Observed disequilibrium is, in about 60% of the cases, due to heterozygote excess. We suggest that most excess disequilibrium can be explained by sequencing problems, and hypothesize mechanisms that can explain exceptional heterozygosities. We report higher rates of disequilibrium for the MHC region on chromosome 6, regions flanking centromeres and p-arms of acrocentric chromosomes. We also detected long-range haplotypes and areas with incidental high disequilibrium. We report disequilibrium to be related to read depth, with variants having extreme read depths being more likely to be out of equilibrium. Disequilibrium rates were found to be 11 times higher in segmental duplications and simple tandem repeat regions. The variants with significant disequilibrium are seen to be concentrated in these areas. For next generation sequence data, Hardy–Weinberg disequilibrium seems to be a major indicator for copy number variation.

SeqCNV: a novel method for identification of copy number variations in targeted next-generation sequencing data

Background

Targeted next-generation sequencing (NGS) has been widely used as a cost-effective way to identify the genetic basis of human disorders. Copy number variations (CNVs) contribute significantly to human genomic variability, some of which can lead to disease. However, effective detection of CNVs from targeted capture sequencing data remains challenging.

Results

Here we present SeqCNV, a novel CNV calling method designed to use capture NGS data. SeqCNV extracts the read depth information and utilizes the maximum penalized likelihood estimation (MPLE) model to identify the copy number ratio and CNV boundary. We applied SeqCNV to both bacterial artificial clone (BAC) and human patient NGS data to identify CNVs. These CNVs were validated by array comparative genomic hybridization (aCGH).

Conclusions

SeqCNV is able to robustly identify CNVs of different size using capture NGS data. Compared with other CNV-calling methods, SeqCNV shows a significant improvement in both sensitivity and specificity.

Electronic supplementary material

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

High mutation rates explain low population genetic divergence at copy-number-variable loci in Homo sapiens

Copy-number-variable (CNV) loci differ from single nucleotide polymorphic (SNP) sites in size, mutation rate, and mechanisms of maintenance in natural populations. It is therefore hypothesized that population genetic divergence at CNV loci will differ from that found at SNP sites. Here, we test this hypothesis by analysing 856 CNV loci from the genomes of 1184 healthy individuals from 11 HapMap populations with a wide range of ancestry. The results show that population genetic divergence at the CNV loci is generally more than three times lower than at genome-wide SNP sites. Populations generally exhibit very small genetic divergence (G_st = 0.05 ± 0.049). The smallest divergence is among African populations (G_st = 0.0081 ± 0.0025), with increased divergence among non-African populations (G_st = 0.0217 ± 0.0109) and then among African and non-African populations (G_st = 0.0324 ± 0.0064). Genetic diversity is high in African populations (~0.13), low in Asian populations (~0.11), and intermediate in the remaining 11 populations. Few significant linkage disequilibria (LDs) occur between the genome-wide CNV loci. Patterns of gametic and zygotic LDs indicate the absence of epistasis among CNV loci. Mutation rate is about twice as large as the migration rate in the non-African populations, suggesting that the high mutation rates play dominant roles in producing the low population genetic divergence at CNV loci.

Copy Number Variations in Tilapia Genomes

Discovering the nature and pattern of genome variation is fundamental in understanding phenotypic diversity among populations. Although several millions of single nucleotide polymorphisms (SNPs) have been discovered in tilapia, the genome-wide characterization of larger structural variants, such as copy number variation (CNV) regions has not been carried out yet. We conducted a genome-wide scan for CNVs in 47 individuals from three tilapia populations. Based on 254 Gb of high-quality paired-end sequencing reads, we identified 4642 distinct high-confidence CNVs. These CNVs account for 1.9% (12.411 Mb) of the used Nile tilapia reference genome. A total of 1100 predicted CNVs were found overlapping with exon regions of protein genes. Further association analysis based on linear model regression found 85 CNVs ranging between 300 and 27,000 base pairs significantly associated to population types (R ² > 0.9 and P > 0.001). Our study sheds first insights on genome-wide CNVs in tilapia. These CNVs among and within tilapia populations may have functional effects on phenotypes and specific adaptation to particular environments.

Microduplication of 7q36.3 encompassing the SHH long‑range regulator (ZRS) in a patient with triphalangeal thumb‑polysyndactyly syndrome and congenital heart disease

Triphalangeal thumb‑polysyndactyly syndrome (TPT‑PS) is an autosomal dominant disorder with complete penetrance and a variable expression consisting of opposable triphalangeal thumbs, duplication of the distal thumb phalanx, pre‑axial polydactyly and duplication of the big toes (hallux). The causative gene of TPT‑PS has been mapped to 7q36.3. Sonic hedgehog (SHH) expressed in the zone of polarizing activity (ZPA) has an important role in defining the anterior‑posterior axis and numbers of digits in limb bud development. Point mutation or duplication in the ZPA regulatory sequence (ZRS), a cis‑regulator of SHH, will lead to TPT‑PS. The present study describes a 1‑year‑old female congenital heart disease (CHD) patient with TPT‑PS phenotype. In this Han Chinese family with TPT‑PS, high resolution single nucleotide polymorphism array technology identified a novel 0.29 Mb duplication comprising ZRS at 7q36.3 where LMBR1 is located. Additionally, a novel deletion of 22q11.21 was detected in the proband with Tetralogy of Fallot. However, the parents and other relatives of the patient did not harbor this genomic lesion nor CHD. The findings supported the hypothesis that an increased copy number variation of ZRS is the genetic mechanism underlying the phenotype of TPT‑PS, and corroborated that 22q11.21 deletion is a genetic cause of CHD.

Application of array-comparative genomic hybridization in tetralogy of Fallot

To explore the underlying pathogenesis and provide references for genetic counseling and prenatal gene diagnosis, we analyzed the chromosome karyotypes and genome-wide copy number variations (CNVs) in 86 patients with tetralogy of Fallot (TOF) by G-banding karyotype analysis and array-comparative genomic hybridization (aCGH), respectively. And then quantitative polymerase chain reaction was used to validate these candidate CNVs. Based on their different properties, CNVs were categorized into benign CNVs, suspiciously pathogenic CNVs, and indefinite CNVs. Data analysis was based on public databases such as UCSC, DECIPHER, DGV, ISCA, and OMIM.The karyotype was normal in all the 86 patients with TOF. CNVs were detected in 11 patients by aCGH and quantitative polymerase chain reaction. Patient no. 0001, 0010, and 0029 had 2.52-Mb deletion in the chromosome 22q11.21 region; patient no. 0008 had both 595- and 428-kb duplications, respectively, in 12p12.3p12.2 and 14q23.2q23.3 regions; patient no. 0009 had 1.46-Mb duplication in the 1q21.1q21.2 region; patient no. 0016 had 513-kb duplication in the 1q42.13 region; patient no. 0024 had 292-kb duplication in the 16q11.2 region; patient no. 0026 had 270-kb duplication in the 16q24.1 region; patient no. 0028 had 222-kb deletion in the 7q31.1 region; patient no. 0033 had 1.73-Mb duplication in the 17q12 region; and patient no. 0061 had 5.79-Mb deletion in the 1p36.33p36.31 region.aCGH can accurately detect CNVs in the patients with TOF. This is conducive to genetic counseling and prenatal diagnosis for TOF and provides a new clue and theoretical basis for exploring the pathogenesis of congenital heart disease.

UGT2B28 genomic variation is associated with hepatitis B e-antigen seroconversion in response to antiviral therapy

Seroconversion of hepatitis B virus (HBV) e-antigen (HBeAg) is a critical but often-missed therapeutic goal in standard antiviral treatments. An extreme-phenotype genome-wide association study was performed, comparing untreated spontaneous recoverers (with seroconversion of HBV surface antigen) versus entecavir-treated patients failing to achieve HBeAg seroconversion. A single-nucleotide-polymorphism rs2132039 on the UGT2B28 gene, alongside an adjacent copy number polymorphism (CNP605), manifested the strongest clinical associations (P = 3.4 × 10-8 and 0.001, respectively). Multivariate analysis showed that rs2132039-TT genotypes, but not CNP605 copy numbers, remained associated to spontaneous recoverers (P = 0.009). The clinical association of rs2132039 was validated successfully in an independent cohort (n = 302; P = 0.002). Longitudinal case-only analyses revealed that the rs2132039-TT genotype predicted shorter time-to-HBeAg-seroconversion in all antiviral-treated patients (n = 380, P = 0.012), as well as the peginterferon-treated subgroup (n = 123; P = 0.024, Hazard ratio [HR] = 2.104, Confidence interval [CI] = 1.105-4.007). In the entecavir-treated subgroup, the predictive effect was restricted by pretreatment alanine aminotransferase (ALT) and aspartate aminotransferase (AST) levels, with effective prediction observed in patients with ALT < 200 IU/ml and ALT/AST ratio <2 (n = 132; P = 0.013, HR = 10.538, CI = 1.420-78.196).

Dealing with Gene-Dosage Imbalance during S Phase

DNA replication perturbs the dosage balance between genes that replicate early during S phase and those that replicate late. If propagated to influence protein content, this dosage imbalance could influence cellular functions. In bacteria, mechanisms have evolved to use this imbalance to tune certain processes with the rate of cell growth. By contrast, eukaryotes buffer this dosage imbalance to ensure gene expression homeostasis also during S phase. Here, we outline classical and more recent studies describing how different organisms deal with this replication-dependent dosage imbalance, and describe recent results linking the eukaryotic buffering mechanism to replication-dependent histone acetylation. Finally, we discuss the possible implications of this buffering mechanism and speculate why it is specific to eukaryote cells.

Testing neutrality at copy-number-variable loci under the finite-allele and finite-site models

Copy-number variation (CNV) is an important form of DNA structural variation because a certain proportion of genomes in many eukaryotic species can contribute to such variations. Owing to the differences between CNVs and single nucleotide polymorphisms (SNPs) in size, mutation rate and maintaining mechanism, it is more realistic to characterize CNV evolution under the finite-allele and finite-site models. Here, we propose a method to test multiple CNVs neutrality under the finite-allele and finite-site models and the assumption of mutation-drift process. The statistical property of the method is evaluated through Monte Carlo simulations under the effects of the sample size, the scaled mutation rates, the number of CNVs, the population demographic change, and selection. Different from Tajima's D test, a bootstrap or a permutation approach is suggested to conduct a neutrality test. Application of this method is illustrated using the diploid CNV genotypes measured in discrete copy numbers in 11 HapMap phase III populations. The results show that the mutation-drift process can explain the variation of genome-wide CNVs among 1184 individuals (856 CNVs, ∼0.02Mb on average in size), irrespective of the historical demographic changes. Patterns from allele-frequency-spectrum analysis also support the hypothesis of neutral CNVs. Our results suggest that most human chromosomal changes in healthy individuals via unbalanced rearrangements of the segments with certain sizes are neutral.

Systematic reanalysis of partial trisomy 21 cases with or without Down syndrome suggests a small region on 21q22.13 as critical to the phenotype

A 'Down Syndrome critical region' (DSCR) sufficient to induce the most constant phenotypes of Down syndrome (DS) had been identified by studying partial (segmental) trisomy 21 (PT21) as an interval of 0.6-8.3 Mb within human chromosome 21 (Hsa21), although its existence was later questioned. We propose an innovative, systematic reanalysis of all described PT21 cases (from 1973 to 2015). In particular, we built an integrated, comparative map from 125 cases with or without DS fulfilling stringent cytogenetic and clinical criteria. The map allowed to define or exclude as candidates for DS fine Hsa21 sequence intervals, also integrating duplication copy number variants (CNVs) data. A highly restricted DSCR (HR-DSCR) of only 34 kb on distal 21q22.13 has been identified as the minimal region whose duplication is shared by all DS subjects and is absent in all non-DS subjects. Also being spared by any duplication CNV in healthy subjects, HR-DSCR is proposed as a candidate for the typical DS features, the intellectual disability and some facial phenotypes. HR-DSCR contains no known gene and has relevant homology only to the chimpanzee genome. Searching for HR-DSCR functional loci might become a priority for understanding the fundamental genotype-phenotype relationships in DS.

Fast Bayesian Inference of Copy Number Variants using Hidden Markov Models with Wavelet Compression

By integrating Haar wavelets with Hidden Markov Models, we achieve drastically reduced running times for Bayesian inference using Forward-Backward Gibbs sampling. We show that this improves detection of genomic copy number variants (CNV) in array CGH experiments compared to the state-of-the-art, including standard Gibbs sampling. The method concentrates computational effort on chromosomal segments which are difficult to call, by dynamically and adaptively recomputing consecutive blocks of observations likely to share a copy number. This makes routine diagnostic use and re-analysis of legacy data collections feasible; to this end, we also propose an effective automatic prior. An open source software implementation of our method is available at http://schlieplab.org/Software/HaMMLET/ (DOI: 10.5281/zenodo.46262). This paper was selected for oral presentation at RECOMB 2016, and an abstract is published in the conference proceedings.

Copy number variation identification and analysis of the chicken genome using a 60K SNP BeadChip

Copy number variation (CNV) is an important source of genetic variation in organisms and a main factor that affects phenotypic variation. A comprehensive study of chicken CNV can provide valuable information on genetic diversity and facilitate future analyses of associations between CNV and economically important traits in chickens. In the present study, an F2 full-sib chicken population (554 individuals), established from a cross between Xinghua and White Recessive Rock chickens, was used to explore CNV in the chicken genome. Genotyping was performed using a chicken 60K SNP BeadChip. A total of 1,875 CNV were detected with the PennCNV algorithm, and the average number of CNV was 3.42 per individual. The CNV were distributed across 383 independent CNV regions (CNVR) and covered 41 megabases (3.97%) of the chicken genome. Seven CNVR in 108 individuals were validated by quantitative real-time PCR, and 81 of these individuals (75%) also were detected with the PennCNV algorithm. In total, 274 CNVR (71.54%) identified in the current study were previously reported. Of these, 147 (38.38%) were reported in at least 2 studies. Additionally, 109 of the CNVR (28.46%) discovered here are novel. A total of 709 genes within or overlapping with the CNVR was retrieved. Out of the 2,742 quantitative trait loci (QTL) collected in the chicken QTL database, 43 QTL had confidence intervals overlapping with the CNVR, and 32 CNVR encompassed one or more functional genes. The functional genes located in the CNVR are likely to be the QTG that are associated with underlying economic traits. This study considerably expands our insight into the structural variation in the genome of chickens and provides an important resource for genomic variation, especially for genomic structural variation related to economic traits in chickens.

An integrated analysis tool for analyzing hybridization intensities and genotypes using new-generation population-optimized human arrays

BACKGROUND

Affymetrix Axiom single nucleotide polymorphism (SNP) arrays provide a cost-effective, high-density, and high-throughput genotyping solution for population-optimized analyses. However, no public software is available for the integrated genomic analysis of hybridization intensities and genotypes for this new-generation population-optimized genotyping platform.

RESULTS

A set of statistical methods was developed for an integrated analysis of allele frequency (AF), allelic imbalance (AI), loss of heterozygosity (LOH), long contiguous stretch of homozygosity (LCSH), and copy number variation or alteration (CNV/CNA) on the basis of SNP probe hybridization intensities and genotypes. This study analyzed 3,236 samples that were genotyped using different SNP platforms. The proposed AF adjustment method considerably increased the accuracy of AF estimation. The proposed quick circular binary segmentation algorithm for segmenting copy number reduced the computation time of the original segmentation method by 30-67 %. The proposed CNV/CNA detection, which integrates AI and LOH/LCSH detection, had a promising true positive rate and well-controlled false positive rate in simulation studies. Moreover, our real-time quantitative polymerase chain reaction experiments successfully validated the CNVs/CNAs that were identified in the Axiom data analyses using the proposed methods; some of the validated CNVs/CNAs were not detected in the Affymetrix Array 6.0 data analysis using the Affymetrix Genotyping Console. All the analysis functions are packaged into the ALICE (AF/LOH/LCSH/AI/CNV/CNA Enterprise) software.

CONCLUSIONS

ALICE and the used genomic reference databases, which can be downloaded from http://hcyang.stat.sinica.edu.tw/software/ALICE.html , are useful resources for analyzing genomic data from the Axiom and other SNP arrays.

Analysis of genomic copy number variation in equine recurrent airway obstruction (heaves)

We explored the involvement of genomic copy number variants (CNVs) in susceptibility to recurrent airway obstruction (RAO), or heaves-an asthmalike inflammatory disease in horses. Analysis of 16 RAO-susceptible (cases) and six RAO-resistant (control) horses on a custom-made whole-genome 400K equine tiling array identified 245 CNV regions (CNVRs), 197 previously known and 48 new, distributed on all horse autosomes and the X chromosome. Among the new CNVRs, 30 were exclusively found in RAO cases and were further analyzed by quantitative PCR, including additional cases and controls. Suggestive association (P = 0.03; corrected P = 0.06) was found between RAO and a loss on chromosome 5 involving NME7, a gene necessary for ciliary functions in lungs and involved in primary ciliary dyskinesia in humans. The CNVR could be a potential marker for RAO susceptibility but needs further study in additional RAO cohorts. Other CNVRs were not associated with RAO, although several involved genes of interest, such as SPI2/SERPINA1 from the serpin gene family, which are associated with chronic obstructive pulmonary disease and asthma in humans. The SPI2/SERPINA1 CNVR showed striking variation among horses, but it was not significantly different between RAO cases and controls. The findings provide baseline information on the relationship between CNVs and RAO susceptibility. Discovery of new CNVs and the use of a larger population of RAO-affected and control horses are needed to shed more light on their significance in modulating this complex and heterogeneous disease.

Enhancing ergosterol production in Pichia pastoris GS115 by overexpressing squalene synthase gene from Glycyrrhiza uralensis

The present study was designed to determine the effects of copy number variations (CNVs) of squalene synthase 1(SQS1) gene on the mevalonate (MVA) pathway. SQS1 gene from G. uralensis (GuSQS1) was cloned and over-expressed in Pichia pastoris GS115. Six recombinant P. pastoris strains containing different copy number of GuSQS1 were constructed. HPLC was used to assay the level of ergosterol in all transgenic P. pastoris strains containing GuSQS1. HPLC analysis showed that the contents of ergosterol in all of the transgenic P. pastoris containing GuSQS1 were higher than that in the negative control. And with the increase of copy number of GuSQS1, the content of ergosterol showed an increasing-decreasing-increasing pattern. The contents of ergosterol in 10-copy-GuSQS1 P. pastoris and 47-copy-GuSQS1 P. pastoris were significantly higher than that in the rest recombinant P. pastoris strains. In conclusion, the CNVs of GuSQS1 influence the content of secondary metabolites in the MVA pathway. The present study provides a basis for over-expressing GuSQS1 and increasing the content of glycyrrhizin in G. uralensis cultivars.

Can we afford to sequence every newborn baby's genome?

Whole-exome sequencing and whole-genome sequencing are gradually entering into the clinical arena. Drops in sequencing prices have led some to suggest that these analyses could be extended to the screening of whole populations or subsets thereof. Herein, we argue that this optimism is presently still unfounded. While cost estimates take into account the generation of sequence data, they fail to properly evaluate both the price of accurate and efficient interpretation and of the proper return of genomic information to the consulting individuals. Thus, short of inventing new, cost-effective ways of achieving these goals, the latter are likely to ruin our healthcare systems. We posit that due to lack of available resources, generalization of this practice remains, for the time being, unrealistic.

Noninvasive detection of fetal subchromosomal abnormalities by semiconductor sequencing of maternal plasma DNA

Noninvasive prenatal testing (NIPT) using sequencing of fetal cell-free DNA from maternal plasma has enabled accurate prenatal diagnosis of aneuploidy and become increasingly accepted in clinical practice. We investigated whether NIPT using semiconductor sequencing platform (SSP) could reliably detect subchromosomal deletions/duplications in women carrying high-risk fetuses. We first showed that increasing concentration of abnormal DNA and sequencing depth improved detection. Subsequently, we analyzed plasma from 1,456 pregnant women to develop a method for estimating fetal DNA concentration based on the size distribution of DNA fragments. Finally, we collected plasma from 1,476 pregnant women with fetal structural abnormalities detected on ultrasound who also underwent an invasive diagnostic procedure. We used SSP of maternal plasma DNA to detect subchromosomal abnormalities and validated our results with array comparative genomic hybridization (aCGH). With 3.5 million reads, SSP detected 56 of 78 (71.8%) subchromosomal abnormalities detected by aCGH. With increased sequencing depth up to 10 million reads and restriction of the size of abnormalities to more than 1 Mb, sensitivity improved to 69 of 73 (94.5%). Of 55 false-positive samples, 35 were caused by deletions/duplications present in maternal DNA, indicating the necessity of a validation test to exclude maternal karyotype abnormalities. This study shows that detection of fetal subchromosomal abnormalities is a viable extension of NIPT based on SSP. Although we focused on the application of cell-free DNA sequencing for NIPT, we believe that this method has broader applications for genetic diagnosis, such as analysis of circulating tumor DNA for detection of cancer.

Detection of CAPN10 copy number variation in Thai patients with type 2 diabetes by denaturing high performance liquid chromatography and real-time quantitative polymerase chain reaction

Aims/Introduction

A combination of multiple genetic and environmental factors contribute to the pathogenesis of type 2 diabetes. Copy number variations (CNVs) are associated with complex human diseases. However, CNVs can cause genotype deviation from the Hardy–Weinberg equilibrium (HWE). A genetic case–control association study in 216 Thai diabetic patients and 192 non-diabetic controls found that, after excluding genotyping errors, genotype distribution of calpain 10 (CAPN10) SNP44 (rs2975760) deviated from HWE. Here, we aimed to detect CNV within the CAPN10 SNP44 region.

Materials and Methods

CNV within the CAPN10 SNP44 region was detected using denaturing high-performance liquid chromatography, and the results confirmed by real-time quantitative polymerase chain reaction with SYBR Green I.

Results

Both methods successfully identified CNV in the CAPN10 SNP44 region, obtaining concordant results. Correction of genotype calling based on the status of identified CNVs showed that the CAPN10 SNP44 genotype is in good agreement with HWE (P > 0.05). However, no association between CNV genotypes and risk of type 2 diabetes was observed.

Conclusions

Identified CNVs for CAPN10 SNP44 genotypes lead to deviation from HWE. Furthermore, both denaturing high-performance liquid chromatography and real-time quantitative polymerase chain reaction are useful for detecting CNVs.

Copy number variation in exportin-4 (XPO4) gene and its association with histological severity of non-alcoholic fatty liver disease

A recent genome-wide copy number (CNV) scan identified a 13q12.11 duplication in the exportin-4 (XPO4) gene to be associated with non-alcoholic steatohepatitis (NASH). We sought to confirm the finding in a larger cohort and to assess the serum XPO4 pattern in a broad spectrum of non-alcoholic fatty liver disease (NAFLD) cases. We analysed 249 NAFLD patients and 232 matched controls using TaqMan assay and serum XPO4 was measured. Copy number distribution was as follows: copy number neutral (NAFLD: 53.8%, controls: 68.6%), copy number losses (NAFLD: 13.3%, controls: 12.9%), copy number gains (NAFLD: 32.9%, controls: 18.5%). CNV gain was significantly associated with a greater risk of NAFLD (adjusted OR 2.22, 95% CI 1.42–3.46, P = 0.0004) and NASH (adjusted OR 2.33, 95% CI 1.47–3.68, P = 0.0003). Interestingly, subjects carrying extra copy number showed significantly higher serum ALT and triglyceride (P < 0.05). Serum XPO4 levels progressively declined (P = 0.043) from controls (24.6 ng/mL) to simple steatosis (20.8 ng/mL) to NASH (13.8 ng/mL). In conclusion, XPO4 CNV duplication was associated with histological severity of NAFLD, and accompanied by changes in serum XPO4 levels providing insights into NAFLD pathogenesis, and has the potential for biomarker development.

Copy number variations and gene polymorphisms of complement components in ocular Behcet's disease and Vogt-Koyanagi-Harada syndrome

Complement is involved in many immune-mediated diseases. However, the association of its copy number variations (CNVs) and polymorphisms with Behcet’s disease (BD) and Vogt-Koyanagi-Harada syndrome (VKH) is unknown. We examined copy number and mRNA expression by real-time PCR. Cytokine production by stimulated peripheral blood mononuclear cells (PBMCs) in genotyped individuals was measured by ELISA. The frequencies of having more than two copies of C3 were significantly increased in BD and VKH, whereas CNV of C5 was only associated with BD. Increased frequencies of the GG genotype of C3 rs408290 and C5 rs2269067 were found in BD. No association was observed between C3 and C5 SNPs and VKH. mRNA expression in the high CNV group and GG cases of C3 and C5 was significantly higher compared to other genotypes. Increased interleukin-17 and IFN-γ was observed in the high CNV group and GG genotype cases of C3. Interleukin-17 but not IFN-γ was increased in the high CNV group and GG genotype cases of C5. No effect of C3 or C5 genetic variants was seen on the production of TNF-α, IL-10, IL-1β, MCP-1, IL-6 and IL-8. Our study thus provides further evidence for a role of complement in the pathogenesis of uveitis.

Production of Extrachromosomal MicroDNAs Is Linked to Mismatch Repair Pathways and Transcriptional Activity

MicroDNAs are <400-base extrachromosomal circles found in mammalian cells. Tens of thousands of microDNAs have been found in all tissue types, including sperm. MicroDNAs arise preferentially from areas with high gene density, GC content, and exon density from promoters with activating chromatin modifications and in sperm from the 5'-UTR of full-length LINE-1 elements, but are depleted from lamin-associated heterochromatin. Analysis of microDNAs from a set of human cancer cell lines revealed lineage-specific patterns of microDNA origins. A survey of microDNAs from chicken cells defective in various DNA repair proteins reveals that homologous recombination and non-homologous end joining repair pathways are not required for microDNA production. Deletion of the MSH3 DNA mismatch repair protein results in a significant decrease in microDNA abundance, specifically from non-CpG genomic regions. Thus, microDNAs arise as part of normal cellular physiology—either from DNA breaks associated with RNA metabolism or from replication slippage followed by mismatch repair.

Identification of an Alu-repeat-mediated deletion of OPTN upstream region in a patient with a complex ocular phenotype

Genetic causes of ocular conditions remain largely unknown. To reveal the molecular basis for a congenital ocular phenotype associated with glaucoma we performed whole‐exome sequencing (WES) and whole‐genome copy number analyses of patient DNA. WES did not identify a causative variant. Copy number variation analysis identified a deletion of 10p13 in the patient and his unaffected father; the deletion breakpoint contained a single 37‐bp sequence that is normally present in two distinct Alu repeats separated by ~181 kb. The deletion removed part of the upstream region of optineurin (OPTN) as well as the upstream sequence and two coding exons of coiled‐coil domain containing 3 (CCDC3); analysis of the patient's second allele showed normal OPTN and CCDC3 sequences. Studies of zebrafish orthologs identified expression in the developing eye for both genes. OPTN is a known factor in dominant adult‐onset glaucoma and Amyotrophic Lateral Sclerosis (ALS). The deletion eliminates 98 kb of the OPTN upstream sequence leaving only ~1 kb of the proximal promoter region. Comparison of transcriptional activation capability of the 3 kb normal and the rearranged del(10)(p13) OPTN promoter sequences demonstrated a statistically significant decrease for the deleted allele; sequence analysis of the entire deleted region identified multiple conserved elements with possible cis‐regulatory activity. Additional screening of CCDC3 indicated that heterozygous loss‐of‐function alleles are unlikely to cause congenital ocular disease. In summary, we report the first regulatory region deletion involving OPTN, caused by Alu‐mediated nonallelic homologous recombination and possibly contributing to the patient's ocular phenotype. In addition, our data indicate that Alu‐mediated rearrangements of the OPTN upstream region may represent a new source of affected alleles in human conditions. Evaluation of the upstream OPTN sequences in additional ocular and ALS patients may help to determine the role of this region, if any, in human disease.

Modeling a model: Mouse genetics, 22q11.2 Deletion Syndrome, and disorders of cortical circuit development

Understanding the developmental etiology of autistic spectrum disorders, attention deficit/hyperactivity disorder and schizophrenia remains a major challenge for establishing new diagnostic and therapeutic approaches to these common, difficult-to-treat diseases that compromise neural circuits in the cerebral cortex. One aspect of this challenge is the breadth and overlap of ASD, ADHD, and SCZ deficits; another is the complexity of mutations associated with each, and a third is the difficulty of analyzing disrupted development in at-risk or affected human fetuses. The identification of distinct genetic syndromes that include behavioral deficits similar to those in ASD, ADHC and SCZ provides a critical starting point for meeting this challenge. We summarize clinical and behavioral impairments in children and adults with one such genetic syndrome, the 22q11.2 Deletion Syndrome, routinely called 22q11DS, caused by micro-deletions of between 1.5 and 3.0 MB on human chromosome 22. Among many syndromic features, including cardiovascular and craniofacial anomalies, 22q11DS patients have a high incidence of brain structural, functional, and behavioral deficits that reflect cerebral cortical dysfunction and fall within the spectrum that defines ASD, ADHD, and SCZ. We show that developmental pathogenesis underlying this apparent genetic "model" syndrome in patients can be defined and analyzed mechanistically using genomically accurate mouse models of the deletion that causes 22q11DS. We conclude that "modeling a model", in this case 22q11DS as a model for idiopathic ASD, ADHD and SCZ, as well as other behavioral disorders like anxiety frequently seen in 22q11DS patients, in genetically engineered mice provides a foundation for understanding the causes and improving diagnosis and therapy for these disorders of cortical circuit development.

Decoding the complex genetic causes of heart diseases using systems biology

The pace of disease gene discovery is still much slower than expected, even with the use of cost-effective DNA sequencing and genotyping technologies. It is increasingly clear that many inherited heart diseases have a more complex polygenic aetiology than previously thought. Understanding the role of gene-gene interactions, epigenetics, and non-coding regulatory regions is becoming increasingly critical in predicting the functional consequences of genetic mutations identified by genome-wide association studies and whole-genome or exome sequencing. A systems biology approach is now being widely employed to systematically discover genes that are involved in heart diseases in humans or relevant animal models through bioinformatics. The overarching premise is that the integration of high-quality causal gene regulatory networks (GRNs), genomics, epigenomics, transcriptomics and other genome-wide data will greatly accelerate the discovery of the complex genetic causes of congenital and complex heart diseases. This review summarises state-of-the-art genomic and bioinformatics techniques that are used in accelerating the pace of disease gene discovery in heart diseases. Accompanying this review, we provide an interactive web-resource for systems biology analysis of mammalian heart development and diseases, CardiacCode ( http://CardiacCode.victorchang.edu.au/ ). CardiacCode features a dataset of over 700 pieces of manually curated genetic or molecular perturbation data, which enables the inference of a cardiac-specific GRN of 280 regulatory relationships between 33 regulator genes and 129 target genes. We believe this growing resource will fill an urgent unmet need to fully realise the true potential of predictive and personalised genomic medicine in tackling human heart disease.