What's in your next-generation sequence data? An exploration of unmapped DNA and RNA sequence reads from the bovine reference individual.

BackgroundNext-generation sequencing projects commonly commence by aligning reads to a reference genome assembly. While improvements in alignment algorithms and computational hardware have greatly enhanced the efficiency and accuracy of alignments, a significant percentage of reads often remain unmapped.ResultsWe generated de novo assemblies of unmapped reads from the DNA and RNA sequencing of the Bos taurus reference individual and identified the closest matching sequence to each contig by alignment to the NCBI non-redundant nucleotide database using BLAST. As expected, many of these contigs represent vertebrate sequence that is absent, incomplete, or misassembled in the UMD3.1 reference assembly. However, numerous additional contigs represent invertebrate species. Most prominent were several species of Spirurid nematodes and a blood-borne parasite, Babesia bigemina. These species are either not present in the US or are not known to infect taurine cattle and the reference animal appears to have been host to unsequenced sister species.ConclusionsWe demonstrate the importance of exploring unmapped reads to ascertain sequences that are either absent or misassembled in the reference assembly and for detecting sequences indicative of parasitic or commensal organisms

Greenland ice sheet surface mass loss: recent developments in observation and modeling

Surface processes currently dominate Greenland ice sheet (GrIS) mass loss. We review recent developments in the observation and modelling of GrIS surface mass balance (SMB), published after the July 2012 deadline for the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR5). Since IPCC AR5 our understanding of GrIS SMB has further improved, but new observational and model studies have also revealed that temporal and spatial variability of many processes are still poorly quantified and understood, e.g. bio-albedo, the formation of ice lenses and their impact on lateral meltwater transport, heterogeneous vertical meltwater transport (‘piping’), the impact of atmospheric circulation changes and mixed-phase clouds on the surface energy balance and the magnitude of turbulent heat exchange over rough ice surfaces. As a result, these processes are only schematically or not at all included in models that are currently used to assess and predict future GrIS surface mass loss

Gene number determination and genetic polymorphism of the gamma delta T cell co-receptor WC1 genes

Background WC1 co-receptors belong to the scavenger receptor cysteine-rich (SRCR) superfamily and are encoded by a multi-gene family. Expression of particular WC1 genes defines functional subpopulations of WC1+ γδ T cells. We have previously identified partial or complete genomic sequences for thirteen different WC1 genes through annotation of the bovine genome Btau_3.1 build. We also identified two WC1 cDNA sequences from other cattle that did not correspond to sequences in the Btau_3.1 build. Their absence in the Btau_3.1 build may have reflected gaps in the genome assembly or polymorphisms among animals. Since the response of γδ T cells to bacterial challenge is determined by WC1 gene expression, it was critical to understand whether individual cattle or breeds differ in the number of WC1 genes or display polymorphisms. Results Real-time quantitative PCR using DNA from the animal whose genome was sequenced (“Dominette”) and sixteen other animals representing ten breeds of cattle, showed that the number of genes coding for WC1 co-receptors is thirteen. The complete coding sequences of those thirteen WC1 genes is presented, including the correction of an error in the WC1-2 gene due to mis-assembly in the Btau_3.1 build. All other cDNA sequences were found to agree with the previous annotation of complete or partial WC1 genes. PCR amplification and sequencing of the most variable N-terminal SRCR domain (domain 1 which has the SRCR “a” pattern) of each of the thirteen WC1 genes showed that the sequences are highly conserved among individuals and breeds. Of 160 sequences of domain 1 from three breeds of cattle, no additional sequences beyond the thirteen described WC1 genes were found. Analysis of the complete WC1 cDNA sequences indicated that the thirteen WC1 genes code for three distinct WC1 molecular forms. Conclusion The bovine WC1 multi-gene family is composed of thirteen genes coding for three structural forms whose sequences are highly conserved among individual cattle and breeds. The sequence diversity necessary for WC1 genes to function as a multi-genic pattern recognition receptor array is encoded in the genome, rather than generated by recombinatorial diversity or hypermutation

Data from: Origins of cattle on Chirikof Island, Alaska, elucidated from genome-wide SNP genotypes

Feral livestock may harbor genetic variation of commercial, scientific, historical or esthetic value. The origins and uniqueness of feral cattle on Chirikof Island, Alaska, are uncertain. The island is now part of the Alaska Maritime Wildlife Refuge and Federal wildlife managers want grazing to cease, presumably leading to demise of the cattle. Here we characterize the cattle of Chirikof Island relative to extant breeds and discern their origins. Our analyses support the inference that Yakut cattle from Russia arrived first on Chirikof Island, then ~120 years ago the first European taurine cattle were introduced to the island, and finally a large wave of Hereford cattle were introduced on average 40 years ago. In addition, this mixture of European and East-Asian cattle is unique compared with other North American breeds and we find evidence that natural selection in the relatively harsh environment of Chirikof Island has further impacted their genetic architecture. These results provide an objective basis for decisions regarding conservation of the Chirikof Island cattle

Gene number determination and genetic polymorphism of the gamma delta T cell co-receptor WC1 genes

Abstract Background WC1 co-receptors belong to the scavenger receptor cysteine-rich (SRCR) superfamily and are encoded by a multi-gene family. Expression of particular WC1 genes defines functional subpopulations of WC1+ γδ T cells. We have previously identified partial or complete genomic sequences for thirteen different WC1 genes through annotation of the bovine genome Btau_3.1 build. We also identified two WC1 cDNA sequences from other cattle that did not correspond to sequences in the Btau_3.1 build. Their absence in the Btau_3.1 build may have reflected gaps in the genome assembly or polymorphisms among animals. Since the response of γδ T cells to bacterial challenge is determined by WC1 gene expression, it was critical to understand whether individual cattle or breeds differ in the number of WC1 genes or display polymorphisms. Results Real-time quantitative PCR using DNA from the animal whose genome was sequenced (“Dominette”) and sixteen other animals representing ten breeds of cattle, showed that the number of genes coding for WC1 co-receptors is thirteen. The complete coding sequences of those thirteen WC1 genes is presented, including the correction of an error in the WC1-2 gene due to mis-assembly in the Btau_3.1 build. All other cDNA sequences were found to agree with the previous annotation of complete or partial WC1 genes. PCR amplification and sequencing of the most variable N-terminal SRCR domain (domain 1 which has the SRCR “a” pattern) of each of the thirteen WC1 genes showed that the sequences are highly conserved among individuals and breeds. Of 160 sequences of domain 1 from three breeds of cattle, no additional sequences beyond the thirteen described WC1 genes were found. Analysis of the complete WC1 cDNA sequences indicated that the thirteen WC1 genes code for three distinct WC1 molecular forms. Conclusion The bovine WC1 multi-gene family is composed of thirteen genes coding for three structural forms whose sequences are highly conserved among individual cattle and breeds. The sequence diversity necessary for WC1 genes to function as a multi-genic pattern recognition receptor array is encoded in the genome, rather than generated by recombinatorial diversity or hypermutation.</p

Discovery of novel genetic networks associated with 19 economically important traits in beef cattle

<p>Quantitative or complex traits are determined by the combined effects of many loci, and are affected by genetic networks or molecular pathways. In the present study, we genotyped a total of 138 mutations, mainly single nucleotide polymorphisms derived from 71 functional genes on a Wagyu x Limousin reference population. Two hundred forty six F₂ animals were measured for 5 carcass, 6 eating quality and 8 fatty acid composition traits. A total of 2,280 single marker-trait association runs with 120 tagged mutations selected based on the HAPLOVIEW analysis revealed 144 significant associations (P &#60; 0.05), but 50 of them were removed from the analysis due to the small number of animals (&#8804; 9) in one genotype group or absence of one genotype among three genotypes. The remaining 94 single-trait associations were then placed into three groups of quantitative trait modes (QTMs) with additive, dominant and overdominant effects. All significant markers and their QTMs associated with each of these 19 traits were involved in a linear regression model analysis, which confirmed single-gene associations for 4 traits, but revealed two-gene networks for 8 traits and three-gene networks for 5 traits. Such genetic networks involving both genotypes and QTMs resulted in high correlations between predicted and actual values of performance, thus providing evidence that the classical Mendelian principles of inheritance can be applied in understanding genetic complexity of complex phenotypes. Our present study also indicated that carcass, eating quality and fatty acid composition traits rarely share genetic networks. Therefore, marker-assisted selection for improvement of one category of these traits would not interfere with improvement of another.</p