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.     

Endometriosis Is Associated with Rare Copy Number Variants

Endometriosis is a complex gynecological condition that affects 6–10% of women in their reproductive years and is defined by the presence of endometrial glands and stroma outside the uterus. Twin, family, and genome-wide association (GWA) studies have confirmed a genetic role, yet only a small part of the genetic risk can be explained by SNP variation. Copy number variants (CNVs) account for a greater portion of human genetic variation than SNPs and include more recent mutations of large effect. CNVs, likely to be prominent in conditions with decreased reproductive fitness, have not previously been examined as a genetic contributor to endometriosis. Here we employ a high-density genotyping microarray in a genome-wide survey of CNVs in a case-control population that includes 2,126 surgically confirmed endometriosis cases and 17,974 population controls of European ancestry. We apply stringent quality filters to reduce the false positive rate common to many CNV-detection algorithms from 77.7% to 7.3% without noticeable reduction in the true positive rate. We detected no differences in the CNV landscape between cases and controls on the global level which showed an average of 1.92 CNVs per individual with an average size of 142.3 kb. On the local level we identify 22 CNV-regions at the nominal significance threshold (P<0.05), which is greater than the 8.15 CNV-regions expected based on permutation analysis (P<0.001). Three CNV''s passed a genome-wide P-value threshold of 9.3×10⁻⁴; a deletion at SGCZ on 8p22 (P = 7.3×10⁻⁴, OR = 8.5, Cl = 2.3–31.7), a deletion in MALRD1 on 10p12.31 (P = 5.6×10⁻⁴, OR = 14.1, Cl = 2.7–90.9), and a deletion at 11q14.1 (P = 5.7×10⁻⁴, OR = 33.8, Cl = 3.3–1651). Two SNPs within the 22 CNVRs show significant genotypic association with endometriosis after adjusting for multiple testing; rs758316 in DPP6 on 7q36.2 (P = 0.0045) and rs4837864 in ASTN2 on 9q33.1 (P = 0.0002). Together, the CNV-loci are detected in 6.9% of affected women compared to 2.1% in the general population.     

In Ovo Injection of Betaine Affects Hepatic Cholesterol Metabolism through Epigenetic Gene Regulation in Newly Hatched Chicks

Betaine is reported to regulate hepatic cholesterol metabolism in mammals. Chicken eggs contain considerable amount of betaine, yet it remains unknown whether and how betaine in the egg affects hepatic cholesterol metabolism in chicks. In this study, eggs were injected with betaine at 2.5 mg/egg and the hepatic cholesterol metabolism was investigated in newly hatched chicks. Betaine did not affect body weight or liver weight, but significantly increased the serum concentration (P < 0.05) and the hepatic content (P < 0.01) of cholesterol. Accordingly, the cholesterol biosynthetic enzyme HMGCR was up-regulated (P < 0.05 for both mRNA and protein), while CYP7A1 which converts cholesterol to bile acids was down-regulated (P < 0.05 for mRNA and P = 0.07 for protein). Moreover, hepatic protein content of the sterol-regulatory element binding protein 1 which regulates cholesterol and lipid biosynthesis, and the mRNA abundance of ATP binding cassette sub-family A member 1 (ABCA1) which mediates cholesterol counter transport were significantly (P < 0.05) increased in betaine-treated chicks. Meanwhile, hepatic protein contents of DNA methyltransferases 1 and adenosylhomocysteinase-like 1 were increased (P < 0.05), which was associated with global genomic DNA hypermethylation (P < 0.05) and diminished gene repression mark histone H3 lysine 27 trimethylation (P < 0.05). Furthermore, CpG methylation level on gene promoters was found to be increased (P < 0.05) for CYP7A1 yet decreased (P < 0.05) for ABCA1. These results indicate that in ovo betaine injection regulates hepatic cholesterol metabolism in chicks through epigenetic mechanisms including DNA and histone methylations.     

Genome-Wide Identification of Copy Number Variations in Chinese Holstein

Recent studies of mammalian genomes have uncovered the vast extent of copy number variations (CNVs) that contribute to phenotypic diversity. Compared to SNP, a CNV can cover a wider chromosome region, which may potentially incur substantial sequence changes and induce more significant effects on phenotypes. CNV has been becoming an alternative promising genetic marker in the field of genetic analyses. Here we firstly report an account of CNV regions in the cattle genome in Chinese Holstein population. The Illumina Bovine SNP50K Beadchips were used for screening 2047 Holstein individuals. Three different programes (PennCNV, cnvPartition and GADA) were implemented to detect potential CNVs. After a strict CNV calling pipeline, a total of 99 CNV regions were identified in cattle genome. These CNV regions cover 23.24 Mb in total with an average size of 151.69 Kb. 52 out of these CNV regions have frequencies of above 1%. 51 out of these CNV regions completely or partially overlap with 138 cattle genes, which are significantly enriched for specific biological functions, such as signaling pathway, sensory perception response and cellular processes. The results provide valuable information for constructing a more comprehensive CNV map in the cattle genome and offer an important resource for investigation of genome structure and genomic variation underlying traits of interest in cattle.     

Mitochondrial DNA Copy Number Is Associated with Breast Cancer Risk

Mitochondrial DNA (mtDNA) copy number in peripheral blood is associated with increased risk of several cancers. However, data from prospective studies on mtDNA copy number and breast cancer risk are lacking. We evaluated the association between mtDNA copy number in peripheral blood and breast cancer risk in a nested case-control study of 183 breast cancer cases with pre-diagnostic blood samples and 529 individually matched controls among participants of the Singapore Chinese Health Study. The mtDNA copy number was measured using real time PCR. Conditional logistic regression analyses showed that there was an overall positive association between mtDNA copy number and breast cancer risk (P_trend = 0.01). The elevated risk for higher mtDNA copy numbers was primarily seen for women with <3 years between blood draw and cancer diagnosis; ORs (95% CIs) for 2^nd, 3^rd, 4^th, and 5^th quintile of mtDNA copy number were 1.52 (0.61, 3.82), 2.52 (1.03, 6.12), 3.12 (1.31, 7.43), and 3.06 (1.25, 7.47), respectively, compared with the 1^st quintile (P_trend = 0.004). There was no association between mtDNA copy number and breast cancer risk among women who donated a blood sample ≥3 years before breast cancer diagnosis (P_trend = 0.41). This study supports a prospective association between increased mtDNA copy number and breast cancer risk that is dependent on the time interval between blood collection and breast cancer diagnosis. Future studies are warranted to confirm these findings and to elucidate the biological role of mtDNA copy number in breast cancer risk.     

Identification of Copy Number Variations in Xiang and Kele Pigs

Xiang and Kele pigs are two well-known local Chinese pig breeds that possess rich genetic resources and have enormous economic and scientific value. We performed a comprehensive genomic analysis of the copy number variations (CNVs) in these breeds. CNVs are one of the most important forms of genomic variation and have profound effects on phenotypic variation. In this study, PorcineSNP60 genotyping data from 98 Xiang pigs and 22 Kele pigs were used to identify CNVs. In total, 172 candidate CNV regions (CNVRs) were identified, ranging from 3.19 kb to 8175.26 kb and covering 80.41 Mb of the pig genome. Approximately 56.40% (97/172) of the CNVRs overlapped with those identified in seven previous studies, and 43.60% (75/172) of the identified CNVRs were novel. Of the identified CNVRs, 82 (47 gain, 33 loss, and two gain-loss events that covered 4.58 Mb of the pig genome) were found only in a Xiang population with a large litter size. In contrast, 13 CNVRs (8 gain and 5 loss events) were unique to a Xiang population with small litter sizes, and 30 CNVRs (14 loss and 16 gain events) were unique to Kele pigs. The CNVRs span approximately 660 annotated Sus scrofa genes that are significantly enriched for specific biological functions, such as sensory perception, cognition, reproduction, ATP biosynthetic processes, and neurological processes. Many CNVR-associated genes, particularly the genes involved in reproductive traits, differed between the Xiang populations with large and small litter sizes, and these genes warrant further investigation due to their importance in determining the reproductive performance of Xiang pigs. Our results provide meaningful information about genomic variation, which may be useful in future assessments of the associations between CNVs and important phenotypes in Xiang and Kele pigs to ultimately help protect these rare breeds.     

Proteomic Changes Resulting from Gene Copy Number Variations in Cancer Cells

Along the transformation process, cells accumulate DNA aberrations, including mutations, translocations, amplifications, and deletions. Despite numerous studies, the overall effects of amplifications and deletions on the end point of gene expression—the level of proteins—is generally unknown. Here we use large-scale and high-resolution proteomics combined with gene copy number analysis to investigate in a global manner to what extent these genomic changes have a proteomic output and therefore the ability to affect cellular transformation. We accurately measure expression levels of 6,735 proteins and directly compare them to the gene copy number. We find that the average effect of these alterations on the protein expression is only a few percent. Nevertheless, by using a novel algorithm, we find the combined impact that many of these regional chromosomal aberrations have at the protein level. We show that proteins encoded by amplified oncogenes are often overexpressed, while adjacent amplified genes, which presumably do not promote growth and survival, are attenuated. Furthermore, regulation of biological processes and molecular complexes is independent of general copy number changes. By connecting the primary genome alteration to their proteomic consequences, this approach helps to interpret the data from large-scale cancer genomics efforts.     

Enhancer Chip: Detecting Human Copy Number Variations in Regulatory Elements

Critical functional properties are embedded in the non-coding portion of the human genome. Recent successful studies have shown that variations in distant-acting gene enhancer sequences can contribute to disease. In fact, various disorders, such as thalassaemias, preaxial polydactyly or susceptibility to Hirschsprung’s disease, may be the result of rearrangements of enhancer elements. We have analyzed the distribution of enhancer loci in the genome and compared their localization to that of previously described copy-number variations (CNVs). These data suggest a negative selection of copy number variable enhancers. To identify CNVs covering enhancer elements, we have developed a simple and cost-effective test. Here we describe the gene selection, design strategy and experimental validation of a customized oligonucleotide Array-Based Comparative Genomic Hybridization (aCGH), designated Enhancer Chip. It has been designed to investigate CNVs, allowing the analysis of all the genome with a 300 Kb resolution and specific disease regions (telomeres, centromeres and selected disease loci) at a tenfold higher resolution. Moreover, this is the first aCGH able to test over 1,250 enhancers, in order to investigate their potential pathogenic role. Validation experiments have demonstrated that Enhancer Chip efficiently detects duplications and deletions covering enhancer loci, demonstrating that it is a powerful instrument to detect and characterize copy number variable enhancers.     

Leukocyte Mitochondrial DNA Copy Number Is Associated with Chronic Obstructive Pulmonary Disease

Background

Oxidative stress is known to be involved in the pathogenesis of chronic obstructive pulmonary disease (COPD). Evidence suggests that leukocytes mitochondria DNA (mtDNA) is susceptible to undergo mutations, insertions, or depletion in response to reactive oxidative stress (ROS). We hypothesize that mtDNA copy number is associated with the development of COPD.

Methodology/Principal Findings

Relative mtDNA copy number was measured by a quantitative real-time PCR assay using DNA extracted from peripheral leukocytes. MtDNA copy number of peripheral leukocytes in the COPD group (n = 86) is significantly decreased compared with non-smoker group (n = 77) (250.3± 21.5 VS. 464.2± 49.9, P<0.001). MtDNA copy number in the COPD group was less than that in the healthy smoking group, but P value nearly achieved significance (250.3± 21.5 VS. 404.0± 76.7, P = 0.08) MtDNA copy number has no significance with age, gender, body mass index, current smoking, and pack-years in COPD group, healthy smoker group and no smoker group, respectively. Serum glutathione level in the COPD group is significantly decreased compared with healthy smoker and non-smoker groups (4.5± 1.3 VS. 6.2± 1.9 and 4.5± 1.3 VS. 7.1±1.1 mU/mL; P<0.001 respectively). Pearson correlation test shows a significant liner correlation between mtDNA copy number and serum glutathione level (R = 0.2, P = 0.009).

Conclusions/Significance

COPD is associated with decreased leukocyte mtDNA copy number and serum glutathione. COPD is a regulatory disorder of leukocytes mitochondria. However, further studies are needed to determine the real mechanisms about the gene and the function of mitochondria.     

Transgene Expression Is Associated with Copy Number and Cytomegalovirus Promoter Methylation in Transgenic Pigs

Transgenic animals have been used for years to study gene function, produce important proteins, and generate models for the study of human diseases. However, inheritance and expression instability of the transgene in transgenic animals is a major limitation. Copy number and promoter methylation are known to regulate gene expression, but no report has systematically examined their effect on transgene expression. In the study, we generated two transgenic pigs by somatic cell nuclear transfer (SCNT) that express green fluorescent protein (GFP) driven by cytomegalovirus (CMV). Absolute quantitative real-time PCR and bisulfite sequencing were performed to determine transgene copy number and promoter methylation level. The correlation of transgene expression with copy number and promoter methylation was analyzed in individual development, fibroblast cells, various tissues, and offspring of the transgenic pigs. Our results demonstrate that transgene expression is associated with copy number and CMV promoter methylation in transgenic pigs.     

Somatic Copy Number Alterations Associated with Japanese or Endometriosis in Ovarian Clear Cell Adenocarcinoma

When compared with other epithelial ovarian cancers, the clinical characteristics of ovarian clear cell adenocarcinoma (CCC) include 1) a higher incidence among Japanese, 2) an association with endometriosis, 3) poor prognosis in advanced stages, and 4) a higher incidence of thrombosis as a complication. We used high resolution comparative genomic hybridization (CGH) to identify somatic copy number alterations (SCNAs) associated with each of these clinical characteristics of CCC. The Human Genome CGH 244A Oligo Microarray was used to examine 144 samples obtained from 120 Japanese, 15 Korean, and nine German patients with CCC. The entire 8q chromosome (minimum corrected p-value: q = 0.0001) and chromosome 20q13.2 including the ZNF217 locus (q = 0.0078) were amplified significantly more in Japanese than in Korean or German samples. This copy number amplification of the ZNF217 gene was confirmed by quantitative real-time polymerase chain reaction (Q-PCR). ZNF217 RNA levels were also higher in Japanese tumor samples than in non-Japanese samples (P = 0.027). Moreover, endometriosis was associated with amplification of EGFR gene (q = 0.047), which was again confirmed by Q-PCR and correlated with EGFR RNA expression. However, no SCNAs were significantly associated with prognosis or thrombosis. These results indicated that there may be an association between CCC and ZNF217 amplification among Japanese patients as well as between endometriosis and EGFR gene amplifications.     

Copy Number Variations Associated With Obesity‐Related Traits in African Americans: A Joint Analysis Between GENOA and HyperGEN

Obesity is a highly heritable trait and a growing public health problem. African Americans (AAs) are a genetically diverse, yet understudied population with a high prevalence of obesity (BMI >30 kg/m²). Recent studies based upon single‐nucleotide polymorphisms (SNPs) have identified genetic markers associated with obesity. However, a large proportion of the heritability of obesity remains unexplained. Copy number variation (CNV) has been cited as a possible source of missing heritability in common diseases such as obesity. We conducted a CNV genome‐wide association study of BMI in two African‐American cohorts from Genetic Epidemiology Network of Arteriopathy (GENOA) and Hypertension Genetic Epidemiology Network (HyperGEN). We performed independent and identical association analyses in each study, then combined the results in a meta‐analysis. We identified three CNVs associated with BMI, obesity, and other obesity‐related traits after adjusting for multiple testing. These CNVs overlap the PARK2, GYPA, and SGCZ genes. Our results suggest that CNV may play a role in the etiology of obesity in AAs.     

Apolipoprotein(a) Kringle-IV Type 2 Copy Number Variation Is Associated with Venous Thromboembolism

In addition to the established association between high lipoprotein(a) [Lp(a)] concentrations and coronary artery disease, an association between Lp(a) and venous thromboembolism (VTE) has also been described. Lp(a) is controlled by genetic variants in LPA gene, coding for apolipoprotein(a), including the kringle-IV type 2 (KIV-2) size polymorphism. Aim of the study was to investigate the role of LPA gene KIV-2 size polymorphism and single nucleotide polymorphisms (SNPs) (rs1853021, rs1800769, rs3798220, rs10455872) in modulating VTE susceptibility. Five hundred and sixteen patients with VTE without hereditary and acquired thrombophilia and 1117 healthy control subjects, comparable for age and sex, were investigated. LPA KIV-2 polymorphism, rs3798220 and rs10455872 SNPs were genotyped by TaqMan technology. Concerning rs1853021 and rs1800769 SNPs, PCR-RFLP assay was used. LPA KIV-2 repeat number was significantly lower in patients than in controls [median (interquartile range) 11(6–17) vs 15(9–25), p<0.0001]. A significantly higher prevalence of KIV-2 repeat number ≤7 was observed in patients than in controls (33.5% vs 15.5%, p<0.0001). KIV-2 repeat number was independently associated with VTE (p = 4.36 x10^-9), as evidenced by the general linear model analysis adjusted for transient risk factors. No significant difference in allele frequency for all SNPs investigated was observed. Haplotype analysis showed that LPA haplotypes rather than individual SNPs influenced disease susceptibility. Receiver operating characteristic curves analysis showed that a combined risk prediction model, including KIV-2 size polymorphism and clinical variables, had a higher performance in identifying subjects at VTE risk than a clinical-only model, also separately in men and women.     

Peripheral Blood Mitochondrial DNA Copy Number Is Associated with Prostate Cancer Risk and Tumor Burden

Alterations of mitochondrial DNA (mtDNA) have been associated with the risk of a number of human cancers; however, the relationship between mtDNA copy number in peripheral blood leukocytes (PBLs) and the risk of prostate cancer (PCa) has not been investigated. In a case-control study of 196 PCa patients and 196 age-paired healthy controls in a Chinese Han population, the association between mtDNA copy number in PBLs and PCa risk was evaluated. The relative mtDNA copy number was measured using quantitative real-time PCR; samples from three cases and two controls could not be assayed, leaving 193 cases and 194 controls for analysis. PCa patients had significantly higher mtDNA copy numbers than controls (medians 0.91 and 0.82, respectively; P<0.001). Dichotomized at the median value of mtDNA copy number in the controls, high mtDNA copy number was significantly associated with an increased risk of PCa (adjusted odds ratio  = 1.85, 95% confidence interval: 1.21–2.83). A significant dose-response relationship was observed between mtDNA copy number and risk of PCa in quartile analysis (P_trend = 0.011). Clinicopathological analysis showed that high mtDNA copy numbers in PCa patients were significantly associated with high Gleason score and advanced tumor stage, but not serum prostate-specific antigen level (P = 0.002, 0.012 and 0.544, respectively). These findings of the present study indicate that increased mtDNA copy number in PBLs is significantly associated with an increased risk of PCa and may be a reflection of tumor burden.     

Leukocyte Mitochondrial DNA Copy Number in Blood Is Not Associated with Major Depressive Disorder in Young Adults

Background

Major depressive disorder (MDD) is the leading cause of disability worldwide, and has significant genetic predisposition. Mitochondria may have a role in MDD and so mitochondrial DNA (mtDNA) has been suggested as a possible biomarker for this disease. We aimed to test whether the mtDNA copy number of peripheral blood leukocytes is related to MDD in young adults.

Methods

A case-control study was conducted with 210 MDD patients and 217 healthy controls (HC). The mtDNA copy number was measured by quantitative polymerase chain reaction (qPCR) method. Depression severity was assessed by the Hamilton-17 Depression Rating Scale (HDRS-17).

Results

We found no significant differences in mtDNA copy number between MDD patients and HC, though the power analysis showed that our sample size has enough power to detect the difference. There were also no significant correlations between mtDNA copy number and the clinical characteristics (such as age, age of onset, episodes, Hamilton Depression Rating Scale (HDRS) score and Global Assessment of Function Scale (GAF) score) in MDD patients.

Conclusion

Our study suggests that leukocyte mtDNA copy number is unlikely to contribute to MDD, but it doesn’t mean that we can exclude the possibility of involvement of mitochondria in the disease. Further studies are required to elucidate whether mtDNA can be a biomarker of MDD.