Palaeoproteomics identifies beaver fur in Danish high-status Viking Age burials - direct evidence of fur trade.

Fur is known from contemporary written sources to have been a key commodity in the Viking Age. Nevertheless, the fur trade has been notoriously difficult to study archaeologically as fur rarely survives in the archaeological record. In Denmark, fur finds are rare and fur in clothing has been limited to a few reports and not recorded systematically. We were therefore given access to fur from six Danish high status graves dated to the Viking Age. The fur was analysed by aDNA and palaeoproteomics methods to identify the species of origin in order to explore the Viking Age fur trade. Endogenous aDNA was not recovered, but fur proteins (keratins) were analysed by MALDI-TOF-MS and LC-MS/MS. We show that Viking Age skin clothing were often composites of several species, showing highly developed manufacturing and material knowledge. For example, fur was produced from wild animals while leather was made of domesticates. Several examples of beaver fur were identified, a species which is not native to Denmark, and therefore indicative of trade. We argue that beaver fur was a luxury commodity, limited to the elite and worn as an easily recognisable indicator of social status.

Differences in gene expression despite identical histomorphology in sinonasal intestinal-type adenocarcinoma and metastases from colorectal adenocarcinoma.

Sinonasal intestinal-type adenocarcinoma (sITAC) is histomorphologically indistinguishable from colorectal adenocarcinoma (CRC) leading to diagnostic challenges. Metastases from CRCs to the sinonasal tract have been reported. The aim of the study was to identify a biomarker making it possible to distinguish between sITAC and metastases of colorectal origin. Formalin-fixated paraffin-embedded (FFPE) tissue from 20 consecutive patients with sITAC treated at Rigshospitalet, Denmark from 2005 to 2017, 20 patients with CRC, and second patients with both sinonasal and colorectal carcinomas were included, and RNA-sequencing was performed on all samples. Moreover, a series of 26 samples from metastasizing CRC were included (in-house data). 3139 differentially expressed genes were identified, of these several were deemed as possible biomarkers, including CSDE1, for which immunohistochemical staining was performed. sITAC and CRC differ in genomic expression. CSDE1, previously found upregulated in CRC, was significantly differentially expressed. Using immunohistochemical staining, no sITACs displayed strong and diffuse staining for CSDE1, which represents a potential marker to use in distinguishing sITAC from a metastasis of colorectal origin. This knowledge could improve the diagnostic process and hopefully the outcome in patients with this rare tumor.

Viral and Genomic Drivers of Squamous Cell Neoplasms Arising in the Lacrimal Drainage System.

The pathogenesis of squamous cell neoplasms arising in the lacrimal drainage system is poorly understood, and the underlying genomic drivers for disease development remain unexplored. We aimed to investigate the genomic aberrations in carcinomas arising in the LDS and correlate the findings to human papillomavirus (HPV) status. The HPV analysis was performed using HPV DNA PCR, HPV E6/E7 mRNA in-situ hybridization, and p16 immunohistochemistry. The genomic characterization was performed by targeted DNA sequencing of 523 cancer-relevant genes. Patients with LDS papilloma (n = 17) and LDS carcinoma (n = 15) were included. There was a male predominance (68%) and a median age at diagnosis of 46.0 years (range 27.5-65.5 years) in patients with papilloma and 63.8 years (range 34.0-87.2 years) in patients with carcinoma. Transcriptional activity of the HPV E6/E7 oncogenes was detected in the whole tumor thickness in 12/15 (80%) papillomas (HPV6, 11, 16) and 10/15 (67%) squamous cell carcinomas (SCC) (HPV11: 3/15 (20%) and HPV16: 7/15 (47%)). Pathogenic variants in PIK3CA, FGFR3, AKT1, and PIK3R1, wildtype TP53, p16 overexpression, and deregulated high-risk E6/E7 transcription characterized the HPV16-positive SCC. The deregulated pattern of HPV E6/E7 expression, correlating with HPV DNA presence and p16 positivity, supports a causal role of HPV in a subset of LDS papillomas and carcinomas. The viral and molecular profile of LDS SCC resembles that of other HPV-driven SCC.

A Comparison of Tools for Copy-Number Variation Detection in Germline Whole Exome and Whole Genome Sequencing Data.

Copy-number variations (CNVs) have important clinical implications for several diseases and cancers. Relevant CNVs are hard to detect because common structural variations define large parts of the human genome. CNV calling from short-read sequencing would allow single protocol full genomic profiling. We reviewed 50 popular CNV calling tools and included 11 tools for benchmarking in a reference cohort encompassing 39 whole genome sequencing (WGS) samples paired current clinical standard-SNP-array based CNV calling. Additionally, for nine samples we also performed whole exome sequencing (WES), to address the effect of sequencing protocol on CNV calling. Furthermore, we included Gold Standard reference sample NA12878, and tested 12 samples with CNVs confirmed by multiplex ligation-dependent probe amplification (MLPA). Tool performance varied greatly in the number of called CNVs and bias for CNV lengths. Some tools had near-perfect recall of CNVs from arrays for some samples, but poor precision. Several tools had better performance for NA12878, which could be a result of overfitting. We suggest combining the best tools also based on different methodologies: GATK gCNV, Lumpy, DELLY, and cn.MOPS. Reducing the total number of called variants could potentially be assisted by the use of background panels for filtering of frequently called variants.

Genomic Alterations in Human Papillomavirus-Positive and -Negative Conjunctival Squamous Cell Carcinomas.

Purpose: The genomic alterations contributing to the pathogenesis of conjunctival squamous cell carcinomas (SCCs) and their precursor lesions are poorly understood and hamper our ability to develop molecular therapies to reduce the recurrence rates and treatment-related morbidities of this disease. We aimed to characterize the somatic DNA alterations in human papillomavirus (HPV)-positive and HPV-negative conjunctival SCC. Methods: Patients diagnosed with conjunctival SCC in situ or SCC treated in ocular oncology referral centers in Denmark were included. HPV detection (HPV DNA PCR, p16 immunohistochemistry, and mRNA in situ hybridization) and targeted capture-based next-generation sequencing of 523 genes frequently involved in cancer were performed to describe the mutational profile based on HPV status. Results: Tumor tissue was available in 33 cases (n = 8 conjunctival SCCs in situ, n = 25 conjunctival SCCs), constituting 25 male and 8 female patients. Nine cases were HPV positive. The HPV-positive SCCs in situ and SCCs were characterized by transcriptionally active high-risk HPV (types 16 and 39) within the tumor cells, frequent mutations in PIK3CA (n = 5/9), and wild-type TP53, CDKN2A, and RB1, while the HPV-negative counterparts harbored frequent mutations in TP53 (n = 21/24), CDKN2A (n = 7/24), and RB1 (n = 6/24). Conclusions: Our findings have delineated two potentially distinct distributions of somatic mutations in conjunctival SCC based on HPV status-pointing to different biological mechanisms of carcinogenesis. The present findings support a causal role of HPV in a subset of conjunctival SCC.

Major driver mutations are shared between sinonasal intestinal-type adenocarcinoma and the morphologically identical colorectal adenocarcinoma.

PURPOSE: The purpose of our study was to compare genomic changes in sinonasal intestinal-type adenocarcinoma (sITAC) and colorectal adenocarcinoma (CRC), as they are histomorphologically indistinguishable. This can cause diagnostic difficulties as sinonasal tumours initially diagnosed as sITAC may represent metastasis from CRC, a frequent cancer. Previous studies have not uncovered the underlying mechanism behind the histomorphological resemblance. METHODS/PATIENTS: Tissue samples from all consecutive patients with sITAC at our facility (20 patients) were compared to samples from 20 patients with CRC as well as samples from 2 patients with both CRC and sinonasal tumours. DNA sequencing was performed using Illumina TruSight Oncology 500 panel consisting of 523 cancer-associated genes. Frequent mutations were inspected manually using the Integrative Genomics Viewer. RESULTS: Several well-known cancer-associated genes were mutated in the CRC group, but also in the sinonasal ITAC group. These genes included APC mutated in 65% of the CRC group and 37% of the sinonasal ITAC group, and TP53 mutated in 65% of CRC samples and 58% of ITAC samples. These shared mutations may explain the histomorphological similarities. Successful DNA sequencing was performed on the colorectal sample from one of the two patients with both CRC and sinonasal tumour. Comparing mutations in these samples from one patient we have shown that the sinonasal tumour in all probability was a CRC metastasis. CONCLUSION: We have identified several genetic similarities between sITAC and CRC. This discovery brings us closer to understanding mechanisms behind the development of sITAC-and hopefully in the future targeted therapy.

A genomic exploration of the early evolution of extant cats and their sabre-toothed relatives.

Background: The evolutionary relationships of Felidae during their Early-Middle Miocene radiation is contentious. Although the early common ancestors have been subsumed under the grade-group Pseudaelurus, this group is thought to be paraphyletic, including the early ancestors of both modern cats and extinct sabretooths. Methods: Here, we sequenced a draft nuclear genome of Smilodon populator, dated to 13,182 ± 90 cal BP, making this the oldest palaeogenome from South America to date, a region known to be problematic for ancient DNA preservation. We analysed this genome, together with genomes from other extinct and extant cats to investigate their phylogenetic relationships. Results: We confirm a deep divergence (~20.65 Ma) within sabre-toothed cats. Through the analysis of both simulated and empirical data, we show a lack of gene flow between Smilodon and contemporary Felidae. Conclusions: Given that some species traditionally assigned to Pseudaelurus originated in the Early Miocene ~20 Ma, this indicates that some species of Pseudaelurus may be younger than the lineages they purportedly gave rise to, further supporting the hypothesis that Pseudaelurus was paraphyletic.

Genomic Adaptations and Evolutionary History of the Extinct Scimitar-Toothed Cat, Homotherium latidens.

Homotherium was a genus of large-bodied scimitar-toothed cats, morphologically distinct from any extant felid species, that went extinct at the end of the Pleistocene [1-4]. They possessed large, saber-form serrated canine teeth, powerful forelimbs, a sloping back, and an enlarged optic bulb, all of which were key characteristics for predation on Pleistocene megafauna [5]. Previous mitochondrial DNA phylogenies suggested that it was a highly divergent sister lineage to all extant cat species [6-8]. However, mitochondrial phylogenies can be misled by hybridization [9], incomplete lineage sorting (ILS), or sex-biased dispersal patterns [10], which might be especially relevant for Homotherium since widespread mito-nuclear discrepancies have been uncovered in modern cats [10]. To examine the evolutionary history of Homotherium, we generated a ∼7x nuclear genome and a ∼38x exome from H. latidens using shotgun and target-capture sequencing approaches. Phylogenetic analyses reveal Homotherium as highly divergent (∼22.5 Ma) from living cat species, with no detectable signs of gene flow. Comparative genomic analyses found signatures of positive selection in several genes, including those involved in vision, cognitive function, and energy consumption, putatively consistent with diurnal activity, well-developed social behavior, and cursorial hunting [5]. Finally, we uncover relatively high levels of genetic diversity, suggesting that Homotherium may have been more abundant than the limited fossil record suggests [3, 4, 11-14]. Our findings complement and extend previous inferences from both the fossil record and initial molecular studies, enhancing our understanding of the evolution and ecology of this remarkable lineage.

Darwin's Fancy Revised: An Updated Understanding of the Genomic Constitution of Pigeon Breeds.

Through its long history of artificial selection, the rock pigeon (Columba livia Gmelin 1789) was forged into a large number of domestic breeds. The incredible amount of phenotypic diversity exhibited in these breeds has long held the fascination of scholars, particularly those interested in biological inheritance and evolution. However, exploiting them as a model system is challenging, as unlike with many other domestic species, few reliable records exist about the origins of, and relationships between, each of the breeds. Therefore, in order to broaden our understanding of the complex evolutionary relationships among pigeon breeds, we generated genome-wide data by performing the genotyping-by-sequencing (GBS) method on close to 200 domestic individuals representing over 60 breeds. We analyzed these GBS data alongside previously published whole-genome sequencing data, and this combined analysis allowed us to conduct the most extensive phylogenetic analysis of the group, including two feral pigeons and one outgroup. We improve previous phylogenies, find considerable population structure across the different breeds, and identify unreported interbreed admixture events. Despite the reduced number of loci relative to whole-genome sequencing, we demonstrate that GBS data provide sufficient analytical power to investigate intertwined evolutionary relationships, such as those that are characteristic of animal domestic breeds. Thus, we argue that future studies should consider sequencing methods akin to the GBS approach as an optimal cost-effective approach for addressing complex phylogenies.

Identification of transcription factor genes involved in anthocyanin biosynthesis in carrot (Daucus carota L.) using RNA-Seq.

BACKGROUND: Anthocyanins are water-soluble colored flavonoids present in multiple organs of various plant species including flowers, fruits, leaves, stems and roots. DNA-binding R2R3-MYB transcription factors, basic helix-loop-helix (bHLH) transcription factors, and WD40 repeat proteins are known to form MYB-bHLH-WD repeat (MBW) complexes, which activates the transcription of structural genes in the anthocyanin pathway. Although black cultivars of carrots (Daucus carota L.) can accumulate large quantities of anthocyanin in their storage roots, the regulatory genes responsible for their biosynthesis are not well characterized. The current study aimed to analyze global transcription profiles based on RNA sequencing (RNA-Seq), and mine MYB, bHLH and WD40 genes that may function as positive or negative regulators in the carrot anthocyanin biosynthesis pathways. RESULTS: RNA was isolated from differently colored calli, as well as tissue samples from taproots of various black carrot cultivars across the course of development, and gene expression levels of colored and non-colored tissue and callus samples were compared. The expression of 32 MYB, bHLH and WD40 genes were significantly correlated with anthocyanin content in black carrot taproot. Of those, 11 genes were consistently up- or downregulated in a purple color-specific manner across various calli and cultivar comparisons. The expression of 10 out of these 11 genes was validated using real-time quantitative reverse transcriptase polymerase chain reaction (qRT-PCR). CONCLUSIONS: The results of this study provide insights into regulatory genes that may be responsible for carrot anthocyanin biosynthesis, and suggest that future focus on them may help improve our overall understanding of the anthocyanin synthesis pathway.

The Power of Inbreeding: NGS-Based GWAS of Rice Reveals Convergent Evolution during Rice Domestication.

Low-coverage whole-genome sequencing is an effective strategy for genome-wide association studies in humans, due to the availability of large reference panels for genotype imputation. However, it is unclear whether this strategy can be utilized in other species without reference panels. Using simulations, we show that this approach is even more relevant in inbred species such as rice (Oryza sativa L.), which are effectively haploid, allowing easy haplotype construction and imputation-based genotype calling, even without the availability of large reference panels. We sequenced 203 rice varieties with well-characterized phenotypes from the United States Department of Agriculture Rice Mini-Core Collection at an average depth of 1.5× and used the data for mapping three traits. For the first two traits, amylose content and seed length, our approach leads to direct identification of the previously identified causal SNPs in the major-effect loci. For the third trait, pericarp color, an important trait underwent selection during domestication, we identified a new major-effect locus. Although known loci can explain color variation in the varieties of two main subspecies of Asian domesticated rice, japonica and indica, the new locus identified is unique to another domesticated rice subgroup, aus, and together with existing loci, can fully explain the major variation in pericarp color in aus. Our discovery of a unique genetic basis of white pericarp in aus provides an example of convergent evolution during rice domestication and suggests that aus may have a domestication history independent of japonica and indica.

Genomic population structure of freshwater-resident and anadromous ide (Leuciscus idus) in north-western Europe.

Climate change experts largely agree that future climate change and associated rises in oceanic water levels over the upcoming decades, will affect marine salinity levels. The subsequent effects on fish communities in estuarine ecosystems however, are less clear. One species that is likely to become increasingly affected by changes in salinity is the ide (Leuciscus idus). The ide is a stenohaline freshwater fish that primarily inhabits rivers, with frequent anadromous behavior when sea salinity does not exceed 15%. Unlike most other anadromous Baltic Sea fish species, the ide has yet to be subjected to large-scale stocking programs, and thus provides an excellent opportunity for studying the natural population structure across the current salinity gradient in the Danish Belts. To explore this, we used Genotyping-by-Sequencing to determine genomic population structure of both freshwater resident and anadromous ide populations in the western Baltic Sea region, and relate the results to the current salinity gradient and the demographic history of ide in the region. The sample sites separate into four clusters, with all anadromous populations in one cluster and the freshwater resident populations in the remaining three. Results demonstrate high level of differentiation between sites hosting freshwater resident populations, but little differentiation among anadromous populations. Thus ide exhibit the genomic population structure of both a typical freshwater species, and a typical anadromous species. In addition to providing a first insight into the population structure of north-western European ide, our data also (1) provide indications of a single illegal introduction by man; (2) suggest limited genetic effects of heavy pollution in the past; and (3) indicate possible historical anadromous behavior in a now isolated freshwater population.

The origin and evolution of maize in the Southwestern United States.

The origin of maize (Zea mays mays) in the US Southwest remains contentious, with conflicting archaeological data supporting either coastal(1-4) or highland(5,6) routes of diffusion of maize into the United States. Furthermore, the genetics of adaptation to the new environmental and cultural context of the Southwest is largely uncharacterized(7). To address these issues, we compared nuclear DNA from 32 archaeological maize samples spanning 6,000 years of evolution to modern landraces. We found that the initial diffusion of maize into the Southwest about 4,000 years ago is likely to have occurred along a highland route, followed by gene flow from a lowland coastal maize beginning at least 2,000 years ago. Our population genetic analysis also enabled us to differentiate selection during domestication for adaptation to the climatic and cultural environment of the Southwest, identifying adaptation loci relevant to drought tolerance and sugar content.