Long-read assembly of a Great Dane genome highlights the contribution of GC-rich sequence and mobile elements to canine genomes.

Technological advances have allowed improvements in genome reference sequence assemblies. Here, we combined long- and short-read sequence resources to assemble the genome of a female Great Dane dog. This assembly has improved continuity compared to the existing Boxer-derived (CanFam3.1) reference genome. Annotation of the Great Dane assembly identified 22,182 protein-coding gene models and 7,049 long noncoding RNAs, including 49 protein-coding genes not present in the CanFam3.1 reference. The Great Dane assembly spans the majority of sequence gaps in the CanFam3.1 reference and illustrates that 2,151 gaps overlap the transcription start site of a predicted protein-coding gene. Moreover, a subset of the resolved gaps, which have an 80.95% median GC content, localize to transcription start sites and recombination hotspots more often than expected by chance, suggesting the stable canine recombinational landscape has shaped genome architecture. Alignment of the Great Dane and CanFam3.1 assemblies identified 16,834 deletions and 15,621 insertions, as well as 2,665 deletions and 3,493 insertions located on secondary contigs. These structural variants are dominated by retrotransposon insertion/deletion polymorphisms and include 16,221 dimorphic canine short interspersed elements (SINECs) and 1,121 dimorphic long interspersed element-1 sequences (LINE-1_Cfs). Analysis of sequences flanking the 3' end of LINE-1_Cfs (i.e., LINE-1_Cf 3'-transductions) suggests multiple retrotransposition-competent LINE-1_Cfs segregate among dog populations. Consistent with this conclusion, we demonstrate that a canine LINE-1_Cf element with intact open reading frames can retrotranspose its own RNA and that of a SINEC_Cf consensus sequence in cultured human cells, implicating ongoing retrotransposon activity as a driver of canine genetic variation.

Generating Controlled, Dynamic Chemical Landscapes to Study Microbial Behavior.

We demonstrate a method for the generation of controlled, dynamic chemical pulses-where localized chemoattractant becomes suddenly available at the microscale-to create micro-environments for microbial chemotaxis experiments. To create chemical pulses, we developed a system to introduce amino acid sources near-instantaneously by photolysis of caged amino acids within a polydimethylsiloxane (PDMS) microfluidic chamber containing a bacterial suspension. We applied this method to the chemotactic bacterium, Vibrio ordalii, which can actively climb these dynamic chemical gradients while being tracked by video microscopy. Amino acids, rendered biologically inert ('caged') by chemical modification with a photoremovable protecting group, are uniformly present in the suspension but not available for consumption until their sudden release, which occurs at user-defined points in time and space by means of a near-UV-A focused LED beam. The number of molecules released in the pulse can be determined by a calibration relationship between exposure time and uncaging fraction, where the absorption spectrum after photolysis is characterized by using UV-Vis spectroscopy. A nanoporous polycarbonate (PCTE) membrane can be integrated into the microfluidic device to allow the continuous removal by flow of the uncaged compounds and the spent media. A strong, irreversible bond between the PCTE membrane and the PDMS microfluidic structure is achieved by coating the membrane with a solution of 3-aminopropyltriethoxysilane (APTES) followed by plasma activation of the surfaces to be bonded. A computer-controlled system can generate user-defined sequences of pulses at different locations and with different intensities, so as to create resource landscapes with prescribed spatial and temporal variability. In each chemical landscape, the dynamics of bacterial movement at the individual scale and their accumulation at the population level can be obtained, thereby allowing the quantification of chemotactic performance and its effects on bacterial aggregations in ecologically relevant environments.

Strain-level diversity drives alternative community types in millimetre-scale granular biofilms.

Microbial communities are often highly diverse in their composition, both at a coarse-grained taxonomic level, such as genus, and at a highly resolved level, such as strains, within species. This variability can be driven by either extrinsic factors such as temperature and or by intrinsic ones, for example demographic fluctuations or ecological interactions. The relative contributions of these factors and the taxonomic level at which they influence community composition remain poorly understood, in part because of the difficulty in identifying true community replicates assembled under the same environmental parameters. Here, we address this problem using an activated granular sludge reactor in which millimetre-scale biofilm granules represent true community replicates. Differences in composition are then expected to be driven primarily by biotic factors. Using 142 shotgun metagenomes of single biofilm granules we found that, at the commonly used genus-level resolution, community replicates varied much more in their composition than would be expected from neutral assembly processes. This variation did not translate into any clear partitioning into discrete community types, that is, distinct compositional states, such as enterotypes in the human gut. However, a strong partition into community types did emerge at the strain level for the dominant organism: genotypes of Candidatus Accumulibacter that coexisted in the metacommunity (the reactor) excluded each other within community replicates (granules). Individual granule communities maintained a significant lineage structure, whereby the strain phylogeny of Accumulibacter correlated with the overall composition of the community, indicating a high potential for co-diversification among species and communities. Our results suggest that due to the high functional redundancy and competition between close relatives, alternative community types are most probably observed at the level of recently differentiated genotypes but not at higher orders of genetic resolution.

Modeling Human Population Separation History Using Physically Phased Genomes.

Phased haplotype sequences are a key component in many population genetic analyses since variation in haplotypes reflects the action of recombination, selection, and changes in population size. In humans, haplotypes are typically estimated from unphased sequence or genotyping data using statistical models applied to large reference panels. To assess the importance of correct haplotype phase on population history inference, we performed fosmid pool sequencing and resolved phased haplotypes of five individuals from diverse African populations (including Yoruba, Esan, Gambia, Maasai, and Mende). We physically phased 98% of heterozygous SNPs into haplotype-resolved blocks, obtaining a block N50 of 1 Mbp. We combined these data with additional phased genomes from San, Mbuti, Gujarati, and Centre de'Etude du Polymorphism Humain European populations and analyzed population size and separation history using the pairwise sequentially Markovian coalescent and multiple sequentially Markovian coalescent models. We find that statistically phased haplotypes yield a more recent split-time estimation compared with experimentally phased haplotypes. To better interpret patterns of cross-population coalescence, we implemented an approximate Bayesian computation approach to estimate population split times and migration rates by fitting the distribution of coalescent times inferred between two haplotypes, one from each population, to a standard isolation-with-migration model. We inferred that the separation between hunter-gatherer populations and other populations happened ∼120-140 KYA, with gene flow continuing until 30-40 KYA; separation between west-African and out-of-African populations happened ∼70-80 KYA; while the separation between Maasai and out-of-African populations happened ∼50 KYA.

Microbial interactions lead to rapid micro-scale successions on model marine particles.

In the ocean, organic particles harbour diverse bacterial communities, which collectively digest and recycle essential nutrients. Traits like motility and exo-enzyme production allow individual taxa to colonize and exploit particle resources, but it remains unclear how community dynamics emerge from these individual traits. Here we track the taxon and trait dynamics of bacteria attached to model marine particles and demonstrate that particle-attached communities undergo rapid, reproducible successions driven by ecological interactions. Motile, particle-degrading taxa are selected for during early successional stages. However, this selective pressure is later relaxed when secondary consumers invade, which are unable to use the particle resource but, instead, rely on carbon from primary degraders. This creates a trophic chain that shifts community metabolism away from the particle substrate. These results suggest that primary successions may shape particle-attached bacterial communities in the ocean and that rapid community-wide metabolic shifts could limit rates of marine particle degradation.

Sequencing Y chromosomes resolves discrepancy in time to common ancestor of males versus females.

The Y chromosome and the mitochondrial genome have been used to estimate when the common patrilineal and matrilineal ancestors of humans lived. We sequenced the genomes of 69 males from nine populations, including two in which we find basal branches of the Y-chromosome tree. We identify ancient phylogenetic structure within African haplogroups and resolve a long-standing ambiguity deep within the tree. Applying equivalent methodologies to the Y chromosome and the mitochondrial genome, we estimate the time to the most recent common ancestor (T(MRCA)) of the Y chromosome to be 120 to 156 thousand years and the mitochondrial genome T(MRCA) to be 99 to 148 thousand years. Our findings suggest that, contrary to previous claims, male lineages do not coalesce significantly more recently than female lineages.

BDNF Val(66)Met genotype is associated with drug-seeking phenotypes in heroin-dependent individuals: a pilot study.

Brain-derived neurotrophic factor (BDNF) Val(66)Met genotype has been associated with neurobehavioral deficits. To examine its relevance for addiction, we examined BDNF genotype differences in drug-seeking behavior. Heroin-dependent volunteers (n = 128) completed an interview that assessed past-month naturalistic drug-seeking/use behaviors. In African Americans (n = 74), the Met allele was uncommon (carrier frequency 6.8%); thus, analyses focused on European Americans (n = 54), in whom the Met allele was common (carrier frequency 37.0%). In their natural setting, Met carriers (n = 20) reported more time- and cost-intensive heroin-seeking and more cigarette use than Val homozygotes (n = 34). BDNF Val(66)Met genotype predicted 18.4% of variance in 'weekly heroin investment' (purchasing time × amount × frequency). These data suggest that the BDNF Met allele may confer a 'preferred drug-invested' phenotype, resistant to moderating effects of higher drug prices and non-drug reinforcement. These preliminary hypothesis-generating findings require replication, but are consistent with pre-clinical data that demonstrate neurotrophic influence in drug reinforcement. Whether this genotype is relevant to other abused substances besides opioids or nicotine, or treatment response, remains to be determined.

PER3 polymorphism and insomnia severity in alcohol dependence.

STUDY OBJECTIVES: Insomnia is common, persistent, and associated with relapse in alcohol-dependent (AD) patients. Although the underlying mechanisms are mostly unstudied, AD patients have impaired circadian rhythms and sleep drive, which may be genetically influenced. A polymorphism in the PER3 gene (PER3(4/4), PER3(4/5), PER3(5/5)) has previously been associated with circadian preference and sleep homeostasis, and the PER3(4/4)genotype has been characterized by evening preference and decreased sleep drive. The purpose of this study was to examine the influence of this polymorphism on insomnia severity in AD patients. We hypothesized that the PER3 polymorphism would be an independent predictor of insomnia severity with greatest severity observed in those with the PER3(4/4)genotype. DESIGN: Cross-sectional association of patient characteristics, genotype, and insomnia severity. Significant (P < 0.05) bivariate correlates were further analyzed by hierarchical, forced entry multiple linear regression. SETTING: Alcohol treatment programs in Warsaw, Poland. PATIENTS: Diagnosed with alcohol dependence (n = 285), according to the Diagnostic and Statistical Manual of Mental Disorders, 4(th) edition. MEASUREMENTS AND RESULTS: Drinking frequency, mental and physical health status, childhood abuse, and PER3 genotype were independent predictors of insomnia severity, as measured by a 7-item subscale of the Sleep Disorders Questionnaire, explaining 28.9% of the variance. Addition of the genotype in the final step significantly increased the amount of variance explained by 1.1% (P = 0.027). Those with the PER3(4/4)genotype had the greatest severity of insomnia symptoms. CONCLUSIONS: PER3 genotype contributed unique variance in predicting insomnia severity in AD patients. These results are consistent with genetically influenced impairment in sleep regulation mechanisms in AD patients with insomnia.

The CC genotype in HTR2A T102C polymorphism is associated with behavioral impulsivity in alcohol-dependent patients.

High levels of impulsivity can increase the vulnerability for development of alcohol dependence. Moreover, impulsivity is considered to be a predictor of poor treatment outcomes. Few studies, however, have directly examined the genetics of impulsivity in alcohol-dependent patients. We analyzed the relationships between a well-recognized genetic marker of serotonin activity and levels of impulsivity as measured by both the Barratt Impulsiveness Scale (BIS-11) and the stop-signal task among 304 alcohol-dependent patients. The stop-signal task was used as an independent, objective method of estimating the level of behavioral impulsivity, and the BIS-11 as a self-report measure of global impulsivity. Blood was collected and analyzed for the T102C (rs6313) polymorphism in the serotonin type 2A receptor gene (HTR2A). Our results indicate a significant association between high levels of behavioral impulsivity and the C/C genotype of rs6313 in alcohol-dependent patients. The CC genotype has been previously found to be associated with a reduction in 5HT2A receptors in the central nervous system. These results support the hypothesis that genetic factors are important determinants of behavioral impulsivity in alcohol-dependent patients, and that the serotonin system plays an important role in establishing its level.

Family-based association analysis to finemap bipolar linkage peak on chromosome 8q24 using 2,500 genotyped SNPs and 15,000 imputed SNPs.

OBJECTIVE: Multiple linkage and association studies have suggested chromosome 8q24 as a promising candidate region for bipolar disorder (BP). We performed a detailed association analysis assessing the contribution of common genetic variation in this region to the risk of BP. METHODS: We analyzed 2,756 single nucleotide polymorphism (SNP) markers in the chromosome 8q24 region of 3,512 individuals from 737 families. In addition, we extended genotype imputation methods to family-based data and imputed 22,725 HapMap SNPs in the same region on 8q24. We applied a family-based method to test 15,552 high-quality genotyped or imputed SNPs for association with BP. RESULTS: Our association analysis identified the most significant marker (p=4.80 × 10(-5) ), near the gene encoding potassium voltage-gated channel KQT-like protein (KCNQ3). Other marginally significant markers were located near adenylate cyclase 8 (ADCY8) and ST3 beta-galactoside alpha-2,3-sialyltransferase 1 (ST3GAL1). CONCLUSIONS: We developed an approach to apply MACH imputation to family-based data, which can increase the power to detect association signals. Our association results showed suggestive evidence of association of BP with loci near KCNQ3, ADCY8, and ST3GAL1. Consistent with genes identified by genome-wide association studies for BP, our results suggest the involvement of ion channelopathy in BP pathogenesis. However, common variants are insufficient to explain linkage findings in 8q24; other genetic variation should be explored.

Increased activity of Diaphanous homolog 3 (DIAPH3)/diaphanous causes hearing defects in humans with auditory neuropathy and in Drosophila.

Auditory neuropathy is a rare form of deafness characterized by an absent or abnormal auditory brainstem response with preservation of outer hair cell function. We have identified Diaphanous homolog 3 (DIAPH3) as the gene responsible for autosomal dominant nonsyndromic auditory neuropathy (AUNA1), which we previously mapped to chromosome 13q21-q24. Genotyping of additional family members narrowed the interval to an 11-Mb, 3.28-cM gene-poor region containing only four genes, including DIAPH3. DNA sequencing of DIAPH3 revealed a c.-172G>A, g. 48G>A mutation in a highly conserved region of the 5' UTR. The c.-172G>A mutation occurs within a GC box sequence element and was not found in 379 controls. Using genome-wide expression arrays and quantitative RT-PCR, we demonstrate a 2- to 3-fold overexpression of DIAPH3 mRNA in lymphoblastoid cell lines from affected individuals. Likewise, a significant increase (approximately 1.5-fold) in DIAPH3 protein was found by quantitative immunoblotting of lysates from lymphoblastoid cell lines derived from affected individuals in comparison with controls. In addition, the c.-172G>A mutation is sufficient to drive overexpression of a luciferase reporter. Finally, the expression of a constitutively active form of diaphanous protein in the auditory organ of Drosophila melanogaster recapitulates the phenotype of impaired response to sound. To date, only two genes, the otoferlin gene OTOF and the pejvakin gene PJVK, are known to underlie nonsyndromic auditory neuropathy. Genetic testing for DIAPH3 may be useful for individuals with recessive as well as dominant inheritance of nonsyndromic auditory neuropathy.

SSRI response in depression may be influenced by SNPs in HTR1B and HTR1A.

OBJECTIVES: Desensitization of serotonin 1A (HTR1A) and 1B (HTR1B) autoreceptors has been proposed to be involved in the delayed onset of response to selective serotonin reuptake inhibitors (SSRIs). Variations in gene expression in these genes may thus affect SSRI response. METHODS: Here, we test this hypothesis in two samples from the Sequenced Treatment Alternatives to Relieve Depression (STAR*D), and show evidence for involvement of several genetic variants alone and in interaction. Initially, three functional single nucleotide polymorphisms (SNPs) in the HTR1B gene and in the HTR1A gene were analyzed in 153 depressed patients treated with citalopram. The 16-item Quick Inventory of Depressive Symptomatology Clinician scores were evaluated over time with respect to genetic variation. RESULTS: Individuals homozygous for the -1019 G allele (rs6295) in HTR1A showed the higher baseline 16-item Quick Inventory of Depressive Symptomatology Clinician scores (P=0.033), and by 12 weeks had a significantly lower response rate (P=0.005). HTR1B haplotypes were estimated according to the previously reported in-vitro expression levels. Individuals who were homozygous for the high-expression haplotype showed significantly slower response to citalopram (P=0.034). We then analyzed more SNPs in the extended overall STAR*D sample. Although we could not directly test the same functional SNPs, we found that homozygotes for the G allele at rs1364043 in HTR1A (P=0.045) and the C allele of rs6298 in HTR1B showed better response to citalopram over time (P=0.022). Test for interaction between rs6298 in HTR1B and rs1364043 in HTR1A was significant (overall P=0.032). CONCLUSION: Our data suggest that an enhanced capacity of HTR1B or HTR1A transcriptional activity may impair desensitization of the autoreceptors during SSRI treatment.

Association between Val66Met brain-derived neurotrophic factor (BDNF) gene polymorphism and post-treatment relapse in alcohol dependence.

BACKGROUND: The purpose of this study was to examine relationships between genetic markers of central serotonin (5-HT) and dopamine function, and risk for post-treatment relapse, in a sample of alcohol-dependent patients. METHODS: The study included 154 patients from addiction treatment programs in Poland, who met DSM-IV criteria for alcohol dependence. After assessing demographics, severity of alcohol use, suicidality, impulsivity, depression, hopelessness, and severity of alcohol use at baseline, patients were followed for approximately 1 year to evaluate treatment outcomes. Genetic polymorphisms in several genes (TPH2, SLC6A4, HTR1A, HTR2A, COMT, and BDNF) were tested as predictors of relapse (defined as any drinking during follow-up) while controlling for baseline measures. RESULTS: Of 154 eligible patients, 123 (80%) completed follow-up and 48% (n = 59) of these individuals relapsed. Patients with the Val allele in the Val66Met BDNF polymorphism and the Met allele in the Val158Met COMT polymorphism were more likely to relapse. Only the BDNF Val/Val genotype predicted post-treatment relapse [odds ratio (OR) = 2.62; p = 0.019], and time to relapse (OR = 2.57; p = 0.002), after adjusting for baseline measures and other significant genetic markers. When the analysis was restricted to patients with a family history of alcohol dependence (n = 73), the associations between the BDNF Val/Val genotype and relapse (OR = 5.76, p = 0.0045) and time to relapse (hazard ratio = 4.93, p = 0.001) were even stronger. CONCLUSIONS: The Val66Met BDNF gene polymorphism was associated with a higher risk and earlier occurrence of relapse among patients treated for alcohol dependence. The study suggests a relationship between genetic markers and treatment outcomes in alcohol dependence. Because a large number of statistical tests were conducted for this study and the literature on genetics and relapse is so novel, the results should be considered as hypothesis generating and need to be replicated in independent studies.

SNPs on chips: the hidden genetic code in expression arrays.

Gene expression microarray analysis in postmortem brains is one of the fastest growing fields of psychiatric research. Here we show that common polymorphisms (SNPs) present on probe sets can masquerade as significant "gene expression" differences. After first observing this artifact in the Catechol-O-methyl transferase (COMT) gene, we replicate the finding in two additional genes predicted to show this artifact. Many Affymetrix chips contain thousands of SNPs that are both common and in the central probe region affecting hybridization, and thus have the potential to confound expression analysis.

[Influence of impulsivity, suicidality and serotonin genes on treatment outcomes in alcohol dependence]. / Wplyw impulsywnosci, sklonnosci samobójczych oraz genów ukladu serotoninowego na wyniki leczenia uzaleinienia od alkoholu.

AIM: The aim of the study was to identify risk factors of relapse by investigating relationships among suicidality, impulsivity, genetic markers of serotonin activity, and relapse in alcohol-dependent patients. METHODS: 90 alcohol dependent patients were followed for 12 months after the baseline assessment, which entailed evaluation of suicidality and impulsivity as well as collection of DNA samples. Polymorphisms of genes involved in the synthesis and activity of the serotonin system were analyzed. After 12 months from the first visit, the patients were re-contacted and interviewed for relapse. RESULTS: Relapse rates were significantly higher among patients with the history of suicidal attempts recorded at the baseline assessment. The genetic analysis showed that patients with the G/G genotype in the 5HTR1A gene were more likely to relapse, whereas patients with the C/C genotype were more likely to abstain. Moreover, there was a strong trend for an association between the G/G genotype and a history of suicide attempts. CONCLUSIONS: High level of suicidality may predict relapse in alcoholic patients. Altered serotonergic function increases the risk of a suicide attempt and may contribute to higher risk of relapse in alcohol dependent patients.

Mutation incidence in folate metabolism genes and regulatory genes in Polish families with neural tube defects.

Neural tube defects (NTDs) are a common cause of disability or death of new-borns, but the aetiology and genetic background of this disease are still poorly understood. Therefore, it was decided to determine the conditions for the identification of several polymorphisms and to perform a preliminary study on Polish NTD patients and their parents. According to the results of this study, the genetic predisposition to NTD can be correlated with the 677TT genotype in the MTHFR gene, 677CT/1298AC haplotype (the MTHFR gene), 2756G allele in the MTR gene, 66AG variant and minisatellite sequence with 5 or 10 repeats in intron 6 of the MTRR gene. The 530GG and TIVS7-2/TIVS7-2 genotypes in the T gene could also be considered as a risk factor for NTD. The analysis also revealed no correlation between neurulation disturbances and A4956G and A1186G mutations in the BRCA1 gene and the 844ins68bp in CBS gene. Although a correlation was found of some molecular markers with NTD, an additional examination should be conducted on more numerous groups to obtain statistically significant results.

[Mutations of MTHFR, MTR, MTRR genes as high risk factors for neural tube defects]. / Analiza mutacji w genach MTHFR, MTR, MTRR w rodzinach ryzyka wad cewy nerwowej.

Neural tube defects (NTDs) have a polygenic background. There are numerous genes known to be high-risk genetic factors for NTDs. Ones of them are mutations of foliate metabolisms pathways genes. This paper shows the results of analysis of common mutations of MTHFR, MTR and MTRR genes. Results of screening mutations 2756A-->G and 66A-->G in MTR and MTRR genes respectively show that are might have an effect on NTDs incidence among the examined population. Analysis of data for the studied population does not prove the influence of mutations 677C-->T and 1298A-->C of MTHFR gene on NTDs.

[Genetic risk factors of neural tube defects]. / Genetyczne czynniki ryzyka wad cewy nerwowej.

Neural tube defects (NTDs) are a group of diseases caused by a failure of closure of the neural tube. Its aetiology contains both environmental and genetic factors. NTDs have a polygenic background. Genes, which are linked with NTDs occurrence, are both directly and indirectly connected with controlling the process of closure of the neural tube. Ones of those are genes of metabolism of folic acid as MTHFR, MTR, MTRR, CBS, MTHFD, folic acid receptors (FR) regulator genes from PAX family, T, PDGFRA and BRCA1 genes.