Recurrent Rearrangements of Human Amylase Genes Create Multiple Independent CNV Series

Identification of the genetic characteristics of copy number variations in experimental specific pathogen-free ducks using whole-genome resequencing

BACKGROUND

Specific pathogen-free ducks are a valuable laboratory resource for waterfowl disease research and poultry vaccine development. High throughput sequencing allows the systematic identification of structural variants in genomes. Copy number variation (CNV) can explain the variation of important duck genetic traits. Herein, the genome-wide CNVs of the three experimental duck species in China (Jinding ducks (JD), Shaoxing ducks (SX), and Fujian Shanma ducks (SM)) were characterized using resequencing to determine their genetic characteristics and selection signatures.

RESULTS

We obtained 4,810 CNV regions (CNVRs) by merging 73,012 CNVs, covering 4.2% of the duck genome. Functional analysis revealed that the shared CNVR-harbored genes were significantly enriched for 31 gene ontology terms and 16 Kyoto Encyclopedia of Genes and Genomes pathways (e.g., olfactory transduction and immune system). Based on the genome-wide fixation index for each CNVR, growth (SPAG17 and PTH1R), disease resistance (CATHL3 and DMBT1), and thermoregulation (TRPC4 and SLIT3) candidate genes were identified in strongly selected signatures specific to JD, SM, and SX, respectively.

CONCLUSIONS

In conclusion, we investigated the genome-wide distribution of experimental duck CNVs, providing a reference to establish the genetic basis of different phenotypic traits, thus contributing to the management of experimental animal genetic resources.

Analysis of the Selection Signal of the Tibetan Black Chicken Genome Based on Whole-Genome Sequencing

BACKGROUND

The Tibetan chicken has adapted well to high altitudes genetically after its long-term habitation in the plateau. In this study, we analyzed the selection signal of Tibetan black chickens (TBCs) and discovered genes associated with the characteristics of germplasm.

METHODS

Whole-genome sequencing (WGS) was used to identify the single-nucleotide polymorphism (SNP) markers and genetic structures in the genome of Tibetan black chickens. Further, we performed a comparative population genomics analysis between the genomic data obtained in this present study and the genomic data for five wild red jungle fowls (RJFs) accessed from the NCBI database (GenBank accession number PRJNA241474). Thereafter, the Fst and Pi selections were used to identify genes under positive selection in the Tibetan black chicken genome.

RESULTS

A total of 9,490,690 SNPs were identified in the Tibetan black chickens. In addition, the results from the gene ontology (GO) analysis showed that 732 genes of TBCs were enriched in a total of 210 GO terms with specific molecular functions such as regulation of cellular catabolic process, the MAPK signaling pathway, regulation of ion transport, growth, morphogenesis and lung alveolus development which may provide a better mechanism to facilitate oxygen transport and utilization in TBCs. Moreover, the results from the KEGG analysis showed that 732 genes of the TBCs were significantly enriched in the calcium signaling pathway, circadian entrainment (ADCY1, GNG7 and PER3), oxytocin signaling pathway and pathways of multiple neurodegeneration diseases. In addition, the CD86 antigen (CD86) was identified as a gene associated with the immune response in chickens. It was also revealed that genes such as TRIT1, HPCAL4, NT5C1A and HEYL were discovered under selection in Tibetan black chickens on chromosome 23. These genes may be related to the local adaptive characteristics of Tibetan black chickens, for instance, NT5C1A and HEYL may be involved in the high-altitude adaption of oxygen delivery in Tibetan black chickens.

CONCLUSIONS

In summary, we found that selection mainly affects the disease resistance and cold acclimatization of Tibetan black chickens. Hence, these results may provide important genetic information for the evolution and breeding of Tibetan black chickens.

Segregation analysis identifies specific alpha-defensin (DEFA1A3) SNP-CNV haplotypes in predisposition to IgA nephropathy

BACKGROUND

Immunoglobulin A (IgA) nephropathy is a disorder of the immune system affecting kidney function, and genome-wide association studies (GWAS) have defined numerous loci with associated variation, all implicating components of innate or adaptive immunity. Among these, single nucleotide polymorphisms (SNPs) in a region including the multiallelic copy number variation (CNV) of DEFA1A3 are associated with IgA nephropathy in both European and Asian populations. At present, the precise factors underlying the observed associations at DEFA1A3 have not been defined, although the key alleles differ between Asian and European populations, and multiple independent factors may be involved even within a single population.

METHODS

In this study, we measured DEFA1A3 copy number in UK family trios with an offspring affected by IgA nephropathy, used the population distributions of joint SNP-CNV haplotypes to infer the likely segregation in trios, and applied transmission disequilibrium tests (TDT) to examine joint SNP-CNV haplotypes for over- or undertransmission into affected offspring from heterozygous parents.

RESULTS AND CONCLUSIONS

We observed overtransmission of 3-copy class 2 haplotypes (raw p = 0.029) and some evidence for under-transmission of 3-copy class 1 haplotypes (raw p = 0.051), although these apparent effects were not statistically significant after correction for testing of multiple haplotypes.

Molecular combing and its application in clinical settings

Molecular combing technology (MCT) is an effective means for stretching DNA molecules and making them thus accessible for in situ studies. MCT uses the force exerted in the process of liquid flow via surface tension to stretch DNA molecules and spread them on solid surfaces, i.e. glass cover slips. Many DNA molecules can be stretched at the same time in parallel and neatly arranged side-by-side, making the approach convenient for statistical analysis. Accordingly, DNA replication and transcription can be studied at the single molecule level. In this paper, the principle, experimental methods, important applications, advantages and shortcuts of MCT in medical field are presented and discussed.

Influence of AMY1A copy number variations on obesity and other cardiometabolic risk factors: A review of the evidence

The rising incidence of obesity and type 2 diabetes is contributing to the escalating burden of disease globally. These metabolic disorders are closely linked with diet and in particular with carbohydrate consumption; hence, it is important to understand the underlying mechanisms that influence carbohydrate metabolism. Amylase, the enzyme responsible for the digestion of starch, is coded by the genes AMY1A, AMY1B, and AMY1C (salivary amylase) and AMY2A and AMY2B (pancreatic amylase). Previous studies demonstrate wide variations in AMY1A copy numbers, which can be attributed to several genetic, nutritional, and geographical diversities seen in populations globally. Current literature suggests that AMY1A copy number variations are important in obesity and other cardiometabolic disorders through their effects on glucose and lipid homeostasis, inflammatory markers, and the gut microbiome. This review synthesizes the available evidence to improve understanding of the role of AMY1A in obesity and related cardiometabolic risk factors and disorders including insulin resistance and type 2 diabetes, cardiovascular risk and inflammation, and the gut microbiome.

Population Structure, and Selection Signatures Underlying High-Altitude Adaptation Inferred From Genome-Wide Copy Number Variations in Chinese Indigenous Cattle

Copy number variations (CNVs) have been demonstrated as crucial substrates for evolution, adaptation and breed formation. Chinese indigenous cattle breeds exhibit a broad geographical distribution and diverse environmental adaptability. Here, we analyzed the population structure and adaptation to high altitude of Chinese indigenous cattle based on genome-wide CNVs derived from the high-density BovineHD SNP array. We successfully detected the genome-wide CNVs of 318 individuals from 24 Chinese indigenous cattle breeds and 37 yaks as outgroups. A total of 5,818 autosomal CNV regions (683 bp-4,477,860 bp in size), covering ~14.34% of the bovine genome (UMD3.1), were identified, showing abundant CNV resources. Neighbor-joining clustering, principal component analysis (PCA), and population admixture analysis based on these CNVs support that most Chinese cattle breeds are hybrids of Bos taurus taurus (hereinafter to be referred as Bos taurus) and Bos taurus indicus (Bos indicus). The distribution patterns of the CNVs could to some extent be related to the geographical backgrounds of the habitat of the breeds, and admixture among cattle breeds from different districts. We analyzed the selective signatures of CNVs positively involved in high-altitude adaptation using pairwise Fst analysis within breeds with a strong Bos taurus background (taurine-type breeds) and within Bos taurus×Bos indicus hybrids, respectively. CNV-overlapping genes with strong selection signatures (at top 0.5% of Fst value), including LETM1 (Fst = 0.490), TXNRD2 (Fst = 0.440), and STUB1 (Fst = 0.420) within taurine-type breeds, and NOXA1 (Fst = 0.233), RUVBL1 (Fst = 0.222), and SLC4A3 (Fst=0.154) within hybrids, were potentially involved in the adaptation to hypoxia. Thus, we provide a new profile of population structure from the CNV aspects of Chinese indigenous cattle and new insights into high-altitude adaptation in cattle.

No Evidence for Association of BMI with Salivary Amylase Gene Copy Number in the UK 1958 Birth Cohort

OBJECTIVE

In a 2014 publication, evidence was presented supporting the association of BMI with the copy number of the salivary amylase 1 (AMY1) gene, with an unprecedented effect size of -0.15 kg/m² (SE 0.02) per copy of AMY1. Most well-powered attempts to reproduce these findings have not been successful. However, because of different study designs, a significant association may still apply under restricted conditions such as in particular age groups. This study specifically tested the BMI-AMY1 association at different age points in the same individuals using longitudinal BMI information from participants in the UK 1958 Birth Cohort study.

METHODS

This study measured the AMY1 copy number by paralogue ratio tests in genomic DNA and by using array comparative genomic hybridization data. BMI data from 1958 Birth Cohort participants were available from eight different age points between 7 and 50 years.

RESULTS

No evidence, even at nominal significance, was found for association of the AMY1 copy number with BMI at any age point in approximately 1,400 members of the 1958 Birth Cohort or in 2,835 people from two disease cohorts from the Wellcome Trust Case Control Consortium.

CONCLUSIONS

The results do not support an association between BMI and AMY1 copy number at any age point.

Salivary α-amylase copy number is not associated with weight trajectories and glycemic improvements following clinical weight loss: results from a 2-phase dietary intervention study

BACKGROUND

Several studies recently reported contradicting results regarding the link between amylase 1 (AMY1) copy numbers (CNs), obesity, and type 2 diabetes.

OBJECTIVE

The aim of this study was to assess the impact of AMY1 CN on anthropometrics and glycemic outcomes in obese individuals following a 2-phase dietary weight loss intervention.

METHODS

Using the paralog ratio test, AMY1 CNs were accurately measured in 761 obese individuals from the DiOGenes study. Subjects first underwent an 8-wk low-calorie diet (LCD, at 800 kcal/d) and then were randomly assigned to a 6-mo weight maintenance dietary (WMD) intervention with arms having different glycemic loads.

RESULTS

At baseline, a modest association between AMY1 CN and BMI (P = 0.04) was observed. AMY1 CN was not associated with baseline glycemic variables. In addition, AMY1 CN was not associated with anthropometric or glycemic outcomes following either LCD or WMD. Interaction analyses between AMY1 CN and nutrient intake did not reveal any significant association with clinical parameters (at baseline and following LCD or WMD) or when testing gene × WMD interactions during the WMD phase.

CONCLUSION

In the absence of association with weight trajectories or glycemic improvements, the AMY1 CN cannot be considered as an important biomarker for response to a clinical weight loss and weight maintenance programs in overweight/obese subjects. This trial was registered at www.clinicaltrials.gov as NCT00390637.

A global analysis of CNVs in diverse yak populations using whole-genome resequencing

Background

Genomic structural variation represents a source for genetic and phenotypic variation, which may be subject to selection during the environmental adaptation and population differentiation. Here, we described a genome-wide analysis of copy number variations (CNVs) in 16 populations of yak based on genome resequencing data and CNV-based cluster analyses of these populations.

Results

In total, we identified 51,461 CNV events and defined 3174 copy number variation regions (CNVRs) that covered 163.8 Mb (6.2%) of yak genome with more “loss” events than both “gain” and “both” events, and we confirmed 31 CNVRs in 36 selected yaks using quantitative PCR. Of the total 163.8 Mb CNVR coverage, a 10.8 Mb region of high-confidence CNVRs directly overlapped with the 52.9 Mb of segmental duplications, and we confirmed their uneven distributions across chromosomes. Furthermore, functional annotation indicated that the CNVR-harbored genes have a considerable variety of molecular functions, including immune response, glucose metabolism, and sensory perception. Notably, some of the identified CNVR-harbored genes associated with adaptation to hypoxia (e.g., DCC, MRPS28, GSTCD, MOGAT2, DEXI, CIITA, and SMYD1). Additionally, cluster analysis, based on either individuals or populations, showed that the CNV clustering was divided into two origins, indicating that some yak CNVs are likely to arisen independently in different populations and contribute to population difference.

Conclusions

Collectively, the results of the present study advanced our understanding of CNV as an important type of genomic structural variation in yak, and provide a useful genomic resource to facilitate further research on yak evolution and breeding.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5451-5) contains supplementary material, which is available to authorized users.

The Malaria-Protective Human Glycophorin Structural Variant DUP4 Shows Somatic Mosaicism and Association with Hemoglobin Levels

Glycophorin A and glycophorin B are red blood cell surface proteins and are both receptors for the parasite Plasmodium falciparum, which is the principal cause of malaria in sub-Saharan Africa. DUP4 is a complex structural genomic variant that carries extra copies of a glycophorin A-glycophorin B fusion gene and has a dramatic effect on malaria risk by reducing the risk of severe malaria by up to 40%. Using fiber-FISH and Illumina sequencing, we validate the structural arrangement of the glycophorin locus in the DUP4 variant and reveal somatic variation in copy number of the glycophorin B-glycophorin A fusion gene. By developing a simple, specific, PCR-based assay for DUP4, we show that the DUP4 variant reaches a frequency of 13% in the population of a malaria-endemic village in south-eastern Tanzania. We genotype a substantial proportion of that village and demonstrate an association of DUP4 genotype with hemoglobin levels, a phenotype related to malaria, using a family-based association test. Taken together, we show that DUP4 is a complex structural variant that may be susceptible to somatic variation and show that DUP4 is associated with a malarial-related phenotype in a longitudinally followed population.

Human amylase gene copy number variation as a determinant of metabolic state

INTRODUCTION

Humans have multiple genes encoding amylase that are broadly divided into salivary (AMY1) and pancreatic (AMY2) genes. They exhibit some of the greatest copy numbers of any human gene, an expansion possibly driven by increased dietary starch intake. Within the population, amylase gene copy number is highly variable and there is evidence of an inverse association between AMY1 copy number and BMI.

AREAS COVERED

We examine the evidence for the link between AMY1 and BMI, its potential mechanisms, and the metabolic effects of salivary and pancreatic amylase, both in the gastrointestinal tract and the blood

EXPERT COMMENTARY

Salivary amylase may influence postprandial 'cephalic phase' insulin release, which improves glucose tolerance, while serum amylase may have insulin-sensitizing properties. This could explain the favorable metabolic status associated with higher AMY1 copy number. The association with BMI is harder to explain and is potentially mediated by increased flux of undigested starch into the ileum, with resultant effects on short-chain fatty acids (SCFAs), changes in gut microbiota and effects on appetite and energy expenditure in those with low copy number. Future research on the role of amylase as a determinant of metabolic health and BMI may lead to novel therapies to target obesity.