Genomic regions showing copy number variations associate with resistance or susceptibility to gastrointestinal nematodes in Angus cattle

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. We previously reported an initial analysis of copy number variations (CNVs) in Angus cattle selected for resistance or susceptibility to gastrointestinal nematodes. In this study, we performed a large-scale analysis of CNVs using SNP genotyping data from 472 animals of the same population. We detected 811 candidate CNV regions, which represent 141.8 Mb (~4.7%) of the genome. To investigate the functional impacts of CNVs, we created 2 groups of 100 individual animals with extremely low or high estimated breeding values of eggs per gram of feces and referred to these groups as parasite resistant (PR) or parasite susceptible (PS), respectively. We identified 297 (~51 Mb) and 282 (~48 Mb) CNV regions from PR and PS groups, respectively. Approximately 60% of the CNV regions were specific to the PS group or PR group of animals. Selected PR- or PS-specific CNVs were further experimentally validated by quantitative PCR. A total of 297 PR CNV regions overlapped with 437 Ensembl genes enriched in immunity and defense, like WC1 gene which uniquely expresses on gamma/delta T cells in cattle. Network analyses indicated that the PR-specific genes were predominantly involved in gastrointestinal disease, immunological disease, inflammatory response, cell-to-cell signaling and interaction, lymphoid tissue development, and cell death. By contrast, the 282 PS CNV regions contained 473 Ensembl genes which are overrepresented in environmental interactions. Network analyses indicated that the PS-specific genes were particularly enriched for inflammatory response, immune cell trafficking, metabolic disease, cell cycle, and cellular organization and movement.     

What a difference copy number variation makes

DNA copy number variation (CNV) represents a considerable source of human genetic diversity. Recently,1 a global map of copy number variation in the human genome has been drawn up which reveals not only the ubiquity but also the complexity of this type of variation. Thus, two human genomes may differ by more than 20 Mb and it is likely that the full extent of CNV still remains to be discovered. Nearly 3000 genes are associated with CNV. This high degree of variability with regard to gene copy number between two individuals challenges definitions of normality. Many CNVs are located in regions of complex genomic structure and this currently limits the extent to which these variants can be genotyped by using tagging SNPs. However, some CNVs are already amenable to genome-wide association studies so that their influence on human phenotypic diversity and disease susceptibility may soon be determined.     

Novel origins of copy number variation in the dog genome

ABSTRACT: BACKGROUND: Copy number variants (CNVs) account for substantial variation between genomes and are a major source of normal and pathogenic phenotypic differences. The dog is an ideal model to investigate mutational mechanisms that generate CNVs as its genome lacks a functional ortholog of the PRDM9 gene implicated in recombination and CNV formation in humans. Here we comprehensively assay CNVs using high-density array comparative genomic hybridization in 50 dogs from 17 dog breeds and 3 gray wolves. RESULTS: We use a stringent new method to identify a total of 430 high-confidence CNV loci, that range in size from 9 kb to 1.6 Mb and span 26.4 Mb, or 1.08%, of the assayed dog genome, overlapping 413 annotated genes. 98% of CNVs observed in each breed are also observed in multiple breeds. CNVs predicted to disrupt gene function are significantly less common than expected by chance. We identify a significant overrepresentation of peaks of GC content, previously shown to be enriched in dog recombination hotspots, in the vicinity of CNV breakpoints. CONCLUSIONS: A number of the CNVs identified by this study are candidates for generating breed-specific phenotypes. Purifying selection seems to be a major factor shaping structural variation in the dog genome, suggesting that many CNVs are deleterious. Localized peaks of GC content appear to be novel sites of CNV formation in the dog genome by non-allelic homologous recombination, potentially activated by the loss of PRDM9. These sequence features may have driven genome instability and chromosomal rearrangements throughout canid evolution.     

Gene Expression Profiling of Lymphoblasts from Autistic and Nonaffected Sib Pairs: Altered Pathways in Neuronal Development and Steroid Biosynthesis

Despite the identification of numerous autism susceptibility genes, the pathobiology of autism remains unknown. The present “case-control” study takes a global approach to understanding the molecular basis of autism spectrum disorders based upon large-scale gene expression profiling. DNA microarray analyses were conducted on lymphoblastoid cell lines from over 20 sib pairs in which one sibling had a diagnosis of autism and the other was not affected in order to identify biochemical and signaling pathways which are differentially regulated in cells from autistic and nonautistic siblings. Bioinformatics and gene ontological analyses of the data implicate genes which are involved in nervous system development, inflammation, and cytoskeletal organization, in addition to genes which may be relevant to gastrointestinal or other physiological symptoms often associated with autism. Moreover, the data further suggests that these processes may be modulated by cholesterol/steroid metabolism, especially at the level of androgenic hormones. Elevation of male hormones, in turn, has been suggested as a possible factor influencing susceptibility to autism, which affects ∼4 times as many males as females. Preliminary metabolic profiling of steroid hormones in lymphoblastoid cell lines from several pairs of siblings reveals higher levels of testosterone in the autistic sibling, which is consistent with the increased expression of two genes involved in the steroidogenesis pathway. Global gene expression profiling of cultured cells from ASD probands thus serves as a window to underlying metabolic and signaling deficits that may be relevant to the pathobiology of autism.     

Identification of a novel gene on chromosome 7q11.2 interrupted by a translocation breakpoint in a pair of autistic twins

We report here the identification and characterization of a novel gene (AUTS2) that spans the 7q11.2 breakpoint in a monozygotic twin pair concordant for autism and a t(7;20) (q11.2; p11.2) translocation. AUTS2 is 1.2 Mb and has 19 exons. The predicted protein is 1295 amino acids and does not correspond to any known protein. DNA sequence analysis of autism subjects and controls revealed 22 biallelic polymorphic sites. For all sites, both alleles were observed in both cases and controls. Thus no autism-specific mutation was observed. Association analysis with two exonic polymorphic sites and linkage analysis of four dinucleotide repeat markers, two within and two flanking AUTS2, was negative. Thus, although it is unlikely that AUTS2 is an autism susceptibility gene for idiopathic autism, it may be the gene responsible for the disorder in the twins studied here.     

The physical map of wheat chromosome 5DS revealed gene duplications and small rearrangements

Background

The substantially large bread wheat genome, organized into highly similar three sub-genomes, renders genomic research challenging. The construction of BAC-based physical maps of individual chromosomes reduces the complexity of this allohexaploid genome, enables elucidation of gene space and evolutionary relationships, provides tools for map-based cloning, and serves as a framework for reference sequencing efforts. In this study, we constructed the first comprehensive physical map of wheat chromosome arm 5DS, thereby exploring its gene space organization and evolution.

Results

The physical map of 5DS was comprised of 164 contigs, of which 45 were organized into 21 supercontigs, covering 176 Mb with an N50 value of 2,173 kb. Fifty-eight of the contigs were larger than 1 Mb, with the largest contig spanning 6,649 kb. A total of 1,864 molecular markers were assigned to the map at a density of 10.5 markers/Mb, anchoring 100 of the 120 contigs (>5 clones) that constitute ~95 % of the cumulative length of the map. Ordering of 80 contigs along the deletion bins of chromosome arm 5DS revealed small-scale breaks in syntenic blocks. Analysis of the gene space of 5DS suggested an increasing gradient of genes organized in islands towards the telomere, with the highest gene density of 5.17 genes/Mb in the 0.67-0.78 deletion bin, 1.4 to 1.6 times that of all other bins.

Conclusions

Here, we provide a chromosome-specific view into the organization and evolution of the D genome of bread wheat, in comparison to one of its ancestors, revealing recent genome rearrangements. The high-quality physical map constructed in this study paves the way for the assembly of a reference sequence, from which breeding efforts will greatly benefit.

Electronic supplementary material

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

Analysis of copy number variations in Holstein cows identify potential mechanisms contributing to differences in residual feed intake

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. In this study, we performed an initial analysis of copy number variations (CNVs) using BovineHD SNP genotyping data from 147 Holstein cows identified as having high or low feed efficiency as estimated by residual feed intake (RFI). We detected 443 candidate CNV regions (CNVRs) that represent 18.4?Mb (0.6?%) of the genome. To investigate the functional impacts of CNVs, we created two groups of 30 individual animals with extremely low or high estimated breeding values (EBVs) for RFI, and referred to these groups as low intake (LI; more efficient) or high intake (HI; less efficient), respectively. We identified 240 (~9.0?Mb) and 274 (~10.2?Mb) CNVRs from LI and HI groups, respectively. Approximately 30–40?% of the CNVRs were specific to the LI group or HI group of animals. The 240 LI CNVRs overlapped with 137 Ensembl genes. Network analyses indicated that the LI-specific genes were predominantly enriched for those functioning in the inflammatory response and immunity. By contrast, the 274 HI CNVRs contained 177 Ensembl genes. Network analyses indicated that the HI-specific genes were particularly involved in the cell cycle, and organ and bone development. These results relate CNVs to two key variables, namely immune response and organ and bone development. The data indicate that greater feed efficiency relates more closely to immune response, whereas cattle with reduced feed efficiency may have a greater capacity for organ and bone development.     

The gene distribution in the genomes of pea, tomato and date palm

The vast majority of genes of maize, rice, barley and wheat are contained in long gene-rich regions (collectively called the 'gene space') separated by long gene-empty regions. The gene space covers a narrow, 0.8-1.6%, GC range, possibly because of the presence of abundant transposons. Here we report that the gene space is not an exclusive property of Gramineae, because it also exists in the large genome of pea (5000 Mb). Moreover, the gene space is not just dependent upon genome size, since a gene space is found in rice (415 Mb), but not in Arabidopsis (120 Mb), nor in two other plants investigated in the present work, date palm (250 Mb) and tomato (1000 Mb).     

Genome-wide patterns of copy number variation in the diversified chicken genomes using next-generation sequencing

Background

Copy number variation (CNV) is important and widespread in the genome, and is a major cause of disease and phenotypic diversity. Herein, we performed a genome-wide CNV analysis in 12 diversified chicken genomes based on whole genome sequencing.

Results

A total of 8,840 CNV regions (CNVRs) covering 98.2 Mb and representing 9.4% of the chicken genome were identified, ranging in size from 1.1 to 268.8 kb with an average of 11.1 kb. Sequencing-based predictions were confirmed at a high validation rate by two independent approaches, including array comparative genomic hybridization (aCGH) and quantitative PCR (qPCR). The Pearson’s correlation coefficients between sequencing and aCGH results ranged from 0.435 to 0.755, and qPCR experiments revealed a positive validation rate of 91.71% and a false negative rate of 22.43%. In total, 2,214 (25.0%) predicted CNVRs span 2,216 (36.4%) RefSeq genes associated with specific biological functions. Besides two previously reported copy number variable genes EDN3 and PRLR, we also found some promising genes with potential in phenotypic variation. Two genes, FZD6 and LIMS1, related to disease susceptibility/resistance are covered by CNVRs. The highly duplicated SOCS2 may lead to higher bone mineral density. Entire or partial duplication of some genes like POPDC3 may have great economic importance in poultry breeding.

Conclusions

Our results based on extensive genetic diversity provide a more refined chicken CNV map and genome-wide gene copy number estimates, and warrant future CNV association studies for important traits in chickens.

Electronic supplementary material

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

Copy number variants at Williams–Beuren syndrome 7q11.23 region

Copy number variants (CNVs) of the Williams–Beuren syndrome (WBS) 7q11.23 region are responsible for neurodevelopmental disorders with multi-system involvement and variable expressivity. Typical features of WBS microdeletion comprise a recognizable pattern of facial dysmorphisms, supravalvular aortic stenosis, connective tissue abnormalities, hypercalcemia, and a distinctive neurobehavioral phenotype. Conversely, the phenotype of patients carrying the 7q11.23 reciprocal duplications includes less distinctive facial dysmorphisms and prominent speech delay. The common deletion/duplication ranges in size from 1.5 to 1.8 Mb and encompasses approximately 28 genes. This region is flanked by low copy repeats (LCRs) with greater than ~97% identity, which can mediate non-allelic homologous recombination resulting from misalignment of LCRs during meiosis. A clear genotype–phenotype correlation has been established in WBS only for the elastin gene, which is responsible for the vascular and connective tissue abnormalities. The molecular substrates underlying the other clinical features of 7q11.23 CNVs, including the neurocognitive phenotypes, are still debated. Recent studies suggest that besides the role of the genes in the deleted/duplicated interval, multiple factors such as regulatory sequences, epigenetic mechanisms, parental origin of the CNV, and nucleotide variations in the non-deleted/duplicated allele may be important in determining the variable expressivity of 7q11.23 CNV phenotypes. Here, we review the clinical and molecular findings and the recent insights on genomic disorders associated with CNVs involving the 7q11.23 region.     

Social approach in genetically engineered mouse lines relevant to autism

Profound impairment in social interaction is a core symptom of autism, a severe neurodevelopmental disorder. Deficits can include a lack of interest in social contact and low levels of approach and proximity to other children. In this study, a three-chambered choice task was used to evaluate sociability and social novelty preference in five lines of mice with mutations in genes implicated in autism spectrum disorders. Fmr1^tm1Cgr/Y ( Fmr1^−/y ) mice represent a model for fragile X, a mental retardation syndrome that is partially comorbid with autism. We tested Fmr1^−/y mice on two genetic backgrounds, C57BL/6J and FVB/N-129/OlaHsd (FVB/129). Targeted disruption of Fmr1 resulted in low sociability on one measure, but only when the mutation was expressed on FVB/129. Autism has been associated with altered serotonin levels and polymorphisms in SLC6A4 (SERT) , the serotonin transporter gene. Male mice with targeted disruption of Slc6a4 displayed significantly less sociability than wild-type controls. Mice with conditional overexpression of Igf-1 (insulin-like growth factor-1) offered a model for brain overgrowth associated with autism. Igf-1 transgenic mice engaged in levels of social approach similar to wild-type controls. Targeted disruption in other genes of interest, En2 (engrailed-2) and Dhcr7 , was carried on genetic backgrounds that showed low levels of exploration in the choice task, precluding meaningful interpretations of social behavior scores. Overall, results show that loss of Fmr1 or Slc6a4 gene function can lead to deficits in sociability. Findings from the fragile X model suggest that the FVB/129 background confers enhanced susceptibility to consequences of Fmr1 mutation on social approach.     

Molecular cytogenetic analysis and resequencing of contactin associated protein-like 2 in autism spectrum disorders

Autism spectrum disorders (ASD) are a group of related neurodevelopmental syndromes with complex genetic etiology. We identified a de novo chromosome 7q inversion disrupting Autism susceptibility candidate 2 (AUTS2) and Contactin Associated Protein-Like 2 (CNTNAP2) in a child with cognitive and social delay. We focused our initial analysis on CNTNAP2 based on our demonstration of disruption of Contactin 4 (CNTN4) in a patient with ASD; the recent finding of rare homozygous mutations in CNTNAP2 leading to intractable seizures and autism; and in situ and biochemical analyses reported herein that confirm expression in relevant brain regions and demonstrate the presence of CNTNAP2 in the synaptic plasma membrane fraction of rat forebrain lysates. We comprehensively resequenced CNTNAP2 in 635 patients and 942 controls. Among patients, we identified a total of 27 nonsynonymous changes; 13 were rare and unique to patients and 8 of these were predicted to be deleterious by bioinformatic approaches and/or altered residues conserved across all species. One variant at a highly conserved position, I869T, was inherited by four affected children in three unrelated families, but was not found in 4010 control chromosomes (p = 0.014). Overall, this resequencing data demonstrated a modest nonsignificant increase in the burden of rare variants in cases versus controls. Nonetheless, when viewed in light of two independent studies published in this issue of AJHG showing a relationship between ASD and common CNTNAP2 alleles, the cytogenetic and mutation screening data suggest that rare variants may also contribute to the pathophysiology of ASD, but place limits on the magnitude of this contribution.     

Estimation of association of <Emphasis Type="Italic">CNTN6</Emphasis> copy number variation with idiopathic intellectual disability

Analysis of the prevalence of copy number variations of the CNTN6 gene, recently selected as a new candidate gene for intellectual disorders, was performed. Real-time PCR did not detect any change in the number of CNTN6 gene copies in a group of 200 patients with impaired intellectual development. However, taking into account our data from the previous aCGH analysis and published data, the overall frequency of microdeletions and microduplications of CNTN6 was estimated as 1: 265 (0.4%). The common phenotypic features of 40 patients with microdeletions and microduplications of CNTN6 appeared to be the autism spectrum disorders, developmental delay, intellectual disability, seizures, cognitive impairment, cardiological defects, and behavioral problems.     

Deletion 2q37.3 and autism: molecular cytogenetic mapping of the candidate region for autistic disorder

Fine mapping of deletion regions in autistic patients represents a valuable screening tool for identifying candidate genes for autism. A number of studies have ascertained associations between autism and terminal 2q deletion with the breakpoint within 2q37. Here we describe a 12-year-old female patient with terminal 2q37.3 cryptic deletion and autistic behaviour. Her clinical features included hypotonia and feeding difficulties during infancy, coarse face with notably prominent forehead, prominent eyebrows, broad flat nasal bridge and round cheeks, small hands and feet with bilateral brachymetaphalangism, proximal implantation of the thumbs and short toenails, mild mental retardation and autistic behaviour. Recorded autistic features included early lack of eye contact and, during infancy, little social interactions, propensity to be stereotypically busy and to get anxious. In order to more closely delineate the linkage region for autism within 2q37, the findings in this patient were combined to those in 2 previously reported siblings with a well documented 2q37.3 deletion, but without autistic disorder. The exact size of the deleted segment was determined by mapping the deleted region in each group with a series of specific BAC clones linearly ordered on the 2q37 region. The deletion in the autistic patient appeared to be larger [breakpoint flanked by more centromeric clones RP11-680016 (236.9 Mb) and 201F21 (237.4 Mb)] than in the non autistic siblings [more telomeric clones RP11-205L13 (237.8 Mb) and 346114 (238.2 Mb)], revealing a distance of maximum 1.3 Mb between the breakpoints. Accordingly, the extent of the candidate region for susceptibility genes for autism on distal 2q is reduced to maximum 1.3 Mb. Comparison with another well documented autistic patient from the literature results in the same conclusion. These findings represent thus a further step towards identifying genes predisposing to autism.     

Genome-wide association study of copy number variants suggests LTBP1 and FGD4 are important for alcohol drinking

Alcohol dependence (AD) is a complex disorder characterized by psychiatric and physiological dependence on alcohol. AD is reflected by regular alcohol drinking, which is highly inheritable. In this study, to identify susceptibility genes associated with alcohol drinking, we performed a genome-wide association study of copy number variants (CNVs) in 2,286 Caucasian subjects with Affymetrix SNP6.0 genotyping array. We replicated our findings in 1,627 Chinese subjects with the same genotyping array. We identified two CNVs, CNV207 (combined p-value 1.91E-03) and CNV1836 (combined p-value 3.05E-03) that were associated with alcohol drinking. CNV207 and CNV1836 are located at the downstream of genes LTBP1 (870 kb) and FGD4 (400 kb), respectively. LTBP1, by interacting TGFB1, may down-regulate enzymes directly participating in alcohol metabolism. FGD4 plays a role in clustering and trafficking GABA(A) receptor and subsequently influence alcohol drinking through activating CDC42. Our results provide suggestive evidence that the newly identified CNV regions and relevant genes may contribute to the genetic mechanism of alcohol dependence.     

Defining the Contribution of CNTNAP2 to Autism Susceptibility

Multiple lines of genetic evidence suggest a role for CNTNAP2 in autism. To assess its population impact we studied 2148 common single nucleotide polymorphisms (SNPs) using transmission disequilibrium test (TDT) across the entire ~3.3 Mb CNTNAP2 locus in 186 (408 trios) multiplex and 323 simplex families with autistic spectrum disorder (ASD). This analysis yielded two SNPs with nominal statistical significance (rs17170073, p = 2.0 x 10^-4; rs2215798, p = 1.6 x 10^-4) that did not survive multiple testing. In a combined analysis of all families, two highly correlated (r ² = 0.99) SNPs in intron 14 showed significant association with autism (rs2710093, p = 9.0 x 10^-6; rs2253031, p = 2.5 x 10^-5). To validate these findings and associations at SNPs from previous autism studies (rs7794745, rs2710102 and rs17236239) we genotyped 2051 additional families (572 multiplex and 1479 simplex). None of these variants were significantly associated with ASD after corrections for multiple testing. The analysis of Mendelian errors within each family did not indicate any segregating deletions. Nevertheless, a study of CNTNAP2 gene expression in brains of autistic patients and of normal controls, demonstrated altered expression in a subset of patients (p = 1.9 x10^-5). Consequently, this study suggests that although CNTNAP2 dysregulation plays a role in some cases, its population contribution to autism susceptibility is limited.     

A first comparative map of copy number variations in the sheep genome

We carried out a cross species cattle-sheep array comparative genome hybridization experiment to identify copy number variations (CNVs) in the sheep genome analysing ewes of Italian dairy or dual-purpose breeds (Bagnolese, Comisana, Laticauda, Massese, Sarda, and Valle del Belice) using a tiling oligonucleotide array with ~385,000 probes designed on the bovine genome. We identified 135 CNV regions (CNVRs; 24 reported in more than one animal) covering ~10.5 Mb of the virtual sheep genome referred to the bovine genome (0.398%) with a mean and a median equal to 77.6 and 55.9 kb, respectively. A comparative analysis between the identified sheep CNVRs and those reported in cattle and goat genomes indicated that overlaps between sheep and both other species CNVRs are highly significant (P<0.0001), suggesting that several chromosome regions might contain recurrent interspecies CNVRs. Many sheep CNVRs include genes with important biological functions. Further studies are needed to evaluate their functional relevance.