Human genes involved in copy number variation: mechanisms of origin, functional effects and implications for disease

Analysis of genomic copy number variations through whole-genome scan in Yunling cattle

Yunling cattle is a new breed of beef cattle bred in Yunnan Province, China, which has the advantages of fast growth, excellent meat quality, improved tolerance ability, and important landscape value. Copy number variation (CNV) is a significant source of gene structural variation and plays a crucial role in evolution and phenotypic diversity. Based on the latest reference genome ARS-UCD2.0, this study analyzed the genome-wide distribution of CNVs in Yunling cattle using short-read whole-genome sequencing data (n = 129) and single-molecule long-read sequencing data (n = 1), and a total of 16,507 CNVs were detected. After merging CNVs with overlapping genomic positions, 3,728 CNV regions (CNVRs) were obtained, accounting for 0.61% of the reference genome. The functional analysis indicated significant enrichment of CNVRs in 96 GO terms and 57 KEGG pathways, primarily related to cell adhesion, signal transduction, neuromodulation, and nutritional metabolism. Additionally, 111 CNVRs overlapped with 76 quantitative trait loci (QTLs), including Subcutaneous fat thickness QTL, Longissimus muscle area QTL, and Marbling score QTL. Several CNVR-overlapping genes, including BZW1, AOX1, and LOC100138449, overlap with regions associated with meat color and quality QTLs. Furthermore, Vst analysis showed that PSMB4, ERICH1, SMC2, and PPP4R3A were highly divergent between Yunling and Brahman cattle. In summary, we have constructed the genomic CNV map of Yunling cattle for the first time using whole-genome resequencing. This provides valuable genetic variation resources for the study of the Yunling cattle genome and contributes to the study of economic traits in Yunling cattle.

Prostaglandin synthase activity of sigma- and mu-class glutathione transferases in a parasitic trematode, Clonorchis sinensis

Sigma-class glutathione transferase (GST) proteins with dual GST and prostaglandin synthase (PGS) activities play a crucial role in the establishment of Clonorchis sinensis infection. Herein, we analyzed the structural and enzymatic properties of sigma-class GST (CsGST-σ) proteins to obtain insight into their antioxidant and immunomodulatory functions in comparison with mu-class GST (CsGST-μ) proteins. CsGST-σ proteins conserved characteristic structures, which had been described in mammalian hematopoietic prostaglandin D2 synthases. Recombinant forms of these CsGST-σ and CsGST-μ proteins expressed in Escherichia coli exhibited considerable degrees of GST and PGS activities with substantially different specific activities. All recombinant proteins displayed higher affinities toward prostaglandin H2 (PGS substrate; average Km of 30.7 and 3.0 μm for prostaglandin D2 [PGDS] and E2 synthase [PGES], respectively) than those toward CDNB (GST substrate; average Km of 1,205.1 μm). Furthermore, the catalytic efficiency (Kcat/Km) of the PGDS/PGES activity was higher than that of GST activity (average Kcat/Km of 3.1, 0.7, and 7.0×10-3 s-1μm-1 for PGDS, PGES, and GST, respectively). Our data strongly suggest that the C. sinensis sigma- and mu-class GST proteins are deeply involved in regulating host immune responses by generating PGD2 and PGE2 in addition to their roles in general detoxification.

The Expansion of House Mouse Major Urinary Protein Genes Likely Did Not Facilitate Commensalism with Humans

Mouse wild-derived strains (WDSs) combine the advantages of classical laboratory stocks and wild animals, and thus appear to be promising tools for diverse biomedical and evolutionary studies. We employed 18 WDSs representing three non-synanthropic species (Mus spretus, Mus spicilegus, and M. macedonicus) and three house mouse subspecies (Mus musculus musculus, M. m. domesticus, M. m. castaneus), which are all important human commensals to explore whether the number of major urinary protein (MUP) genes and their final protein levels in urine are correlated with the level of commensalism. Contrary to expectations, the MUP copy number (CN) and protein excretion in the strains derived from M. m. castaneus, which is supposed to be the strongest commensal, were not significantly different from the non-commensal species. Regardless of an overall tendency for higher MUP amounts in taxa with a higher CN, there was no significant correlation at the strain level. Our study thus suggests that expansion of the Mup cluster, which appeared before the house mouse diversification, is unlikely to facilitate commensalism with humans in three house mouse subspecies. Finally, we found considerable variation among con(sub)specific WDSs, warning against generalisations of results based on a few strains.

Genome-wide detection of CNV regions between Anqing six-end-white and Duroc pigs

BACKGROUND

Anqing six-end-white pig is a native breed in Anhui Province. The pigs have the disadvantages of a slow growth rate, low proportion of lean meat, and thick back fat, but feature the advantages of strong stress resistance and excellent meat quality. Duroc pig is an introduced pig breed with a fast growth rate and high proportion of lean meat. With the latter breed featuring superior growth characteristics but inferior meat quality traits, the underlying molecular mechanism that causes these phenotypic differences between Chinese and foreign pigs is still unclear.

RESULTS

In this study, copy number variation (CNV) detection was performed using the re-sequencing data of Anqing Six-end-white pigs and Duroc pigs, A total of 65,701 CNVs were obtained. After merging the CNVs with overlapping genomic positions, 881 CNV regions (CNVRs) were obtained. Based on the obtained CNVR information combined with their positions on the 18 chromosomes, a whole-genome map of the pig CNVs was drawn. GO analysis of the genes in the CNVRs showed that they were primarily involved in the cellular processes of proliferation, differentiation, and adhesion, and primarily involved in the biological processes of fat metabolism, reproductive traits, and immune processes.

CONCLUSION

The difference analysis of the CNVs between the Chinese and foreign pig breeds showed that the CNV of the Anqing six-end-white pig genome was higher than that of the introduced pig breed Duroc. Six genes related to fat metabolism, reproductive performance, and stress resistance were found in genome-wide CNVRs (DPF3, LEPR, MAP2K6, PPARA, TRAF6, NLRP4).

Resolving single-cell copy number profiling for large datasets

The advances of single-cell DNA sequencing (scDNA-seq) enable us to characterize the genetic heterogeneity of cancer cells. However, the high noise and low coverage of scDNA-seq impede the estimation of copy number variations (CNVs). In addition, existing tools suffer from intensive execution time and often fail on large datasets. Here, we propose SeCNV, an efficient method that leverages structural entropy, to profile the copy numbers. SeCNV adopts a local Gaussian kernel to construct a matrix, depth congruent map (DCM), capturing the similarities between any two bins along the genome. Then, SeCNV partitions the genome into segments by minimizing the structural entropy from the DCM. With the partition, SeCNV estimates the copy numbers within each segment for cells. We simulate nine datasets with various breakpoint distributions and amplitudes of noise to benchmark SeCNV. SeCNV achieves a robust performance, i.e. the F1-scores are higher than 0.95 for breakpoint detections, significantly outperforming state-of-the-art methods. SeCNV successfully processes large datasets (>50 000 cells) within 4 min, while other tools fail to finish within the time limit, i.e. 120 h. We apply SeCNV to single-nucleus sequencing datasets from two breast cancer patients and acoustic cell tagmentation sequencing datasets from eight breast cancer patients. SeCNV successfully reproduces the distinct subclones and infers tumor heterogeneity. SeCNV is available at https://github.com/deepomicslab/SeCNV.

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

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

The genetic informational network: how DNA conveys semantic information

The question of whether "genetic information" is a merely causal factor in development or can be made sense of semantically, in a way analogous to a language or other type of representation, has generated a long debate in the philosophy of biology. It is intimately connected with another intense debate, concerning the limits of genetic determinism. In this paper I argue that widespread attempts to draw analogies between genetic information and information contained in books, blueprints or computer programs, are fundamentally inadequate. In development, gene exons are the central part of an intricate and densely ramified semantic Genetic Informational Network. DNA in the entire genome is in a state of continuous positive and negative feedback with itself and with its 'environment', and is 'read' and acted upon by the cell in various alternative and complementary ways. The linear combinatorial coding relation between codons and amino acids is but one aspect of semantic genetic information, which is, when considered in its entirety, a far wider and richer concept.

Quantitative neurogenetics: applications in understanding disease

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

Identification of Copy Number Variation Among Nonsyndromic Cleft Lip and or Without Cleft Palate With Hypodontia: A Genome-Wide Association Study

Nonsyndromic cleft lip and or without cleft palate (NSCL/P) with the hypodontia is a common developmental abnormality in humans and animals. This study identified the genetic aberration involved in both NSCL/P and hypodontia pathogenesis. A cross-sectional study using genome-wide study copy number variation-targeted CytoScan 750K array carried out on salivary samples from 61 NSCL/P and 20 noncleft with and without hypodontia Malay subjects aged 7-13 years old. Copy number variations (CNVs) of SKI and fragile histidine triad (FHIT) were identified in NSCL/P and noncleft children using quantitative polymerase chain reaction (qPCR) as a validation analysis. Copy number calculated (CNC) for each gene determined with Applied Biosystems CopyCaller Software v2.0. The six significant CNVs included gains (12q14.3, 15q26.3, 1p36.32, and 1p36.33) and losses (3p14.2 and 4q13.2) in NSCL/P with hypodontia patients compared with the NSCL/P only. The genes located in these regions encoded LEMD3, IGF1R, TP73, SKI, FHIT, and UGT2β15. There were a significant gain and loss of both SKI and FHIT copy number in NSCL/P with hypodontia compared with the noncleft group (p < 0.05). The results supported that CNVs significantly furnish to the development of NSCL/P with hypodontia.

Whole-exome sequencing of de novo genetic variants in a Chinese family with a sporadic case of congenital nonsyndromic hearing loss

Background: We examined the genetic variants of a Chinese family with a 22-month-old infant with sporadic non-syndromic sensorineural hearing loss (NSHL). Methods: The whole-exome sequence data in the family, especially the de novo variants presented in the patient, were analyzed and the effect of the disease-causing genetic variants on the protein expression level and cellular localization were examined by cell-based functional assay. Results: The infant had no known NSHL-causing variants, except two compound heterozygous variants in connexin26 gene GJB2; one was the c.79G>A, c.341A>G haplotype from the asymptomatic mother who was benign, and the other was a de novo pathogenic c.262G>C (p.A88P). In vitro, GJB2 with c.262G>C was weakly expressed and displayed a punctate distribution in the cytoplasm and cytomembrane, while wild type GJB2 was robustly expressed in the cytomembrane. We deduced that the de novo pathogenic GJB2 c.262G>C exacerbated loss-of-function in the context of leaky variants c.79G>A, c.341A>G in the patient. Interestingly, further analysis of exome sequences revealed that the occurrence of de novo pathogenic variants in the infant was frequent. Among the total~47,000 variants, 143 were de novo in the patient, whereas among all 74 variants predicted to be pathogenic/likely pathogenic, 21 were heterozygous and two were homozygous de novo. The occurrence rate of de novo deleterious variants was much higher (31.1%, 23/74) than that in total (0.34%, 143/47,000). It is notable that most genes with de novo deleterious variants were environment-sensitive, such as GJB2, MNK1, MNK2, MUC4, RAD21 and DNA copy number variations. Conclusions: The full picture of genetic variants in the exome might help us to interpret the NSHL-causing variants. More research is needed into the causes of de novo deleterious variants and gene-environment interactions in congenital NSHL.

Whole-exome sequencing of de novo genetic variants in a Chinese family with a sporadic case of congenital nonsyndromic hearing loss

Background: We examined the genetic variants of a Chinese family with a 22-month-old infant with sporadic non-syndromic sensorineural hearing loss (NSHL). Methods: The whole-exome sequence data in the family, especially the de novo variants presented in the patient, were analyzed and the effect of the disease-causing genetic variants on the protein expression level and cellular localization were examined by cell-based functional assay. Results: The infant had no known NSHL-causing variants, except two compound heterozygous variants in connexin26 gene GJB2; one was the c.79G>A, c.341A>G haplotype from the asymptomatic mother which was benign, and the other was a de novo pathogenic c.262G>C (p.A88P). In vitro, GJB2 with c.262G>C was weakly expressed and displayed a punctate distribution in the cytoplasm and cytomembrane, while wild type GJB2 was robustly expressed in the cytomembrane. We deduced that the de novo pathogenic GJB2 c.262G>C exacerbated loss-of-function in the context of leaky variants c.79G>A, c.341A>G in the patient. Interestingly, further analysis of exome sequences revealed that the occurrence of de novo pathogenic variants in the infant was frequent. Among the total~47,000 variants, 143 were de novo in the patient, whereas among all 74 variants predicted to be pathogenic/likely pathogenic, 21 were heterozygous and two were homozygous de novo. The occurrence rate of de novo deleterious variants was much higher (31.1%, 23/74) than that in total (0.34%, 143/47,000). It is notable that most genes with de novo deleterious variants were environment-sensitive, such as GJB2, MNK1, MNK2, MUC4, RAD21 and DNA copy number variations. Conclusions: The full picture of genetic variants in the exome might help us to interpret the NSHL-causing variants. More research is needed into the causes of de novo deleterious variants and gene-environment interactions in congenital NSHL.

Analysis of copy number variation in dogs implicates genomic structural variation in the development of anterior cruciate ligament rupture

Anterior cruciate ligament (ACL) rupture is an important condition of the human knee. Second ruptures are common and societal costs are substantial. Canine cranial cruciate ligament (CCL) rupture closely models the human disease. CCL rupture is common in the Labrador Retriever (5.79% prevalence), ~100-fold more prevalent than in humans. Labrador Retriever CCL rupture is a polygenic complex disease, based on genome-wide association study (GWAS) of single nucleotide polymorphism (SNP) markers. Dissection of genetic variation in complex traits can be enhanced by studying structural variation, including copy number variants (CNVs). Dogs are an ideal model for CNV research because of reduced genetic variability within breeds and extensive phenotypic diversity across breeds. We studied the genetic etiology of CCL rupture by association analysis of CNV regions (CNVRs) using 110 case and 164 control Labrador Retrievers. CNVs were called from SNPs using three different programs (PennCNV, CNVPartition, and QuantiSNP). After quality control, CNV calls were combined to create CNVRs using ParseCNV and an association analysis was performed. We found no strong effect CNVRs but found 46 small effect (max(T) permutation P<0.05) CCL rupture associated CNVRs in 22 autosomes; 25 were deletions and 21 were duplications. Of the 46 CCL rupture associated CNVRs, we identified 39 unique regions. Thirty four were identified by a single calling algorithm, 3 were identified by two calling algorithms, and 2 were identified by all three algorithms. For 42 of the associated CNVRs, frequency in the population was <10% while 4 occurred at a frequency in the population ranging from 10–25%. Average CNVR length was 198,872bp and CNVRs covered 0.11 to 0.15% of the genome. All CNVRs were associated with case status. CNVRs did not overlap previous canine CCL rupture risk loci identified by GWAS. Associated CNVRs contained 152 annotated genes; 12 CNVRs did not have genes mapped to CanFam3.1. Using pathway analysis, a cluster of 19 homeobox domain transcript regulator genes was associated with CCL rupture (P = 6.6E-13). This gene cluster influences cranial-caudal body pattern formation during embryonic limb development. Clustered genes were found in 3 CNVRs on chromosome 14 (HoxA), 28 (NKX6-2), and 36 (HoxD). When analysis was limited to deletion CNVRs, the association was strengthened (P = 8.7E-16). This study suggests a component of the polygenic risk of CCL rupture in Labrador Retrievers is associated with small effect CNVs and may include aspects of stifle morphology regulated by homeobox domain transcript regulator genes.

Molecular characterization of the calcium release channel deficiency syndrome

We identified a potentially novel homozygous duplication involving the promoter region and exons 1–4 of the gene encoding type 2 cardiac ryanodine receptor (RYR2) that is responsible for highly penetrant, exertion-related sudden deaths/cardiac arrests in the Amish community without an overt phenotype to suggest RYR2-mediated catecholaminergic polymorphic ventricular tachycardia (CPVT). Homozygous RYR2 duplication (RYR2-DUP) induced pluripotent stem cell cardiomyocytes (iPSC-CMs) were generated from 2 unrelated patients. There was no difference in baseline Ca²⁺ handling measurements between WT-iPSC-CM and RYR2-DUP-iPSC-CM lines. However, compared with WT-iPSC-CMs, both patient lines demonstrated a dramatic reduction in caffeine-stimulated and isoproterenol-stimulated (ISO-stimulated) Ca²⁺ transient amplitude, suggesting RyR2 loss of function. There was a greater than 50% reduction in RYR2 transcript/RyR2 protein expression in both patient iPSC-CMs compared with WT. Delayed afterdepolarization was observed in the RYR2-DUP-iPSC-CMs but not in the WT-iPSC-CMs. Compared with WT-iPSC-CMs, there was significantly elevated arrhythmic activity in the RYR2-DUP-iPSC-CMs in response to ISO. Nadolol, propranolol, and flecainide reduced erratic activity by 8.5-fold, 6.8-fold, and 2.4-fold, respectively, from ISO challenge. Unlike the gain-of-function mechanism observed in RYR2-mediated CPVT, the homozygous multiexon duplication precipitated a dramatic reduction in RYR2 transcription and RyR2 protein translation, a loss of function in calcium handling, and a calcium-induced calcium release apparatus that is insensitive to catecholamines and caffeine.

Molecular and functional characterization of the calcium release channel deficiency syndrome in patient-specific induced pluripotent stem cell-cardiomyocytes.

Genome-Wide Detection of CNVs and Association With Body Weight in Sheep Based on 600K SNP Arrays

Copy number variations (CNVs) are important genomic structural variations and can give rise to significant phenotypic diversity. Herein, we used high-density 600K SNP arrays to detect CNVs in two synthetic lines of sheep (DS and SHH) and in Hu sheep (a local Chinese breed). A total of 919 CNV regions (CNVRs) were detected with a total length of 48.17 Mb, accounting for 1.96% of the sheep genome. These CNVRs consisted of 730 gains, 102 losses, and 87 complex CNVRs. These CNVRs were significantly enriched in the segmental duplication (SD) region. A CNVR-based cluster analysis of the three breeds revealed that the DS and SHH breeds share a close genetic relationship. Functional analysis revealed that some genes in these CNVRs were also significantly enriched in the olfactory transduction pathway (oas04740), including members of the OR gene family such as OR6C76, OR4Q2, and OR4K14. Using association analyses and previous gene annotations, we determined that a subset of identified genes was likely to be associated with body weight, including FOXF2, MAPK12, MAP3K11, STRBP, and C14orf132. Together, these results offer valuable information that will guide future efforts to explore the genetic basis for body weight in sheep.

Copy number variation of CCL3L1 among three major ethnic groups in Malaysia

Background

C-C motif Chemokine Ligand 3 Like 1 (CCL3L1) is a multiallelic copy number variable, which plays a crucial role in immunoregulatory and hosts defense through the production of macrophage inflammatory protein (MIP)-1α. Variable range of the CCL3L1 copies from 0 to 14 copies have been documented in several different populations. However, there is still lack of report on the range of CCL3L1 copy number exclusively among Malaysians who are a multi-ethnic population. Thus, this study aims to extensively examine the distribution of CCL3L1 copy number in the three major populations from Malaysia namely Malay, Chinese and Indian. A diploid copy number of CCL3L1 for 393 Malaysians (Malay = 178, Indian = 90, and Chinese = 125) was quantified using Paralogue Ratio Tests (PRTs) and then validated with microsatellites analysis.

Results

To our knowledge, this is the first report on the CCL3L1 copy number that has been attempted among Malaysians and the Chinese ethnic group exhibits a diverse pattern of CCL3L1 distribution copy number from the Malay and Indian (p < 0.0001). The CCL3L1 ranged from 0 to 8 copies for both the Malay and Indian ethnic groups while 0 to 10 copies for the Chinese ethnic. Consequently, the CCL3L1 copy number among major ethnic groups in the Malaysian population is found to be significantly varied when compared to the European population (p < 0.0001). The mean/median reported for the Malay, Chinese, Indian, and European are 2.759/2.869, 3.453/3.290, 2.437/1.970 and 2.001/1.940 respectively.

Conclusion

This study reveals the existence of genetic variation of CCL3L1 in the Malaysian population, and suggests by examining genetic diversity on the ethnicity, and specific geographical region could help in reconstructing human evolutionary history and for the prediction of disease risk related to the CCL3L1 copy number.

Midnolin is a confirmed genetic risk factor for Parkinson's disease

Objective

Genetic analysis of patients with familial Parkinson’s disease (PD) identified many causative genes. However, the majority of PD cases are sporadic, and the mechanisms of onset still remain unclear. Previously, we found that Midnolin (MIDN) is associated with PD in a Yamagata (Japan) cohort study and that MIDN regulates neurite outgrowth and Parkin expression in neuronal cells. In the present study, we aimed to replicate the genetic association between MIDN and PD in a large British population cohort.

Methods

In this replication study, we analyzed the copy number variations and single‐nucleotide polymorphisms of the MIDN gene in a large British population on a case–control genome‐wide association study dataset including 2,860 controls and 2,168 PD patients.

Results

There was significant copy number loss in the MIDN gene with an odds ratio of 4.35 (P < 2.2 × 10⁻¹⁶). Furthermore, there were many patients in both the British and Yamagata case groups who have a long spanning deletion. The odds ratio dramatically increased to 22.3 (P = 3.59 × 10⁻¹⁵) when a deletion spanning more than 50,000 bp was defined as the copy number loss. There were no significant differences between the controls and study cases for two relatively frequent single‐nucleotide polymorphisms (rs3746106 and rs3746107).

Interpretation

We showed the strong genetic association of MIDN with PD development in a British population and in a Japanese population, suggesting MIDN is a confirmed and universal genetic risk factor for PD.

Identification of co-evolving temporal networks

BACKGROUND

Biological networks describes the mechanisms which govern cellular functions. Temporal networks show how these networks evolve over time. Studying the temporal progression of network topologies is of utmost importance since it uncovers how a network evolves and how it resists to external stimuli and internal variations. Two temporal networks have co-evolving subnetworks if the evolving topologies of these subnetworks remain similar to each other as the network topology evolves over a period of time. In this paper, we consider the problem of identifying co-evolving subnetworks given a pair of temporal networks, which aim to capture the evolution of molecules and their interactions over time. Although this problem shares some characteristics of the well-known network alignment problems, it differs from existing network alignment formulations as it seeks a mapping of the two network topologies that is invariant to temporal evolution of the given networks. This is a computationally challenging problem as it requires capturing not only similar topologies between two networks but also their similar evolution patterns.

RESULTS

We present an efficient algorithm, Tempo, for solving identifying co-evolving subnetworks with two given temporal networks. We formally prove the correctness of our method. We experimentally demonstrate that Tempo scales efficiently with the size of network as well as the number of time points, and generates statistically significant alignments-even when evolution rates of given networks are high. Our results on a human aging dataset demonstrate that Tempo identifies novel genes contributing to the progression of Alzheimer's, Huntington's and Type II diabetes, while existing methods fail to do so.

CONCLUSIONS

Studying temporal networks in general and human aging specifically using Tempo enables us to identify age related genes from non age related genes successfully. More importantly, Tempo takes the network alignment problem one huge step forward by moving beyond the classical static network models.

Gene editing of the multi-copy H2A.B gene and its importance for fertility

BACKGROUND

Altering the biochemical makeup of chromatin by the incorporation of histone variants during development represents a key mechanism in regulating gene expression. The histone variant H2A.B, H2A.B.3 in mice, appeared late in evolution and is most highly expressed in the testis. In the mouse, it is encoded by three different genes. H2A.B expression is spatially and temporally regulated during spermatogenesis being most highly expressed in the haploid round spermatid stage. Active genes gain H2A.B where it directly interacts with polymerase II and RNA processing factors within splicing speckles. However, the importance of H2A.B for gene expression and fertility are unknown.

RESULTS

Here, we report the first mouse knockout of this histone variant and its effects on fertility, nuclear organization, and gene expression. In view of the controversy related to the generation of off-target mutations by gene editing approaches, we test the specificity of TALENs by disrupting the H2A.B multi-copy gene family using only one pair of TALENs. We show that TALENs do display a high level of specificity since no off-target mutations are detected by bioinformatics analyses of exome sequences obtained from three consecutive generations of knockout mice and by Sanger DNA sequencing. Male H2A.B.3 knockout mice are subfertile and display an increase in the proportion of abnormal sperm and clogged seminiferous tubules. Significantly, a loss of proper RNA Pol II targeting to distinct transcription-splicing territories and changes to pre-mRNA splicing are observed.

CONCLUSION

We have produced the first H2A.B knockout mouse using the TALEN approach.

Testicular Cancer: Genes, Environment, Hormones

Testicular cancer (TC) represents one of the most peculiar clinical challenges at present. In fact, currently treatments are so effective ensuring a 5 years disease-free survival rate in nearly 95% of patients. On the other hand however, TC represents the most frequent newly diagnosed form of cancer in men between the ages of 14 and 44 years, with an incidence ranging from <1 to 9.9 affected individuals per 100,000 males across countries, while the overall incidence is also increasing worldwide. Furthermore, cancer survivors show a 2% risk of developing cancer in the contralateral testis within 15 years of initial diagnosis. This complex and multifaceted scenario requires a great deal of effort to understand the clinical base of available evidence. It is now clear that genetic, environmental and hormonal risk factors concur and mutually influence both the development of the disease and its prognosis, in terms of response to treatment and the risk of recurrence. In this paper, the most recent issues describing the relative contribution of the aforementioned risk factors in TC development are discussed. In addition, particular attention is paid to the exposure to environmental chemical substances and thermal stress, whose role in cancer development and progression has recently been investigated at the molecular level.

Perspective: DNA Copy Number Variations in Cardiovascular Diseases

Human genome contains many variations, often called mutations, which are difficult to detect and have remained a challenge for years. A substantial part of the genome encompasses repeats and when such repeats are in the coding region they may lead to change in the gene expression profile followed by pathological conditions. Structural variants are alterations which change one or more sequence feature in the chromosome such as change in the copy number, rearrangements, and translocations of a sequence and can be balanced or unbalanced. Copy number variants (CNVs) may increase or decrease the copies of a given region and have a pivotal role in the onset of many diseases including cardiovascular disorders. Cardiovascular disorders have a magnitude of well-established risk factors and etiology, but their correlation with CNVs is still being studied. In this article, we have discussed history of CNVs and a summary on the diseases associated with CNVs. To detect such variations, we shed light on the number of techniques introduced so far and their limitations. The lack of studies on cardiovascular diseases to determine the frequency of such variants needs clinical studies with larger cohorts. This review is a compilation of articles suggesting the importance of CNVs in multitude of cardiovascular anomalies. Finally, future perspectives for better understanding of CNVs and cardiovascular disorders have also been discussed.

Gene copy number variations involved in balsam poplar (Populus balsamifera L.) adaptive variations

Gene copy number variations (CNVs) involved in phenotypic variations have already been shown in plants, but genomewide testing of CNVs for adaptive variation was not doable until recent technological developments. Thus, reports of the genomic architecture of adaptation involving CNVs remain scarce to date. Here, we investigated F₁ progenies of an intraprovenance cross (north-north cross, 58th parallel) and an interprovenances cross (north-south cross, 58th/49th parallels) for CNVs using comparative genomic hybridization on arrays of probes targeting gene sequences in balsam poplar (Populus balsamifera L.), a widespread North American forest tree. A total of 1,721 genes were found in varying copy numbers over the set of 19,823 tested genes. These gene CNVs presented an estimated average size of 8.3 kb and were distributed over poplar's 19 chromosomes including 22 hotspot regions. Gene CNVs number was higher for the interprovenance progeny in accordance with an expected higher genetic diversity related to the composite origin of this family. Regression analyses between gene CNVs and seven adaptive trait variations resulted in 23 significant links; among these adaptive gene CNVs, 30% were located in hotspots. One-to-five gene CNVs were found related to each of the measured adaptive traits and annotated for both biotic and abiotic stress responses. These annotations can be related to the occurrence of a higher pathogenic pressure in the southern parts of balsam poplar's distribution, and higher photosynthetic assimilation rates and water-use efficiency at high latitudes. Overall, our findings suggest that gene CNVs typically having higher mutation rates than SNPs may in fact represent efficient adaptive variations against fast-evolving pathogens.

Mapping of genome-wide copy number variations in the Iranian indigenous cattle using a dense SNP data set

The objective of this study was to present the first map of the copy number variations (CNVs) in Iranian indigenous cattle based on a high-density single nucleotide polymorphism (SNP) dataset. A total of 90 individuals were genotyped using the Illumina BovineHD BeadChip containing 777 962 SNPs. The QuantiSNP algorithm was used to perform a genome-wide CNV detection across autosomal genome. After merging the overlapping CNV, a total of 221 CNV regions were identified encompassing 36.4 Mb or 1.44% of the bovine autosomal genome. The length of the CNV regions ranged from 3.5 to 2252.8 Kb with an average of 163.8 Kb. These regions included 147 loss (66.52%) and 74 gain (33.48%) events containing a total of 637 annotated Ensembl genes. Gene ontology analysis revealed that most of genes in the CNV regions were involved in environmental responses, disease susceptibility and immune system functions. Furthermore, 543 of these genes corresponded to the human orthologous genes, which involved in a wide range of biological functions. Altogether, 73% of the 221 CNV regions overlapped either completely or partially with those previously reported in other cattle studies. Moreover, novel CNV regions involved several quantitative trait loci (QTL)-related to adaptative traits of Iranian indigenous cattle. These results provided a basis to conduct future studies on association between CNV regions and phenotypic variations in the Iranian indigenous cattle.

Copy number variations of E2F1: a new genetic risk factor for testicular cancer

Testicular germ cell tumor (TGCT) is one of the most heritable forms of cancer. In last years, many evidence suggested that constitutional genetic factors, mainly single nucleotide polymorphisms, can increase its risk. However, the possible contribution of copy number variations (CNVs) in TGCT susceptibility has not been substantially addressed. Indeed, an increasing number of studies have focused on the effect of CNVs on gene expression and on the role of these structural genetic variations as risk factors for different forms of cancer. E2F1 is a transcription factor that plays an important role in regulating cell growth, differentiation, apoptosis and response to DNA damage. Therefore, deficiency or overexpression of this protein might significantly influence fundamental biological processes involved in cancer development and progression, including TGCT. We analyzed E2F1 CNVs in 261 cases with TGCT and 165 controls. We found no CNVs in controls, but 17/261 (6.5%) cases showed duplications in E2F1 Blot analysis demonstrated higher E2F1 expression in testicular samples of TGCT cases with three copies of the gene. Furthermore, we observed higher phosphorylation of Akt and mTOR in samples with E2F1 duplication. Interestingly, normal, non-tumoral testicular tissue in patient with E2F1 duplication showed lower expression of E2F1 and lower AKT/mTOR phosphorylation with respect to adjacent tumor tissue. Furthermore, increased expression of E2F1 obtained in vitro in NTERA-2 testicular cell line induced increased AKT/mTOR phosphorylation. This study suggests for the first time an involvement of E2F1 CNVs in TGCT susceptibility and supports previous preliminary data on the importance of AKT/mTOR signaling pathway in this cancer.

Copy Number Variants in Alzheimer's Disease

Alzheimer's disease (AD) is a devastating disease mainly afflicting elderly people, characterized by decreased cognition, loss of memory, and eventually death. Although risk and deterministic genes are known, major genetics research programs are underway to gain further insights into the inheritance of AD. In the last years, in particular, new developments in genome-wide scanning methodologies have enabled the association of a number of previously uncharacterized copy number variants (CNVs, gain or loss of DNA) in AD. Because of the exceedingly large number of studies performed, it has become difficult for geneticists as well as clinicians to systematically follow, evaluate, and interpret the growing number of (sometime conflicting) CNVs implicated in AD. In this review, after a brief introduction of this type of structural variation, and a description of available databases, computational analyses, and technologies involved, we provide a systematic review of all published data showing statistical and scientific significance of pathogenic CNVs and discuss the role they might play in AD.

A multifaceted computational report on the variants effect on KIR2DL3 and IFNL3 candidate gene in HCV clearance

HCV infection causes acute and chronic liver diseases including, cirrhosis and hepatocellular carcinoma. Following HCV infection, spontaneous clearance occurs in approximately 20 % of the population dependant upon HCV genotype. In this study, functional and non-functional variant analysis was executed for the classical and the latest HCV clearance candidate genes namely, KIR2DL3 and IFNL3. Initially, the functional effects of non-synonymous SNPs were assigned on exposing to homology based tools, SIFT, PolyPhen-2 and PROVEAN. Further, UTR and splice sites variants were scanned for the gene expression and regulation changes. Subsequently, the haplotype and CNV were also identified. The mutation H77Y of KIR2DL3 and R157Q, H156Y, S63L, R157W, F179V, H128R, T101M, R180C, and F176I of IFNL3 results in conservation, RMSD, total energy, stability, and secondary structures revealed a negative impact on the structural fitness. UTRscan and the splice site result indicate functional change, which may affect gene regulation and expression. The graphical display of selected population shows alleles like rs270779, rs2296370, rs10423751, rs12982559, rs9797797, and rs35987710 of KIR2DL3 and rs12972991, rs12980275, rs4803217, rs8109886, and rs8099917 of IFNL3 are in high LD with a measure of [Formula: see text] broadcasting its protective effect in HCV clearance. Similarly, CNV report suggests major DNA fragment loss that could have a profound impact on the gene expression affecting the overall phenotype. This roundup report specifies the effect of NK cell receptor, KIR2DL3 and IFNL3 variants that can have a better prospect in GWAS and immunogenetic studies leading to better understanding of HCV clearance and progression.

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.

Genetic variants in LEP, LEPR, and MC4R explain 30% of severe obesity in children from a consanguineous population

OBJECTIVE

Single gene mutations leading to severe obesity have so far been identified in 3-5% cases in European populations. However, prevalence of these pathogenic mutations has not systematically been examined in specific consanguineous populations. Here we describe the incidence of obesity-associated mutations through a step-wise sequence analysis, in a cohort of 73 Pakistani children with severe obesity from consanguineous families.

METHODS

Initially, all subjects were screened for mutations in coding regions of leptin (LEP) and melanocortin 4 receptor (MC4R) genes by direct sequencing. Subjects negative for mutation in these genes were screened using microdroplet PCR enrichment and NGS. Genomic structural variation was assessed by genotyping. Serum leptin, insulin, and cortisol were determined by ELISA.

RESULTS

Among 73 children with severe obesity (BMI SDS > 3.0), we identified 22 probands and 5 relatives, carrying 10 different loss-of-function homozygous mutations in LEP, leptin receptor (LEPR), and MC4R genes, including 4 novel variants. Hypercortisolemia was significantly emphasized in LEP mutation carriers.

CONCLUSIONS

The prevalence of pathogenic mutations in genes known to directly influence leptin-melanocortin signaling is 30% in our cohort. The results of this study emphasize the desirability of undertaking systematic and in-depth genetic analysis of cases with severe obesity in specific consanguineous populations.

Improved Detection and Characterization of Copy Number Variations Among Diverse Pig Breeds by Array CGH

As a major component of genomic variation, copy number variations (CNVs) are considered as promising markers for some phenotypic and economically important traits in domestic animals. Using a custom-designed 1M array CGH (aCGH), we performed CNV discovery in 12 pig samples from one Asian wild boar population, six Chinese indigenous breeds, and two European commercial breeds. In total, we identified 758 CNV regions (CNVRs), covering 47.43 Mb of the pig genome sequence. Of the total porcine genes, 1295 genes were completely or partially overlapped with the identified CNVRs, which enriched in the terms related to sensory perception of the environment, neurodevelopmental processes, response to external stimuli, and immunity. Further probing the potential functions of these genes, we also found a suite of genes related important traits, which make them a promising resource for exploring the genetic basis of phenotype differences among diverse pig breeds. Compared with previous relevant studies, the current study highlights that different platforms can complement each other, and the combined implementation of different platforms is beneficial to achieve the most comprehensive CNV calls. CNVs detected in diverse populations herein are essentially complementary to the CNV map in the pig genome, which would be helpful for understanding the pig genome variants and investigating the associations between various phenotypes and CNVs.

Genome-wide analyses of LINE-LINE-mediated nonallelic homologous recombination

Nonallelic homologous recombination (NAHR), occurring between low-copy repeats (LCRs) >10 kb in size and sharing >97% DNA sequence identity, is responsible for the majority of recurrent genomic rearrangements in the human genome. Recent studies have shown that transposable elements (TEs) can also mediate recurrent deletions and translocations, indicating the features of substrates that mediate NAHR may be significantly less stringent than previously believed. Using >4 kb length and >95% sequence identity criteria, we analyzed of the genome-wide distribution of long interspersed element (LINE) retrotransposon and their potential to mediate NAHR. We identified 17 005 directly oriented LINE pairs located <10 Mbp from each other as potential NAHR substrates, placing 82.8% of the human genome at risk of LINE-LINE-mediated instability. Cross-referencing these regions with CNVs in the Baylor College of Medicine clinical chromosomal microarray database of 36 285 patients, we identified 516 CNVs potentially mediated by LINEs. Using long-range PCR of five different genomic regions in a total of 44 patients, we confirmed that the CNV breakpoints in each patient map within the LINE elements. To additionally assess the scale of LINE-LINE/NAHR phenomenon in the human genome, we tested DNA samples from six healthy individuals on a custom aCGH microarray targeting LINE elements predicted to mediate CNVs and identified 25 LINE-LINE rearrangements. Our data indicate that LINE-LINE-mediated NAHR is widespread and under-recognized, and is an important mechanism of structural rearrangement contributing to human genomic variability.

Analysis of genome-wide copy number variations in Chinese indigenous and western pig breeds by 60 K SNP genotyping arrays

Copy number variations (CNVs) represent a substantial source of structural variants in mammals and contribute to both normal phenotypic variability and disease susceptibility. Although low-resolution CNV maps are produced in many domestic animals, and several reports have been published about the CNVs of porcine genome, the differences between Chinese and western pigs still remain to be elucidated. In this study, we used Porcine SNP60 BeadChip and PennCNV algorithm to perform a genome-wide CNV detection in 302 individuals from six Chinese indigenous breeds (Tongcheng, Laiwu, Luchuan, Bama, Wuzhishan and Ningxiang pigs), three western breeds (Yorkshire, Landrace and Duroc) and one hybrid (Tongcheng×Duroc). A total of 348 CNV Regions (CNVRs) across genome were identified, covering 150.49 Mb of the pig genome or 6.14% of the autosomal genome sequence. In these CNVRs, 213 CNVRs were found to exist only in the six Chinese indigenous breeds, and 60 CNVRs only in the three western breeds. The characters of CNVs in four Chinese normal size breeds (Luchuan, Tongcheng and Laiwu pigs) and two minipig breeds (Bama and Wuzhishan pigs) were also analyzed in this study. Functional annotation suggested that these CNVRs possess a great variety of molecular function and may play important roles in phenotypic and production traits between Chinese and western breeds. Our results are important complementary to the CNV map in pig genome, which provide new information about the diversity of Chinese and western pig breeds, and facilitate further research on porcine genome CNVs.

Complexity of gene expression evolution after duplication: protein dosage rebalancing

Ongoing debates about functional importance of gene duplications have been recently intensified by a heated discussion of the “ortholog conjecture” (OC). Under the OC, which is central to functional annotation of genomes, orthologous genes are functionally more similar than paralogous genes at the same level of sequence divergence. However, a recent study challenged the OC by reporting a greater functional similarity, in terms of gene ontology (GO) annotations and expression profiles, among within-species paralogs compared to orthologs. These findings were taken to indicate that functional similarity of homologous genes is primarily determined by the cellular context of the genes, rather than evolutionary history. Subsequent studies suggested that the OC appears to be generally valid when applied to mammalian evolution but the complete picture of evolution of gene expression also has to incorporate lineage-specific aspects of paralogy. The observed complexity of gene expression evolution after duplication can be explained through selection for gene dosage effect combined with the duplication-degeneration-complementation model. This paper discusses expression divergence of recent duplications occurring before functional divergence of proteins encoded by duplicate genes.

Global copy number analyses by next generation sequencing provide insight into pig genome variation

Background

Copy number variations (CNVs) confer significant effects on genetic innovation and phenotypic variation. Previous CNV studies in swine seldom focused on in-depth characterization of global CNVs.

Results

Using whole-genome assembly comparison (WGAC) and whole-genome shotgun sequence detection (WSSD) approaches by next generation sequencing (NGS), we probed formation signatures of both segmental duplications (SDs) and individualized CNVs in an integrated fashion, building the finest resolution CNV and SD maps of pigs so far. We obtained copy number estimates of all protein-coding genes with copy number variation carried by individuals, and further confirmed two genes with high copy numbers in Meishan pigs through an enlarged population. We determined genome-wide CNV hotspots, which were significantly enriched in SD regions, suggesting evolution of CNV hotspots may be affected by ancestral SDs. Through systematically enrichment analyses based on simulations and bioinformatics analyses, we revealed CNV-related genes undergo a different selective constraint from those CNV-unrelated regions, and CNVs may be associated with or affect pig health and production performance under recent selection.

Conclusions

Our studies lay out one way for characterization of CNVs in the pig genome, provide insight into the pig genome variation and prompt CNV mechanisms studies when using pigs as biomedical models for human diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-593) contains supplementary material, which is available to authorized users.

Seasonal variation of urinary microRNA expression in male goats (Capra hircus) as assessed by next generation sequencing

Testosterone plays a key role in preparation of a male domesticated goat (Capra hircus) to breeding season including changes in the urogenital tract of a male goat (buck). microRNAs are important regulators of cellular metabolism, differentiation and function. They are powerful intermediaries of hormonal activity in the body, including the urogenital tract. We investigated seasonal changes in expression of microRNAs in goat buck urine and their potential consequences using next generation sequencing (microRNA-Seq). We determined the location of each microRNA gene in the goat genome. Testosterone was measured by radioimmunoassay and the androgen receptor binding sites (ARBS) in the promoters of the microRNA genes were determined by MatInspector. The overall impact of regulated microRNAs on cellular physiology was assessed by mirPath. We observed high testosterone levels during the breeding season and changes in the expression of forty microRNAs. Nineteen microRNAs were upregulated, while twenty-one were downregulated. We identified several ARBS in the promoters of regulated microRNAs. Notably, the mostly inhibited microRNA, miR-1246, has a unique set of several gene copy variants associated with a cluster of androgen receptor binding sites. Concomitant changes in regulated microRNA expression could promote transcription, proliferation and differentiation of urogenital tract cells. Together, these findings indicate that in a domesticated goat (Capra hircus), there are specific changes in the microRNA expression profile in buck urine during breeding season, which could be attributable to high testosterone levels during breeding, and could help in preparation of the urogenital tract for high metabolic demands of that season.

Enhancing genome-wide copy number variation identification by high density array CGH using diverse resources of pig breeds

Copy number variations (CNVs) are important forms of genomic variation, and have attracted extensive attentions in humans as well as domestic animals. In the study, using a custom-designed 2.1 M array comparative genomic hybridization (aCGH), genome-wide CNVs were identified among 12 individuals from diverse pig breeds, including one Asian wild population, six Chinese indigenous breeds and two modern commercial breeds (Yorkshire and Landrace), with one individual of the other modern commercial breed, Duroc, as the reference. A total of 1,344 CNV regions (CNVRs) were identified, covering 47.79 Mb (∼1.70%) of the pig genome. The length of these CNVRs ranged from 3.37 Kb to 1,319.0 Kb with a mean of 35.56 Kb and a median of 11.11 Kb. Compared with similar studies reported, most of the CNVRs (74.18%) were firstly identified in present study. In order to confirm these CNVRs, 21 CNVRs were randomly chosen to be validated by quantitative real time PCR (qPCR) and a high rate (85.71%) of confirmation was obtained. Functional annotation of CNVRs suggested that the identified CNVRs have important function, and may play an important role in phenotypic and production traits difference among various breeds. Our results are essential complementary to the CNV map in the pig genome, which will provide abundant genetic markers to investigate association studies between various phenotypes and CNVs in pigs.

Identification of genome-wide copy number variations among diverse pig breeds using SNP genotyping arrays

Copy number variations (CNVs) are important forms of genetic variation complementary to SNPs, and can be considered as promising markers for some phenotypic and economically important traits or diseases susceptibility in domestic animals. In the present study, we performed a genome-wide CNV identification in 14 individuals selected from diverse populations, including six types of Chinese indigenous breeds, one Asian wild boar population, as well as three modern commercial foreign breeds. We identified 63 CNVRs in total, which covered 9.98 Mb of polymorphic sequence and corresponded to 0.36% of the genome sequence. The length of these CNVRs ranged from 3.20 to 827.21 kb, with an average of 158.37 kb and a median of 97.85 kb. Functional annotation revealed these identified CNVR have important molecular function, and may play an important role in exploring the genetic basis of phenotypic variability and disease susceptibility among pigs. Additionally, to confirm these potential CNVRs, we performed qPCR for 12 randomly selected CNVRs and 8 of them (66.67%) were confirmed successfully. CNVs detected in diverse populations herein are essential complementary to the CNV map in the pig genome, which provide an important resource for studies of genomic variation and the association between various economically important traits and CNVs.

Analysis of copy number variants by three detection algorithms and their association with body size in horses

Background

Copy number variants (CNVs) have been shown to play an important role in genetic diversity of mammals and in the development of many complex phenotypic traits. The aim of this study was to perform a standard comparative evaluation of CNVs in horses using three different CNV detection programs and to identify genomic regions associated with body size in horses.

Results

Analysis was performed using the Illumina Equine SNP50 genotyping beadchip for 854 horses. CNVs were detected by three different algorithms, CNVPartition, PennCNV and QuantiSNP. Comparative analysis revealed 50 CNVs that affected 153 different genes mainly involved in sensory perception, signal transduction and cellular components. Genome-wide association analysis for body size showed highly significant deleted regions on ECA1, ECA8 and ECA9. Homologous regions to the detected CNVs on ECA1 and ECA9 have also been shown to be correlated with human height.

Conclusions

Comparative analysis of CNV detection algorithms was useful to increase the specificity of CNV detection but had certain limitations dependent on the detection tool. GWAS revealed genome-wide associated CNVs for body size in horses.

Gene duplication as a mechanism of genomic adaptation to a changing environment

A subject of extensive study in evolutionary theory has been the issue of how neutral, redundant copies can be maintained in the genome for long periods of time. Concurrently, examples of adaptive gene duplications to various environmental conditions in different species have been described. At this point, it is too early to tell whether or not a substantial fraction of gene copies have initially achieved fixation by positive selection for increased dosage. Nevertheless, enough examples have accumulated in the literature that such a possibility should be considered. Here, I review the recent examples of adaptive gene duplications and make an attempt to draw generalizations on what types of genes may be particularly prone to be selected for under certain environmental conditions. The identification of copy-number variation in ecological field studies of species adapting to stressful or novel environmental conditions may improve our understanding of gene duplications as a mechanism of adaptation and its relevance to the long-term persistence of gene duplications.

Significance Analysis of Microarrays (SAM) Offers Clues to Differences Between the Genomes of Adult Philadelphia Positive ALL and the Lymphoid Blast Transformation of CML

Philadelphia positive malignant disorders are a clinically divergent group of leukemias. These include chronic myeloid leukemia (CML) and de novo acute Philadelphia positive (Ph(+)) leukemia of both myeloid, and lymphoid origin. Recent whole genome screening of Ph(+)ALL in both children and adults identified an almost obligatory cryptic loss of Ikaros, required for the normal B cell maturation. Although similar losses were found in lymphoid blast crisis the genetic background of the transformation in CML is still poorly defined. We used Significance Analysis of Microarrays (SAM) to analyze comparative genomic hybridization (aCGH) data from 30 CML (10 each of chronic phase, myeloid and lymphoid blast stage), 10 Ph(+)ALL adult patients and 10 disease free controls and were able to: (a) discriminate between the genomes of lymphoid and myeloid blast cells and (b) identify differences in the genome profile of de novo Ph(+)ALL and lymphoid blast transformation of CML (BC/L). Furthermore we were able to distinguish a sub group of Ph(+) ALL characterized by gains in chromosome 9 and recurrent losses at several other genome sites offering genetic evidence for the clinical heterogeneity. The significance of these results is that they not only offer clues regarding the pathogenesis of Ph(+) disorders and highlight the potential clinical implications of a set of probes but also demonstrates what SAM can offer for the analysis of genome data.

Copy number variation in the genomes of domestic animals

Copy number variation (CNV) might be one of the main contributors to phenotypic diversity and evolutionary adaptation in animals and plants, employing a wide variety of mechanisms, such as gene dosage and transcript structure alterations, to modulate organismal plasticity. In the past 4 years, considerable advances have been made in the characterization of the genomic architecture of CNV in domestic species. First, low-resolution CNV maps were produced for cattle, goat, sheep, pig, dog, chicken, duck and turkey, showing that these structural polymorphisms comprise a significant part of these genomes. Furthermore, CNVs have been associated with several pigmentation (white coat in horse, pig and sheep) and morphological (late feathering and pea comb in chicken) traits, as well as with susceptibility to a wide array of diseases and developmental disorders, for example osteopetrosis, anhidrotic ectodermal dysplasia, copper toxicosis, intersexuality, cone degeneration, periodic fever and dermoid sinus, among others. In the future, development of high-resolution tools for CNV detection and typing combined with the implementation of databases integrating CNV, QTL and gene expression data will be essential to identify and measure the impact of this source of structural variation on the many phenotypes that are relevant to animal breeders and veterinary practitioners.

Controversies in modern evolutionary biology: the imperative for error detection and quality control

Background

The data from high throughput genomics technologies provide unique opportunities for studies of complex biological systems, but also pose many new challenges. The shift to the genome scale in evolutionary biology, for example, has led to many interesting, but often controversial studies. It has been suggested that part of the conflict may be due to errors in the initial sequences. Most gene sequences are predicted by bioinformatics programs and a number of quality issues have been raised, concerning DNA sequencing errors or badly predicted coding regions, particularly in eukaryotes.

Results

We investigated the impact of these errors on evolutionary studies and specifically on the identification of important genetic events. We focused on the detection of asymmetric evolution after duplication, which has been the subject of controversy recently. Using the human genome as a reference, we established a reliable set of 688 duplicated genes in 13 complete vertebrate genomes, where significantly different evolutionary rates are observed. We estimated the rates at which protein sequence errors occur and are accumulated in the higher-level analyses. We showed that the majority of the detected events (57%) are in fact artifacts due to the putative erroneous sequences and that these artifacts are sufficient to mask the true functional significance of the events.

Conclusions

Initial errors are accumulated throughout the evolutionary analysis, generating artificially high rates of event predictions and leading to substantial uncertainty in the conclusions. This study emphasizes the urgent need for error detection and quality control strategies in order to efficiently extract knowledge from the new genome data.

Targeted screening and validation of copy number variations

The accessibility of genome-wide screening technologies considerably facilitated the identification and characterization of copy number variations (CNVs). The increasing amount of available data describing these variants, clearly demonstrates their abundance in the human genome. This observation shows that not only SNPs, but also CNVs and other structural variants strongly contribute to genetic variation. Even though not all structural variants have an obvious phenotypic effect, there is evidence that CNVs influence gene dosage and hence can have profound effects on human disease susceptibility, disease manifestation, and disease severity. Therefore, CNV screening and analysis methodologies, specifically focusing on disease-related CNVs are actively progressing. This chapter specifically describes different techniques currently available for the targeted screening and validation of CNVs. We not only provide an overview of all these CNV analysis methods, but also address their strong and weak points. Methods covered include fluorescence in situ hybridization (FISH), quantitative real-time PCR (qPCR), paralogue ratio test (PRT), molecular copy-number counting (MCC), and multiplex PCR-based approaches, such as multiplex amplifiable probe hybridization (MAPH), multiplex ligation-dependent probe amplification (MLPA), multiplex PCR-based real-time invader assay (mPCR-RETINA), quantitative multiplex PCR of short fluorescent fragments (QMPSF), and multiplex amplicon quantification (MAQ). We end with some general remarks and conclusions, furthermore briefly addressing the future perspectives.

Accurate single-nucleotide polymorphism allele assignment in trisomic or duplicated regions by using a single base-extension assay with MALDI-TOF mass spectrometry

BACKGROUND

The accurate assignment of alleles embedded within trisomic or duplicated regions is an essential prerequisite for assessing the combined effects of single-nucleotide polymorphisms (SNPs) and genomic copy number. Such an integrated analysis is challenging because heterozygotes for such a SNP may be one of 2 genotypes-AAB or ABB. Established methods for SNP genotyping, however, can have difficulty discriminating between the 2 heterozygous trisomic genotypes. We developed a method for assigning heterozygous trisomic genotypes that uses the ratio of the height of the 2 allele peaks obtained by mass spectrometry after a single-base extension assay.

METHODS

Eighteen COL6A2 (collagen, type VI, alpha 2) SNPs were analyzed in euploid and trisomic individuals by means of a multiplexed single-base extension assay that generated allele-specific oligonucleotides of differing M(r) values for detection by MALDI-TOF mass spectrometry. Reference data (mean and SD) for the allele peak height ratios were determined from heterozygous euploid samples. The heterozygous trisomic genotypes were assigned by calculating the z score for each trisomic allele peak height ratio and by considering the sign (+/-) of the z score.

RESULTS

Heterozygous trisomic genotypes were assigned in 96.1% (range, 89.9%-100%) of the samples for each SNP analyzed. The genotypes obtained were reproduced in 95 (97.5%) of 97 loci retested in a second assay. Subsequently, the origin of nondisjunction was determined in 108 (82%) of 132 family trios with a Down syndrome child.

CONCLUSIONS

This approach enabled reliable genotyping of heterozygous trisomic samples and the determination of the origin of nondisjunction in Down syndrome family trios.

Functional annotation of genes overlapping copy number variants in autistic patients: focus on axon pathfinding

We have used Gene Ontology (GO) and pathway analyses to uncover the common functions associated to the genes overlapping Copy Number Variants (CNVs) in autistic patients. Our source of data were four published studies [1-4]. We first applied a two-step enrichment strategy for autism-specific genes. We fished out from the four mentioned studies a list of 2928 genes overall overlapping 328 CNVs in patients and we first selected a sub-group of 2044 genes after excluding those ones that are also involved in CNVs reported in the Database of Genomic Variants (enrichment step 1). We then selected from the step 1-enriched list a sub-group of 514 genes each of which was found to be deleted or duplicated in at least two patients (enrichment step 2). The number of statistically significant processes and pathways identified by the Database for Annotation, Visualization and Integrated Discovery and Ingenuity Pathways Analysis softwares with the step 2-enriched list was significantly higher compared to the step 1-enriched list. In addition, statistically significant GO terms, biofunctions and pathways related to nervous system development and function were exclusively identified by the step 2-enriched list of genes. Interestingly, 21 genes were associated to axon growth and pathfinding. The latter genes and other ones associated to nervous system in this study represent a new set of autism candidate genes deserving further investigation. In summary, our results suggest that the autism’s “connectivity genes” in some patients affect very early phases of neurodevelopment, i.e., earlier than synaptogenesis.

Copy number variations in candidate genes in neovascular age-related macular degeneration

PURPOSE

The pathogenesis of age-related macular degeneration (AMD) is strongly influenced by genetic factors, and single nucleotide polymorphisms have been consistently linked to AMD. Copy number variation (CNV), or variation in the number of copies of a particular segment of DNA, may also contribute to AMD pathogenesis. This study evaluated CNVs in candidate genes that have been reported to be linked to AMD.

METHODS

Study participants were 131 patients with neovascular AMD and 103 elderly persons without AMD who were evaluated by retinal specialists at the National Eye Institute. DNA was collected from peripheral whole blood, and duplex RT-PCR based copy number (CN) assays were performed for the genes CCR3, CFH, CX3CR1, ERCC6, HTRA1, and VEGF. Quantitative CNs (CN = 0, 1, 2, or 3+) were determined.

RESULTS

Novel CNVs were discovered in CCR3, CX3CR1, and ERCC6. The unadjusted data suggested that CN = 3+ for CX3CR1 might be mildly protective against AMD, but this trend did not persist after adjustment for age. AMD patients appeared to have an elevated mean CFH CN relative to controls (2.13 [95% confidence interval (CI), 2.05-2.21] vs. 2.01 [95% CI, 1.92-2.09 copies]; P = 0.05). No significant associations between CNV and AMD were observed for the remaining genes.

CONCLUSIONS

The methods described are suitable for quantitative characterization of CNV in candidate genes. The authors identified CNVs in AMD-associated genes but did not find strong evidence for a link with neovascular AMD.

Delineating the Hemostaseome as an aid to individualize the analysis of the hereditary basis of thrombotic and bleeding disorders

Next-generation sequencing and genome-wide association studies represent powerful tools to identify genetic variants that confer disease risk within populations. On their own, however, they cannot provide insight into how these variants contribute to individual risk for diseases that exhibit complex inheritance, or alternatively confer health in a given individual. Even in the case of well-characterized variants that confer a significant disease risk, more healthy individuals carry the variant, with no apparent ill effect, than those who manifest disease. Access to low-cost genome sequence data promises to provide an unprecedentedly detailed view of the nature of the hereditary component of complex diseases, but requires the large-scale comparison of sequence data from individuals with and without disease to deliver a clinical calibration. The provision of informatics support remains problematic as there are currently no means to interpret the data generated. Here, we initiate this process, a prerequisite for such a study, by narrowing the focus from an entire genome to that of a single biological system. To this end, we examine the 'Hemostaseome,' and more specifically focus on DNA sequence changes pertaining to those human genes known to impact upon hemostasis and thrombosis that can be analyzed coordinately, and on an individual basis, to interrogate how specific combinations of variants act to confer disease predisposition. As a first step, we delineate known members of the Hemostaseome and explore the nature of the genetic variants that may cause disease in individuals whose hemostatic balance has become shifted toward either a prothrombotic or anticoagulant phenotype.

Genes, mutations, and human inherited disease at the dawn of the age of personalized genomics

The number of reported germline mutations in human nuclear genes, either underlying or associated with inherited disease, has now exceeded 100,000 in more than 3,700 different genes. The availability of these data has both revolutionized the study of the morbid anatomy of the human genome and facilitated "personalized genomics." With approximately 300 new "inherited disease genes" (and approximately 10,000 new mutations) being identified annually, it is pertinent to ask how many "inherited disease genes" there are in the human genome, how many mutations reside within them, and where such lesions are likely to be located? To address these questions, it is necessary not only to reconsider how we define human genes but also to explore notions of gene "essentiality" and "dispensability."Answers to these questions are now emerging from recent novel insights into genome structure and function and through complete genome sequence information derived from multiple individual human genomes. However, a change in focus toward screening functional genomic elements as opposed to genes sensu stricto will be required if we are to capitalize fully on recent technical and conceptual advances and identify new types of disease-associated mutation within noncoding regions remote from the genes whose function they disrupt.

Molecular genetic analysis of Chd3 and polytene chromosome region 76B-D in Drosophila melanogaster

The Drosophila melanogaster Chd3 gene encodes a member of the CHD group of SNF2/RAD54 ATPases. CHD proteins are conserved from yeast to man and many are subunits of chromatin-remodeling complexes that facilitate transcription. Drosophila CHD3 proteins are not found in protein complexes, but as monomers that remodel chromatin in vitro. CHD3 colocalize with elongating RNA polymerase II on salivary gland polytene chromosomes. Since the role of Chd3 in development was unknown, we isolated and characterized the essential genes within the 640-kb region of the third chromosome (polytene chromosome region 76B-D) that includes Chd3. We recovered mutations in 24 genes that are essential for zygotic viability. We found that transposon-insertion mutants for 46% of the essential genes are included in the Drosophila Gene Disruption Project collection. None of the essential genes that we identified are in a 200-kb region that includes Chd3. We generated a deletion of Chd3 by targeted gene replacement. This deletion had no effect on either viability or fertility.