Whole genome resequencing of black Angus and Holstein cattle for SNP and CNV discovery

Background: One of the goals of livestock genomics research is to identify the genetic differences responsible for variation in phenotypic traits, particularly those of economic importance. Characterizing the genetic variation in livestock species is an important step towards linking genes or genomic regions with phenotypes. The completion of the bovine genome sequence and recent advances in DNA sequencing technology allow for in-depth characterization of the genetic variations present in cattle. Here we describe the whole-genome resequencing of two Bos taurus bulls from distinct breeds for the purpose of identifying and annotating novel forms of genetic variation in cattle.Results: The genomes of a Black Angus bull and a Holstein bull were sequenced to 22-fold and 19-fold coverage, respectively, using the ABI SOLiD system. Comparisons of the sequences with the Btau4.0 reference assembly yielded 7 million single nucleotide polymorphisms (SNPs), 24% of which were identified in both animals. Of the total SNPs found in Holstein, Black Angus, and in both animals, 81%, 81%, and 75% respectively are novel. In-depth annotations of the data identified more than 16 thousand distinct non-synonymous SNPs (85% novel) between the two datasets. Alignments between the SNP-altered proteins and orthologues from numerous species indicate that many of the SNPs alter well-conserved amino acids. Several SNPs predicted to create or remove stop codons were also found. A comparison between the sequencing SNPs and genotyping results from the BovineHD high-density genotyping chip indicates a detection rate of 91% for homozygous SNPs and 81% for heterozygous SNPs. The false positive rate is estimated to be about 2% for both the Black Angus and Holstein SNP sets, based on follow-up genotyping of 422 and 427 SNPs, respectively. Comparisons of read depth between the two bulls along the reference assembly identified 790 putative copy-number variations (CNVs). Ten randomly selected CNVs, five genic and five non-genic, were successfully validated using quantitative real-time PCR. The CNVs are enriched for immune system genes and include genes that may contribute to lactation capacity. The majority of the CNVs (69%) were detected as regions with higher abundance in the Holstein bull.Conclusions: Substantial genetic differences exist between the Black Angus and Holstein animals sequenced in this work and the Hereford reference sequence, and some of this variation is predicted to affect evolutionarily conserved amino acids or gene copy number. The deeply annotated SNPs and CNVs identified in this resequencing study can serve as useful genetic tools, and as candidates in searches for phenotype-altering DNA differences

A Comprehensive DNA Barcode Library for the Looper Moths (Lepidoptera: Geometridae) of British Columbia, Canada

The construction of comprehensive reference libraries is essential to foster the development of DNA barcoding as a tool for monitoring biodiversity and detecting invasive species. The looper moths of British Columbia (BC), Canada present a challenging case for species discrimination via DNA barcoding due to their considerable diversity and limited taxonomic maturity.By analyzing specimens held in national and regional natural history collections, we assemble barcode records from representatives of 400 species from BC and surrounding provinces, territories and states. Sequence variation in the barcode region unambiguously discriminates over 93% of these 400 geometrid species. However, a final estimate of resolution success awaits detailed taxonomic analysis of 48 species where patterns of barcode variation suggest cases of cryptic species, unrecognized synonymy as well as young species.A catalog of these taxa meriting further taxonomic investigation is presented as well as the supplemental information needed to facilitate these investigations

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition

The genomes of two key bumblebee species with primitive eusocial organization

Background: The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats. Results: We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits. Conclusions: These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation

Global gene expression analysis of canine osteosarcoma stem cells reveals a novel role for COX-2 in tumour initiation

Osteosarcoma is the most common primary bone tumour of both children and dogs. It is an aggressive tumour in both species with a rapid clinical course leading ultimately to metastasis. In dogs and children distant metastasis occurs in >80% of individuals treated by surgery alone. Both canine and human osteosarcoma has been shown to contain a sub-population of cancer stem cells (CSCs), which may drive tumour growth, recurrence and metastasis, suggesting that naturally occurring canine osteosarcoma could act as a preclinical model for the human disease. Here we report the successful isolation of CSCs from primary canine osteosarcoma, as well as established cell lines. We show that these cells can form tumourspheres, and demonstrate relative resistance to chemotherapy. We demonstrate similar results for the human osteosarcma cell lines, U2OS and SAOS2. Utilizing the Affymetrix canine microarray, we are able to definitively show that there are significant differences in global gene expression profiles of isolated osteosarcoma stem cells and the daughter adherent cells. We identified 13,221 significant differences (p = 0.05), and significantly, COX-2 was expressed 141-fold more in CSC spheres than daughter adherent cells. To study the role of COX-2 expression in CSCs we utilized the COX-2 inhibitors meloxicam and mavacoxib. We found that COX-2 inhibition had no effect on CSC growth, or resistance to chemotherapy. However inhibition of COX-2 in daughter cells prevented sphere formation, indicating a potential significant role for COX-2 in tumour initiation

Inhibition of Notch pathway prevents osteosarcoma growth by cell cycle regulation

The study shows constitutive activation of the Notch pathway in various types of malignancies. However, it remains unclear how the Notch pathway is involved in the pathogenesis of osteosarcoma. We investigated the expression of the Notch pathway molecules in osteosarcoma biopsy specimens and examined the effect of Notch pathway inhibition. Real-time PCR revealed overexpression of Notch2, Jagged1, HEY1, and HEY2. On the other hand, Notch1 and DLL1 were downregulated in biopsy specimens. Notch pathway inhibition using γ-secretase inhibitor and CBF1 siRNA slowed the growth of osteosarcomas in vitro. In addition, γ-secretase inhibitor-treated xenograft models exhibited significantly slower osteosarcoma growth. Cell cycle analysis revealed that γ-secretase inhibitor promoted G1 arrest. Real-time PCR and western blot revealed that γ-secretase inhibitor reduced the expression of accelerators of the cell cycle, including cyclin D1, cyclin E1, E2, and SKP2. On the other hand, p21cip1 protein, a cell cycle suppressor, was upregulated by γ-secretase inhibitor treatment. These findings suggest that inhibition of Notch pathway suppresses osteosarcoma growth by regulation of cell cycle regulator expression and that the inactivation of the Notch pathway may be a useful approach to the treatment of patients with osteosarcoma

Genetics of Sputum Gene Expression in Chronic Obstructive Pulmonary Disease

Previous expression quantitative trait loci (eQTL) studies have performed genetic association studies for gene expression, but most of these studies examined lymphoblastoid cell lines from non-diseased individuals. We examined the genetics of gene expression in a relevant disease tissue from chronic obstructive pulmonary disease (COPD) patients to identify functional effects of known susceptibility genes and to find novel disease genes. By combining gene expression profiling on induced sputum samples from 131 COPD cases from the ECLIPSE Study with genomewide single nucleotide polymorphism (SNP) data, we found 4315 significant cis-eQTL SNP-probe set associations (3309 unique SNPs). The 3309 SNPs were tested for association with COPD in a genomewide association study (GWAS) dataset, which included 2940 COPD cases and 1380 controls. Adjusting for 3309 tests (p<1.5e-5), the two SNPs which were significantly associated with COPD were located in two separate genes in a known COPD locus on chromosome 15: CHRNA5 and IREB2. Detailed analysis of chromosome 15 demonstrated additional eQTLs for IREB2 mapping to that gene. eQTL SNPs for CHRNA5 mapped to multiple linkage disequilibrium (LD) bins. The eQTLs for IREB2 and CHRNA5 were not in LD. Seventy-four additional eQTL SNPs were associated with COPD at p<0.01. These were genotyped in two COPD populations, finding replicated associations with a SNP in PSORS1C1, in the HLA-C region on chromosome 6. Integrative analysis of GWAS and gene expression data from relevant tissue from diseased subjects has located potential functional variants in two known COPD genes and has identified a novel COPD susceptibility locus

Rapid Selection and Proliferation of CD133(+) Cells from Cancer Cell Lines: Chemotherapeutic Implications

Cancer stem cells (CSCs) are considered a subset of the bulk tumor responsible for initiating and maintaining the disease. Several surface cellular markers have been recently used to identify CSCs. Among those is CD133, which is expressed by hematopoietic progenitor cells as well as embryonic stem cells and various cancers. We have recently isolated and cultured CD133 positive [CD133(+)] cells from various cancer cell lines using a NASA developed Hydrodynamic Focusing Bioreactor (HFB) (Celdyne, Houston, TX). For comparison, another bioreactor, the rotary cell culture system (RCCS) manufactured by Synthecon (Houston, TX) was used. Both the HFB and the RCCS bioreactors simulate aspects of hypogravity. In our study, the HFB increased CD133(+) cell growth from various cell lines compared to the RCCS vessel and to normal gravity control. We observed a (+)15-fold proliferation of the CD133(+) cellular fraction with cancer cells that were cultured for 7-days at optimized conditions. The RCCS vessel instead yielded a (−)4.8-fold decrease in the CD133(+)cellular fraction respect to the HFB after 7-days of culture. Interestingly, we also found that the hypogravity environment of the HFB greatly sensitized the CD133(+) cancer cells, which are normally resistant to chemo treatment, to become susceptible to various chemotherapeutic agents, paving the way to less toxic and more effective chemotherapeutic treatment in patients. To be able to test the efficacy of cytotoxic agents in vitro prior to their use in clinical setting on cancer cells as well as on cancer stem cells may pave the way to more effective chemotherapeutic strategies in patients. This could be an important advancement in the therapeutic options of oncologic patients, allowing for more targeted and personalized chemotherapy regimens as well as for higher response rates

The S phase checkpoint promotes the Smc5/6 complex dependent SUMOylation of Pol2, the catalytic subunit of DNA polymerase ε

Replication fork stalling and accumulation of single-stranded DNA trigger the S phase checkpoint, a signalling cascade that, in budding yeast, leads to the activation of the Rad53 kinase. Rad53 is essential in maintaining cell viability, but its targets of regulation are still partially unknown. Here we show that Rad53 drives the hyper-SUMOylation of Pol2, the catalytic subunit of DNA polymerase ε, principally following replication forks stalling induced by nucleotide depletion. Pol2 is the main target of SUMOylation within the replisome and its modification requires the SUMO-ligase Mms21, a subunit of the Smc5/6 complex. Moreover, the Smc5/6 complex co-purifies with Pol ε, independently of other replisome components. Finally, we map Pol2 SUMOylation to a single site within the N-terminal catalytic domain and identify a SUMO-interacting motif at the C-terminus of Pol2. These data suggest that the S phase checkpoint regulate Pol ε during replication stress through Pol2 SUMOylation and SUMO-binding abilit

Use of different RT-QuIC substrates for detecting CWD prions in the brain of Norwegian cervids

Chronic wasting disease (CWD) is a highly contagious prion disease affecting captive and free-ranging cervid populations. CWD has been detected in United States, Canada, South Korea and, most recently, in Europe (Norway, Finland and Sweden). Animals with CWD release infectious prions in the environment through saliva, urine and feces sustaining disease spreading between cervids but also potentially to other non-cervids ruminants (e.g. sheep, goats and cattle). In the light of these considerations and due to CWD unknown zoonotic potential, it is of utmost importance to follow specific surveillance programs useful to minimize disease spreading and transmission. The European community has already in place specific surveillance measures, but the traditional diagnostic tests performed on nervous or lymphoid tissues lack sensitivity. We have optimized a Real-Time Quaking-Induced Conversion (RT-QuIC) assay for detecting CWD prions with high sensitivity and specificity to try to overcome this problem. In this work, we show that bank vole prion protein (PrP) is an excellent substrate for RT-QuIC reactions, enabling the detection of trace-amounts of CWD prions, regardless of prion strain and cervid species. Beside supporting the traditional diagnostic tests, this technology could be exploited for detecting prions in peripheral tissues from live animals, possibly even at preclinical stages of the disease