Copy number variations at the Prader-Willi syndrome region on chromosome 15 and associations with obesity in whites

Obesity is a serious health problem with strong genetic determination. Copy number variation (CNV) is a common type of genomic variant associated with some complex human diseases. However, it is not clear how CNVs contribute to the etiology of obesity. In this study, we examined 1,000 unrelated US whites to search for CNVs that may predispose to obesity. We focused our analyses on the Prader‐Willi syndrome (PWS) critical region (chromosome 15q11–q13), because the PWS region is a hotspot for CNV generation and obesity is one of the major clinical manifestations for chromosome abnormalities at this region. We constructed a map containing 39 CNVs at the PWS critical region with CNV occurrence rates higher than 1%. Among them, three CNVs were significantly associated with body fat mass (P < 0.05), with a higher copy number (CN) associated with an increase of 5.08–9.77 kg in body fat mass. These three CNVs are close to two known PWS genes, NDN (necdin homolog) and C15orf2 (chromosome 15 open reading frame 2), and partially overlap with another obesity gene PWRN1 (Prader‐Willi region nonprotein‐coding RNA 1). Interestingly, our recently published whole genome association scan study using the same sample by examining single‐nucleotide polymorphisms (SNPs) did not find any significant associations at these CNV regions, suggesting the importance of examining both CNVs and SNPs for better understanding of genetic basis of obesity. Further studies are warranted to validate these CNVs and their importance to obesity.     

Analysis of structural diversity in wolf-like canids reveals post-domestication variants

Background

Although a variety of genetic changes have been implicated in causing phenotypic differences among dogs, the role of copy number variants (CNVs) and their impact on phenotypic variation is still poorly understood. Further, very limited knowledge exists on structural variation in the gray wolf, the ancestor of the dog, or other closely related wild canids. Documenting CNVs variation in wild canids is essential to identify ancestral states and variation that may have appeared after domestication.

Results

In this work, we genotyped 1,611 dog CNVs in 23 wolf-like canids (4 purebred dogs, one dingo, 15 gray wolves, one red wolf, one coyote and one golden jackal) to identify CNVs that may have arisen after domestication. We have found an increase in GC-rich regions close to the breakpoints and around 1 kb away from them suggesting that some common motifs might be associated with the formation of CNVs. Among the CNV regions that showed the largest differentiation between dogs and wild canids we found 12 genes, nine of which are related to two known functions associated with dog domestication; growth (PDE4D, CRTC3 and NEB) and neurological function (PDE4D, EML5, ZNF500, SLC6A11, ELAVL2, RGS7 and CTSB).

Conclusions

Our results provide insight into the evolution of structural variation in canines, where recombination is not regulated by PRDM9 due to the inactivation of this gene. We also identified genes within the most differentiated CNV regions between dogs and wolves, which could reflect selection during the domestication process.

Electronic supplementary material

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

Inherited deletion of 9p22.3-p24.3 and duplication of 18p11.31-p11.32 associated with neurodevelopmental delay: Phenotypic matching of involved genes

We describe a 3.5-year-old Iranian female child and her affected 10-month-old brother with a maternally inherited derivative chromosome 9 [der(9)]. The postnatally detected rearrangement was finely characterized by aCGH analysis, which revealed a 15.056 Mb deletion of 9p22.3-p24.3p22.3 encompassing 14 OMIM morbid genes such as DOCK8, KANK1, DMRT1 and SMARCA2, and a gain of 3.309 Mb on 18p11.31-p11.32 encompassing USP14, THOC1, COLEC12, SMCHD1 and LPIN2. We aligned the genes affected by detected CNVs to clinical and functional phenotypic features using PhenogramViz. In this regard, the patient's phenotype and CNVs data were entered into PhenogramViz. For the 9p deletion CNV, 53 affected genes were identified and 17 of them were matched to 24 HPO terms describing the patient's phenotypes. Also, for CNV of 18p duplication, 22 affected genes were identified and six of them were matched to 13 phenotypes. Moreover, we used DECIPHER for in-depth characterization of involved genes in detected CNVs and also comparison of patient phenotypes with 9p and 18p genomic imbalances. Based on our filtration strategy, in the 9p22.3-p24.3 region, approximately 80 pathogenic/likely pathogenic/uncertain overlapping CNVs were in DECIPHER. The size of these CNVs ranged from 12.01 kb to 18.45 Mb and 52 CNVs were smaller than 1 Mb in size affecting 10 OMIM morbid genes. The 18p11.31-p11.32 region overlapped 19 CNVs in the DECIPHER database with the size ranging from 23.42 kb to 1.82 Mb. These CNVs affect eight haploinsufficient genes.     

Genome wide CNV analysis reveals additional variants associated with milk production traits in Holsteins

Background

Milk production is an economically important sector of global agriculture. Much attention has been paid to the identification of quantitative trait loci (QTL) associated with milk, fat, and protein yield and the genetic and molecular mechanisms underlying them. Copy number variation (CNV) is an emerging class of variants which may be associated with complex traits.

Results

In this study, we performed a genome-wide association between CNVs and milk production traits in 26,362 Holstein bulls and cows. A total of 99 candidate CNVs were identified using Illumina BovineSNP50 array data, and association tests for each production trait were performed using a linear regression analysis with PCA correlation. A total of 34 CNVs on 22 chromosomes were significantly associated with at least one milk production trait after false discovery rate (FDR) correction. Some of those CNVs were located within or near known QTL for milk production traits. We further investigated the relationship between associated CNVs with neighboring SNPs. For all 82 combinations of traits and CNVs (less than 400 kb in length), we found 17 cases where CNVs directly overlapped with tag SNPs and 40 cases where CNVs were adjacent to tag SNPs. In 5 cases, CNVs located were in strong linkage disequilibrium with tag SNPs, either within or adjacent to the same haplotype block. There were an additional 20 cases where CNVs did not have a significant association with SNPs, suggesting that the effects of those CNVs were probably not captured by tag SNPs.

Conclusion

We conclude that combining CNV with SNP analyses reveals more genetic variations underlying milk production traits than those revealed by SNPs alone.

Electronic supplementary material

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

An efficient method for measuring copy number variation applied to improvement of nematode resistance in soybean

Copy number variation (CNV) is implicated in important traits in multiple crop plants, but can be challenging to genotype using conventional methods. The Rhg1 locus of soybean, which confers resistance to soybean cyst nematode (SCN), is a CNV of multiple 31.2‐kb genomic units each containing four genes. Reliable, high‐throughput methods to quantify Rhg1 and other CNVs for selective breeding were developed. The CNV genotyping assay described here uses a homeologous gene copy within the paleopolyploid soybean genome to provide the internal control for a single‐tube TaqMan copy number assay. Using this assay, CNV in breeding populations can be tracked with high precision. We also show that extensive CNV exists within Fayette, a released, inbred SCN‐resistant soybean cultivar with a high copy number at Rhg1 derived from a single donor parent. Copy number at Rhg1 is therefore unstable within a released variety over a relatively small number of generations. Using this assay to select for individuals with altered copy number, plants were obtained with both increased copy number and increased SCN resistance relative to control plants. Thus, CNV genotyping technologies can be used as a new type of marker‐assisted selection to select for desirable traits in breeding populations, and to control for undesirable variation within cultivars.     

A loss-of-function mutant allele of a glycosyl hydrolase gene has been co-opted for seed weight control during soybean domestication

The resultant DNA from loss-of-function mutation can be recruited in biological evolution and development. Here, we present such a rare and potential case of “to gain by loss” as a neomorphic mutation during soybean domestication for increasing seed weight. Using a population derived from a chromosome segment substitution line of Glycine max(SN14) and Glycine soja(ZYD06), a quantitative trait locus(QTL) of 100-seed weight(q HSW) was mapped on chromosome 11, corresponding to a truncated β-1, 3-gl...     

A genome‐wide association study of copy number variations with umbilical hernia in swine

Umbilical hernia (UH) is one of the most common congenital defects in pigs, leading to considerable economic loss and serious animal welfare problems. To test whether copy number variations (CNVs) contribute to pig UH, we performed a case–control genome‐wide CNV association study on 905 pigs from the Duroc, Landrace and Yorkshire breeds using the Porcine SNP60 BeadChip and penncnv algorithm. We first constructed a genomic map comprising 6193 CNVs that pertain to 737 CNV regions. Then, we identified eight CNVs significantly associated with the risk for UH in the three pig breeds. Six of seven significantly associated CNVs were validated using quantitative real‐time PCR. Notably, a rare CNV (CNV14:13030843–13059455) encompassing the NUGGC gene was strongly associated with UH (permutation‐corrected P = 0.0015) in Duroc pigs. This CNV occurred exclusively in seven Duroc UH‐affected individuals. SNPs surrounding the CNV did not show association signals, indicating that rare CNVs may play an important role in complex pig diseases such as UH. The NUGGC gene has been implicated in human omphalocele and inguinal hernia. Our finding supports that CNVs, including the NUGGC CNV, contribute to the pathogenesis of pig UH.     

Nucleotide diversity of the upstream region of the putative MADS-box gene controlling soybean maturity

MADS-box genes are involved in plant reproductive development. However, the role of gene nucleotide diversity in soybean flowering and maturity remains unknown. Therefore, in this study, the distribution of DNA polymorphisms in the putative MADS-box gene located near the quantitative trait loci (QTL) for flowering time and maturity was targeted for association analysis using Glycine max (cultivated soybean) and Glycine soja (wild soybean). Sixteen single nucleotide polymorphisms identified in the upstream region of the putative MADS-box gene around QTL Pod mat 13-7 and Fflr 4-2 on chromosome 7 were found to be highly associated with maturity in soybean. The genetic diversity between cultivated soybeans and the wild relative was comparable, although the early maturity group (EMG) was less diverse than the late maturity group (LMG) of the cultivated soybean. Population size changes of the MADS-box gene in this soybean germplasm appeared to result from non-random selection. A selective pressure seemed to act on this gene in the EMG, while the LMG and G. soja were in genetic equilibrium. Neutrality tests and the constructed neighbor-joining tree indicate that the EMG of G. max has experienced strong artificial selection for its domestication and genetic improvement.     

Genome-wide CNV analysis revealed variants associated with growth traits in African indigenous goats

Using the CaprineSNP50 data generated by the AGIN consortium, we detected common CNVs in 126 samples from four African indigenous goat breeds. A total of 30 CNVs ranging from 30,237 bp to 4,910,757 bp were identified. These CNVs were then associated with six growth traits by a linear regression analysis. Three significant associations were identified between two CNVs and two body traits after false discovery rate (FDR) correction (P < .05). One of them (CNV27) was significantly associated with both chest width and width of pin bones. It overlaps the SNX29 gene, the Gene Ontology (GO) annotations of which indicate CNV27 could be a potential functional candidate for meat production, health and reproduction traits. To our knowledge, this study is the first CNV-based association test of growth traits using SNP chip data in African meat goats.     

A novel multiplex fluorescent competitive PCR for copy number variation detection

The copy number variation (CNV) is an important genetic marker in cancer and other diseases. To detect CNVs of specific genetic loci, the multiplex ligation-dependent probe amplification (MLPA) is an appropriate approach, but the experimental optimization and probe synthesis are still great challenges. The multiplex competitive PCR is an alternative method for CNV detection. However, the construction of internal competitive template and establishment of a stable multiplex PCR system are the main limiting factors for this method. Here, we introduce a novel multiplex fluorescent competitive PCR (NMFC-PCR) for detecting CNVs. In this method, the blunt hairpin primers are used to rapidly establish a stable multiplex PCR system due to the reduction of non-specific amplification, and limited cycles' amplification is used to obtain the internal competitive template instead of artificial synthesis. With this method, we tested 21 clinical samples with potential LIM homeobox 1 (LHX1) or T-box 6 (TBX6) deletion. Every three segments located on the LHX1 and TBX6 were selected as the target regions, while two segments located on X-chromosome and five segments located on autosome were selected as the reference regions for detecting CNVs. The results showed that the gender information of 21 samples can be accurately inferred by the copy number ratio (CNR) of X-chromosomal reference region to autosomal reference region (X/A), and 2 samples had one copy of LHX1 and 9 samples had one copy of TBX6. To evaluate the accuracy of NMFC-PCR, 5 random samples with CNV were also detected by array-based comparative genomic hybridization (aCGH), and the results of aCGH were consistent with the NMFC-PCR results. To further assess the performance of NMFC-PCR, 60 normal samples were simultaneously tested. The results showed that the gender results were exactly the same as known information, and CNVs of LHX1 or TBX6 were not found. In conclusion, the method is a cheap, efficient, accurate, and convenient competitive PCR method for CNV detection.     

Copy number variants in Italian Large White pigs detected using high‐density single nucleotide polymorphisms and their association with back fat thickness

The aim of this study was to identify copy number variants (CNVs) in Italian Large White pigs and test them for association with back fat thickness (BFT). Within a population of 12 000 performance‐tested pigs, two groups of animals with extreme and divergent BFT estimated breeding values (EBVs; 147 with negative and 150 with positive EBVs) were genotyped with the Illumina Porcine SNP60 BeadChip. CNVs were detected with penncnv software. We identified a total of 4146 CNV events in 170 copy number variation regions (CNVRs) located on 15 porcine autosomes. Validation of detected CNVRs was carried out (i) by comparing CNVRs already detected by other studies and (ii) by semiquantitative fluorescent multiplex (SQFM) PCR of a few CNVRs. Most of CNVRs detected in Italian Large White pigs (71.2%) were already reported in other pig breeds/populations, and 82.1% of the CNV events detected by penncnv were confirmed by SQFM PCR. For each CNVR, we compared the occurrence of CNV events between the pigs of the high and low BFT EBV tails. Sixteen regions showed significance at P < 0.10, and seven were significant at P < 0.05 but were not significant after Bonferroni correction (Fisher's exact test). These results indicated that CNVs could explain a limited fraction of the genetic variability of fat deposition in Italian Large White pigs. However, it was interesting to note that one of these CNVRs encompassed the ZPLD1 gene. In humans, a rare CNV event including this gene is associated with obesity. Studies identifying CNVs in pigs could assist in elucidating the genetic mechanisms underlying human obesity.     

Analysis of chromosomal structural variation in patients with congenital left‐sided cardiac lesions

Background: We sought to characterize the landscape of structural variation associated with the subset of congenital cardiac defects characterized by left‐sided obstruction. Methods: Cases with left‐sided cardiac defects (LSCD) and pediatric controls were uniformly genotyped and assessed for copy number variant (CNV) calls. Significance testing was performed to ascertain differences in overall CNV incidence, and for CNV enrichment of specific genes and gene functions in LSCD cases relative to controls. Results: A total of 257 cases of European descent and 962 ethnically matched, disease‐free pediatric controls were included. Although there was no difference in CNV rate between cases and controls, a significant enrichment in rare LSCD CNVs was detected overall (p = 7.30 × 10^?3, case/control ratio = 1.26) and when restricted either to deletions (p = 7.58 × 10^?3, case/control ratio = 1.20) or duplications (3.02 × 10^?3, case/control ratio = 1.43). Neither gene‐based, functional nor knowledge‐based analyses identified genes, loci or pathways that were significantly enriched in cases as compared to controls when appropriate corrections for multiple tests were applied. However, several genes of interest were identified by virtue of their association with cardiac development, known human conditions, or reported disruption by CNVs in other patient cohorts. Conclusion: This study examines the largest cohort to date with LSCD for structural variation. These data suggest that CNVs play a role in disease risk and identify numerous genes disrupted by CNVs of potential disease relevance. These findings further highlight the genetic heterogeneity and complexity of these disorders. Birth Defects Research (Part A) 100:951–964, 2014. © 2014 Wiley Periodicals, Inc.     

Domestication‐related variation in social preferences in chickens is affected by genotype on a growth QTL

A growth‐related QTL on chicken chromosome 1 has previously been shown to influence domestication behaviour in chickens. In this study, we used Red Junglefowl (RJF) and White Leghorn (WL) as well as the intercross between them to investigate whether stress affects the way birds allocate their time between familiar and unfamiliar conspecifics in a social preference test (‘social support seeking’), and how this is related to genotype at specific loci within the growth QTL. Red Junglefowl males spent more time with unfamiliar chickens before the stressful event compared to the other birds, whereas all birds except WL males tended to spend less time with unfamiliar ones after stress. A significant QTL locus was found to influence both social preference under undisturbed circumstances and social support seeking. The WL allele at this QTL was associated not only with a preference for unfamiliar individuals but also with a shift towards familiar ones in response to stress (social support seeking). A second, suggestive QTL also affected social support seeking, but in the opposite direction; the WL allele was associated with increased time spent with unfamiliar individuals. The region contains several possible candidate genes, and gene expression analysis of a number of them showed differential expression between RJF and WL of AVPR2 (receptor for vasotocin), and possibly AVPR1a (another vasotocin receptor) and NRCAM (involved in neural development) in the lower frontal lobes of the brains of RJF and WL animals. These three genes continue to be interesting candidates for the observed behavioural effects .     

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.     

Mapping of QTL associated with chilling tolerance during reproductive growth in soybean

Low temperatures in summer bring about drastic reduction in seed yield of soybean [Glycine max (L.) Merr.]. To identify quantitative trait loci (QTL) associated with chilling tolerance during the reproductive growth in soybean, a recombinant inbred line (RIL) population consisting of 104 F₆-derived lines was created from a cross between two cultivars, chilling-tolerant Hayahikari and chilling-sensitive Toyomusume. The RIL were genotyped with 181 molecular and phenotypic markers and were scored with regard to chilling tolerance, which was evaluated by comparison of seed-yielding abilities in two artificial climatic environments at chilling and usual temperatures. Three QTL were detected for chilling tolerance in seed-yielding ability. Two of them, qCTTSW1 and qCTTSW2, were mapped near QTL for flowering time, and the latter had an epistatic interaction with a marker locus located near another QTL for flowering time, where no significant QTL for chilling tolerance was detected. The analysis of an F₂ population derived from the cross between Hayahikari and an RIL of the Hayahikari genotype at all QTL for flowering time confirmed the effect of the third QTL, qCTTSW3, on chilling tolerance and suggested that qCTTSW1 was basically independent of the QTL for flowering time. The findings and QTL found in this study may provide useful information for marker-assisted selection (MAS) and further genetic studies on soybean chilling tolerance.     

Dynamic QTL analysis of linolenic acid content in different developmental stages of soybean seed

Linolenic acid (LN) in soybean (Glycine max L. Merr.) seed mainly contributes to the undesirable odors and flavors commonly associated with poor oil quality. LN deposition at various stages of soybean seed development had not been reported by 2010. The objects of this study were (1) to identify and measure quantitative trait loci (QTL) underlying LN content and (2) to estimate the QTL effects expressed from earlier seed developmental stages to drying seed of soybean. One hundred and twenty-five F_5:8 and F_5:9 recombinant inbred lines derived from the cross of soybean cultivars ‘Hefeng 25’ and ‘Dongnong L5’ were used for the identification of QTL underlying LN content from the 37 day (D) to 86D stages after flowering, at Harbin in 2008 and 2009. QTL × Environment interactions (QE) effects were evaluated using a mixed genetic model (Zhu in J Zhejiang Univ (Natural Science) 33:327–335, 1999). Twelve unconditional QTL and 12 conditional QTL associated with LN content were identified at different developmental stages. Most of the QTL explained <10% of phenotypic variation of LN content. Unconditional QTL QLNF-1, QLNC2-1, QLND1b-1, QLNA2-1 and QLNH-1 influenced LN content across different development stages and environments. Conditional QTL QLNF-1, QLNC2-1 and QLNH-1 were identified in multiple developmental stages and environments. Conditional and unconditional QTL clustered in neighboring intervals on linkage groups A2, C2 and D1b. Ten QTL with conditional additive main effects (a) and/or conditional additive × environment interaction effects (ae) at specific developmental stage were identified on nine linkage groups. Of them, six QTL only possessed additive main effects and seven QTL had significant ae effects in different developmental stages. A total of 13 epistatic pairwise QTL were identified by conditional mapping in different developmental stages. Two pairs of QTL only showed aa effects and five pairs of QTL only showed aae effects at different developmental stages. QTL with aa effects, as well as their environmental interaction effects, appeared to vary at different developmental stages.     

Identification of copy number variation hotspots in human populations

Copy number variants (CNVs) in the human genome contribute to both Mendelian and complex traits as well as to genomic plasticity in evolution. The investigation of mutational rates of CNVs is critical to understanding genomic instability and the etiology of the copy number variation (CNV)-related traits. However, the evaluation of the CNV mutation rate at the genome level poses an insurmountable practical challenge that requires large samples and accurate typing. In this study, we show that an approximate estimation of the CNV mutation rate could be achieved by using the phylogeny information of flanking SNPs. This allows a genome-wide comparison of mutation rates between CNVs with the use of vast, readily available data of SNP genotyping. A total of 4187 CNV regions (CNVRs) previously identified in HapMap populations were investigated in this study. We showed that the mutation rates for the majority of these CNVRs are at the order of 10⁻⁵ per generation, consistent with experimental observations at individual loci. Notably, the mutation rates of 104 (2.5%) CNVRs were estimated at the order of 10⁻³ per generation; therefore, they were identified as potential hotspots. Additional analyses revealed that genome architecture at CNV loci has a potential role in inciting mutational hotspots in the human genome. Interestingly, 49 (47%) CNV hotspots include human genes, some of which are known to be functional CNV loci (e.g., CNVs of C4 and β-defensin causing autoimmune diseases and CNVs of HYDIN with implication in control of cerebral cortex size), implicating the important role of CNV in human health and evolution, especially in common and complex diseases.     

Copy number variation in the cattle genome

Copy number variations (CNVs) are gains and losses of genomic sequence greater than 50?bp between two individuals of a species. While single nucleotide polymorphisms (SNPs) are more frequent, CNVs impact a higher percentage of genomic sequence and have potentially greater effects, including the changing of gene structure and dosage, altering gene regulation and exposing recessive alleles. In particular, segmental duplications (SDs) were shown to be one of the catalysts and hotspots for CNV formation. Substantial progress has been made in understanding CNVs in mammals, especially in humans and rodents. CNVs have been shown to be important in both normal phenotypic variability and disease susceptibility. Recently, interest in CNV study has extended into domesticated animals, including cattle. Multiple genome-wide cattle CNV studies have been carried out using both microarray and next generation sequencing technologies. Integration of SD and CNV results with SNP and other datasets are beginning to reveal impacts of CNVs on cattle domestication, health, and production traits.