A massively parallel barcoded sequencing pipeline enables generation of the first ORFeome and interactome map for rice.

Systematic mappings of protein interactome networks have provided invaluable functional information for numerous model organisms. Here we develop PCR-mediated Linkage of barcoded Adapters To nucleic acid Elements for sequencing (PLATE-seq) that serves as a general tool to rapidly sequence thousands of DNA elements. We validate its utility by generating the ORFeome for Oryza sativa covering 2,300 genes and constructing a high-quality protein-protein interactome map consisting of 322 interactions between 289 proteins, expanding the known interactions in rice by roughly 50%. Our work paves the way for high-throughput profiling of protein-protein interactions in a wide range of organisms.

A novel iterative mixed model to remap three complex orthopedic traits in dogs.

Hip dysplasia (HD), elbow dysplasia (ED), and rupture of the cranial (anterior) cruciate ligament (RCCL) are the most common complex orthopedic traits of dogs and all result in debilitating osteoarthritis. We reanalyzed previously reported data: the Norberg angle (a quantitative measure of HD) in 921 dogs, ED in 113 cases and 633 controls, and RCCL in 271 cases and 399 controls and their genotypes at ~185,000 single nucleotide polymorphisms. A novel fixed and random model with a circulating probability unification (FarmCPU) function, with marker-based principal components and a kinship matrix to correct for population stratification, was used. A Bonferroni correction at p<0.01 resulted in a P< 6.96 ×10-8. Six loci were identified; three for HD and three for RCCL. An associated locus at CFA28:34,369,342 for HD was described previously in the same dogs using a conventional mixed model. No loci were identified for RCCL in the previous report but the two loci for ED in the previous report did not reach genome-wide significance using the FarmCPU model. These results were supported by simulation which demonstrated that the FarmCPU held no power advantage over the linear mixed model for the ED sample but provided additional power for the HD and RCCL samples. Candidate genes for HD and RCCL are discussed. When using FarmCPU software, we recommend a resampling test, that a positive control be used to determine the optimum pseudo quantitative trait nucleotide-based covariate structure of the model, and a negative control be used consisting of permutation testing and the identical resampling test as for the non-permuted phenotypes.

Multi-platform assessment of transcriptome profiling using RNA-seq in the ABRF next-generation sequencing study.

High-throughput RNA sequencing (RNA-seq) greatly expands the potential for genomics discoveries, but the wide variety of platforms, protocols and performance capabilitites has created the need for comprehensive reference data. Here we describe the Association of Biomolecular Resource Facilities next-generation sequencing (ABRF-NGS) study on RNA-seq. We carried out replicate experiments across 15 laboratory sites using reference RNA standards to test four protocols (poly-A-selected, ribo-depleted, size-selected and degraded) on five sequencing platforms (Illumina HiSeq, Life Technologies PGM and Proton, Pacific Biosciences RS and Roche 454). The results show high intraplatform (Spearman rank R > 0.86) and inter-platform (R > 0.83) concordance for expression measures across the deep-count platforms, but highly variable efficiency and cost for splice junction and variant detection between all platforms. For intact RNA, gene expression profiles from rRNA-depletion and poly-A enrichment are similar. In addition, rRNA depletion enables effective analysis of degraded RNA samples. This study provides a broad foundation for cross-platform standardization, evaluation and improvement of RNA-seq.

IQCB1 and PDE6B mutations cause similar early onset retinal degenerations in two closely related terrier dog breeds.

PURPOSE: To identify the causative mutations in two early-onset canine retinal degenerations, crd1 and crd2, segregating in the American Staffordshire terrier and the Pit Bull Terrier breeds, respectively. METHODS: Retinal morphology of crd1- and crd2-affected dogs was evaluated by light microscopy. DNA was extracted from affected and related unaffected controls. Association analysis was undertaken using the Illumina Canine SNP array and PLINK (crd1 study), or the Affymetrix Version 2 Canine array, the "MAGIC" genotype algorithm, and Fisher's Exact test for association (crd2 study). Positional candidate genes were evaluated for each disease. RESULTS: Structural photoreceptor abnormalities were observed in crd1-affected dogs as young as 11-weeks old. Rod and cone inner segment (IS) and outer segments (OS) were abnormal in size, shape, and number. In crd2-affected dogs, rod and cone IS and OS were abnormal as early as 3 weeks of age, progressing with age to severe loss of the OS, and thinning of the outer nuclear layer (ONL) by 12 weeks of age. Genome-wide association study (GWAS) identified association at the telomeric end of CFA3 in crd1-affected dogs and on CFA33 in crd2-affected dogs. Candidate gene evaluation identified a three bases deletion in exon 21 of PDE6B in crd1-affected dogs, and a cytosine insertion in exon 10 of IQCB1 in crd2-affected dogs. CONCLUSIONS: Identification of the mutations responsible for these two early-onset retinal degenerations provides new large animal models for comparative disease studies and evaluation of potential therapeutic approaches for the homologous human diseases.

Antisense sequencing improves the accuracy and precision of A-to-I editing measurements using the peak height ratio method.

BACKGROUND: A-to-I RNA editing is found in all phyla of animals and contributes to transcript diversity that may have profound impacts on behavior and physiology. Many transcripts of genes involved in axonal conductance, synaptic transmission and modulation are the targets of A-to-I RNA editing. There are a number of methods to measure the extent of A-to-I RNA editing, but they are generally costly and time consuming. One way to determine the frequency of A-to-I RNA editing is the peak height ratio method, which compares the size of peaks on electropherograms that represent unedited and edited sites. FINDINGS: Sequencing of 4 editing sites of the Dα6 nicotinic acetylcholine receptor subunit with an antisense primer (which uses T/C peaks to measure unedited and edited sites, respectively) showed very accurate and precise measurements of A-to-I RNA editing. The accuracy and precision were excellent for all editing sites, including those edited with high or low frequencies. The frequency of A-to-I RNA editing was comparable to the editing frequency as measured by clone counting from the same sample. Sequencing these same sites with the sense primer (which uses A/G peaks) yielded inaccurate and imprecise measurements. CONCLUSIONS: We have validated and improved the accuracy and precision of the peak height ratio method to measure the frequency of A-to-I RNA editing, and shown that results are primer specific. Thus, the correct sequencing primer must be utilized for the most dependable data. When compared to other methods used to measure the frequency of A-to-I RNA editing, the major benefits of the peak height ratio are that this method is inexpensive, fast, non-labor intensive and easily adaptable to many laboratory and field settings.

Evaluation of a fibrillin 2 gene haplotype associated with hip dysplasia and incipient osteoarthritis in dogs.

OBJECTIVE: To determine whether a mutation in the fibrillin 2 gene (FBN2) is associated with canine hip dysplasia (CHD) and osteoarthritis in dogs. ANIMALS: 1,551 dogs. Procedures-Hip conformation was measured radiographically. The FBN2 was sequenced from genomic DNA of 21 Labrador Retrievers and 2 Greyhounds, and a haplotype in intron 30 of FBN2 was sequenced in 90 additional Labrador Retrievers and 143 dogs of 6 other breeds. Steady-state values of FBN2 mRNA and control genes were measured in hip joint tissues of fourteen 8-month-old Labrador Retriever-Greyhound crossbreeds. RESULTS: The Labrador Retrievers homozygous for a 10-bp deletion haplotype in intron 30 of FBN2 had significantly worse CHD as measured via higher distraction index and extended-hip joint radiograph score and a lower Norberg angle and dorsolateral subluxation score. Among 143 dogs of 6 other breeds, those homozygous for the same deletion haplotype also had significantly worse radiographic CHD. Among the 14 crossbred dogs, as the dorsolateral subluxation score decreased, the capsular FBN2 mRNA increased significantly. Those dogs with incipient hip joint osteoarthritis had significantly increased capsular FBN2 mRNA, compared with those dogs without osteoarthritis. Dogs homozygous for the FBN2 deletion haplotype had significantly less FBN2 mRNA in their femoral head articular cartilage. CONCLUSIONS AND CLINICAL RELEVANCE: The FBN2 deletion haplotype was associated with CHD. Capsular gene expression of FBN2 was confounded by incipient secondary osteoarthritis in dysplastic hip joints. Genes influencing complex traits in dogs can be identified by genome-wide screening, fine mapping, and candidate gene screening.

Evaluation of quantitative trait loci for hip dysplasia in Labrador Retrievers.

OBJECTIVE: To identify the quantitative trait loci (QTL) that contribute to hip dysplasia in dogs. ANIMALS: 192 Labrador Retrievers. PROCEDURES: Hip dysplasia was measured by use of the Norberg angle (NA), dorsolateral subluxation (DLS) score, and distraction index (DI). Genome-wide screening was conducted by use of 276 unique microsatellites. Linkage analysis was performed with a variance-based linear model. Logarithm of the odds (LOD) scores were reported when values were > 2.0. RESULTS: Canis familiaris autosomes (CFAs) 01, 02, 10, 20, 22, and 32 harbored significant QTL at LOD scores > 2.0. Among the 6 QTL, the QTL on CFA02 had not been reported to harbor QTL for hip dysplasia. The highest LOD score of 3.32 on CFA20 contributed to the second principal component of the DLS score and NA of the right hip joint. The QTL that was mapped on CFA01 (LOD score of 3.13 at 55 centimorgans) was located on the same chromosome reported to harbor a QTL for hip dysplasia in Portuguese Water Dogs and German Shepherd Dogs. In this study, CFAs 10, 20, 22, and 32 harbored QTL for hip dysplasia that have been identified in a Labrador Retriever-Greyhound pedigree and in German Shepherd Dogs. CONCLUSIONS AND CLINICAL RELEVANCE: Multiple QTL were clearly involved with hip dysplasia. Identification of these QTL will enable fine-resolution mapping and subsequent assessment of candidate genes within the refined intervals to enable researchers to develop genetic screening tests and preventative and novel therapeutic regimens.

Illumina DNA test panel-based genotyping of whole genome amplified-DNA extracted from hair samples: performance and agreement with genotyping results from genomic DNA from buccal cells.

BACKGROUND: Hair is a DNA source that can be collected easily and inexpensively from participants in epidemiological studies. However, there is concern about DNA quality and quantity. Therefore, we assessed genotyping performance of whole genome amplified (WGA)-DNA extracted from hair using the GenomePlex method and evaluated its agreement with genotyping results of buccal cell DNA from the same individuals, using the Illumina GoldenGate platform. METHODS: The Illumina DNA test panel includes 360 highly validated single nucleotide polymorphisms (SNPs) selected from the Linkage IV Panel that are distributed across the entire genome. DNA was extracted from both archived hair and buccal cell samples obtained from 44 randomly selected subjects participating in a large cohort study in Canada. RESULTS: The genotyping success rate was 97.7% for 44 paired samples. However, WGA-DNA from hair failed more during genotyping in comparison to buccal cell DNA. Hair samples with a pre-WGA-DNA>or=1 ng/microL quantified using the PicoGreen assay (n=33) showed an average genotyping completion rate of 98.8% and SNP concordance of 91.2% with genotyping performance of buccal cell DNA. In contrast, samples with a pre-WGA-DNA<1 ng/microL had lower genotyping completion rate (94%) and poor SNP concordance (49%). CONCLUSIONS: Results suggest that WGA-DNA obtained from hair can produce excellent genotyping call rates and show relatively good SNP concordance with results from buccal cell DNA using high-throughput technology. DNA quantity obtained from hair samples is a crucial determinant of genotyping performance. Larger studies are needed to examine the utility of hair DNA with different genotyping platforms.

Development and characterization of a normalized canine retinal cDNA library for genomic and expression studies.

PURPOSE: Identification of causative mutations for retinal blinding disorders is often limited by restricted understanding of gene expression and underlying molecular mechanisms that trigger degenerative processes. This study was conducted to develop a catalog of canine retina-expressed genes that would provide a unique tool to investigate normal and altered function in the adult retina. Because of the conserved syntenies between the dog and human, this approach would identify new potential disease candidate genes for both species. METHODS: A canine normalized retinal cDNA library was produced and analyzed by using a modified PhredPhrap algorithm. Computerized annotation provided gene homology and chromosomal location for individual clones and contigs in a Web-accessible database. RESULTS: From 6316 cDNA clones, 3980 retinal expressed sequence tags (ESTs) were derived. Homology to the canine genome draft sequence was found for more than 99% of all ESTs, but only for 32% when compared with annotated canine cDNAs. Functional analysis suggests an enrichment of this library for genes involved with eye function and development, chaperone, or ribosomal functions when compared with mouse and human National Center for Biotechnology Information (NCBI) RefSeq entries. CONCLUSIONS: A combination of annotation approaches with ongoing mapping and expression studies provide functional data covering at least 27% to 30% of the currently proposed canine catalog of genes expressed in the retina. This is an essential first step toward establishing an integrated network for gene identification and expression patterns suitable for functional genetics, comparative genomics and evolutionary analysis of genes and gene families with respect to the developmental and degenerative processes of the retina.

A cattle-human comparative map built with cattle BAC-ends and human genome sequence.

As a step toward the goal of adding the cattle genome to those available for multispecies comparative genome analysis, 40,224 cattle BAC clones were end-sequenced, yielding 60,547 sequences (BAC end sequences, BESs) after trimming with an average read length of 515 bp. Cattle BACs were anchored to the human and mouse genome sequences by BLASTN search, revealing 29.4% and 10.1% significant hits (E < e-5), respectively. More than 60% of all cattle BES hits in both the human and mouse genomes are located within known genes. In order to confirm in silico predictions of orthologyand their relative position on cattle chromosomes, 84 cattle BESs with similarity to sequences on HSA11 were mapped using a cattle-hamster radiation hybrid (RH) panel. Resulting RH maps of BTA15 and BTA29 cover approximately 85% of HSA11 sequence, revealing a complex patchwork shuffling of segments not explained by a simple translocation followed by internal rearrangements. Overlay of the mouse conserved syntenies onto HSA11 revealed that segmental boundaries appear to be conserved in all three species. The BAC clone-based comparative map provides a foundation for the evolutionary analysis of mammalian karyotypes and for sequencing of the cattle genome.

Mutations at the mouse ichthyosis locus are within the lamin B receptor gene: a single gene model for human Pelger-Huët anomaly.

The nature of the wild-type gene product at the mouse ichthyosis (ic) locus has been of great interest because mutations at this locus cause marked abnormalities in nuclear heterochromatin, similar to those observed in Pelger-Huët anomaly (PHA). We recently found that human PHA is caused by mutations in the gene (LBR) encoding lamin B receptor, an evolutionarily conserved inner nuclear membrane protein involved in nuclear assembly and chromatin binding. Mice homozygous for deleterious alleles at the ichthyosis (ic) locus present with a blood phenotype similar to PHA, and develop other phenotypic abnormalities, including alopecia, variable expression of syndactyly and hydrocephalus. The ic locus on mouse chromosome 1 shares conserved synteny with the chromosomal location of the human LBR locus on human chromosome 1. In this study, we identified one nonsense (815ins) and two frameshift mutations (1088insCC and 1884insGGAA) within the Lbr gene of mice homozygous for either of three independent mutations (ic, ic(J) and ic(4J), respectively) at the ichthyosis locus. These allelic mutations are predicted to result in truncated or severely impaired LBR protein. Our studies of mice homozygous for the ic(J) mutation revealed a complete loss of LBR protein as shown by immunofluorescence microscopy and immunoblotting. The findings provide the molecular basis for the heterochromatin clumping and other distinct phenotypes caused by ic mutations. These spontaneous Lbr mutations confirm the molecular basis of human PHA and provide a small animal model for determination of the precise function of LBR in normal and pathological states.