Disruption of global hypothalamic microRNA (miR) profiles and associated behavioral changes in California mice (Peromyscus californicus) developmentally exposed to endocrine disrupting chemicals.

Developmental exposure to endocrine disrupting chemicals (EDCs), e.g., bisphenol A (BPA) or genistein (GEN), causes longstanding epigenome effects. MicroRNAs (miRs) regulate which mRNAs will be translated to proteins and thereby serve as the final checkpoint in epigenetic control. Scant amount is known, however, whether EDCs affect neural miRNA (miR) patterns. We aimed to test the hypothesis that developmental exposure of California mice (Peromyscus californicus) to GEN, BPA, or both chemicals influences hypothalamic miR/small RNA profiles and ascertain the extent such biomolecular alterations correlate with behavioral and metabolic changes. California mice were developmentally exposed to GEN (250 mg/kg feed weight, FW), GEN (250 mg/kg FW)+BPA (5 mg/kg FW), low dose (LD) BPA (5 mg/kg FW), or upper dose (UD) BPA (50 mg/kg FW). Adult offspring were tested in a battery of behavioral and metabolic tests; whereupon, mice were euthanized, brains were collected and frozen, small RNAs were isolated from hypothalamic punches, and subsequently sequenced. California mice exposed to one or both EDCs engaged in one or more repetitive behaviors. GEN, LD BPA, and UD BPA altered aspects of ultrasonic and audible vocalizations. Each EDC exposure led to sex-dependent differences in differentially expressed miR/small RNAs with miR7-2, miR146, and miR148a being increased in all female and male EDC exposed groups. Current findings reveal that developmental exposure to GEN and/or BPA affects hypothalamic miR/small RNA expression patterns, and such changes correlate with EDC-induced behavioral and metabolic alterations. miR146 is likely an important mediator and biomarker of EDC exposure in mammals, including humans.

Refinement of the canine CD1 locus topology and investigation of antibody binding to recombinant canine CD1 isoforms.

CD1 molecules are antigen-presenting glycoproteins primarily found on dendritic cells (DCs) responsible for lipid antigen presentation to CD1-restricted T cells. Despite their pivotal role in immunity, little is known about CD1 protein expression in dogs, notably due to lack of isoform-specific antibodies. The canine (Canis familiaris) CD1 locus was previously found to contain three functional CD1A genes: canCD1A2, canCD1A6, and canCD1A8, where two variants of canCD1A8, canCD1A8.1 and canCD1A8.2, were assumed to be allelic variants. However, we hypothesized that these rather represented two separate genes. Sequencing of three overlapping bacterial artificial chromosomes (BACs) spanning the entire canine CD1 locus revealed canCD1A8.2 and canCD1A8.1 to be located in tandem between canCD1A7 and canCD1C, and canCD1A8.1 was consequently renamed canCD1A9. Green fluorescent protein (GFP)-fused canine CD1 transcripts were recombinantly expressed in 293T cells. All proteins showed a highly positive GFP expression except for canine CD1d and a splice variant of canine CD1a8 lacking exon 3. Probing with a panel of anti-CD1 monoclonal antibodies (mAbs) showed that Ca13.9H11 and Ca9.AG5 only recognized canine CD1a8 and CD1a9 isoforms, and Fe1.5F4 mAb solely recognized canine CD1a6. Anti-CD1b mAbs recognized the canine CD1b protein, but also bound CD1a2, CD1a8, and CD1a9. Interestingly, Ca9.AG5 showed allele specificity based on a single nucleotide polymorphism (SNP) located at position 321. Our findings have refined the structure of the canine CD1 locus and available antibody specificity against canine CD1 proteins. These are important fundamentals for future investigation of the role of canine CD1 in lipid immunity.

Assemblathon 1: a competitive assessment of de novo short read assembly methods.

Low-cost short read sequencing technology has revolutionized genomics, though it is only just becoming practical for the high-quality de novo assembly of a novel large genome. We describe the Assemblathon 1 competition, which aimed to comprehensively assess the state of the art in de novo assembly methods when applied to current sequencing technologies. In a collaborative effort, teams were asked to assemble a simulated Illumina HiSeq data set of an unknown, simulated diploid genome. A total of 41 assemblies from 17 different groups were received. Novel haplotype aware assessments of coverage, contiguity, structure, base calling, and copy number were made. We establish that within this benchmark: (1) It is possible to assemble the genome to a high level of coverage and accuracy, and that (2) large differences exist between the assemblies, suggesting room for further improvements in current methods. The simulated benchmark, including the correct answer, the assemblies, and the code that was used to evaluate the assemblies is now public and freely available from http://www.assemblathon.org/.

Identifying human-rhesus macaque gene orthologs using heterospecific SNP probes.

We genotyped a Chinese and an Indian-origin rhesus macaque using the Affymetrix Genome-Wide Human SNP Array 6.0 and cataloged 85,473 uniquely mapping heterospecific SNPs. These SNPs were assigned to rhesus chromosomes according to their probe sequence alignments as displayed in the human and rhesus reference sequences. The conserved gene order (synteny) revealed by heterospecific SNP maps is in concordance with that of the published human and rhesus macaque genomes. Using these SNPs' original human rs numbers, we identified 12,328 genes annotated in humans that are associated with these SNPs, 3674 of which were found in at least one of the two rhesus macaques studied. Due to their density, the heterospecific SNPs allow fine-grained comparisons, including approximate boundaries of intra- and extra-chromosomal rearrangements involving gene orthologs, which can be used to distinguish rhesus macaque chromosomes from human chromosomes.

Comparative analysis of immune infiltrates in head and neck cancers across anatomical sites.

BACKGROUND: The response rate to immune checkpoint inhibitors targeting programmed cell death 1 (PD-1) receptor is 13%-18% for patients with recurrent or metastatic head and neck squamous cell carcinoma (HNSCC). Detailed understanding of the tumor immune microenvironment (TIME) is crucial in order to explain and improve this response rate. HNSCCs arise at various anatomical locations including the oral cavity, hypopharynx, larynx and oropharynx. Studies directly comparing immune infiltration between anatomical sites are scarce. Since the distinct locations could drive deviating microenvironments, we questioned whether the immune composition varies across these HNSCC sites. METHODS: Here, we characterized the TIME of 76 fresh tumor specimens using flow cytometry and performed single-cell RNA-sequencing on nine head and neck tumor samples. RESULTS: We found major differences in the composition of the TIME between patients. When comparing anatomical sites: tumors originating from the oral cavity had higher T cell infiltrates than tumors from other anatomical sites. The percentage of tumor-infiltrating T-lymphocytes positive for the immune checkpoint PD-1 varied considerably between patients, with the highest fraction of PD-1+ T cells found in larynx squamous cell carcinomas (SCCs). While we had hypothesized that the anatomical sites of tumor origin would drive sample clustering, our data showed that the type of TIME was more dominant and was particularly driven by the fraction of T cells positive for PD-1. Moreover, a high proportion of PD-1+ CD8+ T cells associated with an improved overall survival. Using single-cell RNA-sequencing, we observed that PD-1 expression was highest in the CD8-ENTPD1 tissue resident memory T cell/exhausted T cell and CD4-CXCL13 type 1 T helper cell clusters. CONCLUSIONS: We found that oral cavity SCCs had the highest frequencies of T cells. We also observed considerable interpatient heterogeneity for PD-1 on T cells, with noticeably higher frequencies of PD-1+ CD4+ T helper cells in larynx SCCs. Within the entire cohort, a higher fraction of CD8+ T cells positive for PD-1 was linked to improved overall survival. Whether the fraction of PD-1+ T cells within the TIME enables immune checkpoint inhibitor response prediction for patients with head and neck cancer remains to be determined.

Probing functional polymorphisms in the dengue vector, Aedes aegypti.

BACKGROUND: Dengue is the most prevalent arboviral disease world-wide and its primary vector is the mosquito Aedes aegypti. The current lack of commercially-available vaccines makes control of vector populations the only effective strategy to prevent dengue transmission. Aedes aegypti geographic populations exhibit great variability in insecticide resistance and susceptibility to dengue infection. The characterization of single nucleotide polymorphisms (SNPs) as molecular markers to study quantitatively this variation is needed greatly because this species has a low abundance of microsatellite markers and limited known restriction fragments length polymorphisms (RFLPs) and single-strand conformation polymorphism (SSCP) markers. RESULTS: We used RNA-seq to characterize SNPs in three Ae. aegypti strains, including the Liverpool (LVP) strain, from which the current genome annotation is derived. We identified 131,764 unique genome locations with at least one alternative nucleotide to what is reported in the reference annotation. These comprised changes in both open-reading frames (ORFs) and untranslated regions (UTRs) of transcripts. An in depth-look at sequence variation in immunity genes revealed that those associated with autophagy, MD2-like receptors and Peptidoglycan Recognition Proteins had more sequence variation in their 3'UTRs than mutations associated with non-synonymous changes. This supports the conclusion that these genes had maintained their functional specificity while being adapted to different regulatory domains. In contrast, a number of peroxidases, serpins and Clip-domain serine proteases exhibited conservation of putative UTR regulatory sequences while displaying diversification of the ORFs. Transcriptome evidence also was found for ~2500 novel transcriptional units (NTUs) not annotated in the reference genome. CONCLUSIONS: The transcriptome-wide assessment of within and inter-strain polymorphisms in Ae. aegypti adds considerably to the number of molecular markers available for genetic studies in this mosquito. Additionally, data supporting NTU discovery emphasizes the need for continuous amendments of the reference genome annotation.

Reference genome-independent assessment of mutation density using restriction enzyme-phased sequencing.

BACKGROUND: The availability of low cost sequencing has spurred its application to discovery and typing of variation, including variation induced by mutagenesis. Mutation discovery is challenging as it requires a substantial amount of sequencing and analysis to detect very rare changes and distinguish them from noise. Also challenging are the cases when the organism of interest has not been sequenced or is highly divergent from the reference. RESULTS: We describe the development of a simple method for reduced representation sequencing. Input DNA was digested with a single restriction enzyme and ligated to Y adapters modified to contain a sequence barcode and to provide a compatible overhang for ligation. We demonstrated the efficiency of this method at SNP discovery using rice and arabidopsis. To test its suitability for the discovery of very rare SNP, one control and three mutagenized rice individuals (1, 5 and 10 mM sodium azide) were used to prepare genomic libraries for Illumina sequencers by ligating barcoded adapters to NlaIII restriction sites. For genome-dependent discovery 15-30 million of 80 base reads per individual were aligned to the reference sequence achieving individual sequencing coverage from 7 to 15×. We identified high-confidence base changes by comparing sequences across individuals and identified instances consistent with mutations, i.e. changes that were found in a single treated individual and were solely GC to AT transitions. For genome-independent discovery 70-mers were extracted from the sequence of the control individual and single-copy sequence was identified by comparing the 70-mers across samples to evaluate copy number and variation. This de novo "genome" was used to align the reads and identify mutations as above. Covering approximately 1/5 of the 380 Mb genome of rice we detected mutation densities ranging from 0.6 to 4 per Mb of diploid DNA depending on the mutagenic treatment. CONCLUSIONS: The combination of a simple and cost-effective library construction method, with Illumina sequencing, and the use of a bioinformatic pipeline allows practical SNP discovery regardless of whether a genomic reference is available.

Genome-wide SNP discovery in walnut with an AGSNP pipeline updated for SNP discovery in allogamous organisms.

BACKGROUND: A genome-wide set of single nucleotide polymorphisms (SNPs) is a valuable resource in genetic research and breeding and is usually developed by re-sequencing a genome. If a genome sequence is not available, an alternative strategy must be used. We previously reported the development of a pipeline (AGSNP) for genome-wide SNP discovery in coding sequences and other single-copy DNA without a complete genome sequence in self-pollinating (autogamous) plants. Here we updated this pipeline for SNP discovery in outcrossing (allogamous) species and demonstrated its efficacy in SNP discovery in walnut (Juglans regia L.). RESULTS: The first step in the original implementation of the AGSNP pipeline was the construction of a reference sequence and the identification of single-copy sequences in it. To identify single-copy sequences, multiple genome equivalents of short SOLiD reads of another individual were mapped to shallow genome coverage of long Sanger or Roche 454 reads making up the reference sequence. The relative depth of SOLiD reads was used to filter out repeated sequences from single-copy sequences in the reference sequence. The second step was a search for SNPs between SOLiD reads and the reference sequence. Polymorphism within the mapped SOLiD reads would have precluded SNP discovery; hence both individuals had to be homozygous. The AGSNP pipeline was updated here for using SOLiD or other type of short reads of a heterozygous individual for these two principal steps. A total of 32.6X walnut genome equivalents of SOLiD reads of vegetatively propagated walnut scion cultivar 'Chandler' were mapped to 48,661 'Chandler' bacterial artificial chromosome (BAC) end sequences (BESs) produced by Sanger sequencing during the construction of a walnut physical map. A total of 22,799 putative SNPs were initially identified. A total of 6,000 Infinium II type SNPs evenly distributed along the walnut physical map were selected for the construction of an Infinium BeadChip, which was used to genotype a walnut mapping population having 'Chandler' as one of the parents. Genotyping results were used to adjust the filtering parameters of the updated AGSNP pipeline. With the adjusted filtering criteria, 69.6% of SNPs discovered with the updated pipeline were real and could be mapped on the walnut genetic map. A total of 13,439 SNPs were discovered by BES re-sequencing. BESs harboring SNPs were in 677 FPC contigs covering 98% of the physical map of the walnut genome. CONCLUSION: The updated AGSNP pipeline is a versatile SNP discovery tool for a high-throughput, genome-wide SNP discovery in both autogamous and allogamous species. With this pipeline, a large set of SNPs were identified in a single walnut cultivar.

Multiple clonal MLL fusions in a patient receiving CHOP-based chemotherapy.

MLL rearrangements were analysed in the blood of a patient receiving chemotherapy for diffuse large B-cell lymphoma using inverse polymerase chain reaction targeting exon 12, parallel sequencing and a custom algorithm design. Of thirteen MLL rearrangements detected, five were capable of generating MLL fusion genes, including MLL-MLLT3, the most common fusion in acute myeloid leukaemia (AML). Other fusions, all previously clinically unobserved, included MLL-NKD1, a fusion to the negative regulator of Wnt/ß-catenin signaling, a pathway linked to leukaemic cell proliferation. The majority of the fusions exhibited clonal persistence from before treatment until 6 months post-chemotherapy, suggesting the fusions may confer a survival advantage to the mutant clone. MLL breakpoints were partly clustered at a specific location, indicating commonality in the process of their formation. Further, the same MLL breakpoint location exhibited a 50-100-fold increase in C to T transitions, consistent with attack by activation-induced cytidine deaminase (AICDA). As is also observed in AML and acute lymphoblastic leukaemia, in this single patient setting, MLL is capable of interacting with multiple fusion partners. This finding defines a discrete site of MLL susceptibility to fragmentation, linked to possible deregulation of AICDA function.

Discovery of rare mutations in populations: TILLING by sequencing.

Discovery of rare mutations in populations requires methods, such as TILLING (for Targeting Induced Local Lesions in Genomes), for processing and analyzing many individuals in parallel. Previous TILLING protocols employed enzymatic or physical discrimination of heteroduplexed from homoduplexed target DNA. Using mutant populations of rice (Oryza sativa) and wheat (Triticum durum), we developed a method based on Illumina sequencing of target genes amplified from multidimensionally pooled templates representing 768 individuals per experiment. Parallel processing of sequencing libraries was aided by unique tracer sequences and barcodes allowing flexibility in the number and pooling arrangement of targeted genes, species, and pooling scheme. Sequencing reads were processed and aligned to the reference to identify possible single-nucleotide changes, which were then evaluated for frequency, sequencing quality, intersection pattern in pools, and statistical relevance to produce a Bayesian score with an associated confidence threshold. Discovery was robust both in rice and wheat using either bidimensional or tridimensional pooling schemes. The method compared favorably with other molecular and computational approaches, providing high sensitivity and specificity.

T Cells Engineered to Express Immunoreceptors Targeting the Frequently Expressed Medullary Thyroid Cancer Antigens Calcitonin, CEA, and RET M918T.

Background: Medullary thyroid cancer (MTC) is a rare malignancy originating from the calcitonin-producing C cells of the thyroid. Despite recent therapeutic advances, metastatic MTC remains incurable. Adoptive cell therapy (ACT) using genetically engineered T cells targeting either tissue-restricted tumor-associated antigens or mutated neoantigens has led to durable remissions in other metastatic solid tumors. The majority of MTC express the tumor-associated antigens calcitonin and carcinoembryonic antigen (CEA), and ∼40% of MTC harbor the RET M918T oncogenic driver mutation. Methods: We developed and characterized three immunoreceptors that recognize extracellular CEA, a calcitonin epitope presented by HLA-A*24:02, or an RET M918T neoepitope restricted by HLA-DPB1*04:01/02. The chimeric antigen receptor (CAR) targeting CEA was synthetically designed, while the T cell receptors (TCRs) targeting calcitonin and RET M918T were isolated from a transgenic mouse and patient with MTC, respectively. These immunoreceptors were genetically engineered into peripheral blood T cells and tested for antigen specificity and antitumor activity. Results: T cells expressing the anti-CEA CAR or the calcitonin-reactive TCR produced effector cytokines and displayed cytotoxicity against cell lines expressing their cognate antigen in vitro. In immunodeficient mice harboring a human MTC cell line, the adoptive transfer of T cells engineered to express the anti-CEA CAR or calcitonin-reactive TCR led to complete tumor regression. T cells expressing the HLA-DPB1*04:01/02-restricted TCR targeting RET M918T, which was cloned from peripheral blood CD4+ T cells of a patient with MTC, demonstrated specific reactivity against cells pulsed with the mutated peptide and MTC tumor cells that expressed HLA-DPB1*04:01 and RET M918T. Conclusion: The preclinical data presented herein demonstrate the potential of using genetically engineered T cells targeting CEA, calcitonin, and/or RET M918T to treat metastatic MTC.

Transcriptomic profiles of neoantigen-reactive T cells in human gastrointestinal cancers.

Tumor-infiltrating neoantigen-reactive T cells can mediate regression of metastatic gastrointestinal cancers yet remain poorly characterized. We performed immunological screening against personalized neoantigens in combination with single-cell RNA sequencing on tumor-infiltrating lymphocytes from bile duct and pancreatic cancer patients to characterize the transcriptomic landscape of neoantigen-reactive T cells. We found that most neoantigen-reactive CD8+ T cells displayed an exhausted state with significant CXCL13 and GZMA co-expression compared with non-neoantigen-reactive bystander cells. Most neoantigen-reactive CD4+ T cells from a patient with bile duct cancer also exhibited an exhausted phenotype but with overexpression of HOPX or ADGRG1 while lacking IL7R expression. Thus, neoantigen-reactive T cells infiltrating gastrointestinal cancers harbor distinct transcriptomic signatures, which may provide new opportunities for harnessing these cells for therapy.

Whole genome sequencing of peach (Prunus persica L.) for SNP identification and selection.

BACKGROUND: The application of next generation sequencing technologies and bioinformatic scripts to identify high frequency SNPs distributed throughout the peach genome is described. Three peach genomes were sequenced using Roche 454 and Illumina/Solexa technologies to obtain long contigs for alignment to the draft 'Lovell' peach sequence as well as sufficient depth of coverage for 'in silico' SNP discovery. DESCRIPTION: The sequences were aligned to the 'Lovell' peach genome released April 01, 2010 by the International Peach Genome Initiative (IPGI). 'Dr. Davis', 'F8, 1-42' and 'Georgia Belle' were sequenced to add SNPs segregating in two breeding populations, Pop DF ('Dr. Davis' × 'F8, 1-42') and Pop DG ('Dr. Davis' × 'Georgia Belle'). Roche 454 sequencing produced 980,000 total reads with 236 Mb sequence for 'Dr. Davis' and 735,000 total reads with 172 Mb sequence for 'F8, 1-42'. 84 bp × 84 bp paired end Illumina/Solexa sequences yielded 25.5, 21.4, 25.5 million sequences for 'Dr. Davis', 'F8, 1-42' and 'Georgia Belle', respectively. BWA/SAMtools were used for alignment of raw reads and SNP detection, with custom PERL scripts for SNP filtering. Velvet's Columbus module was used for sequence assembly. Comparison of aligned and overlapping sequences from both Roche 454 and Illumina/Solexa resulted in the selection of 6654 high quality SNPs for 'Dr. Davis' vs. 'F8, 1-42' and 'Georgia Belle', distributed on eight major peach genome scaffolds as defined from the 'Lovell' assembly. CONCLUSION: The eight scaffolds contained about 215-225 Mb of peach genomic sequences with one SNP/~ 40,000 bases. All sequences from Roche 454 and Illumina/Solexa have been submitted to NCBI for public use in the Short Read Archive database. SNPs have been deposited in the NCBI SNP database.

Exploring the Transcriptome Landscape of Pomegranate Fruit Peel for Natural Product Biosynthetic Gene and SSR Marker Discovery(F).

Pomegranate fruit peel is rich in bioactive plant natural products, such as hydrolyzable tannins and anthocyanins. Despite their documented roles in human nutrition and fruit quality, genes involved in natural product biosynthesis have not been cloned from pomegranate and very little sequence information is available on pomegranate in the public domain. Shotgun transcriptome sequencing of pomegranate fruit peel cDNA was performed using RNA-Seq on the Illumina Genome Analyzer platform. Over 100 million raw sequence reads were obtained and assembled into 9,839 transcriptome assemblies (TAs) (>200 bp). Candidate genes for hydrolyzable tannin, anthocyanin, flavonoid, terpenoid and fatty acid biosynthesis and/or regulation were identified. Three lipid transfer proteins were obtained that may contribute to the previously reported IgE reactivity of pomegranate fruit extracts. In addition, 115 SSR markers were identified from the pomegranate fruit peel transcriptome and primers were designed for 77 SSR markers. The pomegranate fruit peel transcriptome set provides a valuable platform for natural product biosynthetic gene and SSR marker discovery in pomegranate. This work also demonstrates that next-generation transcriptome sequencing is an economical and effective approach for investigating natural product biosynthesis, identifying genes controlling important agronomic traits, and discovering molecular markers in non-model specialty crop species.

A label-free differential quantitative mass spectrometry method for the characterization and identification of protein changes during citrus fruit development.

BACKGROUND: Citrus is one of the most important and widely grown commodity fruit crops. In this study a label-free LC-MS/MS based shot-gun proteomics approach was taken to explore three main stages of citrus fruit development. These approaches were used to identify and evaluate changes occurring in juice sac cells in various metabolic pathways affecting citrus fruit development and quality. RESULTS: Protein changes in citrus juice sac cells were identified and quantified using label-free shotgun methodologies. Two alternative methods, differential mass-spectrometry (dMS) and spectral counting (SC) were used to analyze protein changes occurring during earlier and late stages of fruit development. Both methods were compared in order to develop a proteomics workflow that could be used in a non-model plant lacking a sequenced genome. In order to resolve the bioinformatics limitations of EST databases from species that lack a full sequenced genome, we established iCitrus. iCitrus is a comprehensive sequence database created by merging three major sources of sequences (HarvEST:citrus, NCBI/citrus/unigenes, NCBI/citrus/proteins) and improving the annotation of existing unigenes. iCitrus provided a useful bioinformatics tool for the high-throughput identification of citrus proteins. We have identified approximately 1500 citrus proteins expressed in fruit juice sac cells and quantified the changes of their expression during fruit development. Our results showed that both dMS and SC provided significant information on protein changes, with dMS providing a higher accuracy. CONCLUSION: Our data supports the notion of the complementary use of dMS and SC for label-free comparative proteomics, broadening the identification spectrum and strengthening the identification of trends in protein expression changes during the particular processes being compared.

MYBPC3 Haplotype Linked to Hypertrophic Cardiomyopathy in Rhesus Macaques (Macaca mulatta).

In humans, abnormal thickening of the left ventricle of the heart clinically defines hypertrophic cardiomyopathy (HCM), a common inherited cardiovascular disorder that can precede a sudden cardiac death event. The wide range of clinical presentations in HCM obscures genetic variants that may influence an individual's susceptibility to sudden cardiac death. Although exon sequencing of major sarcomere genes can be used to detect high-impact causal mutations, this strategy is successful in only half of patient cases. The incidence of left ventricular hypertrophy (LVH) in a managed research colony of rhesus macaques provides an excellent comparative model in which to explore the genomic etiology of severe HCM and sudden cardiac death. Because no rhesus HCM-associated mutations have been reported, we used a next-generation genotyping assay that targets 7 sarcomeric rhesus genes within 63 genomic sites that are orthologous to human genomic regions known to harbor HCM disease variants. Amplicon sequencing was performed on 52 macaques with confirmed LVH and 42 unrelated, unaffected animals representing both the Indian and Chinese rhesus macaque subspecies. Bias-reduced logistic regression uncovered a risk haplotype in the rhesus MYBPC3 gene, which is frequently disrupted in both human and feline HCM; this haplotype implicates an intronic variant strongly associated with disease in either homozygous or carrier form. Our results highlight that leveraging evolutionary genomic data provides a unique, practical strategy for minimizing population bias in complex disease studies.

Viral genome integration of canine papillomavirus 16.

Papillomaviruses infect humans and animals, most often causing benign proliferations on skin or mucosal surfaces. Rarely, these infections persist and progress to cancer. In humans, this transformation most often occurs with high-risk papillomaviruses, where viral integration is a critical event in carcinogenesis. The first aim of this study was to sequence the viral genome of canine papillomavirus (CPV) 16 from a pigmented viral plaque that progressed to metastatic squamous cell carcinoma in a dog. The second aim was to characterize multiple viral genomic deletions and translocations as well as host integration sites. The full viral genome was identified using a combination of PCR and high throughput sequencing. CPV16 is most closely related to chipapillomaviruses CPV4, CPV9, and CPV12 and we propose CPV16 be classified as a chipapillomavirus. Assembly of the full viral genome enabled identification of deletion of portions of the E1 and E2/E4 genes and two viral translocations within the squamous cell carcinoma. Genome walking was performed which identified four sites of viral integration into the host genome. This is the first description of integration of a canine papillomavirus into the host genome, raising the possibility that CPV16 may be a potential canine high-risk papillomavirus type.

Stochastic simulation of actin dynamics reveals the role of annealing and fragmentation.

Recent observations of F-actin dynamics call for theoretical models to interpret and understand the quantitative data. A number of existing models rely on simplifications and do not take into account F-actin fragmentation and annealing. We use Gillespie's algorithm for stochastic simulations of the F-actin dynamics including fragmentation and annealing. The simulations vividly illustrate that fragmentation and annealing have little influence on the shape of the polymerization curve and on nucleotide profiles within filaments but drastically affect the F-actin length distribution, making it exponential. We find that recent surprising measurements of high length diffusivity at the critical concentration cannot be explained by fragmentation and annealing events unless both fragmentation rates and frequency of undetected fragmentation and annealing events are greater than previously thought. The simulations compare well with experimentally measured actin polymerization data and lend additional support to a number of existing theoretical models.

Identifying rhesus macaque gene orthologs using heterospecific human CNV probes.

We used the Affymetrix(®) Genome-Wide Human SNP Array 6.0 to identify heterospecific markers and compare copy number and structural genomic variation between humans and rhesus macaques. Over 200,000 human copy number variation (CNV) probes were mapped to a Chinese and an Indian rhesus macaque sample. Observed genomic rearrangements and synteny were in agreement with the results of a previously published genomic comparison between humans and rhesus macaques. Comparisons between each of the two rhesus macaques and humans yielded 206 regions with copy numbers that differed by at least two fold in the Indian rhesus macaque and human, 32 in the Chinese rhesus macaque and human, and 147 in both rhesus macaques. The detailed genomic map and preliminary CNV data are useful for better understanding genetic variation in rhesus macaques, identifying derived changes in human CNVs that may have evolved by selection, and determining the suitability of rhesus macaques as human models for particular biomedical studies.

Complete genome sequence of canine papillomavirus type 16.

Papillomaviruses are epitheliotropic, circular, double-stranded DNA viruses within the family Papillomaviridae that are associated with benign and malignant tumors in humans and animals. We report the complete genome sequence of canine papillomavirus type 16 identified within multiple pigmented cutaneous plaques and squamous cell carcinoma from an intact female Basenji dog.