Exploring the tumor genomic landscape of aggressive prostate cancer by whole-genome sequencing of tissue or liquid biopsies.

Treatment resistance remains a major issue in aggressive prostate cancer (PC), and novel genomic biomarkers may guide better treatment selection. Circulating tumor DNA (ctDNA) can provide minimally invasive information about tumor genomes, but the genomic landscape of aggressive PC based on whole-genome sequencing (WGS) of ctDNA remains incompletely characterized. Thus, we here performed WGS of tumor tissue (n = 31) or plasma ctDNA (n = 10) from a total of 41 aggressive PC patients, including 11 hormone-naïve, 15 hormone-sensitive, and 15 castration-resistant patients. Across all variant types, we found progressively more altered tumor genomic profiles in later stages of aggressive PC. The potential driver genes most frequently affected by single-nucleotide variants or insertions/deletions included the known PC-related genes TP53, CDK12, and PTEN and the novel genes COL13A1, KCNH3, and SENP3. Etiologically, aggressive PC was associated with age-related and DNA repair-related mutational signatures. Copy number variants most frequently affected 14q11.2 and 8p21.2, where no well-recognized PC-related genes are located, and also frequently affected regions near the known PC-related genes MYC, AR, TP53, PTEN, and BRCA1. Structural variants most frequently involved not only the known PC-related genes TMPRSS2 and ERG but also the less extensively studied gene in this context, PTPRD. Finally, clinically actionable variants were detected throughout all stages of aggressive PC and in both plasma and tissue samples, emphasizing the potential clinical applicability of WGS of minimally invasive plasma samples. Overall, our study highlights the feasibility of using liquid biopsies for comprehensive genomic characterization as an alternative to tissue biopsies in advanced/aggressive PC.

Multi-ancestry study of the genetics of problematic alcohol use in over 1 million individuals.

Problematic alcohol use (PAU), a trait that combines alcohol use disorder and alcohol-related problems assessed with a questionnaire, is a leading cause of death and morbidity worldwide. Here we conducted a large cross-ancestry meta-analysis of PAU in 1,079,947 individuals (European, N = 903,147; African, N = 122,571; Latin American, N = 38,962; East Asian, N = 13,551; and South Asian, N = 1,716 ancestries). We observed a high degree of cross-ancestral similarity in the genetic architecture of PAU and identified 110 independent risk variants in within- and cross-ancestry analyses. Cross-ancestry fine mapping improved the identification of likely causal variants. Prioritizing genes through gene expression and chromatin interaction in brain tissues identified multiple genes associated with PAU. We identified existing medications for potential pharmacological studies by a computational drug repurposing analysis. Cross-ancestry polygenic risk scores showed better performance of association in independent samples than single-ancestry polygenic risk scores. Genetic correlations between PAU and other traits were observed in multiple ancestries, with other substance use traits having the highest correlations. This study advances our knowledge of the genetic etiology of PAU, and these findings may bring possible clinical applicability of genetics insights-together with neuroscience, biology and data science-closer.

Depression pathophysiology, risk prediction of recurrence and comorbid psychiatric disorders using genome-wide analyses.

Depression is a common psychiatric disorder and a leading cause of disability worldwide. Here we conducted a genome-wide association study meta-analysis of six datasets, including >1.3 million individuals (371,184 with depression) and identified 243 risk loci. Overall, 64 loci were new, including genes encoding glutamate and GABA receptors, which are targets for antidepressant drugs. Intersection with functional genomics data prioritized likely causal genes and revealed new enrichment of prenatal GABAergic neurons, astrocytes and oligodendrocyte lineages. We found depression to be highly polygenic, with ~11,700 variants explaining 90% of the single-nucleotide polymorphism heritability, estimating that >95% of risk variants for other psychiatric disorders (anxiety, schizophrenia, bipolar disorder and attention deficit hyperactivity disorder) were influencing depression risk when both concordant and discordant variants were considered, and nearly all depression risk variants influenced educational attainment. Additionally, depression genetic risk was associated with impaired complex cognition domains. We dissected the genetic and clinical heterogeneity, revealing distinct polygenic architectures across subgroups of depression and demonstrating significantly increased absolute risks for recurrence and psychiatric comorbidity among cases of depression with the highest polygenic burden, with considerable sex differences. The risks were up to 5- and 32-fold higher than cases with the lowest polygenic burden and the background population, respectively. These results deepen the understanding of the biology underlying depression, its disease progression and inform precision medicine approaches to treatment.

Multi-ancestry study of the genetics of problematic alcohol use in >1 million individuals.

Problematic alcohol use (PAU) is a leading cause of death and disability worldwide. To improve our understanding of the genetics of PAU, we conducted a large cross-ancestry meta-analysis of PAU in 1,079,947 individuals. We observed a high degree of cross-ancestral similarity in the genetic architecture of PAU and identified 110 independent risk variants in within- and cross-ancestry analyses. Cross-ancestry fine-mapping improved the identification of likely causal variants. Prioritizing genes through gene expression and/or chromatin interaction in brain tissues identified multiple genes associated with PAU. We identified existing medications for potential pharmacological studies by drug repurposing analysis. Cross-ancestry polygenic risk scores (PRS) showed better performance in independent sample than single-ancestry PRS. Genetic correlations between PAU and other traits were observed in multiple ancestries, with other substance use traits having the highest correlations. The analysis of diverse ancestries contributed significantly to the findings, and fills an important gap in the literature.

Genome-wide analyses of ADHD identify 27 risk loci, refine the genetic architecture and implicate several cognitive domains.

Attention-deficit hyperactivity disorder (ADHD) is a prevalent neurodevelopmental disorder with a major genetic component. Here, we present a genome-wide association study meta-analysis of ADHD comprising 38,691 individuals with ADHD and 186,843 controls. We identified 27 genome-wide significant loci, highlighting 76 potential risk genes enriched among genes expressed particularly in early brain development. Overall, ADHD genetic risk was associated with several brain-specific neuronal subtypes and midbrain dopaminergic neurons. In exome-sequencing data from 17,896 individuals, we identified an increased load of rare protein-truncating variants in ADHD for a set of risk genes enriched with probable causal common variants, potentially implicating SORCS3 in ADHD by both common and rare variants. Bivariate Gaussian mixture modeling estimated that 84-98% of ADHD-influencing variants are shared with other psychiatric disorders. In addition, common-variant ADHD risk was associated with impaired complex cognition such as verbal reasoning and a range of executive functions, including attention.

Authoritative subspecies diagnosis tool for European honey bees based on ancestry informative SNPs.

BACKGROUND: With numerous endemic subspecies representing four of its five evolutionary lineages, Europe holds a large fraction of Apis mellifera genetic diversity. This diversity and the natural distribution range have been altered by anthropogenic factors. The conservation of this natural heritage relies on the availability of accurate tools for subspecies diagnosis. Based on pool-sequence data from 2145 worker bees representing 22 populations sampled across Europe, we employed two highly discriminative approaches (PCA and FST) to select the most informative SNPs for ancestry inference. RESULTS: Using a supervised machine learning (ML) approach and a set of 3896 genotyped individuals, we could show that the 4094 selected single nucleotide polymorphisms (SNPs) provide an accurate prediction of ancestry inference in European honey bees. The best ML model was Linear Support Vector Classifier (Linear SVC) which correctly assigned most individuals to one of the 14 subspecies or different genetic origins with a mean accuracy of 96.2% ± 0.8 SD. A total of 3.8% of test individuals were misclassified, most probably due to limited differentiation between the subspecies caused by close geographical proximity, or human interference of genetic integrity of reference subspecies, or a combination thereof. CONCLUSIONS: The diagnostic tool presented here will contribute to a sustainable conservation and support breeding activities in order to preserve the genetic heritage of European honey bees.

Molecular characteristics of porcine alpha-synuclein splicing variants.

Alpha-synuclein (α-syn) is a 140 amino acid, intrinsically disordered protein with a potential role in neurotransmitter vesicle release. The protein is natively unfolded under physiological conditions, and is expressed predominantly in neural tissue. α-syn is associated with neuropathological conditions in Parkinson's disease, where the protein misfolds into oligomers and fibrils resulting in aggregates in Lewy bodies. Here we report the molecular cloning of SNCA cDNA encoding porcine α-syn and transcript variants hereof. Six transcripts coding for porcine α-syn are presented in the report, of which three result from exon skipping, generating in-frame splicing of coding exons 3 and 5. The splicing pattern of these alternative spliced variants is conserved between human and pig. All the observed in-frame deletions yield significantly shorter α-syn proteins compared with the 140 amino acid full-length protein. Expression analysis performed by real-time quantitative RT-PCR revealed a differential expression of the six transcript splicing variants in different pig organs and tissues. Common for all splicing variants, a very high transcript expression was detected in brain tissues and in spinal cord and very low or no expression outside the central nervous system. The porcine α-syn protein demonstrated markedly different biophysical characteristics compared with its human counterpart. No fibrillation of porcine α-syn was observed with the pig wild-type α-syn and A30P α-syn, and both variants show significantly reduced ability to bind to lipid vesicles. Overexpression of mutated porcine α-syn might recapitulate the human PD pathogenesis and lead to the identification of genetic modifiers of the disease.

Molecular cloning and characterization of porcine Na⁺/K⁺-ATPase isoform α4.

Na+/K+-ATPase is responsible for maintaining electrochemical gradients of Na+ and K+, which is essential for a variety of cellular functions including neuronal activity. The α-subunit of the Na+/K+-ATPase is composed of four different polypeptides (α1-α4) encoded by different genes. Na,K-ATPase α4, encoded by the ATP1A4 gene, is expressed in testis and in male germ cells of humans, rats and mice. The α4 polypeptide has an important role in sperm motility, and is essential for male fertility. Here we present the RT-PCR cloning and characterization of the porcine ATP1A4 cDNA coding for Na⁺/K⁺-ATPase polypeptide α4. The Na⁺/K⁺-ATPase polypeptide α4, consisting of 1030 amino acids, displays a high homology with its human counterpart (86%). Phylogenetic analysis demonstrated that porcine Na⁺/K⁺-ATPase polypeptide α4 is closely related to other mammalian counterparts. In addition, the genomic structure of the porcine ATP1A4 gene was determined, and the intron-exon organization was found to be similar to that of the human ATP1A4 gene. The promoter sequence for the porcine ATP1A4 gene was also identified. Investigation of the genetic variation in the porcine ATP1A4 gene revealed a missense A/G SNP in exon 18. This A/G polymorphism results in a substitution of a methionine to a glycine residue (M888G). A very high overall DNA methylation rate of the ATP1A4 gene, 70-80%, was observed in both brain and liver. Expression analysis demonstrated that the porcine ATP1A4 gene is predominantly expressed in testis. The sequence of the porcine ATP1A4 cDNA encoding the Na⁺/K⁺-ATPase α4 protein has been submitted to GenBank under the accession number GenBank Accession No. MG587082.

A-to-I RNA editing of the IGFBP7 transcript increases during aging in porcine brain tissues.

The IGFBP7 gene encodes insulin-like growth factor protein 7. IGFPB7 is involved in diverse biological functions including cell growth regulation, senescence and apoptosis, and also acts as a tumor suppressor in multiple cancers. The IGFBP7 mRNA is subject to A-to-I RNA editing mediated by adenosine deaminases acting on RNA 1 and 2 (ADAR1 and ADAR2). In the current study we have examined molecular characteristics of the porcine IGFBP7 gene, and determined the mRNA editing in different tissues. The A-to-I RNA editing of human IGFBP7 in positions Arg78 and Lys95 was shown to be conserved in the porcine homologue. In addition, a novel editing site was discovered in position Lys97 in the porcine IGFBP7 transcript. A differential editing was demonstrated at the three positions in the IGFBP7 transcript with very high degrees of editing in frontal cortex, cerebellum and lung. Interestingly, the degree of editing increased during aging in porcine frontal cortex and cerebellum. The IGFBP7 gene was mapped to pig chromosome 8. The porcine IGFBP7 gene was found to be ubiquitously expressed in examined organs and tissues. The methylation status of the IGFBP gene was examined in brain and liver by bisulfate sequencing and a high degree of methylation was found in the two tissues, 52% and 54%, respectively.

Differential A-to-I RNA editing of the serotonin-2C receptor G-protein-coupled, HTR2C, in porcine brain tissues.

The HTR2C gene encodes the 5-Hydroxytryptamine (serotonin) receptor 2C, G-protein-coupled protein which functions as a serotonin receptor. The HTR2C mRNA is subject to A-to-I RNA editing mediated by adenosine deaminases acting on RNA 1 and 2 (ADAR1 and ADAR2). In the current study we examined the molecular characteristics of the porcine HTR2C gene and determined the mRNA editing of the HTR2C transcript in different tissues. The A-to-I RNA editing of HTR2C was shown to be conserved in the porcine homologue with five nucleotides edited in exon 5. A differential editing was demonstrated with a high editing frequency in the frontal cortex, parietal cortex, occipital cortex, hypothalamus, brain stem and spinal cord and significantly lower in the cerebellum. No editing was seen in the liver and kidney. The porcine HTR2C gene was found to be exclusively expressed in brain tissues. The HTR2C gene was mapped to pig chromosome X. The methylation status of the HTR2C gene was examined in brain and liver by bisulfate sequencing and a high degree of methylation was found in the two tissues, at 89 and 72%, respectively. Our data describe differences in RNA editing in various regions of the porcine brain. The differences might reflect functional differences. Similarities between pigs and humans in differential RNA editing support the use of the pig as a model organism for the study of neurological diseases.

Molecular characterization and analysis of the porcine NURR1 gene.

Orphan receptor NURR1 (also termed NR4A2) belongs to the nuclear receptor superfamily and functions as a regulatory factor of differentiation, migration, maturation and maintenance of mesencephalic dopaminergic neurons. NURR1 plays an important role in nigrostriatal dopamine neuron development and is therefore implicated in the pathogenesis of neurodegenerative diseases linked to the dopamine system of the midbrain. Here we report the isolation and characterization of porcine NURR1 cDNA. The NURR1 cDNA was RT-PCR cloned using NURR1-specific oligonucleotide primers derived from in silico sequences. The porcine NURR1 cDNA encodes a polypeptide of 598 amino acids, displaying a very high similarity with bovine, human and mouse (99%) NURR1 protein. Expression analysis revealed a differential NURR1 mRNA expression in various organs and tissues. NURR1 transcripts could be detected as early as at 60 days of embryo development in different brain tissues. A significant increase in NURR1 transcript in the cerebellum and a decrease in NURR1 transcript in the basal ganglia was observed during embryo development. The porcine NURR1 gene was mapped to chromosome 15. Two missense mutations were found in exon 3, the first coding exon of NURR1. Methylation analysis of the porcine NURR1 gene body revealed a high methylation degree in brain tissue, whereas methylation of the promoter was very low. A decrease in DNA methylation in a discrete region of the NURR1 promoter was observed in pig frontal cortex during pig embryo development. This observation correlated with an increase in NURR1 transcripts. Therefore, methylation might be a determinant of NURR1 expression at certain time points in embryo development.

Splicing variants of porcine synphilin-1.

Parkinson's disease (PD), idiopathic and familial, is characterized by degradation of dopaminergic neurons and the presence of Lewy bodies (LB) in the substantia nigra. LBs contain aggregated proteins of which α-synuclein is the major component. The protein synphilin-1 interacts and colocalizes with α-synuclein in LBs. The aim of this study was to isolate and characterize porcine synphilin-1 and isoforms hereof with the future perspective to use the pig as a model for Parkinson's disease. The porcine SNCAIP cDNA was cloned by reverse transcriptase PCR. The spatial expression of SNCAIP mRNA was investigated by RNAseq. The presented work reports the molecular cloning and characterization of the porcine (Sus scrofa) synphilin-1 cDNA (SNCAIP) and three splice variants hereof. The porcine SNCAIP cDNA codes for a protein (synphilin-1) of 919 amino acids which shows a high similarity to human (90%) and to mouse (84%) synphilin-1. Three shorter transcript variants of the synphilin-1 gene were identified, all lacking one or more exons. SNCAIP transcripts were detected in most examined organs and tissues and the highest expression was found in brain tissues and lung. Conserved splicing variants and a novel splice form of synhilin-1 were found in this study. All synphilin-1 isoforms encoded by the identified transcript variants lack functional domains important for protein degradation.

Cloning and characterization of the porcine DBC1 gene encoding deleted in bladder cancer.

Deleted in bladder cancer 1 (DBC1) is a tumour suppressor which is involved in the regulation of cell growth and programmed cell death. In this study we report the cloning and characterization of porcine DBC1 cDNA. RT-PCR cloning produced a cDNA with an open reading frame of 2,283 bp encoding a polypeptide of 761 amino acids with a predicted molecular mass of 88.6 kDa and estimated isoelectric point of 9.1. The encoded pig DBC1 protein shows a very high amino acid similarity to human (99 %) and to mouse (98 %) DBC1. The porcine DBC1 gene was mapped to chromosome 1. The nucleotide sequence of the promoter displayed a high degree of conservation of elements responsible for neuron-specific expression. The porcine DBC1 gene was found to be highly expressed in brain tissues. The methylation status of the porcine DBC1 gene was examined in brain and liver by bisulfite sequencing. Methylation percentages of 53-61 were observed for the gene body whereas significantly lower values (1-4 %) were found in exon 1 and the promoter sequence of DBC1. The sequences of the porcine DBC1 cDNA and the DBC1 promoter and exon 1 sequence have been submitted to DDBJ/EMBL/GenBank under the accession numbers KF733442 and KJ396193, respectively.

Porcine SLITRK1: Molecular cloning and characterization.

The membrane protein SLITRK1 functions as a developmentally regulated stimulator of neurite outgrowth and variants in this gene have been implicated in Tourette syndrome. In the current study we have cloned and characterized the porcine SLITRK1 gene. The genomic organization of SLITRK1 lacks introns, as does its human and mouse counterparts. RT-PCR cloning revealed two SLITRK1 transcripts: a full-length mRNA and a transcript variant that results in a truncated protein. The encoded SLITRK1 protein, consisting of 695 amino acids, displays a very high homology to human SLITRK1 (99%). The porcine SLITRK1 gene is expressed exclusively in brain tissues.

Pairwise comparisons of ten porcine tissues identify differential transcriptional regulation at the gene, isoform, promoter and transcription start site level.

The transcriptome is the absolute set of transcripts in a tissue or cell at the time of sampling. In this study RNA-Seq is employed to enable the differential analysis of the transcriptome profile for ten porcine tissues in order to evaluate differences between the tissues at the gene and isoform expression level, together with an analysis of variation in transcription start sites, promoter usage, and splicing. Totally, 223 million RNA fragments were sequenced leading to the identification of 59,930 transcribed gene locations and 290,936 transcript variants using Cufflinks with similarity to approximately 13,899 annotated human genes. Pairwise analysis of tissues for differential expression at the gene level showed that the smallest differences were between tissues originating from the porcine brain. Interestingly, the relative level of differential expression at the isoform level did generally not vary between tissue contrasts. Furthermore, analysis of differential promoter usage between tissues, revealed a proportionally higher variation between cerebellum (CBE) versus frontal cortex and cerebellum versus hypothalamus (HYP) than in the remaining comparisons. In addition, the comparison of differential transcription start sites showed that the number of these sites is generally increased in comparisons including hypothalamus in contrast to other pairwise assessments. A comprehensive analysis of one of the tissue contrasts, i.e. cerebellum versus heart for differential variation at the gene, isoform, and transcription start site (TSS), and promoter level showed that several of the genes differed at all four levels. Interestingly, these genes were mainly annotated to the "electron transport chain" and neuronal differentiation, emphasizing that "tissue important" genes are regulated at several levels. Furthermore, our analysis shows that the "across tissue approach" has a promising potential when screening for possible explanations for variations, such as those observed at the gene expression levels.