Identification of three genetic variants as novel susceptibility loci for body mass index in a Japanese population

Identification of Heilongjiang crossbred beef cattle pedigrees and reveals functional genes related to economic traits based on whole-genome SNP data

Introduction: To enhance the beef cattle industry, Heilongjiang Province has developed a new Crossbred beef cattle variety through crossbreeding with exotic commercial breeds. This new variety exhibits relatively excellent meat quality, and efficient reproductive performance, catering to market demands. Method: This study employed whole genome resequencing technology to analyze the genetic pedigree and diversity of 19 Heilongjiang Crossbred beef cattle, alongside 59 published genomes from East Asian, Eurasian, and European taurine cattle as controls. In addition, genes related to production traits were also searched by identifying Runs of Homozygosity (ROH) islands and important fragments from ancestors. Results: A total of 14,427,729 biallelic SNPs were discovered, with the majority located in intergenic and intron regions and a small percentage in exon regions, impacting protein function. Population genetic analyses including Principal Component Analysis (PCA), Neighbor-Joining (NJ) tree, and ADMIXTURE identified Angus, Holstein, and Mishima as the main ancestors of Crossbred beef cattle. In genetic diversity analysis, nucleotide diversity, linkage disequilibrium, and inbreeding coefficient analysis reveal that the genetic diversity of Crossbred beef cattle is at a moderate level, and a higher inbreeding coefficient indicates the need for careful breeding management. In addition, some genes related to economic traits are identified through the identification of Runs of Homozygosity (ROH) islands and important fragments from ancestors. Conclusion: This comprehensive genomic characterization supports the targeted improvement of economically important traits in Crossbred beef cattle, facilitating advanced breeding strategies.

Single-cell RNA sequencing reveals the heterogeneity of epithelial cell and fibroblast cells from non- to metastatic lymph node OTSCC

Lymph node metastasis (LNM) is one of the common features of oral tongue squamous cell carcinoma (OTSCC). LNM is also taken as a sign of advanced OTSCC and poor survival rate. Recently, single-cell RNA sequencing has been applied in investigating the heterogeneity of tumor microenvironment and discovering the potential biomarkers for helping the diagnosis and prognosticating. Pathogenesis of LNM in OTSCC remains unknown. Specifically, cancer-associated fibroblasts (CAFs) and epithelial tumor cells could foster the progression of tumors. Thus, in this study, we aimed to comprehensively analyze the roles of subpopulations of CAFs and epithelial tumor cells in lymph node metastatic OTSCC using the integration of OTSCC single-cell RNA sequencing datasets. Four distinct subtypes of CAFs, namely vascular CAFs, myofibroblast CAFs, inflammatory CAFs, and growth arrest CAFs were successfully discovered in LNM tumor and confirmed the roles of GAS and PTN pathways in the progression of tumor metastasis. In addition, NKAIN2+ epithelial cells and FN1+ epithelial cells specifically exhibited an upregulation of PTN, NRG, MIF, and SPP1 signaling pathways in the metastatic OTSCC. In doing so, we put forth some potential biomarkers that could be utilized for the purpose of diagnosing and prognosticating OTSCC during its metastatic phase and tried to confirm by immunofluorescence assays.

Whole genome sequence analysis of apparent treatment resistant hypertension status in participants from the Trans-Omics for Precision Medicine program

Introduction: Apparent treatment-resistant hypertension (aTRH) is characterized by the use of four or more antihypertensive (AHT) classes to achieve blood pressure (BP) control. In the current study, we conducted single-variant and gene-based analyses of aTRH among individuals from 12 Trans-Omics for Precision Medicine cohorts with whole-genome sequencing data. Methods: Cases were defined as individuals treated for hypertension (HTN) taking three different AHT classes, with average systolic BP ≥ 140 or diastolic BP ≥ 90 mmHg, or four or more medications regardless of BP (n = 1,705). A normotensive control group was defined as individuals with BP < 140/90 mmHg (n = 22,079), not on AHT medication. A second control group comprised individuals who were treatment responsive on one AHT medication with BP < 140/ 90 mmHg (n = 5,424). Logistic regression with kinship adjustment using the Scalable and Accurate Implementation of Generalized mixed models (SAIGE) was performed, adjusting for age, sex, and genetic ancestry. We assessed variants using SKAT-O in rare-variant analyses. Single-variant and gene-based tests were conducted in a pooled multi-ethnicity stratum, as well as self-reported ethnic/racial strata (European and African American). Results: One variant in the known HTN locus, KCNK3, was a top finding in the multi-ethnic analysis (p = 8.23E-07) for the normotensive control group [rs12476527, odds ratio (95% confidence interval) = 0.80 (0.74-0.88)]. This variant was replicated in the Vanderbilt University Medical Center's DNA repository data. Aggregate gene-based signals included the genes AGTPBP, MYL4, PDCD4, BBS9, ERG, and IER3. Discussion: Additional work validating these loci in larger, more diverse populations, is warranted to determine whether these regions influence the pathobiology of aTRH.

Identify known and novel candidate genes associated with backfat thickness in Duroc pigs by large-scale genome-wide association analysis

Backfat thickness (BFT) is complex and economically important traits in the pig industry, since it reflects fat deposition and can be used to measure the carcass lean meat percentage in pigs. In this study, all 6,550 pigs were genotyped using the Geneseek Porcine 50K SNP Chip to identify SNPs related to BFT and to search for candidate genes through genome-wide association analysis in two Duroc populations. In total, 80 SNPs, including 39 significant and 41 suggestive SNPs, and 6 QTLs were identified significantly associated with the BFT. In addition, 9 candidate genes, including a proven major gene MC4R, 3 important candidate genes (RYR1, HMGA1, and NUDT3) which were previously described as related to BFT, and 5 novel candidate genes (SIRT2, NKAIN2, AMH, SORCS1, and SORCS3) were found based on their potential functional roles in BFT. The functions of candidate genes and gene set enrichment analysis indicate that most important pathways are related to energy homeostasis and adipogenesis. Finally, our data suggest that most of the candidate genes can be directly used for genetic improvement through molecular markers, except that the MC4R gene has an antagonistic effect on growth rate and carcass lean meat percentage in breeding. Our results will advance our understanding of the complex genetic architecture of BFT traits and laid the foundation for additional genetic studies to increase carcass lean meat percentage of pig through marker-assisted selection and/or genomic selection.

Comparative population genomic analysis uncovers novel genomic footprints and genes associated with small body size in Chinese pony

Background

Body size is considered as one of the most fundamental properties of an organism. Due to intensive breeding and artificial selection throughout the domestication history, horses exhibit striking variations for heights at withers and body sizes. Debao pony (DBP), a famous Chinese horse, is known for its small body size and lives in Guangxi mountains of southern China. In this study, we employed comparative population genomics to study the genetic basis underlying the small body size of DBP breed based on the whole genome sequencing data. To detect genomic signatures of positive selection, we applied three methods based on population comparison, fixation index (F_ST), cross population composite likelihood ratio (XP-CLR) and nucleotide diversity (θπ), and further analyzed the results to find genomic regions under selection for body size-related traits.

Results

A number of protein-coding genes in windows with the top 1% values of F_ST (367 genes), XP-CLR (681 genes), and log₂ (θπ ratio) (332 genes) were identified. The most significant signal of positive selection was mapped to the NELL1 gene, probably underlies the body size and development traits, and may also have been selected for short stature in the DBP population. In addition, some other loci on different chromosomes were identified to be potentially involved in the development of body size.

Conclusions

Results of our study identified some positively selected genes across the horse genome, which are possibly involved in body size traits. These novel candidate genes may be useful targets for clarifying our understanding of the molecular basis of body size and as such they should be of great interest for future research into the genetic architecture of relevant traits in horse breeding program.

Protein and fat intake interacts with the haplotype of PTPN11_rs11066325, RPH3A_rs886477, and OAS3_rs2072134 to modulate serum HDL concentrations in middle-aged people

BACKGROUND & AIMS

Low serum HDL cholesterol (HDL-C) concentration is a risk factor for cardiovascular diseases and it is influenced by genetic and environmental factors. We hypothesized that genetic variants that decrease serum HDL-C concentrations may interact with nutrient intakes in ways that increase or decrease the risk of cardiovascular disease.

METHODS

Candidate genetic variants that can lower serum HDL-C concentrations were explored by genome-wide association studies (GWAS), after adjusting for covariates, in the Ansan/Ansung cohort (n = 8842) from KoGES. The best genetic variants were selected and used to form a haplotype. According to the haplotype frequencies of SNPs, they were divided into major allele, heterozygote allele, and minor allele. The association of haplotype with serum HDL-C levels was determined using logistic regression after adjusting for confounding factors. Interaction of the haplotype with nutrient intake was also determined.

RESULTS

PTPN11_rs11066325, RPH3A_rs886477 and OAS3_rs2072134 were selected to modulate serum HDL-C levels from GWAS(P = 1.09E-09, 7.04E-10, and 1.27E-09, respectively). The adjusted odds ratios (ORs) for a decrease in serum HDL-C concentration in the minor-allele group of the haplotype were elevated by 1.534 fold, compared to the major-allele group of the haplotype. Furthermore, the adjusted ORs for serum LDL cholesterol and levels increased by 1.645 in the minor-alleles compared to the major-alleles of the haplotype without a significant change of serum cholesterol levels. Interestingly, the adjusted ORs for serum triglyceride were lower in the minor-alleles than in the major-alleles. The haplotype had a significant interaction with the intake of protein, fat, saturated fatty acids (SAF) and polyunsaturated fatty acids (PUFA; P < 0.05). In particular, the minor alleles of the haplotype decreased serum HDL-C levels compared to the major-alleles in the high intake of protein, fat, SFA, and PUFA, not in the low intake.

CONCLUSIONS

People carrying the minor-allele of haplotypes should avoid diets that are high in protein and fat, especially rich in SFA and PUFA.

Evolutionary history of disease-susceptibility loci identified in longitudinal exome-wide association studies

Background

Our longitudinal exome‐wide association studies previously detected various genetic determinants of complex disorders using ~26,000 single‐nucleotide polymorphisms (SNPs) that passed quality control and longitudinal medical examination data (mean follow‐up period, 5 years) in 4884–6022 Japanese subjects. We found that allele frequencies of several identified SNPs were remarkably different among four ethnic groups. Elucidating the evolutionary history of disease‐susceptibility loci may help us uncover the pathogenesis of the related complex disorders.

Methods

In the present study, we conducted evolutionary analyses such as extended haplotype homozygosity, focusing on genomic regions containing disease‐susceptibility loci and based on genotyping data of our previous studies and datasets from the 1000 Genomes Project.

Results

Our evolutionary analyses suggest that derived alleles of rs78338345 of GGA3, rs7656604 at 4q13.3, rs34902660 of SLC17A3, and six SNPs closely located at 12q24.1 associated with type 2 diabetes mellitus, obesity, dyslipidemia, and three complex disorders (hypertension, hyperuricemia, and dyslipidemia), respectively, rapidly expanded after the human dispersion from Africa (Out‐of‐Africa). Allele frequencies of GGA3 and six SNPs at 12q24.1 appeared to have remarkably changed in East Asians, whereas the derived alleles of rs34902660 of SLC17A3 and rs7656604 at 4q13.3 might have spread across Japanese and non‐Africans, respectively, although we cannot completely exclude the possibility that allele frequencies of disease‐associated loci may be affected by demographic events.

Conclusion

Our findings indicate that derived allele frequencies of nine disease‐associated SNPs (rs78338345 of GGA3, rs7656604 at 4q13.3, rs34902660 of SLC17A3, and six SNPs at 12q24.1) identified in the longitudinal exome‐wide association studies largely increased in non‐Africans after Out‐of‐Africa.

Single-nucleotide polymorphisms (rs342275, rs342293, rs7694379, rs11789898, and rs17824620) showed significant association with lobaplatin-induced thrombocytopenia

This study aimed to investigate single-nucleotide polymorphisms (SNPs) associated with lobaplatin-induced thrombocytopenia in patients with advanced lung cancer in China. Thirty-nine patients who received lobaplatin-based chemotherapy in the 307 Hospitals of Chinese People's Liberation Army from April 2017 to March 2018 were enrolled as study subjects. Peripheral blood DNA was extracted, and 79 candidate SNP positions were selected. A Sanger sequencing platform was employed to measure genotypes for locating the SNP positions associated with lobaplatin-induced thrombocytopenia. Of the 79 candidate genes, SNPs rs342275 and rs7694379 were significantly associated with lobaplatin-induced decrease in platelet (PLT) count (P < 0.05). SNPs rs342275, rs342293, rs11789898, and rs17824620 showed significant association with lobaplatin-induced lowest PLT counts (P < 0.05). SNPs rs342275, rs342293, rs11789898, rs17824620, and rs7694379 can be used as predictors of thrombocytopenia induced by lobaplatin-based chemotherapy in patients with advanced lung cancer in China.