Candidate Loci Research Articles

Genome-wide association analysis unveils candidate genes and loci associated with aplasia cutis congenita in pigs

BackgroundAplasia cutis congenita (ACC) is a rare genetic disorder characterized by the localized or widespread absence of skin in humans and animals. Individuals with ACC may experience developmental abnormalities in the skeletal and muscular systems, as well as potential complications. Localized and isolated cases of ACC can be treated through surgical and medical interventions, while extensive cases of ACC may result in neonatal mortality. The presence of ACC in pigs has implications for animal welfare. It contributes to an elevated mortality rate among piglets at birth, leading to substantial economic losses in the pig farming industry. In order to elucidate candidate genetic loci associated with ACC, we performed a Genome-Wide Association Study analysis on 216 Duroc pigs. The primary goal of this study was to identify candidate genes that associated with ACC.ResultsThis study identified nine significant SNPs associated with ACC. Further analysis revealed the presence of two quantitative trait loci, 483 kb (5:18,196,971–18,680,098) on SSC 5 and 159 kb (13:20,713,440–207294431 bp) on SSC13. By annotating candidate genes within a 1 Mb region surrounding the significant SNPs, a total of 11 candidate genes were identified on SSC5 and SSC13, including KRT71, KRT1, KRT4, ITGB7, CSAD, RARG, SP7, PFKL, TRPM2, SUMO3, and TSPEAR.ConclusionsThe results of this study further elucidate the potential mechanisms underlying and genetic architecture of ACC and identify reliable candidate genes. These results lay the foundation for treating and understanding ACC in humans.

Association of MMP3, MMP14, and MMP25 gene polymorphisms with cerebral stroke risk: a case-control study

BackgroundCerebral stroke (CS) is the leading cause of death in China, and a complex disease caused by both alterable risk factors and genetic factors. This study intended to investigate the association of MMP3, MMP14, and MMP25 single nucleotide polymorphisms (SNPs) with CS risk in a Chinese Han population.MethodsA total of 1,348 Han Chinese were recruited in this case-control study. Four candidate loci including rs520540 A/G and rs679620 T/C of MMP3, rs2236302 G/C of MMP14, and rs10431961 T/C of MMP25 were successfully screened. The correlation between the four SNPs and CS risk was assessed by logistic regression analysis. The results were analyzed by false-positive report probability (FPRP) for chance or significance. The interactions between four SNPs associated with CS risk were assessed by multifactor dimensionality reduction (MDR).Resultsrs520540 A/G and rs679620 C/T SNP in MMP3 were associated with risk of CS in allele, codominant, dominant and log-additive models. Ischemic stroke risk were significantly lower in carriers with rs520540-A allele and rs679620-T allele than those with G/G or C/C genotypes. However, rs520540-A allele and rs679620-T allele were associated with higher risk of hemorrhagic stroke. Stratified analysis showed that these two SNPs were associated with reduced risk of CS in aged < 55 years, non-smoking and non-drinking participants, and rs679620 SNP also reduced CS risk in male participants. The levels of uric acid, high-density lipoprotein cholesterol, and eosinophil were different among patients with different genotypes of rs520540 and rs679620. No statistically significant association was found between MMP14 rs2236302 G/C or MMP25 rs10431961 T/C with CS even after stratification by stroke subtypes, age, gender as well as smoking and drinking conditions in all the genetic models.ConclusionMMP3 rs520540 A/G and rs679620 C/T polymorphisms were associated with CS risk in the Chinese Han population, which provides useful information for the prevention and diagnosis of CS.

Fine-mapping of a candidate gene for web blotch resistance in Arachis hypogaea L.

Peanut (Arachis hypogaea L.) is a globally important oil crop. Web blotch is one of the most important foliar diseases affecting peanut, which results in serious yield losses worldwide. Breeding web blotch-resistant peanut varieties is the most effective and economically viable method for minimizing yield losses due to web blotch. In the current study, a bulked segregant analysis with next-generation sequencing was used to analyze an F2:3 segregating population and identify candidate loci related to web blotch resistance. Based on the fine-mapping of the candidate genomic interval using kompetitive allele-specific PCR (KASP) markers, we identified a novel web blotch resistance-related locus spanning approximately 169 kb on chromosome 16. This region included four annotated genes, of which only Arahy.35VVQ3 had a non-synonymous single nucleotide polymorphism in the coding region between the two parents. Two markers (Chr.16.12872635 and Chr.16.12966357) linked to this gene were shown to be co-segregated with the resistance of peanut web blotch by 72 randomly selected recombinant inbred lines (RIL), which could be used in marker-assisted breeding of resistant peanut varieties.

Genome-wide dissection and haplotype analysis identified candidate loci for nitrogen use efficiency under drought conditions in winter wheat.

Climate change causes extreme conditions like prolonged drought, which results in yield reductions due to its effects on nutrient balances such as nitrogen uptake and utilization by plants. Nitrogen (N) is a crucial nutrient element for plant growth and productivity. Understanding the mechanistic basis of nitrogen use efficiency (NUE) under drought conditions is essential to improve wheat (Triticum aestivum L.) yield. Here, we evaluated the genetic variation of NUE-related traits and photosynthesis response in a diversity panel of 200 wheat genotypes under drought and nitrogen stress conditions to uncover the inherent genetic variation and identify quantitative trait loci (QTLs) underlying these traits. The results revealed significant genetic variations among the genotypes in response to drought stress and nitrogen deprivation. Drought impacted plant performance more than N deprivation due to its effect on water and nutrient uptake. GWAS identified a total of 27 QTLs with a significant main effect on the drought-related traits, while 10 QTLs were strongly associated with the NUE traits. Haplotype analysis revealed two different haplotype blocks within the associated region on chromosomes 1B and 5A. The two haplotypes showed contrasting effects on N uptake and use efficiency traits. The in silico and transcript analyses implicated candidate gene coding for cold shock protein. This gene was the most highly expressed gene under several stress conditions, including drought stress. Upon validation, these QTLs on 1B and 5A could be used as a diagnostic marker for NUE and drought tolerance screening in wheat.

Genetic control and prospects of predictive breeding for European winter wheat’s Zeleny sedimentation values and Hagberg-Perten falling number

Key messageSedimentation values and falling number in the last decades have helped maintain high baking quality despite rigorous selection for grain yield in wheat. Allelic combinations of major loci sustained the bread-making quality while improving grain yield. Glu-D1, Pinb-D1, and non-gluten proteins are associated with sedimentation values and falling number in European wheat.Zeleny sedimentation values (ZSV) and Hagberg-Perten falling number (HFN) are among the most important parameters that help determine the baking quality classes of wheat and, thus, influence the monetary benefits for growers. We used a published data set of 372 European wheat varieties evaluated in replicated field trials in multiple environments. ZSV and HFN traits hold a wide and significant genotypic variation and high broad-sense heritability. The genetic correlations revealed positive and significant associations of ZSV and HFN with each other, grain protein content (GPC) and grain hardness; however, they were all significantly negatively correlated with grain yield. Besides, GPC appeared to be the major predictor for ZSV and HFN. Our genome-wide association analyses based on high-quality SSR, SNP, and candidate gene markers revealed a strong quantitative genetic nature of ZSV and HFN by explaining their total genotypic variance as 41.49% and 38.06%, respectively. The association of known Glutenin (Glu-1) and Puroindoline (Pin-1) with ZSV provided positive analytic proof of our studies. We report novel candidate loci associated with globulins and albumins—the non-gluten monomeric proteins in wheat. In addition, predictive breeding analyses for ZSV and HFN suggest using genomic selection in the early stages of breeding programs with an average prediction accuracy of 81 and 59%, respectively.

Helicobacter pylori infection associated DNA methylation in primary gastric cancer significantly correlates with specific molecular and clinicopathological features.

Helicobacter pyloriinduces DNA methylation in gastric mucosa, which links to gastric cancer (GC) risk. In contrast, CpG island methylator phenotype (CIMP) is defined as high levels of cancer-specific methylation and provides distinct molecular and clinicopathological features of GC. The association between those two types of methylation in GC remains unclear. We examined DNA methylation of well-validated H. pylori infection associated genes in GC and its adjacent mucosa and investigated its association with CIMP, various molecular subtypes and clinical features.We studied 50 candidate loci in 24 gastric samples to identify H. pylori infection associated genes. Identified loci were further examined in 624 gastric tissue from 217 primary GC, 217 adjacent mucosa, and 190 mucosae from cancer-free subjects.We identified fivegenes (IGF2, SLC16A2, SOX11, P2RX7, and MYOD1) as hypermethylated in H. pylori infected gastric mucosa. In non-neoplastic mucosa, methylation of H. pylori infection associated genes was higher in patients with GC than those without. In primary GC tissues, higher methylation of H. pylori infection associated genes correlated with CIMP-positive and its related features, such as MLH1 methylated cases. On the other hand, GC with lower methylation of these genes presented aggressive clinicopathological features including undifferentiated histopathology, advanced stage at diagnosis.H. pylori infection associated DNA methylation is correlated with CIMP, specific molecular and clinicopathological features in GC, supporting its utility as promising biomarker in this tumor type.

QTLs and Candidate Loci Associated with Drought Tolerance Traits of Kaybonnet x ZHE733 Recombinant Inbred Lines Rice Population.

Rice is the most important staple crop for the sustenance of the world's population, and drought is a major factor limiting rice production. Quantitative trait locus (QTL) analysis of drought-resistance-related traits was conducted on a recombinant inbred line (RIL) population derived from the self-fed progeny of a cross between the drought-resistant tropical japonica U.S. adapted cultivar Kaybonnet and the drought-sensitive indica cultivar ZHE733. K/Z RIL population of 198 lines was screened in the field at Fayetteville (AR) for three consecutive years under controlled drought stress (DS) and well-watered (WW) treatment during the reproductive stage. The effects of DS were quantified by measuring morphological traits, grain yield components, and root architectural traits. A QTL analysis using a set of 4133 single nucleotide polymorphism (SNP) markers and the QTL IciMapping identified 41 QTLs and 184 candidate genes for drought-related traits within the DR-QTL regions. RT-qPCR in parental lines was used to confirm the putative candidate genes. The comparison between the drought-resistant parent (Kaybonnet) and the drought-sensitive parent (ZHE733) under DS conditions revealed that the gene expression of 15 candidate DR genes with known annotations and two candidate DR genes with unknown annotations within the DR-QTL regions was up-regulated in the drought-resistant parent (Kaybonnet). The outcomes of this research provide essential information that can be utilized in developing drought-resistant rice cultivars that have higher productivity when DS conditions are prevalent.

Global transcriptome analysis reveals fungal disease responsive core gene regulatory landscape in tea

Fungal infections are the inevitable limiting factor for productivity of tea. Transcriptome reprogramming recruits multiple regulatory pathways during pathogen infection. A comprehensive meta-analysis was performed utilizing previously reported, well-replicated transcriptomic datasets from seven fungal diseases of tea. The study identified a cumulative set of 18,517 differentially expressed genes (DEGs) in tea, implicated in several functional clusters, including the MAPK signaling pathway, transcriptional regulation, and the biosynthesis of phenylpropanoids. Gene set enrichment analyses under each pathogen stress elucidated that DEGs were involved in ethylene metabolism, secondary metabolism, receptor kinase activity, and various reactive oxygen species detoxification enzyme activities. Expressional fold change of combined datasets highlighting 2258 meta-DEGs shared a common transcriptomic response upon fungal stress in tea. Pervasive duplication events caused biotic stress-responsive core DEGs to appear in multiple copies throughout the tea genome. The co-expression network of meta-DEGs in multiple modules demonstrated the coordination of appropriate pathways, most of which involved cell wall organization. The functional coordination was controlled by a number of hub genes and miRNAs, leading to pathogenic resistance or susceptibility. This first-of-its-kind meta-analysis of host–pathogen interaction generated consensus candidate loci as molecular signatures, which can be associated with future resistance breeding programs in tea.

The holocentric chromosome microevolution: From phylogeographic patterns to genomic associations with environmental gradients.

Geographic isolation and chromosome evolution are two of the major drivers of diversification in eukaryotes in general, and specifically, in plants. On one hand, range shifts induced by Pleistocene glacial oscillations deeply shaped the evolutionary trajectories of species in the Northern Hemisphere. On the other hand, karyotype variability within species or species complexes may have adaptive potential as different karyotypes may represent different recombination rates and linkage groups that may be associated with locally adapted genes or supergenes. Organisms with holocentric chromosomes are ideal to study the link between local adaptation and chromosome evolution, due to their high cytogenetic variability, especially when it seems to be related to environmental variation. Here, we integrate the study of the phylogeography, chromosomal evolution and ecological requirements of a plant species complex distributed in the Western Euro-Mediterranean region (Carex gr. laevigata, Cyperaceae). We aim to clarify the relative influence of these factors on population differentiation and ultimately on speciation. We obtained a well-resolved RADseq phylogeny that sheds light on the phylogeographic patterns of molecular and chromosome number variation, which are compatible with south-to-north postglacial migration. In addition, landscape genomics analyses identified candidate loci for local adaptation, and also strong significant associations between the karyotype and the environment. We conclude that karyotype distribution in C. gr. laevigata has been constrained by both range shift dynamics and local adaptation. Our study demonstrates that chromosome evolution may be responsible, at least partially, for microevolutionary patterns of population differentiation and adaptation in Carex.

Natural Selection Signatures in the Hondo and Ryukyu Japanese Subpopulations.

Natural selection signatures across Japanese subpopulations are under-explored. Here we conducted genome-wide selection scans with 622,926 single nucleotide polymorphisms for 20,366 Japanese individuals, who were recruited from the main-islands of Japanese Archipelago (Hondo) and the Ryukyu Archipelago (Ryukyu), representing two major Japanese subpopulations. The integrated haplotype score (iHS) analysis identified several signals in one or both subpopulations. We found a novel candidate locus at IKZF2, especially in Ryukyu. Significant signals were observed in the major histocompatibility complex region in both subpopulations. The lead variants differed and demonstrated substantial allele frequency differences between Hondo and Ryukyu. The lead variant in Hondo tags HLA-A*33:03-C*14:03-B*44:03-DRB1*13:02-DQB1*06:04-DPB1*04:01, a haplotype specific to Japanese and Korean. While in Ryukyu, the lead variant tags DRB1*15:01-DQB1*06:02, which had been recognized as a genetic risk factor for narcolepsy. In contrast, it is reported to confer protective effects against type 1 diabetes and human T lymphotropic virus type 1-associated myelopathy/tropical spastic paraparesis. The FastSMC analysis identified 8 loci potentially affected by selection within the past 20-150 generations, including 2 novel candidate loci. The analysis also showed differences in selection patterns of ALDH2 between Hondo and Ryukyu, a gene recognized to be specifically targeted by selection in East Asian. In summary, our study provided insights into the selection signatures within the Japanese and nominated potential sources of selection pressure.

Genome-wide association study to identify novel loci and genes for Fusarium root rot resistance in sweet potato using genotyping-by-sequencing.

Fusarium root rot, caused by Fusarium solani, is a major post-harvest disease in sweet potatoes (Ipomoea batatas (L.) Lam.). An effective strategy for controlling this disease is the development of resistant varieties. In this study, a genome-wide association study (GWAS) was conducted on 96 sweet potato genotypes to identify novel candidate loci and dissect the genetic basis of Fusarium root rot resistance. Genotyping was performed using genotyping-by-sequencing (GBS), and 44,255 SNPs were identified after filtering. The genotypes (n = 96) were evaluated through resistance tests in 2021 and 2022, separately and combined. The GWAS identified two significant SNP markers (LG3_22903756 and LG4_2449919) on chromosomes 3 and 4 associated with Fusarium root rot resistance, respectively. Lesion length showed significant differences between homozygous A and G alleles of LG3_22903756, which can potentially be used to develop molecular markers for selecting accessions resistant to Fusarium root rot. Expression analysis of 11 putative genes flanking the significant SNPs revealed the alteration in the expression of nine genes, indicating their possible involvement in Fusarium root rot resistance. The results of this study will aid in the marker-assisted selection and functional analysis of candidate genes for Fusarium root rot resistance in sweet potatoes.

RyR2 Binding of an Antiarrhythmic Cyclic Depsipeptide Mapped Using Confocal Fluorescence Lifetime Detection of FRET.

Hyperactivity of cardiac sarcoplasmic reticulum (SR) ryanodine receptor (RyR2) Ca2+-release channels contributes to heart failure and arrhythmias. Reducing the RyR2 activity, particularly during cardiac relaxation (diastole), is a desirable therapeutic goal. We previously reported that the unnatural enantiomer (ent) of an insect-RyR activator, verticilide, inhibits porcine and mouse RyR2 at diastolic (nanomolar) Ca2+ and has in vivo efficacy against atrial and ventricular arrhythmia. To determine the ent-verticilide structural mode of action on RyR2 and guide its further development via medicinal chemistry structure-activity relationship studies, here, we used fluorescence lifetime (FLT)-measurements of Förster resonance energy transfer (FRET) in HEK293 cells expressing human RyR2. For these studies, we used an RyR-specific FRET molecular-toolkit and computational methods for trilateration (i.e., using distances to locate a point of interest). Multiexponential analysis of FLT-FRET measurements between four donor-labeled FKBP12.6 variants and acceptor-labeled ent-verticilide yielded distance relationships placing the acceptor probe at two candidate loci within the RyR2 cryo-EM map. One locus is within the Ry12 domain (at the corner periphery of the RyR2 tetrameric complex). The other locus is sandwiched at the interface between helical domain 1 and the SPRY3 domain. These findings document RyR2-target engagement by ent-verticilide, reveal new insight into the mechanism of action of this new class of RyR2-targeting drug candidate, and can serve as input in future computational determinations of the ent-verticilide binding site on RyR2 that will inform structure-activity studies for lead optimization.

Testing a candidate meiotic drive locus identified by pool sequencing.

Meiotic drive biases the transmission of alleles in heterozygous individuals, such that Mendel's law of equal segregation is violated. Most examples of meiotic drive have been discovered over the past century based on causing sex ratio distortion or the biased transmission of easily scoreable genetic markers that were linked to drive alleles. More recently, several approaches have been developed that attempt to identify distortions of Mendelian segregation genome wide. Here, we test a candidate female meiotic drive locus in Drosophila melanogaster, identified previously as causing a ∼54:46 distortion ratio using sequencing of large pools of backcross progeny. We inserted fluorescent visible markers near the candidate locus and scored transmission in thousands of individual progeny. We observed a small but significant deviation from the Mendelian expectation; however, it was in the opposite direction to that predicted based on the original experiments. We discuss several possible causes of the discrepancy between the 2 approaches, noting that subtle viability effects are particularly challenging to disentangle from potential small-effect meiotic drive loci. We conclude that pool sequencing approaches remain a powerful method to identify candidate meiotic drive loci but that genotyping of individual progeny at early developmental stages may be required for robust confirmation.

A genomic association study revealing subphenotypes of childhood steroid-sensitive nephrotic syndrome in a larger genomic sequencing cohort

Dissecting the genetic components that contribute to the two main subphenotypes of steroid-sensitive nephrotic syndrome (SSNS) using genome-wide association studies (GWAS) strategy is important for understanding the disease. We conducted a multicenter cohort study (360 patients and 1835 controls) combined with a GWAS strategy to identify susceptibility variants associated with the following two subphenotypes of SSNS: steroid-sensitive nephrotic syndrome without relapse (SSNSWR, 181 patients) and steroid-dependent/frequent relapse nephrotic syndrome (SDNS/FRNS, 179 patients). The distribution of two single-nucleotide polymorphisms (SNPs) in ANKRD36 and ALPG was significant between SSNSWR and healthy controls, and that of two SNPs in GAD1 and HLA-DQA1 was significant between SDNS/FRNS and healthy controls. Interestingly, rs1047989 in HLA-DQA1 was a candidate locus for SDNS/FRNS but not for SSNSWR. No significant SNPs were observed between SSNSWR and SDNS/FRNS. Meanwhile, chromosome 2:171713702 in GAD1 was associated with a greater steroid dose (>0.75 mg/kg/d) upon relapse to first remission in patients with SDNS/FRNS (odds ratio = 3.14; 95% confidence interval, 0.97–9.87; P = 0.034). rs117014418 in APOL4 was significantly associated with a decrease in eGFR of greater than 20% compared with the baseline in SDNS/FRNS patients (P = 0.0001). Protein–protein intersection network construction suggested that HLA-DQA1 and HLA-DQB1 function together through GSDMA. Thus, SSNSWR belongs to non-HLA region-dependent nephropathy, and the HLA-DQA/DQB region is likely strongly associated with disease relapse, especially in SDNS/FRNS. The study provides a novel approach for the GWAS strategy of SSNS and contributes to our understanding of the pathological mechanisms of SSNSWR and SDNS/FRNS.

Genome-Wide Association Study of Starch Properties in Local Thai Rice.

Rice (Oryza sativa L.) is the main source of energy for humans and a staple food of high cultural significance for much of the world's population. Rice with highly resistant starch (RS) is beneficial for health and can reduce the risk of disease, especially type II diabetes. The identification of loci affecting starch properties will facilitate breeding of high-quality and health-supportive rice. A genome-wide association study (GWAS) of 230 rice cultivars was used to identify candidate loci affecting starch properties. The apparent amylose content (AAC) among rice cultivars ranged from 7.04 to 33.06%, and the AAC was positively correlated with RS (R2 = 0.94) and negatively correlated with rapidly available glucose (RAG) (R2 = -0.73). Three loci responsible for starch properties were detected on chromosomes 1, 6, and 11. On chromosome 6, the most significant SNP corresponded to LOC_Os06g04200 which encodes granule-bound starch synthase I (GBSSI) or starch synthase. Two novel loci associated with starch traits were LOC_Os01g65810 and LOC_Os11g01580, which encode an unknown protein and a sodium/calcium exchanger, respectively. The markers associated with GBSSI and LOC_Os11g01580 were tested in two independent sets of rice populations to confirm their effect on starch properties. The identification of genes associated with starch traits will further the understanding of the molecular mechanisms affecting starch in rice and may be useful in the selection of rice varieties with improved starch.

Evaluating restriction enzyme selection for reduced representation sequencing in conservation genomics.

Conservation genomic studies in non-model organisms generally rely on reduced representation sequencing techniques based on restriction enzymes to identify population structure as well as candidate loci for local adaptation. While the expectation is that the reduced representation of the genome is randomly distributed, the proportion of the genome sampled might depend on the GC content of the recognition site of the restriction enzyme used. Here, we evaluated the distribution and functional composition of loci obtained after a reduced representation approach using Genotyping-by-Sequencing (GBS). To do so, we compared experimental data from two endemic fish species (Symphodus ocellatus and Symphodus tinca, EcoT22I enzyme) and two ecosystem engineer sea urchins (Paracentrotus lividus and Arbacia lixula, ApeKI enzyme). In brief, we mapped the sequenced loci to the phylogenetically closest reference genome available (Labrus bergylta in the fish and Strongylocentrotus purpuratus in the sea urchin datasets), classified them as exonic, intronic and intergenic, and studied their function by using Gene Ontology (GO) terms. We also simulated the effect of using both enzymes in the two reference genomes. In both simulated and experimental data, we detected an enrichment towards exonic or intergenic regions depending on the restriction enzyme used and failed to detect differences between total loci and candidate loci for adaptation in the empirical dataset. Most of the functions assigned to the mapped loci were shared between the four species and involved a myriad of general functions. Our results highlight the importance of restriction enzyme selection and the need for high-quality annotated genomes in conservation genomic studies.

Identification of novel candidate loci and genes for seed vigor-related traits in upland cotton (Gossypium hirsutum L.) via GWAS.

Seed vigor (SV) is a crucial trait determining the quality of crop seeds. Currently, over 80% of China's cotton-planting area is in Xinjiang Province, where a fully mechanized planting model is adopted, accounting for more than 90% of the total fiber production. Therefore, identifying SV-related loci and genes is crucial for improving cotton yield in Xinjiang. In this study, three seed vigor-related traits, including germination potential, germination rate, and germination index, were investigated across three environments in a panel of 355 diverse accessions based on 2,261,854 high-quality single-nucleotide polymorphisms (SNPs). A total of 26 significant SNPs were detected and divided into six quantitative trait locus regions, including 121 predicted candidate genes. By combining gene expression, gene annotation, and haplotype analysis, two novel candidate genes (Ghir_A09G002730 and Ghir_D03G009280) within qGR-A09-1 and qGI/GP/GR-D03-3 were associated with vigor-related traits, and Ghir_A09G002730 was found to be involved in artificial selection during cotton breeding by population genetic analysis. Thus, understanding the genetic mechanisms underlying seed vigor-related traits in cotton could help increase the efficiency of direct seeding by molecular marker-assisted selection breeding.

The Genomic Basis of Adaptation to High Elevations in Africanized Honey Bees.

A range of different genetic architectures underpin local adaptation in nature. Honey bees (Apis mellifera) in the Eastern African Mountains harbor high frequencies of two chromosomal inversions that likely govern adaptation to this high-elevation habitat. In the Americas, honey bees are hybrids of European and African ancestries and adaptation to latitudinal variation in climate correlates with the proportion of these ancestries across the genome. It is unknown which, if either, of these forms of genetic variation governs adaptation in honey bees living at high elevations in the Americas. Here, we performed whole-genome sequencing of 29 honey bees from both high- and low-elevation populations in Colombia. Analysis of genetic ancestry indicated that both populations were predominantly of African ancestry, but the East African inversions were not detected. However, individuals in the higher elevation population had significantly higher proportions of European ancestry, likely reflecting local adaptation. Several genomic regions exhibited particularly high differentiation between highland and lowland bees, containing candidate loci for local adaptation. Genes that were highly differentiated between highland and lowland populations were enriched for functions related to reproduction and sperm competition. Furthermore, variation in levels of European ancestry across the genome was correlated between populations of honey bees in the highland population and populations at higher latitudes in South America. The results are consistent with the hypothesis that adaptation to both latitude and elevation in these hybrid honey bees are mediated by variation in ancestry at many loci across the genome.

A Case-Only Genome-Wide Interaction Study of Smoking and Bladder Cancer Risk: Results from the COBLAnCE Cohort.

Bladder cancer (BC) is the 6th most common cancer worldwide, with tobacco smoking considered as its main risk factor. Accumulating evidence has found associations between genetic variants and the risk of BC. Candidate gene-environment interaction studies have suggested interactions between cigarette smoking and NAT2/GSTM1 gene variants. Our objective was to perform a genome-wide association case-only study using the French national prospective COBLAnCE cohort (COhort to study BLAdder CancEr), focusing on smoking behavior. The COBLAnCE cohort comprises 1800 BC patients enrolled between 2012 and 2018. Peripheral blood samples collected at enrolment were genotyped using the Illumina Global Screening Array with a Multi-Disease drop-in panel. Genotyping data (9,719,614 single nucleotide polymorphisms (SNP)) of 1674, 1283, and 1342 patients were analyzed for smoking status, average tobacco consumption, and age at smoking initiation, respectively. A genome-wide association study (GWAS) was conducted adjusting for gender, age, and genetic principal components. The results suggest new candidate loci (4q22.1, 12p13.1, 16p13.3) interacting with smoking behavior for the risk of BC. Our results need to be validated in other case-control or cohort studies.

The MARS PETCARE BIOBANK protocol: establishing a longitudinal study of health and disease in dogs and cats

BackgroundThe veterinary care of cats and dogs is increasingly embracing innovations first applied to human health, including an increased emphasis on preventative care and precision medicine. Large scale human population biobanks have advanced research in these areas; however, few have been established in veterinary medicine. The MARS PETCARE BIOBANK™ (MPB) is a prospective study that aims to build a longitudinal bank of biological samples, with paired medical and lifestyle data, from 20,000 initially healthy cats and dogs (10,000 / species), recruited through veterinary hospitals over a ten-year period. Here, we describe the MPB protocol and discuss its potential as a platform to increase understanding of why and how diseases develop and how to advance personalised veterinary healthcare.MethodsAt regular intervals, extensive diet, health and lifestyle information, electronic medical records, clinicopathology and activity data are collected, genotypes, whole genome sequences and faecal metagenomes analysed, and blood, plasma, serum, and faecal samples stored for future research.DiscussionProposed areas for research include the early detection and progression of age-related disease, risk factors for common conditions, the influence of the microbiome on health and disease and, through genome wide association studies, the identification of candidate loci for disease associated genetic variants. Genomic data will be open access and research proposals for access to data and samples will be considered. Over the coming years, the MPB will provide the longitudinal data and systematically collected biological samples required to generate important insights into companion animal health, identifying biomarkers of disease, supporting earlier identification of risk, and enabling individually tailored interventions to manage disease.