Fine-mapping and mutation analysis of TRPM1: a candidate gene for leopard complex (LP) spotting and congenital stationary night blindness in horses.

Leopard Complex spotting occurs in several breeds of horses and is caused by an incompletely dominant allele (LP). Homozygosity for LP is also associated with congenital stationary night blindness (CSNB) in Appaloosa horses. Previously, LP was mapped to a 6 cm region on ECA1 containing the candidate gene TRPM1 (Transient Receptor Potential Cation Channel, Subfamily M, Member 1) and decreased expression of this gene, measured by qRT-PCR, was identified as the likely cause of both spotting and ocular phenotypes. This study describes investigations for a mutation causing or associated with the Leopard Complex and CSNB phenotype in horses. Re-sequencing of the gene and associated splice sites within the 105 624 bp genomic region of TRPM1 led to the discovery of 18 SNPs. Most of the SNPs did not have a predictive value for the presence of LP. However, one SNP (ECA1:108,249,293 C>T) found within intron 11 had a strong (P < 0.0005), but not complete, association with LP and CSNB and thus is a good marker but unlikely to be causative. To further localize the association, 70 SNPs spanning over two Mb including the TRPM1 gene were genotyped in 192 horses from three different breeds segregating for LP. A single 173 kb haplotype associated with LP and CSNB (ECA1: 108,197,355- 108,370,150) was identified. Illumina sequencing of 300 kb surrounding this haplotype revealed 57 SNP variants. Based on their localization within expressed sequences or regions of high sequence conservation across mammals, six of these SNPs were considered to be the most likely candidate mutations. While the precise function of TRPM1 remains to be elucidated, this work solidifies its functional role in both pigmentation and night vision. Further, this work has identified several potential regulatory elements of the TRPM1 gene that should be investigated further in this and other species.

KLIKK proteases of Tannerella forsythia: putative virulence factors with a unique domain structure.

Comparative genomics of virulent Tannerella forsythia ATCC 43037 and a close health-associated relative, Tannerella BU063, revealed, in the latter, the absence of an entire array of genes encoding putative secretory proteases that possess a nearly identical C-terminal domain (CTD) that ends with a -Lys-Leu-Ile-Lys-Lys motif. This observation suggests that these proteins, referred to as KLIKK proteases, may function as virulence factors. Re-sequencing of the loci of the KLIKK proteases found only six genes grouped in two clusters. All six genes were expressed by T. forsythia in routine culture conditions, although at different levels. More importantly, a transcript of each gene was detected in gingival crevicular fluid (GCF) from periodontitis sites infected with T. forsythia indicating that the proteases are expressed in vivo. In each protein, a protease domain was flanked by a unique N-terminal profragment and a C-terminal extension ending with the CTD. Partially purified recombinant proteases showed variable levels of proteolytic activity in zymography gels and toward protein substrates, including collagen, gelatin, elastin, and casein. Taken together, these results indicate that the pathogenic strain of T. forsythia secretes active proteases capable of degrading an array of host proteins, which likely represents an important pathogenic feature of this bacterium.

TECPR2 Associated Neuroaxonal Dystrophy in Spanish Water Dogs.

Clinical, pathological and genetic examination revealed an as yet uncharacterized juvenile-onset neuroaxonal dystrophy (NAD) in Spanish water dogs. Affected dogs presented with various neurological deficits including gait abnormalities and behavioral deficits. Histopathology demonstrated spheroid formation accentuated in the grey matter of the cerebral hemispheres, the cerebellum, the brain stem and in the sensory pathways of the spinal cord. Iron accumulation was absent. Ultrastructurally spheroids contained predominantly closely packed vesicles with a double-layered membrane, which were characterized as autophagosomes using immunohistochemistry. The family history of the four affected dogs suggested an autosomal recessive inheritance. SNP genotyping showed a single genomic region of extended homozygosity of 4.5 Mb in the four cases on CFA 8. Linkage analysis revealed a maximal parametric LOD score of 2.5 at this region. By whole genome re-sequencing of one affected dog, a perfectly associated, single, non-synonymous coding variant in the canine tectonin beta-propeller repeat-containing protein 2 (TECPR2) gene affecting a highly conserved region was detected (c.4009C>T or p.R1337W). This canine NAD form displays etiologic parallels to an inherited TECPR2 associated type of human hereditary spastic paraparesis (HSP). In contrast to the canine NAD, the spinal cord lesions in most types of human HSP involve the sensory and the motor pathways. Furthermore, the canine NAD form reveals similarities to cases of human NAD defined by widespread spheroid formation without iron accumulation in the basal ganglia. Thus TECPR2 should also be considered as candidate gene for human NAD. Immunohistochemistry and the ultrastructural findings further support the assumption, that TECPR2 regulates autophagosome accumulation in the autophagic pathways. Consequently, this report provides the first genetic characterization of juvenile canine NAD, describes the histopathological features associated with the TECPR2 mutation and provides evidence to emphasize the association between failure of autophagy and neurodegeneration.

Genome-wide association studies based on sequence-derived genotypes reveal new QTL associated with conformation and performance traits in the Franches-Montagnes horse breed

To identify novel quantitative trait loci (QTL) within horses, we performed genome-wide association studies (GWAS) based on sequence-level genotypes for conformation and performance traits in the Franches-Montagnes (FM) horse breed. Sequence-level genotypes of FM horses were derived by re-sequencing 30 key founders and imputing 50K data of genotyped horses. In total, we included 1077 FM horses genotyped for ~4 million SNPs and their respective de-regressed breeding values of the traits in the analysis. Based on this dataset, we identified a total of 14 QTL associated with 18 conformation traits and one performance trait. Therefore, our results suggest that the application of sequence-derived genotypes increases the power to identify novel QTL which were not identified previously based on 50K SNP chip data.

Comparison of two statistical methods for rare variant association analysis

Studies have shown that rare genetic variants have stronger effects in predisposing common diseases, and several statistical methods have been developed for association studies involving rare variants. In order to better understand how these statistical methods perform, we seek to compare two recently developed rare variant statistical methods (VT and C-alpha) on 10,000 simulated re-sequencing data sets with disease status and the corresponding 10,000 simulated null data sets. The SLC1A1 gene has been suggested to be associated with diastolic blood pressure (DBP) in previous studies. In the current study, we applied VT and C-alpha methods to the empirical re-sequencing data for the SLC1A1 gene from 300 whites and 200 blacks. We found that VT method obtains higher power and performs better than C-alpha method with the simulated data we used. The type I errors were well-controlled for both methods. In addition, both VT and C-alpha methods suggested no statistical evidence for the association between the SLC1A1 gene and DBP. Overall, our findings provided an important comparison of the two statistical methods for future reference and provided preliminary and pioneer findings on the association between the SLC1A1 gene and blood pressure.^

Déterminer le rôle de C1orf106, un gène associé aux maladies inflammatoires de l’intestin

Les maladies inflammatoires de l’intestin (MIIs, [MIM 266600]) sont caractérisées par une inflammation chronique au niveau du tube gastro-intestinal. Les deux principales formes sont la maladie de Crohn (MC) et la colite ulcéreuse (CU). Les MIIs résulteraient d’un défaut du système immunitaire et de l’épithélium intestinal. Ce dernier forme une barrière physique et biochimique qui sépare notre système immunitaire des microorganismes commensaux et pathogènes de la microflore intestinale. Un défaut dans la barrière épithéliale intestinale pourrait donc mener à une réponse immunitaire soutenue contre notre microflore intestinale. Les études d’association pangénomiques (GWAS) ont permis d’identifier 201 régions de susceptibilité aux MIIs. Parmi celles-ci, la région 1q32 associée à la MC (p<2x10-11) et à la CU (p<6x10-7) contient 4 gènes, dont C1orf106, un gène codant pour une protéine de fonction inconnue. Le re-séquençage de la région 1q32 a permis d’identifier une variante génétique rare de C1orf106 (MAF˂1%) associée aux MIIs (p=0,009), Y333F. Nous avons démontré que la substitution de la tyr333 par une phénylalanine semble avoir un effet sur la stabilité protéique de C1orf106 tel que démontré lors de l’inhibition de la synthèse protéique induite par le cycloheximide. Nous avons déterminé que C1orf106 est exprimé dans le côlon et l’intestin grêle. De plus, son expression est augmentée lors de la différenciation des cellules épithéliales Caco-2 en épithélium intestinal polarisé. Son profil d’expression correspond aux types cellulaires et tissulaires affectés dans les MIIs. De plus, C1orf106 est partiellement co-localisée avec le marqueur des jonctions serrées, ZO-1. Toutefois, son marquage reproduit parfaitement celui du marqueur des jonctions adhérentes, E-cadhérine. Les jonctions serrées et adhérentes sont localisées du côté apical de la jonction intercellulaire et sont toutes deux impliquées dans l’établissement de la barrière épithéliale. Nous avons donc testé l’impact de C1orf106 sur la perméabilité de l’épithélium intestinal. Nous avons observé une augmentation de la perméabilité épithéliale chez un épithélium intestinal formé par des cellules Caco-2 sous-exprimant C1orf106. Nos résultats suggèrent que C1orf106 pourrait être le gène causal de la région 1q32.

Déterminer le rôle de C1orf106, un gène associé aux maladies inflammatoires de l’intestin

Les maladies inflammatoires de l’intestin (MIIs, [MIM 266600]) sont caractérisées par une inflammation chronique au niveau du tube gastro-intestinal. Les deux principales formes sont la maladie de Crohn (MC) et la colite ulcéreuse (CU). Les MIIs résulteraient d’un défaut du système immunitaire et de l’épithélium intestinal. Ce dernier forme une barrière physique et biochimique qui sépare notre système immunitaire des microorganismes commensaux et pathogènes de la microflore intestinale. Un défaut dans la barrière épithéliale intestinale pourrait donc mener à une réponse immunitaire soutenue contre notre microflore intestinale. Les études d’association pangénomiques (GWAS) ont permis d’identifier 201 régions de susceptibilité aux MIIs. Parmi celles-ci, la région 1q32 associée à la MC (p<2x10-11) et à la CU (p<6x10-7) contient 4 gènes, dont C1orf106, un gène codant pour une protéine de fonction inconnue. Le re-séquençage de la région 1q32 a permis d’identifier une variante génétique rare de C1orf106 (MAF˂1%) associée aux MIIs (p=0,009), Y333F. Nous avons démontré que la substitution de la tyr333 par une phénylalanine semble avoir un effet sur la stabilité protéique de C1orf106 tel que démontré lors de l’inhibition de la synthèse protéique induite par le cycloheximide. Nous avons déterminé que C1orf106 est exprimé dans le côlon et l’intestin grêle. De plus, son expression est augmentée lors de la différenciation des cellules épithéliales Caco-2 en épithélium intestinal polarisé. Son profil d’expression correspond aux types cellulaires et tissulaires affectés dans les MIIs. De plus, C1orf106 est partiellement co-localisée avec le marqueur des jonctions serrées, ZO-1. Toutefois, son marquage reproduit parfaitement celui du marqueur des jonctions adhérentes, E-cadhérine. Les jonctions serrées et adhérentes sont localisées du côté apical de la jonction intercellulaire et sont toutes deux impliquées dans l’établissement de la barrière épithéliale. Nous avons donc testé l’impact de C1orf106 sur la perméabilité de l’épithélium intestinal. Nous avons observé une augmentation de la perméabilité épithéliale chez un épithélium intestinal formé par des cellules Caco-2 sous-exprimant C1orf106. Nos résultats suggèrent que C1orf106 pourrait être le gène causal de la région 1q32.

Lexical access, story re-telling and sequencing skills in adults who clutter and those who do not

Cluttering is a rate-based disorder of fluency, the scope of whose diagnostic criteria currently remains unclear. This paper reports preliminary findings from a larger study which aims to determine whether cluttering can be associated with language disturbances as well as motor and rate based ones. Subtests from the Mt Wilga High Level Language Test (MWHLLT) were used to determine whether people who clutter (PWC) have word finding difficulties, and use significantly more maze behaviours compared to controls, during story re-telling and simple sequencing tasks. Independent t tests showed that PWC were significantly slower than control participants in lexical access and sentence completion tasks, but returned mixed findings when PWCs were required to name items within a semantic category. PWC produced significantly more maze behaviour than controls in a task where participants were required to explain how to undertake commonly performed actions, but no difference in use of maze behaviour was found between the two groups when retelling a story from memory. The implications of these findings are discussed

miRDeep*: an integrated application tool for miRNA identification from RNA sequencing data

miRDeep and its varieties are widely used to quantify known and novel micro RNA (miRNA) from small RNA sequencing (RNAseq). This article describes miRDeep*, our integrated miRNA identification tool, which is modeled off miRDeep, but the precision of detecting novel miRNAs is improved by introducing new strategies to identify precursor miRNAs. miRDeep* has a user-friendly graphic interface and accepts raw data in FastQ and Sequence Alignment Map (SAM) or the binary equivalent (BAM) format. Known and novel miRNA expression levels, as measured by the number of reads, are displayed in an interface, which shows each RNAseq read relative to the pre-miRNA hairpin. The secondary pre-miRNA structure and read locations for each predicted miRNA are shown and kept in a separate figure file. Moreover, the target genes of known and novel miRNAs are predicted using the TargetScan algorithm, and the targets are ranked according to the confidence score. miRDeep* is an integrated standalone application where sequence alignment, pre-miRNA secondary structure calculation and graphical display are purely Java coded. This application tool can be executed using a normal personal computer with 1.5 GB of memory. Further, we show that miRDeep* outperformed existing miRNA prediction tools using our LNCaP and other small RNAseq datasets. miRDeep* is freely available online at http://www.australianprostatecentre.org/research/software/mirdeep-star

H5N1 avian influenza re-emergence of Lake Communication Qinghai: phylogenetic and antigenic analyses of the newly isolated viruses and roles of migratory birds in virus circulation

Highly pathogenic avian influenza H5N1 virus has swept west across the globe and caused serious debates on the roles of migratory birds in virus circulation since the first large-scale outbreak in migratory birds of Lake Qinghai, 2005. In May 2006, another outbreak struck Lake Qinghai and six novel strains were isolated. To elucidate these QH06 viruses, the six isolates were subjected to whole-genome sequencing. Phylogenetic analyses show that QH06 viruses are derived from the lineages of Lake Qinghai, 2005. Five of the six novel isolates are adjacent to the strain A/Cygnus olor/Croatia/1/05, and the last one is related to the strain A/duck/Novosibirsk/ 02/05, an isolate of the flyway. Antigenic analyses suggest that QH06 and QH05 viruses are similar to each other. These findings implicate that QH06 viruses of Lake Qinghai may travel back via migratory birds, though not ruling out the possibility of local circulation of viruses of Lake Qinghai.

Sequencing three crocodilian genomes to illuminate the evolution of archosaurs and amniotes.

The International Crocodilian Genomes Working Group (ICGWG) will sequence and assemble the American alligator (Alligator mississippiensis), saltwater crocodile (Crocodylus porosus) and Indian gharial (Gavialis gangeticus) genomes. The status of these projects and our planned analyses are described.

A heuristic algorithm for two-machine re-entrant shop scheduling

This paper considers the problem of sequencing n jobs in a two‐machine re‐entrant shopwith the objective of minimizing the maximum completion time. The shop consists of twomachines, M1 and M2 , and each job has the processing route (M1 , M2 , M1 ). An O(n log n)time heuristic is presented which generates a schedule with length at most 4/3 times that ofan optimal schedule, thereby improving the best previously available worst‐case performanceratio of 3/2.

Connectivity Mapping for Candidate Therapeutics Identification Using Next Generation Sequencing RNA-Seq Data

The advent of next generation sequencing technologies (NGS) has expanded the area of genomic research, offering high coverage and increased sensitivity over older microarray platforms. Although the current cost of next generation sequencing is still exceeding that of microarray approaches, the rapid advances in NGS will likely make it the platform of choice for future research in differential gene expression. Connectivity mapping is a procedure for examining the connections among diseases, genes and drugs by differential gene expression initially based on microarray technology, with which a large collection of compound-induced reference gene expression profiles have been accumulated. In this work, we aim to test the feasibility of incorporating NGS RNA-Seq data into the current connectivity mapping framework by utilizing the microarray based reference profiles and the construction of a differentially expressed gene signature from a NGS dataset. This would allow for the establishment of connections between the NGS gene signature and those microarray reference profiles, alleviating the associated incurring cost of re-creating drug profiles with NGS technology. We examined the connectivity mapping approach on a publicly available NGS dataset with androgen stimulation of LNCaP cells in order to extract candidate compounds that could inhibit the proliferative phenotype of LNCaP cells and to elucidate their potential in a laboratory setting. In addition, we also analyzed an independent microarray dataset of similar experimental settings. We found a high level of concordance between the top compounds identified using the gene signatures from the two datasets. The nicotine derivative cotinine was returned as the top candidate among the overlapping compounds with potential to suppress this proliferative phenotype. Subsequent lab experiments validated this connectivity mapping hit, showing that cotinine inhibits cell proliferation in an androgen dependent manner. Thus the results in this study suggest a promising prospect of integrating NGS data with connectivity mapping. © 2013 McArt et al.

Characterisation of a novel panel of polymorphic microsatellite loci for the liver fluke, Fasciola hepatica, using a next generation sequencing approach

The liver fluke, Fasciola hepatica is an economically important pathogen of sheep and cattle and has been described by the WHO as a re-emerging zoonosis. Control is heavily reliant on the use of drugs, particularly triclabendazole and as a result resistance has now emerged. The population structure of F. hepatica is not well known, yet it can impact on host-parasite interactions and parasite control with drugs, particularly regarding the spread of triclabendazole resistance. We have identified 2448 potential microsatellites from 83Mb of F. hepatica genome sequence using msatfinder. Thirty-five loci were developed and optimised for microsatellite PCR, resulting in a panel of 15 polymorphic loci, with a range of three to 15 alleles. This panel was validated on genomic DNA from 46 adult F. hepatica; 38 liver flukes sourced from a Northwest abattoir, UK and 8 liver flukes from an established isolate (Shrewsbury; Ridgeway Research). Evidence for null alleles was found at four loci (Fh_1, Fh_8, Fh_13 and Fh_14), which showed markedly higher levels of homozygosity than the remaining 11 loci. Of the 38 liver flukes isolated from cattle livers (n=10) at the abattoir, 37 genotypes were identified. Using a multiplex approach all 15 loci could be amplified from several life cycle stages that typically yield low amounts of DNA, including metacercariae, the infective life cycle stage present on pasture, highlighting the utility of this multiplex microsatellite panel. This study reports the largest panel of microsatellite markers available to date for population studies of F. hepatica and the first multiplex panel of microsatellite markers that can be used for several life cycle stages.

Gene activation and re-expression of a trypanosoma-brucei variant surface glycoptotein

The expression of the Trypanosoma brucei variant surface glycoprotein AnTat 1.1 proceeds by a mechanism that transfers a duplicated gene copy into a new genomic environment, the so-called expression site, where it will be expressed. We have isolated a genomic fragment containing the region spanning the expression site-transposon junction, and the 5' half of the coding sequence. Comparing this DNA segment with its template copy (basic copy) allowed us to identify the exact breaking point and indicated a base sequence which could be involved in initiating the transposition event. Sequencing data also indicated that the co-transposed segment 5' to the coding sequence is 430 bp in length. The extreme 5' end of the mRNA is derived from a region in the expression site not immediately adjacent to the transposed DNA segment. This particular sequence exists in multiple copies in the genome and is common to the mRNA of all variant surface glycoproteins so far analysed.