A non-coding genomic duplication at the HMX1 locus is associated with crop ears in highland cattle.

Highland cattle with congenital crop ears have notches of variable size on the tips of both ears. In some cases, cartilage deformation can be seen and occasionally the external ears are shortened. We collected 40 cases and 80 controls across Switzerland. Pedigree data analysis confirmed a monogenic autosomal dominant mode of inheritance with variable expressivity. All affected animals could be traced back to a single common ancestor. A genome-wide association study was performed and the causative mutation was mapped to a 4 Mb interval on bovine chromosome 6. The H6 family homeobox 1 (HMX1) gene was selected as a positional and functional candidate gene. By whole genome re-sequencing of an affected Highland cattle, we detected 6 non-synonymous coding sequence variants and two variants in an ultra-conserved element at the HMX1 locus with respect to the reference genome. Of these 8 variants, only a non-coding 76 bp genomic duplication (g.106720058_106720133dup) located in the conserved region was perfectly associated with crop ears. The identified copy number variation probably results in HMX1 misregulation and possible gain-of-function. Our findings confirm the role of HMX1 during the development of the external ear. As it is sometimes difficult to phenotypically diagnose Highland cattle with slight ear notches, genetic testing can now be used to improve selection against this undesired trait.

Genomic signatures of relaxed disruptive selection associated with speciation reversal in whitefish

Background: Speciation reversal: the erosion of species differentiation via an increase in introgressive hybridization due to the weakening of previously divergent selection regimes, is thought to be an important, yet poorly understood, driver of biodiversity loss. Our study system, the Alpine whitefish (Coregonus spp.) species complex is a classic example of a recent postglacial adaptive radiation: forming an array of endemic lake flocks, with the independent origination of similar ecotypes among flocks. However, many of the lakes of the Alpine radiation have been seriously impacted by anthropogenic nutrient enrichment, resulting in a collapse in neutral genetic and phenotypic differentiation within the most polluted lakes. Here we investigate the effects of eutrophication on the selective forces that have shaped this radiation, using population genomics. We studied eight sympatric species assemblages belonging to five independent parallel adaptive radiations, and one species pair in secondary contact. We used AFLP markers, and applied FST outlier (BAYESCAN, DFDIST) and logistic regression analyses (MATSAM), to identify candidate regions for disruptive selection in the genome and their associations with adaptive traits within each lake flock. The number of outlier and adaptive trait associated loci identified per lake were then regressed against two variables (historical phosphorus concentration and contemporary oxygen concentration) representing the strength of eutrophication. Results: Whilst we identify disruptive selection candidate regions in all lake flocks, we find similar trends, across analysis methods, towards fewer disruptive selection candidate regions and fewer adaptive trait/candidate loci associations in the more polluted lakes. Conclusions: Weakened disruptive selection and a concomitant breakdown in reproductive isolating mechanisms in more polluted lakes has lead to increased gene flow between coexisting Alpine whitefish species. We hypothesize that the resulting higher rates of interspecific recombination reduce either the number or extent of genomic islands of divergence surrounding loci evolving under disruptive natural selection. This produces the negative trend seen in the number of selection candidate loci recovered during genome scans of whitefish species flocks, with increasing levels of anthropogenic eutrophication: as the likelihood decreases that AFLP restriction sites will fall within regions of heightened genomic divergence and therefore be classified as FST outlier loci. This study explores for the first time the potential effects of human-mediated relaxation of disruptive selection on heterogeneous genomic divergence between coexisting species.

Sarcoid-derived fibroblasts: links between genomic instability, energy metabolism and senescence.

Bovine papillomavirus 1 (BPV-1) is a well recognized etiopathogenetic factor in a cancer-like state in horses, namely equine sarcoid disease. Nevertheless, little is known about BPV-1-mediated cell transforming effects. It was shown that BPV-1 triggers genomic instability through DNA hypomethylation and oxidative stress. In the present study, we further characterized BPV-1-positive fibroblasts derived from sarcoid tumors. The focus was on cancer-like features of sarcoid-derived fibroblasts, including cell cycle perturbation, comprehensive DNA damage analysis, end-replication problem, energy metabolism and oncogene-induced premature senescence. The S phase of the cell cycle, polyploidy events, DNA double strand breaks (DSBs) and DNA single strand breaks (SSBs) were increased in BPV-1-positive cells compared to control fibroblasts. BPV-1-mediated oxidative stress may contribute to telomere dysfunction in sarcoid-derived fibroblasts. Loss of mitochondrial membrane potential and concurrent elevation in intracellular ATP production may be a consequence of changes in energy-supplying pathways in BPV-1-positive cells which is also typical for cancer cells. Shifts in energy metabolism may support rapid proliferation in cells infected by BPV-1. Nevertheless, sarcoid-derived fibroblasts representing a heterogeneous cell fraction vary in some aspects of metabolic phenotype due to a dual role of BPV-1 in cell transformation and oncogene-induced premature senescence. This was shown with increased senescence-associated β-galactosidase (SA-β-gal) activity. Taken together, metabolic phenotypes in sarcoid-derived fibroblasts are plastic, which are similar to greater plasticity of cancer tissues than normal tissues.

Does Sex Speed Up Evolutionary Rate and Increase Biodiversity?

Most empirical and theoretical studies have shown that sex increases the rate of evolution, although evidence of sex constraining genomic and epigenetic variation and slowing down evolution also exists. Faster rates with sex have been attributed to new gene combinations, removal of deleterious mutations, and adaptation to heterogeneous environments. Slower rates with sex have been attributed to removal of major genetic rearrangements, the cost of finding a mate, vulnerability to predation, and exposure to sexually transmitted diseases. Whether sex speeds or slows evolution, the connection between reproductive mode, the evolutionary rate, and species diversity remains largely unexplored. Here we present a spatially explicit model of ecological and evolutionary dynamics based on DNA sequence change to study the connection between mutation, speciation, and the resulting biodiversity in sexual and asexual populations. We show that faster speciation can decrease the abundance of newly formed species and thus decrease long-term biodiversity. In this way, sex can reduce diversity relative to asexual populations, because it leads to a higher rate of production of new species, but with lower abundances. Our results show that reproductive mode and the mechanisms underlying it can alter the link between mutation, evolutionary rate, speciation and biodiversity and we suggest that a high rate of evolution may not be required to yield high biodiversity.

The Oxytricha trifallax macronuclear genome: a complex eukaryotic genome with 16,000 tiny chromosomes.

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor "silent" germline micronuclear genome by a process of "unscrambling" and fragmentation. The tiny macronuclear "nanochromosomes" typically encode single, protein-coding genes (a small portion, 10%, encode 2-8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.

New generation genomic platforms in investigation of complex diseases and BEN.

(Full text is available at http://www.manu.edu.mk/prilozi). New generation genomic platforms enable us to decipher the complex genetic basis of complex diseases and Balkan Endemic Nephropathy (BEN) at a high-throughput basis. They give valuable information about predisposing Single Nucleotide Polymorphisms (SNPs), Copy Number Variations (CNVs) or Loss of Heterozygosity (LOH) (using SNP-array) and about disease-causing mutations along the whole sequence of candidate-genes (using Next Generation Sequencing). This information could be used for screening of individuals in risk families and moving the main medicine stream to the prevention. They also might have an impact on more effective treatment. Here we discuss these genomic platforms and report some applications of SNP-array technology in a case with familial nephrotic syndrome. Key words: complex diseases, genome wide association studies, SNP, genomic arrays, next generation sequ-encing.

The Draft Assembly of the Radically Organized Stylonychia lemnae Macronuclear Genome.

Stylonychia lemnae is a classical model single-celled eukaryote, and a quintessential ciliate typified by dimorphic nuclei: A small, germline micronucleus and a massive, vegetative macronucleus. The genome within Stylonychia's macronucleus has a very unusual architecture, comprised variably and highly amplified "nanochromosomes," each usually encoding a single gene with a minimal amount of surrounding noncoding DNA. As only a tiny fraction of the Stylonychia genes has been sequenced, and to promote research using this organism, we sequenced its macronuclear genome. We report the analysis of the 50.2-Mb draft S. lemnae macronuclear genome assembly, containing in excess of 16,000 complete nanochromosomes, assembled as less than 20,000 contigs. We found considerable conservation of fundamental genomic properties between S. lemnae and its close relative, Oxytricha trifallax, including nanochromosomal gene synteny, alternative fragmentation, and copy number. Protein domain searches in Stylonychia revealed two new telomere-binding protein homologs and the presence of linker histones. Among the diverse histone variants of S. lemnae and O. trifallax, we found divergent, coexpressed variants corresponding to four of the five core nucleosomal proteins (H1.2, H2A.6, H2B.4, and H3.7) suggesting that these ciliates may possess specialized nucleosomes involved in genome processing during nuclear differentiation. The assembly of the S. lemnae macronuclear genome demonstrates that largely complete, well-assembled highly fragmented genomes of similar size and complexity may be produced from one library and lane of Illumina HiSeq 2000 shotgun sequencing. The provision of the S. lemnae macronuclear genome sets the stage for future detailed experimental studies of chromatin-mediated, RNA-guided developmental genome rearrangements.

Genomic cloning and expression of three murine UDP-galactose: beta-N-acetylglucosamine beta1,3-galactosyltransferase genes.

Based on the detection of expressed sequence tags that are similar to known galactosyltransferase sequences, we have isolated three novel UDP-galactose:beta-N-acetylglucosamine beta1, 3-galactosyltransferase (beta3GalT) genes from a mouse genomic library. The three genes, named beta3GalT-I, -II, and -III, encode type II transmembrane proteins of 326, 422, and 331 amino acids, respectively. The three proteins constitute a distinct subfamily as they do not share any sequence identity with other eucaryotic galactosyltransferases. Also, the entire protein-coding region of the three beta3GalT genes was contained in a single exon, which contrasts with the genomic organization of the beta1,4- and alpha1, 3-galactosyltransferase genes. The three beta3GalT genes were mainly expressed in brain tissue. The expression of the full-length murine genes as recombinant baculoviruses in insect cells revealed that the beta3GalT enzymes share the same acceptor specificity for beta-linked GlcNAc, although they differ in their Km for this acceptor and the donor UDP-Gal. The identification of beta3GalT genes emphasizes the structural diversity present in the galactosyltransferase gene family.

Tumour genomic and microenvironmental heterogeneity for integrated prediction of 5-year biochemical recurrence of prostate cancer: a retrospective cohort study.

BACKGROUND Clinical prognostic groupings for localised prostate cancers are imprecise, with 30-50% of patients recurring after image-guided radiotherapy or radical prostatectomy. We aimed to test combined genomic and microenvironmental indices in prostate cancer to improve risk stratification and complement clinical prognostic factors. METHODS We used DNA-based indices alone or in combination with intra-prostatic hypoxia measurements to develop four prognostic indices in 126 low-risk to intermediate-risk patients (Toronto cohort) who will receive image-guided radiotherapy. We validated these indices in two independent cohorts of 154 (Memorial Sloan Kettering Cancer Center cohort [MSKCC] cohort) and 117 (Cambridge cohort) radical prostatectomy specimens from low-risk to high-risk patients. We applied unsupervised and supervised machine learning techniques to the copy-number profiles of 126 pre-image-guided radiotherapy diagnostic biopsies to develop prognostic signatures. Our primary endpoint was the development of a set of prognostic measures capable of stratifying patients for risk of biochemical relapse 5 years after primary treatment. FINDINGS Biochemical relapse was associated with indices of tumour hypoxia, genomic instability, and genomic subtypes based on multivariate analyses. We identified four genomic subtypes for prostate cancer, which had different 5-year biochemical relapse-free survival. Genomic instability is prognostic for relapse in both image-guided radiotherapy (multivariate analysis hazard ratio [HR] 4·5 [95% CI 2·1-9·8]; p=0·00013; area under the receiver operator curve [AUC] 0·70 [95% CI 0·65-0·76]) and radical prostatectomy (4·0 [1·6-9·7]; p=0·0024; AUC 0·57 [0·52-0·61]) patients with prostate cancer, and its effect is magnified by intratumoral hypoxia (3·8 [1·2-12]; p=0·019; AUC 0·67 [0·61-0·73]). A novel 100-loci DNA signature accurately classified treatment outcome in the MSKCC low-risk to intermediate-risk cohort (multivariate analysis HR 6·1 [95% CI 2·0-19]; p=0·0015; AUC 0·74 [95% CI 0·65-0·83]). In the independent MSKCC and Cambridge cohorts, this signature identified low-risk to high-risk patients who were most likely to fail treatment within 18 months (combined cohorts multivariate analysis HR 2·9 [95% CI 1·4-6·0]; p=0·0039; AUC 0·68 [95% CI 0·63-0·73]), and was better at predicting biochemical relapse than 23 previously published RNA signatures. INTERPRETATION This is the first study of cancer outcome to integrate DNA-based and microenvironment-based failure indices to predict patient outcome. Patients exhibiting these aggressive features after biopsy should be entered into treatment intensification trials. FUNDING Movember Foundation, Prostate Cancer Canada, Ontario Institute for Cancer Research, Canadian Institute for Health Research, NIHR Cambridge Biomedical Research Centre, The University of Cambridge, Cancer Research UK, Cambridge Cancer Charity, Prostate Cancer UK, Hutchison Whampoa Limited, Terry Fox Research Institute, Princess Margaret Cancer Centre Foundation, PMH-Radiation Medicine Program Academic Enrichment Fund, Motorcycle Ride for Dad (Durham), Canadian Cancer Society.

Genomic heterogeneity of small ruminant lentiviruses detected by PCR

In order to detect a large spectrum of small ruminant lentiviruses, primers for PCR were chosen in conserved parts of the LTR and GAG genes of Icelandic Visna virus 1514 and of the POL gene of caprine arthritis-encephalitis virus. This set of primers was tested in six different caprine arthritis-encephalitis virus (CAEV)- and Maedi-Visna virus isolates of Dutch, American and Swiss origin. The LTR primers allowed the detection of the corresponding fragments of all isolates. The GAG primers allowed amplification of the corresponding fragments of all but the Swiss Maedi-Visna virus strain OLV. Using the POL primers, one Maedi-Visna- and two caprine arthritis-encephalitis virus strains were detected after one round of amplification. Sequencing of the GAG and POL amplification products and comparison to Icelandic Visna virus and CAEV strain CO revealed total heterogeneity of 38% for the GAG- and 28% for the POL fragment. The virus strains studied fall into two groups which are more closely related to one another than to Icelandic Visna virus.

A genomic perspective on a new bacterial genus and species from the Alcaligenaceae family, Basilea psittacipulmonis

BACKGROUND A novel Gram-negative, non-haemolytic, non-motile, rod-shaped bacterium was discovered in the lungs of a dead parakeet (Melopsittacus undulatus) that was kept in captivity in a petshop in Basel, Switzerland. The organism is described with a chemotaxonomic profile and the nearly complete genome sequence obtained through the assembly of short sequence reads. RESULTS Genome sequence analysis and characterization of respiratory quinones, fatty acids, polar lipids, and biochemical phenotype is presented here. Comparison of gene sequences revealed that the most similar species is Pelistega europaea, with BLAST identities of only 93% to the 16S rDNA gene, 76% identity to the rpoB gene, and a similar GC content (~43%) as the organism isolated from the parakeet, DSM 24701 (40%). The closest full genome sequences are those of Bordetella spp. and Taylorella spp. High-throughput sequencing reads from the Illumina-Solexa platform were assembled with the Edena de novo assembler to form 195 contigs comprising the ~2 Mb genome. Genome annotation with RAST, construction of phylogenetic trees with the 16S rDNA (rrs) gene sequence and the rpoB gene, and phylogenetic placement using other highly conserved marker genes with ML Tree all suggest that the bacterial species belongs to the Alcaligenaceae family. Analysis of samples from cages with healthy parakeets suggested that the newly discovered bacterial species is not widespread in parakeet living quarters. CONCLUSIONS Classification of this organism in the current taxonomy system requires the formation of a new genus and species. We designate the new genus Basilea and the new species psittacipulmonis. The type strain of Basilea psittacipulmonis is DSM 24701 (= CIP 110308 T, 16S rDNA gene sequence Genbank accession number JX412111 and GI 406042063).

An arranged marriage for precision medicine: hypoxia and genomic assays in localized prostate cancer radiotherapy

Prostate cancer (CaP) is the most commonly diagnosed malignancy in males in the Western world with one in six males diagnosed in their lifetime. Current clinical prognostication groupings use pathologic Gleason score, pre-treatment prostatic-specific antigen and Union for International Cancer Control-TNM staging to place patients with localized CaP into low-, intermediate- and high-risk categories. These categories represent an increasing risk of biochemical failure and CaP-specific mortality rates, they also reflect the need for increasing treatment intensity and justification for increased side effects. In this article, we point out that 30-50% of patients will still fail image-guided radiotherapy or surgery despite the judicious use of clinical risk categories owing to interpatient heterogeneity in treatment response. To improve treatment individualization, better predictors of prognosis and radiotherapy treatment response are needed to triage patients to bespoke and intensified CaP treatment protocols. These should include the use of pre-treatment genomic tests based on DNA or RNA indices and/or assays that reflect cancer metabolism, such as hypoxia assays, to define patient-specific CaP progression and aggression. More importantly, it is argued that these novel prognostic assays could be even more useful if combined together to drive forward precision cancer medicine for localized CaP.

DNA hypomethylation and oxidative stress-mediated increase in genomic instability in equine sarcoid-derived fibroblasts

It is widely accepted that equine sarcoid disease, the most common skin associated neoplasm in equids, is induced by bovine papillomavirus (BPV-1). Although BPV-1 DNA has been found in almost all examined sarcoids so far, its detailed impact on the horse's host cell metabolism is largely unknown. We used equine fibroblast cell lines originating from sarcoid biopsies to study BPV-1-associated changes on DNA methylation status and oxidative stress parameters. Sarcoid-derived fibroblasts manifested increased proliferation in vitro, transcriptional rDNA activity (NORs expression) and DNA hypomethylation compared to control cells. Cells isolated from equine sarcoids suffered from oxidative stress: the expression of antioxidant enzymes was decreased and the superoxide production was increased. Moreover, increased ploidy, oxidative DNA damage and micronuclei formation was monitored in sarcoid cells. We postulate that both altered DNA methylation status and redox milieu may affect genomic stability in BPV-1-infected cells and in turn contribute to sarcoid pathology.