Identify Candidate Single Nucleotide Polymorphisms Research Articles

Deciphering the population structure and genetic basis of growth traits from whole-genome resequencing of the leopard coral grouper ( Plectropomus leopardus).

The leopard coral grouper ( Plectropomus leopardus) is a species of significant economic importance. Although artificial cultivation of P. leopardus has thrived in recent decades, the advancement of selective breeding has been hindered by the lack of comprehensive population genomic data. In this study, we identified over 8.73 million single nucleotide polymorphisms (SNPs) through whole-genome resequencing of 326 individuals spanning six distinct groups. Furthermore, we categorized 226 individuals with high-coverage sequencing depth (≥14×) into eight clusters based on their genetic profiles and phylogenetic relationships. Notably, four of these clusters exhibited pronounced genetic differentiation compared with the other populations. To identify potentially advantageous loci for P. leopardus, we examined genomic regions exhibiting selective sweeps by analyzing the nucleotide diversity ( θπ) and fixation index ( F ST) in these four clusters. Using these high-coverage resequencing data, we successfully constructed the first haplotype reference panel specific to P. leopardus. This achievement holds promise for enabling high-quality, cost-effective imputation methods. Additionally, we combined low-coverage sequencing data with imputation techniques for a genome-wide association study, aiming to identify candidate SNP loci and genes associated with growth traits. A significant concentration of these genes was observed on chromosome 17, which is primarily involved in skeletal muscle and embryonic development and cell proliferation. Notably, our detailed investigation of growth-related SNPs across the eight clusters revealed that cluster 5 harbored the most promising candidate SNPs, showing potential for genetic selective breeding efforts. These findings provide a robust toolkit and valuable insights into the management of germplasm resources and genome-driven breeding initiatives targeting P. leopardus.

Identification of single nucleotide polymorphisms (SNPs) associated with chronic graft-versus-host disease in patients undergoing allogeneic hematopoietic cell transplantation

IntroductionChronic graft-versus-host disease (cGVHD) is a debilitating side effect of allogeneic hematopoietic cell transplantation (HCT), affecting the quality of life of patients. We used whole exome sequencing to identify candidate SNPs and complete a multi-marker gene-level analysis using a cohort of cGVHD( +) (N = 16) and cGVHD( −) (N = 66) HCT patients.MethodsSaliva samples were collected from HCT patients (N = 82) pre-conditioning in a multi-center study from March 2011 to May 2018. Exome sequencing was performed and FASTQ files were processed for sequence alignments. Significant SNPs were identified by logistic regression using PLINK2v3.7 and Fisher’s exact test. One cGVHD( −) patient sample was excluded from further analysis since no SNP was present in at least 10% of the sample population. The FUMA platform’s SNP2GENE was utilized to annotate SNPs and generate a MAGMA output. Chromatin state visualization of lead SNPs was completed using Epilogos tool. FUMA’s GENE2FUNC was used to obtain gene function and tissue expression from lead genomic loci.ResultsLogistic regression classified 986 SNPs associated with cGVHD( +). SNP2GENE returned three genomic risk loci, four lead SNPs, 48 candidate SNPs, seven candidate GWAS tagged SNPs, and four mapped genes. Fisher’s exact test identified significant homozygous genotypes of four lead SNPs (p < 0.05). GENE2FUNC analysis of multi-marker SNP sets identified one positional gene set including lead SNPs for KANK1 and KDM4C and two curated gene sets including lead SNPs for PTPRD, KDM4C, and/or KANK1.ConclusionsOur data suggest that SNPs in three genes located on chromosome 9 confer genetic susceptibility to cGVHD in HCT patients. These genes modulate STAT3 expression and phosphorylation in cancer pathogenesis. The findings may have implications in the modulation of pathways currently targeted by JAK inhibitors in cGVHD clinical trials.

Identification of genetic loci associated with renal dysfunction after lung transplantation using an ethnic-specific single-nucleotide polymorphism array

Renal dysfunction is a long-term complication associated with an increased mortality after lung transplantation (LT). We investigated the association of single-nucleotide polymorphisms (SNPs) with the development of renal dysfunction after LT using a Japanese-specific SNP array. First, eligible samples of 34 LT recipients were genotyped using the SNP array and divided into two groups, according to the presence of homozygous and heterozygous combinations of mutant alleles of the 126 renal-related SNPs. To identify candidate SNPs, the renal function tests were compared between the two groups for each SNP. Next, we investigated the association between the candidate SNPs and the time course of changes of the estimated glomerular filtration rate (eGFR) in the 99 recipients until 10 years after the LT. ΔeGFR was defined as the difference between the postoperative and preoperative eGFR values. Eight SNPs were identified as the candidate SNPs in the 34 recipients. Validation analysis of these 8 candidate SNPs in all the 99 recipients showed that three SNPs, namely, rs10277115, rs4690095, and rs792064, were associated with significant changes of the ΔeGFR. Pre-transplant identification of high-risk patients for the development of renal dysfunction after LT based on the presence of these SNPs might contribute to providing personalized medicine.

Genome-Wide Association Study Using Genotyping by Sequencing for Bacterial Leaf Blight Resistance Loci in Local Thai Indica Rice

Bacterial leaf blight (BLB) is a devastating disease caused by Xanthomonas oryzae pv. oryzae (Xoo), which poses a significant threat to global rice production. In this study, a genome-wide association study (GWAS) was conducted using the genotyping-by-sequencing (GBS) approach to identify candidate single nucleotide polymorphisms (SNPs) associated with BLB resistance genes. The study utilized 200 indica rice accessions inoculated with seven distinct Xoo isolates and filtered highly significant SNPs using a minor allele frequency (MAF) of &gt;5% and a call rate of 75%. Four statistical models were used to explore potential SNPs associated with BLB resistance, resulting in the identification of 32 significant SNPs on chromosomes 1–8 and 12 in the rice genome. Additionally, 179 genes were located within ±100 kb of the SNP region, of which 49 were selected as candidate genes based on their known functions in plant defense mechanisms. Several candidate genes were identified, including two genes in the same linkage disequilibrium (LD) decay as the well-known BLB resistance gene (Xa1). These findings represent a valuable resource for conducting further functional studies and developing novel breeding strategies to enhance the crop’s resistance to this disease.

MicroRNA sequence variation can impact interactions with target mRNA in cattle

BackgroundSmall non-coding microRNAs (miRNAs) are important modulators at post-transcriptional levels. Single-nucleotide polymorphisms (SNPs) located in miRNA genes can alter the secondary structure of pre-miRNA to either impair or promote the miRNA maturation processes. Furthermore, SNPs located in the miRNA seed regions can stabilize or disturb miRNA-target interactions, thereby, quantitatively influence the expression of target genes. Therefore, the main objective of this study was to detect SNPs in bovine miRNAs using the whole-genome re-sequencing datasets of 1632 cattle of five breeds from the 1000 bull genomes project. ResultsIn total, our study identified 1109, 334, and 130 SNPs in the miRNA precursor, mature, and seed regions, respectively. The heterozygosity values were generally less than 0.3, and the minor allele frequencies (MAFs) were mainly less than 0.1. Most SNPs were in Hardy-Weinberg equilibrium (HWE) (HWE-P > 0.05). Furthermore, we found that the majority of SNPs (MAF > 0.1 and HWE-P > 0.05) in the miRNA seed regions altered the repertoire of target genes, which in turn were enriched in different KEGG pathways or GO terms. Thus target prediction for bta-miR-2888 revealed loss of 309 target genes and gain of 691 target genes. The 691 gained target genes were significantly enriched in 60 KEGG pathways and 21 GO terms. ConclusionIn summary, our study identified candidate SNPs in miRNA precursor, mature, and seed regions that are likely to affect target RNA interactions, thereby potentially influencing cattle phenotypic traits.

Genetic polymorphism of pleiotrophin is associated with pain experience in Japanese adults: Case-control study.

Genetic factors play a role in individual differences in pain experience. Here, we performed a genome-wide association study (GWAS) to identify novel loci regulating pain processing. We conducted a 2-stage GWAS and the candidate single-nucleotide polymorphisms (SNPs) association study on pain experience using an exploratory cohort of patients with cancer pain. The confirmatory cohort comprised of participants from the general population with and without habitual use of analgesic medication. In the exploratory cohort, we evaluated pain intensity using a numerical rating scale, recorded daily opioid dosages, and calculated pain reduction rate. In the confirmatory cohort, pain experience was defined as habitual nonsteroidal anti-inflammatory drug usage. Using linear regression models, we identified candidate SNP in the exploratory samples, and tested the association between phenotype and experienced pain in the confirmatory samples. We found 1 novel SNP (rs11764598)—located on the gene encoding for pleiotrophin on chromosome 7—that passed the genome-wide suggestive significance at 20% false discovery rate (FDR) correction in the exploratory samples of patients with cancer pain (P = 1.31 × 10-7, FDR = 0.101). We confirmed its significant association with daily analgesic usage in the confirmatory cohort (P = .028), although the minor allele affected pain experience in an opposite manner. We identified a novel genetic variant associated with pain experience. Further studies are required to validate the role of pleiotrophin in pain processing.

3DFAACTS-SNP: using regulatory T cell-specific epigenomics data to uncover candidate mechanisms of type 1 diabetes (T1D) risk

BackgroundGenome-wide association studies (GWAS) have enabled the discovery of single nucleotide polymorphisms (SNPs) that are significantly associated with many autoimmune diseases including type 1 diabetes (T1D). However, many of the identified variants lie in non-coding regions, limiting the identification of mechanisms that contribute to autoimmune disease progression. To address this problem, we developed a variant filtering workflow called 3DFAACTS-SNP to link genetic variants to target genes in a cell-specific manner. Here, we use 3DFAACTS-SNP to identify candidate SNPs and target genes associated with the loss of immune tolerance in regulatory T cells (Treg) in T1D.ResultsUsing 3DFAACTS-SNP, we identified from a list of 1228 previously fine-mapped variants, 36 SNPs with plausible Treg-specific mechanisms of action. The integration of cell type-specific chromosome conformation capture data in 3DFAACTS-SNP identified 266 regulatory regions and 47 candidate target genes that interact with these variant-containing regions in Treg cells. We further demonstrated the utility of the workflow by applying it to three other SNP autoimmune datasets, identifying 16 Treg-centric candidate variants and 60 interacting genes. Finally, we demonstrate the broad utility of 3DFAACTS-SNP for functional annotation of all known common (> 10% allele frequency) variants from the Genome Aggregation Database (gnomAD). We identified 9376 candidate variants and 4968 candidate target genes, generating a list of potential sites for future T1D or other autoimmune disease research.ConclusionsWe demonstrate that it is possible to further prioritise variants that contribute to T1D based on regulatory function, and illustrate the power of using cell type-specific multi-omics datasets to determine disease mechanisms. Our workflow can be customised to any cell type for which the individual datasets for functional annotation have been generated, giving broad applicability and utility.

Genome-wide association study identified candidate SNPs and genes associated with hypoxia tolerance in large yellow croaker (Larimichthys crocea)

The large yellow croaker (Larimichthys crocea) is an economically important fish farmed along the southeast coasts of China. In recent years, however, hypoxia has caused substantial economic losses in the L. crocea aquaculture industry. At present, the numbers of genes and molecular markers available for the selection and breeding of hypoxia tolerance traits in L. crocea remain very limited, and thus there is an urgent need to conduct breeding studies on the tolerance of this fish to hypoxia. Genome-wide association studies (GWAS) have been widely used to study the genetic basis of economically important traits in aquatic animals, and in this study, we adopted a GWAS approach to genetically characterize hypoxia tolerance in the L. crocea. A total of 1982 fish were subjected to hypoxia experiments, among which 398 individuals characterized by the most sensitive and tolerant phenotypes were genotyped. Following filtering, we obtained a total of 54,224 high-quality single-nucleotide polymorphisms (SNPs), of which four were identified as significant based on mixed linear (MLM) and Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (Blink) model analyses. Using these four SNPs, we identified ten candidate genes, including one gene, egln2, directly associated with hypoxia and other involved in important regulatory pathways, such as oxidative stress, ion regulation, and energy metabolism. These findings provide new insights into the genetic basis of hypoxia tolerance in L. crocea and will facilitate future breeding for marker-assisted selection and the breeding of hypoxia-tolerant L. crocea for commercial production.

Genome-Wide Association Analyses of Osteochondrosis in Belgian Warmbloods Reveal Candidate Genes Associated With Chondrocyte Development

Osteochondrosis (OC) is an important skeletal disease causing profound welfare concerns in horses. Although numerous studies have explored the genetics underlying OC in various breeds, the Belgian Warmblood (BW) remains unstudied despite having a concerning prevalence of 32.0%. As a result, this study aimed to conduct genome-wide association (GWA) analyses to identify candidate variants associated with OC in BWs. To achieve this, blood samples and radiographs were collected from 407 Belgian Warmbloods registered to one of two BW studbooks (Belgisch Warmbloedpaard and Zangersheide), and genotyping was performed using the 670K Axiom Equine Genotyping Array. GWA analyses using a principle component approach were then performed on OC status (OCS; presence or absence of OC at any joint), hock OC status (HOC) and stifle OC status (SOC). These analyses yielded significantly associated (P < .01) SNPs on Equus caballus chromosome (ECA) 3, ECA 12, and ECA 18 for OCS; however, no single nucleotide polymorphisms (SNPs) reached significance for HOC or SOC. Subsequent analysis of candidate genes within 500 kilobases of the significant SNPs revealed functions broadly related to cell differentiation and chondrocyte development. While this study represents another step forward in uncovering variants and biological pathways associated with OC, additional studies are needed to validate the newly identified candidate SNPs for OC in BWs. Further studies of OC in BWs, as well as other breeds, are critical in our efforts to fully understand the disease's etiopathogenesis and ultimately provide breeding programs better equipped to improve horse health and well-being.

Competitive endogenous RNA network and pathway-based analysis of LncRNA single-nucleotide polymorphism in myasthenia gravis

Myasthenia gravis (MG) is a complex neurological autoimmune disease with a pathogenetic mechanism that has yet to be elucidated. Emerging evidence has revealed that genes, non-coding RNAs and genetic variants play significant roles in the pathogenesis of MG. However, the molecular mechanisms of single nucleotide polymorphisms (SNPs) located on lncRNAs could disturb lncRNA-mediated ceRNA regulatory functions still unclear in MG. In this study, we collated 276 experimentally confirmed MG risk genes and 192 MG risk miRNAs. We then constructed a lncRNA-mediated ceRNA network for MG based on multi-step computational strategies. Next, we systematically integrated risk pathways and identified candidate SNPs in lncRNAs for MG based on data acquired from public databases. In addition, we constructed a pathway-based lncRNA-SNP mediated network (LSPN) that contained 128 lncRNAs targeting 8 MG risk pathways. By analyzing network, we propose a latent mechanism for how the “lncRNA-SNP-mRNA-pathway” axis affects the pathogenesis of MG. Moreover, 25 lncRNAs and 51 SNPs on lncRNAs were extracted from the “lncRNA-SNP-mRNA-pathway” axis. Finally, functional analyses demonstrated lncRNA-SNPs mediated ceRNA regulation pairs associated with MG participated in the MAPK signaling pathway. In summary, we constructed MG-specific lncRNA-SNPs mediated ceRNA regulatory networks based on pathway in the present study, which was helpful to elucidate the roles of lncRNA-SNPs in the pathogenesis of MG and provide novel insights into mechanism of lncRNA-SNPs as potential genetic risk biomarkers of MG.

Terc Gene Cluster Variants Predict Liver Telomere Length in Mice.

Variants in a gene cluster upstream-adjacent to TERC on human chromosome 3, which includes genes APRM, LRRC31, LRRC34 and MYNN, have been associated with telomere length in several human populations. Currently, the mechanism by which variants in the TERC gene cluster influence telomere length in humans is unknown. Given the proximity between the TERC gene cluster and TERC (~0.05 Mb) in humans, it is speculated that cluster variants are in linkage disequilibrium with a TERC causal variant. In mice, the Terc gene/Terc gene cluster are also located on chromosome 3; however, the Terc gene cluster is located distantly downstream of Terc (~60 Mb). Here, we initially aim to investigate the interactions between genotype and nicotine exposure on absolute liver telomere length (aTL) in a panel of eight inbred mouse strains. Although we found no significant impact of nicotine on liver aTL, this first experiment identified candidate single nucleotide polymorphisms (SNPs) in the murine Terc gene cluster (within genes Lrrc31, Lrriq4 and Mynn) co-varying with aTL in our panel. In a second experiment, we tested the association of these Terc gene cluster variants with liver aTL in an independent panel of eight inbred mice selected based on candidate SNP genotype. This supported our initial finding that Terc gene cluster polymorphisms impact aTL in mice, consistent with data in human populations. This provides support for mice as a model for telomere dynamics, especially for studying mechanisms underlying the association between Terc cluster variants and telomere length. Finally, these data suggest that mechanisms independent of linkage disequilibrium between the Terc/TERC gene cluster and the Terc/TERC gene mediate the cluster’s regulation of telomere length.

GWAS Identifies a Region Containing the SALL1 Gene in Variation of Pigmentation Intensity Within the Chestnut Coat Color of Horses.

Chestnut coat color in horses is determined by a missense mutation within the MC1R gene. However, the intensity of the chestnut color can vary widely within individuals possessing this genotype. Here, we investigated this variation using standardized photographs of 96 horses. Each horse was ranked lightest to darkest within the cohort for phenotype by 3 blinded observers. A genome-wide association study utilizing the relative shade ranking as the phenotype and using 268 487 single-nucleotide polymorphisms (SNPs) genotyped using the Affymetrix Equine 670k array identified a single significantly associated region on chromosome 3 (P = 2.934 × 10-8). Analysis of whole-genome sequences for horses spanning the diverse range of chestnut color identified candidate SNPs within the coding sequence of the only gene in the region: SALL1. The function of SALL1 is largely unknown, though it is predicted to interact with the Hermansky-Pudlak Syndrome type 1 (HPS1) protein, which causes partial albinism in humans. However, with only one study suggesting a circumstantial influence of the SALL1 protein on pigmentation, additional work is needed to confirm this new coat color locus in larger populations and investigate the function of this protein for impacts on equine health.

Gene expression pattern analysis using dual-color RT-MLPA and integrative genome-wide association studies of eQTL for tuberculosis suscepitibility

BackgroundWhen infected with Mycobacterium tuberculosis, only a small proportion of the population will develop active TB, and the role of host genetic factors in different TB infection status was not fully understood.MethodsForty-three patients with active tuberculosis and 49 with latent tuberculosis were enrolled in the prospective cohort. Expressing levels of 27 candidate mRNAs, which were previously demonstrated to differentially expressed in latent and active TB, were measured by dual color reverse transcription multiplex ligation dependent probe amplification assay (dcRT-MLPA). Using expression levels of these mRNAs as quantitative traits, associations between expression abundance and genome-wild single nucleotide polymorphisms (SNPs) were calculated. Finally, identified candidate SNPs were further assessed for their associations with TB infection status in a validation cohort with 313 Chinese Han cases.ResultsWe identified 9 differentially expressed mRNAs including il7r, il4, il8, tnfrsf1b, pgm5, ccl19, il2ra, marco and fpr1 in the prospective cohort. Through expression quantitative trait loci mapping, we screened out 8 SNPs associated with these mRNAs. Then, CG genotype of the SNP rs62292160 was finally verified to be significantly associated with higher transcription levels of IL4 in LTBI than in TB patients.ConclusionWe reported that the SNP rs62292160 in Chinese Han population may link to higher expression of il4 in latent tuberculosis. Our findings provided a new genetic variation locus for further exploration of the mechanisms of TB and a possible target for TB genetic susceptibility studies, which might aid the clinical decision to precision treatment of TB.

Identification of Candidate Genes and Pathways Associated with Obesity-Related Traits in Canines via Gene-Set Enrichment and Pathway-Based GWAS Analysis.

Simple SummaryObesity is a serious health issue and is increasing at an alarming rate in several dog breeds, but there is limited information on the genetic mechanism underlying it. Moreover, there have been very few reports on genetic markers associated with canine obesity. These studies were limited to the use of a single breed in the association study. In this study, we have performed a GWAS and supplemented it with gene-set enrichment and pathway-based analyses to identify causative loci and genes associated with canine obesity in 18 different dog breeds. From the GWAS, the significant markers associated with obesity-related traits including body weight (CACNA1B, C22orf39, U6, MYH14, PTPN2, SEH1L) and blood sugar (PRSS55, GRIK2), were identified. Furthermore, the gene-set enrichment and pathway-based analysis (GESA) highlighted five enriched pathways (Wnt signaling pathway, adherens junction, pathways in cancer, axon guidance, and insulin secretion) and seven GO terms (fat cell differentiation, calcium ion binding, cytoplasm, nucleus, phospholipid transport, central nervous system development, and cell surface) which were found to be shared among all the traits.The present study aimed to identify causative loci and genes enriched in pathways associated with canine obesity using a genome-wide association study (GWAS). The GWAS was first performed to identify candidate single-nucleotide polymorphisms (SNPs) associated with obesity and obesity-related traits including body weight and blood sugar in 18 different breeds of 153 dogs. A total of 10 and 2 SNPs were found to be significantly (p < 3.74 × 10−7) associated with body weight and blood sugar, respectively. None of the SNPs were identified to be significantly associated with obesity trait. We subsequently followed up the GWAS analysis with gene-set enrichment and pathway analyses. A gene-set with 1057, 1409, and 1243 SNPs annotated to 449, 933 and 820 genes for obesity, body weight, and blood sugar, respectively was created by sub-setting the GWAS result at a threshold of p < 0.01 for the gene-set enrichment analysis. In total, 84 GO and 21 KEGG pathways for obesity, 114 GO and 44 KEGG pathways for blood sugar, 120 GO and 24 KEGG pathways for body weight were found to be enriched. Among the pathways and GO terms, we highlighted five enriched pathways (Wnt signaling pathway, adherens junction, pathways in cancer, axon guidance, and insulin secretion) and seven GO terms (fat cell differentiation, calcium ion binding, cytoplasm, nucleus, phospholipid transport, central nervous system development, and cell surface) that were found to be shared among all the traits. Our data provide insights into the genes and pathways associated with obesity and obesity-related traits.

Characterization of genetic and phenotypic heterogeneity of obstructive sleep apnea using electronic health records

BackgroundObstructive sleep apnea (OSA) is defined by frequent episodes of reduced or complete cessation of airflow during sleep and is linked to negative health outcomes. Understanding the genetic factors influencing expression of OSA may lead to new treatment strategies. Electronic health records (EHRs) can be leveraged to both validate previously reported OSA-associated genomic variation and detect novel relationships between these variants and comorbidities.MethodsWe identified candidate single nucleotide polymorphisms (SNPs) via systematic literature review of existing research. Using datasets available at Geisinger (n = 39,407) and Vanderbilt University Medical Center (n = 24,084), we evaluated associations between 40 previously implicated SNPs and OSA diagnosis, defined using clinical codes. We also evaluated associations between these SNPs and OSA severity measures obtained from sleep reports at Geisinger (n = 6571). Finally, we used a phenome-wide association study approach to help reveal pleiotropic genetic effects between OSA candidate SNPs and other clinical codes and laboratory values available in the EHR.ResultsMost previously reported OSA candidate SNPs showed minimal to no evidence for associations with OSA diagnosis or severity in the EHR-derived datasets. Three SNPs in LEPR, MMP-9, and GABBR1 validated for an association with OSA diagnosis in European Americans; the SNP in GABBR1 was associated following meta-analysis of results from both clinical populations. The GABBR1 and LEPR SNPs, and one additional SNP, were associated with OSA severity measures in European Americans from Geisinger. Three additional candidate OSA SNPs were not associated with OSA-related traits but instead with hyperlipidemia and autoimmune diseases of the thyroid.ConclusionsTo our knowledge, this is one of the largest candidate gene studies and one of the first phenome-wide association studies of OSA genomic variation. Results validate genetic associates with OSA in the LEPR, MMP-9 and GABBR1 genes, but suggest that the majority of previously identified genetic associations with OSA may be false positives. Phenome-wide analyses provide evidence of mediated pleiotropy. Future well-powered genome-wide association analyses of OSA risk and severity across populations with diverse ancestral backgrounds are needed. The comprehensive nature of the analyses represents a platform for informing future work focused on understanding how genetic data can be useful to informing treatment of OSA and related comorbidities.

Identification of single nucleotide pleomorphisms associated with periodontal disease in head and neck cancer irradiation patients by exome sequencing.

Periodontal disease (PD) is a common oral complication in patients with head and neck cancer (HNC) undergoing radiation therapy (RT). Our objective was to identify candidate single nucleotide polymorphisms (SNPs) associated with PD in radiation-treated patients with HNC. DNA was extracted from the saliva of patients with HNC (n=69) before RT. Clinical attachment loss (CAL) increment greater than 0.2 mm over 24 months after RT was used to define PD progression. After exome sequencing, SNPs associated with post-RT PD progression were identified by using logistic regression and homozygosity analyses. The web tools STRING, the Database for Annotation, Visualization and Integrated Discovery (DAVID), GeneCodis, and Ensembl Variant Effect Predictor were used for functional analysis. Of the 48 patients with HNC with post-RT PD progression, 24 had no tooth with 5 mm or greater pocket depth before RT, whereas of the 21 patients with HNC without progression, 11 had PD initially. A total of 330 SNPs (249 genes) with over-represented homozygous genotype (98.5% variant allele) were found to be associated with post-RT PD. Sixty of these corresponded to PD-related pathways, including previously identified genes. In patients with HNC with post-RT PD progression, SNPs were found in genes (n=10) in contrast to those without progression (n=7). The SNPs of collagen genes were identified, potentially defining susceptibility to PD in patients with HNC, and this could be further investigated to characterize PD drug targets.

Using RNA-seq to characterize responses to 4-hydroxyphenylpyruvate dioxygenase (HPPD) inhibitor herbicide resistance in waterhemp (Amaranthus tuberculatus)

BackgroundWaterhemp (Amaranthus tuberculatus (Moq.) J.D. Sauer) is a problem weed commonly found in the Midwestern United States that can cause crippling yield losses for both maize (Zea mays L.) and soybean (Glycine max L. Merr). In 2011, 4-hydroxyphenylpyruvate-dioxygenase (HPPD, EC 1.13.11.27) inhibitor herbicide resistance was first reported in two waterhemp populations. Since the discovery of HPPD-herbicide resistance, studies have identified the mechanism of resistance and described the inheritance of the herbicide resistance. However, no studies have examined genome-wide gene expression changes in response to herbicide treatment in herbicide resistant and susceptible waterhemp.ResultsWe conducted RNA-sequencing (RNA-seq) analyses of two waterhemp populations (HPPD-herbicide resistant and susceptible), from herbicide-treated and mock-treated leaf samples at three, six, twelve, and twenty-four hours after treatment (HAT). We performed a de novo transcriptome assembly using all sample sequences. Following assessments of our assembly, individual samples were mapped to the de novo transcriptome allowing us to identify transcripts specific to a genotype, herbicide treatment, or time point. Our results indicate that the response of HPPD-herbicide resistant and susceptible waterhemp genotypes to HPPD-inhibiting herbicide is rapid, established as soon as 3 hours after herbicide treatment. Further, there was little overlap in gene expression between resistant and susceptible genotypes, highlighting dynamic differences in response to herbicide treatment. In addition, we used stringent analytical methods to identify candidate single nucleotide polymorphisms (SNPs) that distinguish the resistant and susceptible genotypes.ConclusionsThe waterhemp transcriptome, herbicide-responsive genes, and SNPs generated in this study provide valuable tools for future studies by numerous plant science communities. This collection of resources is essential to study and understand herbicide effects on gene expression in resistant and susceptible weeds. Understanding how herbicides impact gene expression could allow us to develop novel approaches for future herbicide development. Additionally, an increased understanding of the prolific traits intrinsic in weed success could lead to crop improvement.

Genetic Characterization of Maximilian Sunflower for the Development of a Locally Adapted Perennial Grain Oilseed

Maximilian sunflower (Helianthus maximiliani Schrad.), a crop wild relative of sunflower (Helianthus annuus L.), has been identified as a species of interest for the development of a perennial oilseed crop. Knowledge of the diversity, the potential for crop development, and genomic resources of this crop wild relative is limited. To facilitate its use in breeding programs, a baseline characterization of locally adapted germplasm is required to develop informed breeding strategies. Individuals were collected from nine sites in southern Manitoba, Canada, and characterized for phenotypic and genotypic divergence to estimate traits of interest for the implementation of a breeding program in Maximilian sunflower. Genotype‐by‐sequencing was used to characterize population genetic parameters and identify candidate single nucleotide polymorphisms (SNPs) associated with phenotypic divergence and environmental differences between collection sites. Candidate SNPs associated with frost‐free period, temperature during the primary vegetative growth period, elevation, soil CaCO3 equivalent, days to anthesis and capitulum size were identified and may be useful for the improvement of H. maximiliani and crop species related to cultivated sunflower. Associations between temperature, population structure, and overall plant size were also identified, suggesting phenotypic divergence across a local temperature gradient. The sampled Maximilian sunflower populations exhibited a high degree of polymorphism, low levels of inbreeding, and a highly heterozygous genome at the local scale, traits that favor the establishment of locally adapted germplasm pools. There appears to be sufficient variation to make selections for agronomic traits in local germplasm of Maximilian sunflower to support its development as a perennial oilseed.

Genome-wide association pathway analysis to identify candidate single nucleotide polymorphisms and molecular pathways associated with TP53 expression status in HBV-related hepatocellular carcinoma.

BackgroundThe aim of this investigation was to identify candidate single nucleotide polymorphisms (SNPs) and molecular pathways associated with tumor protein p53 (TP53) expression status in hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC), clarify their potential mechanisms, and generate SNP-to-gene to pathway hypothesis.Materials and methodsIdentify candidate Causal SNPs and Pathways (ICSNPathway) was used to perform pathway analysis based on the results of our previous genome-wide association study of TP53 expression status in 387 HBV-related HCC patients.ResultsThrough the ICSNPathway analysis, we identified 18 candidate SNPs and 10 candidate pathways that are associated with TP53 expression status in HBV-related HCC. The strongest mechanism involved the modulation of major histocompatibility complex, class II, DP beta 1 (human leukocyte antigen [HLA]-DPB1-rs1042153), major histocompatibility complex, class II, DQ beta 1 (HLA-DQB1-rs1130399, HLA-DQB1-rs1049056, HLA-DQB1-rs1049059, and HLA-DQB1-rs1049060), and major histocompatibility complex, class II, DR beta 1 (HLA-DRB1-rs35445101). SNPs consequently affected regulatory roles in all the candidate pathways except hematopoietic cell lineage pathways. Association analysis using the GSE14520 data set, Gene Multiple Association Network Integration Algorithm, and Search Tool for the Retrieval of Interacting Genes/Proteins suggests that all genes of the candidate SNPs were associated with TP53. Survival analysis showed that the collagen type VI alpha 3 chain (COL6A3) rs111231885 and COL6A3-rs113155945 and COL6A3 block 4 CC haplotypes with TP53 negative status may have protective effects in HBV-related HCC patients after hepatectomy.ConclusionOur pathway analysis identified 18 candidate SNPs and 10 candidate pathways that were associated with TP53 expression status in HBV-related HCC. Among these candidate SNPs, the genetic variation of COL6A3 may be a potential prognostic biomarker of HBV-related HCC.

Novel Sex Chromosomes in 3 Cichlid Fishes from Lake Tanganyika.

African cichlids are well known for their adaptive radiations, but it is now apparent that they also harbor an extraordinary diversity of sex chromosome systems. In this study, we sequenced pools of males and females from species in 3 different genera of cichlids from Lake Tanganyika. We then searched for regions that were differentiated following the patterns expected for sex chromosomes. We report 2 novel sex chromosomes systems, an XY system on LG19 in Tropheus sp. "black" and a ZW system on LG7 in Hemibates stenosoma. We also identify a ZW system on LG5 in Cyprichromis leptosoma that may be convergent with a system previously described in Lake Malawi cichlids. Our data also identify candidate single nucleotide polymorphisms for the blue/yellow tail color polymorphism observed among male C. leptosoma.