Polymorphisms near EXOC4 and LRGUK on chromosome 7q32 are associated with Type 2 Diabetes and fasting glucose; the NHLBI Family Heart Study

Whole genome sequencing of mouse lines divergently selected for fatness (FLI) and leanness (FHI) revealed several genetic variants as candidates for novel obesity genes

BACKGROUND

Analysing genomes of animal model organisms is widely used for understanding the genetic basis of complex traits and diseases, such as obesity, for which only a few mouse models exist, however, without their lean counterparts.

OBJECTIVE

To analyse genetic differences in the unique mouse models of polygenic obesity (Fat line) and leanness (Lean line) originating from the same base population and established by divergent selection over more than 60 generations.

METHODS

Genetic variability was analysed using WGS. Variants were identified with GATK and annotated with Ensembl VEP. g.Profiler, WebGestalt, and KEGG were used for GO and pathway enrichment analysis. miRNA seed regions were obtained with miRPathDB 2.0, LncRRIsearch was used to predict targets of identified lncRNAs, and genes influencing adipose tissue amount were searched using the IMPC database.

RESULTS

WGS analysis revealed 6.3 million SNPs, 1.3 million were new. Thousands of potentially impactful SNPs were identified, including within 24 genes related to adipose tissue amount. SNP density was highest in pseudogenes and regulatory RNAs. The Lean line carries SNP rs248726381 in the seed region of mmu-miR-3086-3p, which may affect fatty acid metabolism. KEGG analysis showed deleterious missense variants in immune response and diabetes genes, with food perception pathways being most enriched. Gene prioritisation considering SNP GERP scores, variant consequences, and allele comparison with other mouse lines identified seven novel obesity candidate genes: 4930441H08Rik, Aff3, Fam237b, Gm36633, Pced1a, Tecrl, and Zfp536.

CONCLUSION

WGS revealed many genetic differences between the lines that accumulated over the selection period, including variants with potential negative impacts on gene function. Given the increasing availability of mouse strains and genetic polymorphism catalogues, the study is a valuable resource for researchers to study obesity.

A Rare BRAF Fusion in Advanced Rectal Cancer Treated with Anti-Epidermal Growth Factor Receptor Therapy

Recently, v-raf murine sarcoma viral oncogene homologue B (BRAF) fusions have been identified in multiple cancer types using comprehensive genomic profiling (CGP) assays. BRAF fusions are extremely rare, occurring in <0.5% of patients with metastatic colorectal cancer (mCRC). Until now, there is no standard treatment for mCRC with BRAF fusions. Here, we report a recurrent colorectal cancer case that harbored an EXOC4-BRAF fusion. A 40-year-old female patient with a 2-year history of type 2 diabetes was diagnosed with pathologically confirmed stage IV rectal adenocarcinoma with liver metastasis. She underwent R0 resection after neoadjuvant therapy; however, her disease recurred at multiple metastatic sites (lymph nodes, ovary, and peritoneal gland). A rectal cancer surgical specimen was submitted for CGP (Foundation One) to identify potential targets to develop treatment strategies. An EXOC4-BRAF fusion was identified, and she achieved partial response to FOLFOX + panitumumab which is a fully human antibody directed against epidermal growth factor receptor. No EXOC4-BRAF fusions in colorectal cancer cases have been reported to date. Further studies investigating molecular mechanisms and novel targeted therapy approaches are required.

Genomic Characterization of the Istrian Shorthaired Hound

Istrian shorthaired hound is an old indigenous Croatian dog breed with historical traces of its origin, which date back to the 14th century. Due to its intelligence and great hunting abilities, it is considered an excellent hunting dog. Despite its ancient origin, there is no data on genetic diversity, population structure, and degree of inbreeding that could be used for advanced management and conservation of this breed. Our study aimed to provide a high-resolution population structure of the Istrian shorthaired hound using a 220K HD SNP array, to compare the obtained data with the genealogical records and to place the breed in a broader context of world dog populations. Relatively high population size and low inbreeding coefficient estimated from genealogical data indicate a preserved genetic diversity in this breed. The principle component analysis, the NeighborNet network, and TreeMix were used to determine the genetic relationship between the Istrian shorthaired hound and other breeds. The Istrian shorthaired hound was found to be genetically related to Italian hunting dogs sharing the same branch with the Segugio Italiano a Pelo Raso and Segugio Italiano a Pelo Forte. The ADMIXTURE analysis indicated that the Istrian shorthaired hound could be involved in the development of some other hunting dog breeds. The estimated effective population size (Ne) based on SNP data was similar to Ne calculated from genealogical data indicating the absence of bottlenecks and well-balanced use of breeding animals. The low genomic inbreeding coefficient, together with the higher number of short runs of homozygosity, observed in the Istrian shorthaired hound, confirms the ancient origin of the breed based on historical documents. The analysis of selective sweeps identified genomic regions with the strongest selection signals in the vicinity of the genes associated with cognitive performance and behavior. Genome analysis proved to be a useful tool for estimating population parameters and can be implemented in the conservation plan for this indigenous breed.

Novel genes involved in the genetic architecture of temperament in Brahman cattle

Cattle temperament is a complex and economically relevant trait. The objective of this study was to identify genomic regions and genes associated with cattle temperament. From a Brahman cattle population of 1,370 animals evaluated for temperament traits (Exit velocity-EV, Pen Score-PS, Temperament Score-TS), two groups of temperament-contrasting animals were identified based on their EV-average values ±1/2 standard deviation (SD). To be considered in the calm group, the EV of females ranged between 0.16–1.82 m/s (n = 50) and the EV of males ranged between 0.4–1.56 m/s (n = 48). Females were classified as temperamental if their EV ranged between 3.13–7.66 m/s (n = 46) and males were classified as temperamental if their EV ranged between 3.05–10.83 m/s (n = 45). Selected animals were genotyped using a total of 139,376 SNPs (GGP-HD-150K), evaluated for their association with EV. The Genome-Wide Association analysis (GWAS) identified fourteen SNPs: rs135340276, rs134895560, rs110190635, rs42949831, rs135982573, rs109393235, rs109531929, rs135087545, rs41839733, rs42486577, rs136661522, rs110882543, rs110864071, rs109722627, (P<8.1E-05), nine of them were located on intergenic regions, harboring seventeen genes, of which only ACER3, VRK2, FANCL and SLCO3A1 were considered candidate associated with bovine temperament due to their reported biological functions. Five SNPs were located at introns of the NRXN3, EXOC4, CACNG4 and SLC9A4 genes. The indicated candidate genes are implicated in a wide range of behavioural phenotypes and complex cognitive functions. The association of the fourteen SNPs on bovine temperament traits (EV, PS and TS) was evaluated; all these SNPs were significant for EV; only some were associated with PS and TS. Fourteen SNPs were associated with EV which allowed the identification of twenty-one candidate genes for Brahman temperament. From a functional point of view, the five intronic SNPs identified in this study, are candidates to address control of bovine temperament, further investigation will probe their role in expression of this trait.

Genome-Wide Tiling Array Analysis of HPV-Induced Warts Reveals Aberrant Methylation of Protein-Coding and Non-Coding Regions

The human papillomaviruses (HPV) are a group of double-stranded DNA viruses that exhibit an exclusive tropism for squamous epithelia. HPV can either be low- or high-risk depending on its ability to cause benign lesions or cancer, respectively. Unsurprisingly, the majority of epigenetic research has focused on the high-risk HPV types, neglecting the low-risk types in the process. Therefore, the main objective of this study is to better understand the epigenetics of wart formation by investigating the differences in methylation between HPV-induced cutaneous warts and normal skin. A number of clear and very significant differences in methylation patterns were found between cutaneous warts and normal skin. Around 55% of the top-ranking 100 differentially methylated genes in warts were protein coding, including the EXOC4, KCNU, RTN1, LGI1, IRF2, and NRG1 genes. Additionally, non-coding RNA genes, such as the AZIN1-AS1, LINC02008, and MGC27382 genes, constituted 11% of the top-ranking 100 differentially methylated genes. Warts exhibited a unique pattern of methylation that is a possible explanation for their transient nature. Since the genetics of cutaneous wart formation are not completely known, the findings of the present study could contribute to a better understanding of how HPV infection modulates host methylation to give rise to warts in the skin.

Genome-Wide Interaction and Pathway Association Studies for Body Mass Index

Objective: We investigated gene interactions (epistasis) for body mass index (BMI) in a European-American adult female cohort via genome-wide interaction analyses (GWIA) and pathway association analyses. Methods: Genome-wide pairwise interaction analyses were carried out for BMI in 493 extremely obese cases (BMI > 35 kg/m2) and 537 never-overweight controls (BMI < 25 kg/m2). To further validate the results, specific SNPs were selected based on the GWIA results for haplotype-based association studies. Pathway-based association analyses were performed using a modified Gene Set Enrichment Algorithm (GSEA) (GenGen program) to further explore BMI-related pathways using our genome wide association study (GWAS) data set, GIANT, ENGAGE, and DIAGRAM Consortia. Results: The EXOC4-1q23.1 interaction was associated with BMI, with the most significant epistasis between rs7800006 and rs10797020 (P = 2.63 × 10-11). In the pathway-based association analysis, Tob1 pathway showed the most significant association with BMI (empirical P < 0.001, FDR = 0.044, FWER = 0.040). These findings were further validated in different populations. Conclusion: Genome-wide pairwise SNP-SNP interaction and pathway analyses suggest that EXOC4 and TOB1-related pathways may contribute to the development of obesity.

Genome-wide analysis of insomnia in 1,331,010 individuals identifies new risk loci and functional pathways

Insomnia is the second most prevalent mental disorder, with no sufficient treatment available. Despite substantial heritability, insight into the associated genes and neurobiological pathways remains limited. Here, we use a large genetic association sample (n = 1,331,010) to detect novel loci and gain insight into the pathways, tissue and cell types involved in insomnia complaints. We identify 202 loci implicating 956 genes through positional, expression quantitative trait loci, and chromatin mapping. The meta-analysis explained 2.6% of the variance. We show gene set enrichments for the axonal part of neurons, cortical and subcortical tissues, and specific cell types, including striatal, hypothalamic, and claustrum neurons. We found considerable genetic correlations with psychiatric traits and sleep duration, and modest correlations with other sleep-related traits. Mendelian randomization identified the causal effects of insomnia on depression, diabetes, and cardiovascular disease, and the protective effects of educational attainment and intracranial volume. Our findings highlight key brain areas and cell types implicated in insomnia, and provide new treatment targets.

Identification of core genes and prediction of miRNAs associated with osteoporosis using a bioinformatics approach

Osteoporosis (OP) is an age-related disease, and osteoporotic fracture is one of the major causes of disability and mortality in elderly patients (>70 years old). As the pathogenesis and molecular mechanism of OP remain unclear, the identification of disease biomarkers is important for guiding research and providing therapeutic targets. In the present study, core genes and microRNAs (miRNAs) associated with OP were identified. Differentially expressed genes (DEGs) between human mesenchymal stem cell specimens from normal osseous tissues and OP tissues were detected using the GEO2R tool of the Gene Expression Omnibus database and Morpheus. Network topological parameters were determined using NetworkAnalyzer. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed using the Database for Annotation, Visualization and Integrated Discovery, and ClueGO. Cytoscape with the Search Tool for the Retrieval of Interacting Genes and Molecular Complex Detection plug-in was used to visualize protein-protein interactions (PPIs). Additionally, miRNA-gene regulatory modules were predicted using CyTargetLinker in order to guide future research. In total, 915 DEGs were identified, including 774 upregulated and 141 downregulated genes. Enriched GO terms and pathways were determined, including 'nervous system development', 'regulation of molecular function', 'glutamatergic synapse pathway' and 'pathways in cancer'. The node degrees of DEGs followed power-law distributions. A PPI network with 541 nodes and 1,431 edges was obtained. Overall, 3 important modules were identified from the PPI network. The following 10 genes were identified as core genes based on high degrees of connectivity: Albumin, PH domain leucine-rich repeat-containing protein phosphatase 2 (PHLPP2), DNA topoisomerase 2-α, kininogen 1 (KNG1), interleukin 2 (IL2), leucine-rich repeats and guanylate kinase domain containing, phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit γ (PIK3CG), leptin, transferrin and RNA polymerase II subunit A (POLR2A). Additionally, 15 miRNA-target interactions were obtained using CyTargetLinker. Overall, 7 miRNAs co-regulated IL2, 3 regulated PHLPP2, 3 regulated KNG1, 1 regulated PIK3CG and 1 modulated POLR2A. These results indicate potential biomarkers in the pathogenesis of OP and therapeutic targets.

SLC35B4, an Inhibitor of Gluconeogenesis, Responds to Glucose Stimulation and Downregulates Hsp60 among Other Proteins in HepG2 Liver Cell Lines

SLC35B4, solute receptor for UDP-N-acetylglucosamine and UDP-xylose, is associated with diabetes and predisposing conditions. This study investigated the localization of SLC35B4 and compared the differential expression between a knockdown of SLC35B4 and controls in HepG2. Responsiveness to glucose, expression, and localization were assayed using Western blot and immunostaining. Localization was confirmed using a proximity ligation assay. Two-dimensional (2D) gel electrophoresis and MALDI-TOF were used to identify differentially expressed proteins and pathway analysis was performed. SLC35B4 was increased by 60% upon glucose stimulation and localized in Golgi apparatus and endoplasmic reticulum. Presence of SLC35B4 in the Golgi apparatus suggests its involvement in the biosynthesis of glycoconjugate proteins. Four proteins were markedly under-expressed (Hsp60, HspA8, TUBA1A, and ENO1) and linked to the pathogenesis of diabetes or post-translationally modified by O-GlcNAc. Glucose levels activate SLC35B4 expression. This triggers a downstream effect via Hsp60 and other proteins. We hypothesize that the downstream effect on the proteins is mediated via altering the glycosylation pattern inside liver cells. The downstream cascade ultimately alters the ability of cultured liver cells to inhibit endogenous glucose production, and this could play a role in the association of the above-listed genes with the pathogenesis of diabetes.

Genomic structural variations for cardiovascular and metabolic comorbidity

The objective of this study was to identify genes targeted by both copy number and copy-neutral changes in the right coronary arteries in the area of advanced atherosclerotic plaques and intact internal mammary arteries derived from the same individuals with comorbid coronary artery disease and metabolic syndrome. The artery samples from 10 patients were screened for genomic imbalances using array comparative genomic hybridization. Ninety high-confidence, identical copy number variations (CNVs) were detected. We also identified eight copy-neutral changes (cn-LOHs) > 1.5 Mb in paired arterial samples in 4 of 10 individuals. The frequencies of the two gains located in the 10q24.31 (ERLIN1) and 12q24.11 (UNG, ACACB) genomic regions were evaluated in 33 paired arteries and blood samples. Two patients contained the gain in 10q24.31 (ERLIN1) and one patient contained the gain in 12q24.11 (UNG, ACACB) that affected only the blood DNA. An additional two patients harboured these CNVs in both the arteries and blood. In conclusion, we discovered and confirmed a gain of the 10q24.31 (ERLIN1) and 12q24.11 (UNG, ACACB) genomic regions in patients with coronary artery disease and metabolic comorbidity. Analysis of DNA extracted from blood indicated a possible somatic origin for these CNVs.

Integrative Analysis of Metabolomic, Proteomic and Genomic Data to Reveal Functional Pathways and Candidate Genes for Drip Loss in Pigs

The aim of this study was to integrate multi omics data to characterize underlying functional pathways and candidate genes for drip loss in pigs. The consideration of different omics levels allows elucidating the black box of phenotype expression. Metabolite and protein profiling was applied in Musculus longissimus dorsi samples of 97 Duroc × Pietrain pigs. In total, 126 and 35 annotated metabolites and proteins were quantified, respectively. In addition, all animals were genotyped with the porcine 60 k Illumina beadchip. An enrichment analysis resulted in 10 pathways, amongst others, sphingolipid metabolism and glycolysis/gluconeogenesis, with significant influence on drip loss. Drip loss and 22 metabolic components were analyzed as intermediate phenotypes within a genome-wide association study (GWAS). We detected significantly associated genetic markers and candidate genes for drip loss and for most of the metabolic components. On chromosome 18, a region with promising candidate genes was identified based on SNPs associated with drip loss, the protein "phosphoglycerate mutase 2" and the metabolite glycine. We hypothesize that association studies based on intermediate phenotypes are able to provide comprehensive insights in the genetic variation of genes directly involved in the metabolism of performance traits. In this way, the analyses contribute to identify reliable candidate genes.

Genome-wide testing of putative functional exonic variants in relationship with breast and prostate cancer risk in a multiethnic population

Rare variation in protein coding sequence is poorly captured by GWAS arrays and has been hypothesized to contribute to disease heritability. Using the Illumina HumanExome SNP array, we successfully genotyped 191,032 common and rare non-synonymous, splice site, or nonsense variants in a multiethnic sample of 2,984 breast cancer cases, 4,376 prostate cancer cases, and 7,545 controls. In breast cancer, the strongest associations included either SNPs in or gene burden scores for genes LDLRAD1, SLC19A1, FGFBP3, CASP5, MMAB, SLC16A6, and INS-IGF2. In prostate cancer, one of the most associated SNPs was in the gene GPRC6A (rs2274911, Pro91Ser, OR = 0.88, P = 1.3 × 10(-5)) near to a known risk locus for prostate cancer; other suggestive associations were noted in genes such as F13A1, ANXA4, MANSC1, and GP6. For both breast and prostate cancer, several of the most significant associations involving SNPs or gene burden scores (sum of minor alleles) were noted in genes previously reported to be associated with a cancer-related phenotype. However, only one of the associations (rs145889899 in LDLRAD1, p = 2.5 × 10(-7) only seen in African Americans) for overall breast or prostate cancer risk was statistically significant after correcting for multiple comparisons. In addition to breast and prostate cancer, other cancer-related traits were examined (body mass index, PSA level, and alcohol drinking) with a number of known and potentially novel associations described. In general, these findings do not support there being many protein coding variants of moderate to high risk for breast and prostate cancer with odds ratios over a range that is probably required for protein coding variation to play a truly outstanding role in risk heritability. Very large sample sizes will be required to better define the role of rare and less penetrant coding variation in prostate and breast cancer disease genetics.

[Artificial selection for cattle based on high-density SNP markers]

With the implementation of genetic improvement in recent years, artificial selection has greatly improved beef cattle production performance and its genetic basis has been dramatically changed. In this study, based on the Illumina BovineSNP50 (54K) and BovineHD (770K) BeadChip and the FST value, we analyzed the genetic differentiation of cattle and screened the imprints of selection in bovine genome. Finally, we found 47104 OUTLIER SNP loci and 3064 candidate genes, for example, CLIC5, TG, CACNA2D1, and FSHR etc. The biological processes and molecular functions of genes were analyzed through gene annotation.The results of this study established a genome-wide map of selection footprints in beef cattle genome and a clue for in-depth study of artificial selection and understanding of biological evolution.Our results indicate that artificial selection has played an important role in cattle breed genetic improvement.

Dietary saturated fat, gender and genetic variation at the TCF7L2 locus predict the development of metabolic syndrome

Transcription factor 7-like 2 (TCF7L2) is the strongest genetic determinant of type 2 diabetes (T2DM) and insulin-related phenotypes to date. Dietary fat is a key environmental factor which may interact with genotype to affect risk of metabolic syndrome (MetS) and T2DM. This study investigated the relationship between the TCF7L2 rs7903146 polymorphism, insulin sensitivity/resistance and MetS in the LIPGENE-SU.VI.MAX study of MetS cases and matched controls (n=1754) and determined potential interactions with dietary fat intake. Female minor T allele carriers of rs7903146 had increased MetS risk (odds ratio [OR] 1.66, confidence interval [CI] 1.02-2.70, P=.04) and displayed elevated insulin concentrations (P=.005), impaired insulin sensitivity (P=.011), increased abdominal obesity (P=.008) and body mass index (P=.001) and higher blood pressure (P<.05) compared to the CC homozygotes. Metabolic syndrome risk was also modulated by dietary saturated fat (SFA) intake (P=.035 for interaction). High dietary SFA intake (≥15.5% energy) exacerbated MetS risk (OR 2.35, 95% CI 1.29-4.27, P=.005) and was associated with further impaired insulin sensitivity in the T allele carriers relative to the CC homozygotes (P=.025) and particularly to the T allele carriers with the lowest SFA intake (P=.008). No significant genotype effect on MetS risk or insulin sensitivity was evident among low-SFA consumers. In conclusion, the TCF7L2 rs7903146 polymorphism influences MetS risk, which is augmented by both gender and dietary SFA intake, suggesting novel gene-diet-gender interactions.

Deep congenic analysis identifies many strong, context-dependent QTLs, one of which, Slc35b4, regulates obesity and glucose homeostasis

Although central to many studies of phenotypic variation and disease susceptibility, characterizing the genetic architecture of complex traits has been unexpectedly difficult. For example, most of the susceptibility genes that contribute to highly heritable conditions such as obesity and type 2 diabetes (T2D) remain to be identified despite intensive study. We took advantage of mouse models of diet-induced metabolic disease in chromosome substitution strains (CSSs) both to characterize the genetic architecture of diet-induced obesity and glucose homeostasis and to test the feasibility of gene discovery. Beginning with a survey of CSSs, followed with genetic and phenotypic analysis of congenic, subcongenic, and subsubcongenic strains, we identified a remarkable number of closely linked, phenotypically heterogeneous quantitative trait loci (QTLs) on mouse chromosome 6 that have unexpectedly large phenotypic effects. Although fine-mapping reduced the genomic intervals and gene content of these QTLs over 3000-fold, the average phenotypic effect on body weight was reduced less than threefold, highlighting the "fractal" nature of genetic architecture in mice. Despite this genetic complexity, we found evidence for 14 QTLs in only 32 recombination events in less than 3000 mice, and with an average of four genes located within the three body weight QTLs in the subsubcongenic strains. For Obrq2a1, genetic and functional studies collectively identified the solute receptor Slc35b4 as a regulator of obesity, insulin resistance, and gluconeogenesis. This work demonstrated the unique power of CSSs as a platform for studying complex genetic traits and identifying QTLs.

Distinct genetic alterations in colorectal cancer

Background

Colon cancer (CRC) development often includes chromosomal instability (CIN) leading to amplifications and deletions of large DNA segments. Epidemiological, clinical, and cytogenetic studies showed that there are considerable differences between CRC tumors from African Americans (AAs) and Caucasian patients. In this study, we determined genomic copy number aberrations in sporadic CRC tumors from AAs, in order to investigate possible explanations for the observed disparities.

Methodology/Principal Findings

We applied genome-wide array comparative genome hybridization (aCGH) using a 105k chip to identify copy number aberrations in samples from 15 AAs. In addition, we did a population comparative analysis with aCGH data in Caucasians as well as with a widely publicized list of colon cancer genes (CAN genes). There was an average of 20 aberrations per patient with more amplifications than deletions. Analysis of DNA copy number of frequently altered chromosomes revealed that deletions occurred primarily in chromosomes 4, 8 and 18. Chromosomal duplications occurred in more than 50% of cases on chromosomes 7, 8, 13, 20 and X. The CIN profile showed some differences when compared to Caucasian alterations.

Conclusions/Significance

Chromosome X amplification in male patients and chromosomes 4, 8 and 18 deletions were prominent aberrations in AAs. Some CAN genes were altered at high frequencies in AAs with EXOC4, EPHB6, GNAS, MLL3 and TBX22 as the most frequently deleted genes and HAPLN1, ADAM29, SMAD2 and SMAD4 as the most frequently amplified genes. The observed CIN may play a distinctive role in CRC in AAs.