Human CCL3L1 copy number variation, gene expression, and the role of the CCL3L1-CCR5 axis in lung function

Exome-wide analysis of copy number variation shows association of the human leukocyte antigen region with asthma in UK Biobank

BACKGROUND

The role of copy number variants (CNVs) in susceptibility to asthma is not well understood. This is, in part, due to the difficulty of accurately measuring CNVs in large enough sample sizes to detect associations. The recent availability of whole-exome sequencing (WES) in large biobank studies provides an unprecedented opportunity to study the role of CNVs in asthma.

METHODS

We called common CNVs in 49,953 individuals in the first release of UK Biobank WES using ClinCNV software. CNVs were tested for association with asthma in a stage 1 analysis comprising 7098 asthma cases and 36,578 controls from the first release of sequencing data. Nominally-associated CNVs were then meta-analysed in stage 2 with an additional 17,280 asthma cases and 115,562 controls from the second release of UK Biobank exome sequencing, followed by validation and fine-mapping.

RESULTS

Five of 189 CNVs were associated with asthma in stage 2, including a deletion overlapping the HLA-DQA1 and HLA-DQB1 genes, a duplication of CHROMR/PRKRA, deletions within MUC22 and TAP2, and a duplication in FBRSL1. The HLA-DQA1, HLA-DQB1, MUC22 and TAP2 genes all reside within the human leukocyte antigen (HLA) region on chromosome 6. In silico analyses demonstrated that the deletion overlapping HLA-DQA1 and HLA-DQB1 is likely to be an artefact arising from under-mapping of reads from non-reference HLA haplotypes, and that the CHROMR/PRKRA and FBRSL1 duplications represent presence/absence of pseudogenes within the HLA region. Bayesian fine-mapping of the HLA region suggested that there are two independent asthma association signals. The variants with the largest posterior inclusion probability in the two credible sets were an amino acid change in HLA-DQB1 (glutamine to histidine at residue 253) and a multi-allelic amino acid change in HLA-DRB1 (presence/absence of serine, glycine or leucine at residue 11).

CONCLUSIONS

At least two independent loci characterised by amino acid changes in the HLA-DQA1, HLA-DQB1 and HLA-DRB1 genes are likely to account for association of SNPs and CNVs in this region with asthma. The high divergence of haplotypes in the HLA can give rise to spurious CNVs, providing an important, cautionary tale for future large-scale analyses of sequencing data.

CCL3L3-null status is associated with susceptibility to systemic lupus erythematosus

The correlation between copy number variation (CNV) and the susceptibility to systemic lupus erythematosus (SLE) has been reported for various immunity-related genes. However, the contribution of CNVs to SLE susceptibility awaits more investigation. To evaluate the copy numbers in immunity-related genes such as TNFAIP3, TNIP1, IL12B, TBX21 (T-bet), TLR7, C4A, C4B, CCL3L1, and CCL3L3, the modified real competitive polymerase chain reaction (mrcPCR) assay was employed, and the association between the copy numbers and SLE susceptibility was analyzed in 334 SLE patients and 338 controls. CCL3L3-null status was significantly associated with SLE susceptibility (OR > 18, P < 0.0001), which remained significant by Bonferroni's correction (corrected P = 0.0007). However, the significant association between C4B low-copy status and SLE susceptibility (OR = 1.6051, P = 0.0331) became non-significant by Bonferroni's correction (corrected P = 0.3938). Except for these results, no other significant association between SLE susceptibility and copy number status in other genes was observed. The CCL3L3-null status may be a significant factor for SLE susceptibility.

Identification of differentially expressed genes and the role of PDK4 in CD14+ monocytes of coronary artery disease

Background. Coronary artery disease (CAD) is a chronic inflammatory disease caused by development of atherosclerosis (AS), which is the leading cause of mortality and disability. Our study aimed to identify the differentially expressed genes (DEGs) in CD14⁺ monocytes from CAD patients compared with those from non-CAD controls, which might pave the way to diagnosis and treatment for CAD. Methods. The RNA-sequencing (RNA-seq) was performed by BGISEQ-500, followed by analyzing with R package to screening DEGs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed by R package. In addition, we validated the results of RNA-seq using real-time quantitative polymerase chain reaction (RT-qPCR). Furthermore, we explored the function of selected ten genes in LDL-treated CD14⁺ monocytes by RT-qPCR. Results. a total of 2897 DEGs were identified, including 753 up- and 2144 down-regulated genes in CD14⁺ monocytes from CAD patients. These DEGs were mainly enriched in plasma membrane and cell periphery of cell component, immune system process of biological process, NF-κB signaling pathway, cell adhesion molecules signaling pathway and cytokine–cytokine receptor interaction signaling pathway. In LDL-treated CD14⁺ monocytes, the mRNA expression of pyruvate dehydrogenase kinase 4 (PDK4) was significantly up-regulated. Conclusion. In the present study, we suggested that PDK4 might play a role in progression of CAD. The study will provide some pieces of evidence to investigate the role and mechanism of key genes in the pathogenesis of CAD.

Genotyping complex structural variation at the malaria-associated human glycophorin locus using a PCR-based strategy

Structural variation in the human genome can affect risk of disease. An example is a complex structural variant of the human glycophorin gene cluster, called DUP4, which is associated with a clinically significant level of protection against severe malaria. The human glycophorin gene cluster harbours at least 23 distinct structural variants, and accurate genotyping of this complex structural variation remains a challenge. Here, we use a polymerase chain reaction-based strategy to genotype structural variation at the human glycophorin gene cluster, including the alleles responsible for the U- blood group. We validate our approach, based on a triplex paralogue ratio test, on publically available samples from the 1000 Genomes project. We then genotype 574 individuals from a longitudinal birth cohort (Tori-Bossito cohort) using small amounts of DNA at low cost. Our approach readily identifies known deletions and duplications, and can potentially identify novel variants for further analysis. It will allow exploration of genetic variation at the glycophorin locus, and investigation of its relationship with malaria, in large sample sets at minimal cost, using standard molecular biology equipment.

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.