Newcastle disease virus genotype VII gene expression in experimentally infected birds.

Newcastle disease virus genotype VII (NDV-GVII) is a highly contagious pathogen responsible for pandemics that have caused devastating economic losses in the poultry industry. Several features in the transcription of NDV mRNA, including differentially expressed genes across the viral genome, are shared with that for other single, non-segmented, negative-strand viruses. Previous studies measuring viral gene expression using northern blotting indicated that the NDV transcription produced non-equimolar levels of viral mRNAs. However, deep high-throughput sequencing of virus-infected tissues can provide a better insight into the patterns of viral transcription. In this report, the transcription pattern of virulent NDV-GVII was analysed using RNA-seq and qRT-PCR. This study revealed the transcriptional profiling of these highly pathogenic NDV-GVII genes: NP:P:M:F:HN:L, in which there was a slight attenuation at the NP:P and HN:L gene boundaries. Our result also provides a fully comprehensive qPCR protocol for measuring viral transcript abundance that may be more convenient for laboratories where accessing RNA-seq is not feasible.

Genomic Prediction for Whole Weight, Body Shape, Meat Yield, and Color Traits in the Portuguese Oyster Crassostrea angulata.

Genetic improvement for quality traits, especially color and meat yield, has been limited in aquaculture because the assessment of these traits requires that the animals be slaughtered first. Genotyping technologies do, however, provide an opportunity to improve the selection efficiency for these traits. The main purpose of this study is to assess the potential for using genomic information to improve meat yield (soft tissue weight and condition index), body shape (cup and fan ratios), color (shell and mantle), and whole weight traits at harvest in the Portuguese oyster, Crassostrea angulata. The study consisted of 647 oysters: 188 oysters from 57 full-sib families from the first generation and 459 oysters from 33 full-sib families from the second generation. The number per family ranged from two to eight oysters for the first and 12-15 oysters for the second generation. After quality control, a set of 13,048 markers were analyzed to estimate the genetic parameters (heritability and genetic correlation) and predictive accuracy of the genomic selection for these traits. The multi-locus mixed model analysis indicated high estimates of heritability for meat yield traits: 0.43 for soft tissue weight and 0.77 for condition index. The estimated genomic heritabilities were 0.45 for whole weight, 0.24 for cup ratio, and 0.33 for fan ratio and ranged from 0.14 to 0.54 for color traits. The genetic correlations among whole weight, meat yield, and body shape traits were favorably positive, suggesting that the selection for whole weight would have beneficial effects on meat yield and body shape traits. Of paramount importance is the fact that the genomic prediction showed moderate to high accuracy for the traits studied (0.38-0.92). Therefore, there are good prospects to improve whole weight, meat yield, body shape, and color traits using genomic information. A multi-trait selection program using the genomic information can boost the genetic gain and minimize inbreeding in the long-term for this population.

Molecular Diagnosis of Koala Retrovirus (KoRV) in South Australian Koalas (Phascolarctos cinereus).

Koala retrovirus, a recent discovery in Australian koalas, is endogenised in 100% of northern koalas but has lower prevalence in southern populations, with lower proviral and viral loads, and an undetermined level of endogenisation. KoRV has been associated with lymphoid neoplasia, e.g., lymphoma. Recent studies have revealed high complexity in southern koala retroviral infections, with a need to clarify what constitutes positive and negative cases. This study aimed to define KoRV infection status in Mount Lofty Ranges koalas in South Australia using RNA-seq and proviral analysis (n = 216). The basis for positivity of KoRV was deemed the presence of central regions of the KoRV genome (gag 2, pol, env 1, and env 2) and based on this, 41% (89/216) koalas were positive, 57% (124/216) negative, and 2% inconclusive. These genes showed higher expression in lymph node tissue from KoRV positive koalas with lymphoma compared with other KoRV positive koalas, which showed lower, fragmented expression. Terminal regions (LTRs, partial gag, and partial env) were present in SA koalas regardless of KoRV status, with almost all (99.5%, 215/216) koalas positive for gag 1 by proviral PCR. Further investigation is needed to understand the differences in KoRV infection in southern koala populations.

Distinctive gene expression patterns and imprinting signatures revealed in reciprocal crosses between cattle sub-species.

BACKGROUND: There are two genetically distinct subspecies of cattle, Bos taurus taurus and Bos taurus indicus, which arose from independent domestication events. The two types of cattle show substantial phenotypic differences, some of which emerge during fetal development and are reflected in birth outcomes, including birth weight. We explored gene expression profiles in the placenta and four fetal tissues at mid-gestation from one taurine (Bos taurus taurus; Angus) and one indicine (Bos taurus indicus; Brahman) breed and their reciprocal crosses. RESULTS: In total 120 samples were analysed from a pure taurine breed, an indicine breed and their reciprocal cross fetuses, which identified 6456 differentially expressed genes (DEGs) between the two pure breeds in at least one fetal tissue of which 110 genes were differentially expressed in all five tissues examined. DEGs shared across tissues were enriched for pathways related to immune and stress response functions. Only the liver had a substantial number of DEGs when reciprocal crossed were compared among which 310 DEGs were found to be in common with DEGs identified between purebred livers; these DEGs were significantly enriched for metabolic process GO terms. Analysis of DEGs across purebred and crossbred tissues suggested an additive expression pattern for most genes, where both paternal and maternal alleles contributed to variation in gene expression levels. However, expression of 5% of DEGs in each tissue was consistent with parent of origin effects, with both paternal and maternal dominance effects identified. CONCLUSIONS: These data identify candidate genes potentially driving the tissue-specific differences between these taurine and indicine breeds and provide a biological insight into parental genome effects underlying phenotypic differences in bovine fetal development.

Prediction Accuracies of Genomic Selection for Nine Commercially Important Traits in the Portuguese Oyster (Crassostrea angulata) Using DArT-Seq Technology.

Genomic selection has been widely used in terrestrial animals but has had limited application in aquaculture due to relatively high genotyping costs. Genomic information has an important role in improving the prediction accuracy of breeding values, especially for traits that are difficult or expensive to measure. The purposes of this study were to (i) further evaluate the use of genomic information to improve prediction accuracies of breeding values from, (ii) compare different prediction methods (BayesA, BayesCπ and GBLUP) on prediction accuracies in our field data, and (iii) investigate the effects of different SNP marker densities on prediction accuracies of traits in the Portuguese oyster (Crassostrea angulata). The traits studied are all of economic importance and included morphometric traits (shell length, shell width, shell depth, shell weight), edibility traits (tenderness, taste, moisture content), and disease traits (Polydora sp. and Marteilioides chungmuensis). A total of 18,849 single nucleotide polymorphisms were obtained from genotyping by sequencing and used to estimate genetic parameters (heritability and genetic correlation) and the prediction accuracy of genomic selection for these traits. Multi-locus mixed model analysis indicated high estimates of heritability for edibility traits; 0.44 for moisture content, 0.59 for taste, and 0.72 for tenderness. The morphometric traits, shell length, shell width, shell depth and shell weight had estimated genomic heritabilities ranging from 0.28 to 0.55. The genomic heritabilities were relatively low for the disease related traits: Polydora sp. prevalence (0.11) and M. chungmuensis (0.10). Genomic correlations between whole weight and other morphometric traits were from moderate to high and positive (0.58-0.90). However, unfavourably positive genomic correlations were observed between whole weight and the disease traits (0.35-0.37). The genomic best linear unbiased prediction method (GBLUP) showed slightly higher accuracy for the traits studied (0.240-0.794) compared with both BayesA and BayesCπ methods but these differences were not significant. In addition, there is a large potential for using low-density SNP markers for genomic selection in this population at a number of 3000 SNPs. Therefore, there is the prospect to improve morphometric, edibility and disease related traits using genomic information in this species.

The distribution of runs of homozygosity in the genome of river and swamp buffaloes reveals a history of adaptation, migration and crossbred events.

BACKGROUND: Water buffalo is one of the most important livestock species in the world. Two types of water buffalo exist: river buffalo (Bubalus bubalis bubalis) and swamp buffalo (Bubalus bubalis carabanensis). The buffalo genome has been recently sequenced, and thus a new 90 K single nucleotide polymorphism (SNP) bead chip has been developed. In this study, we investigated the genomic population structure and the level of inbreeding of 185 river and 153 swamp buffaloes using runs of homozygosity (ROH). Analyses were carried out jointly and separately for the two buffalo types. RESULTS: The SNP bead chip detected in swamp about one-third of the SNPs identified in the river type. In total, 18,116 ROH were detected in the combined data set (17,784 SNPs), and 16,251 of these were unique. ROH were present in both buffalo types mostly detected (~ 59%) in swamp buffalo. The number of ROH per animal was larger and genomic inbreeding was higher in swamp than river buffalo. In the separated datasets (46,891 and 17,690 SNPs for river and swamp type, respectively), 19,760 and 10,581 ROH were found in river and swamp, respectively. The genes that map to the ROH islands are associated with the adaptation to the environment, fitness traits and reproduction. CONCLUSIONS: Analysis of ROH features in the genome of the two water buffalo types allowed their genomic characterization and highlighted differences between buffalo types and between breeds. A large ROH island on chromosome 2 was shared between river and swamp buffaloes and contained genes that are involved in environmental adaptation and reproduction.

Author Correction: A robust benchmark for detection of germline large deletions and insertions.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Genome Sequences of Newly Emerged Newcastle Disease Virus Strains Isolated from Disease Outbreaks in Indonesia.

Here, we report two genomes of newly emerged strains of Newcastle disease virus (NDV), Chicken/Indonesia/Tangerang/004WJ/14 and Chicken/Indonesia/VD/003WJ/11, from disease outbreaks in chickens in Indonesia. Phylogenetic study results of the fusion (F) protein's gene-coding sequences of different genotypes of NDV revealed that these two strains belong to genotype VII.1 in the class II cluster of avian paramyxoviruses.

Genome Sequences of Newcastle Disease Virus Strains from Two Outbreaks in Indonesia.

The genomes of two newly emerged Newcastle disease virus strains, chicken/Indonesia/Mega/001WJ/2013 and chicken/Indonesia/Cimanglid/002WJ/2015, from disease outbreaks in chickens in Indonesia are reported. Phylogenetic analysis of different genotypes of Newcastle disease virus using the F gene coding sequences suggests that these two strains belong to genotype VII.2, in class II of avian paramyxoviruses.

A robust benchmark for detection of germline large deletions and insertions.

New technologies and analysis methods are enabling genomic structural variants (SVs) to be detected with ever-increasing accuracy, resolution and comprehensiveness. To help translate these methods to routine research and clinical practice, we developed a sequence-resolved benchmark set for identification of both false-negative and false-positive germline large insertions and deletions. To create this benchmark for a broadly consented son in a Personal Genome Project trio with broadly available cells and DNA, the Genome in a Bottle Consortium integrated 19 sequence-resolved variant calling methods from diverse technologies. The final benchmark set contains 12,745 isolated, sequence-resolved insertion (7,281) and deletion (5,464) calls ≥50 base pairs (bp). The Tier 1 benchmark regions, for which any extra calls are putative false positives, cover 2.51 Gbp and 5,262 insertions and 4,095 deletions supported by ≥1 diploid assembly. We demonstrate that the benchmark set reliably identifies false negatives and false positives in high-quality SV callsets from short-, linked- and long-read sequencing and optical mapping.

Haplotype-resolved genomes provide insights into structural variation and gene content in Angus and Brahman cattle.

Inbred animals were historically chosen for genome analysis to circumvent assembly issues caused by haplotype variation but this resulted in a composite of the two genomes. Here we report a haplotype-aware scaffolding and polishing pipeline which was used to create haplotype-resolved, chromosome-level genome assemblies of Angus (taurine) and Brahman (indicine) cattle subspecies from contigs generated by the trio binning method. These assemblies reveal structural and copy number variants that differentiate the subspecies and that variant detection is sensitive to the specific reference genome chosen. Six genes with immune related functions have additional copies in the indicine compared with taurine lineage and an indicus-specific extra copy of fatty acid desaturase is under positive selection. The haplotyped genomes also enable transcripts to be phased to detect allele-specific expression. This work exemplifies the value of haplotype-resolved genomes to better explore evolutionary and functional variations.

Metal resistant bacteria on gold particles: Implications of how anthropogenic contaminants could affect natural gold biogeochemical cycling.

In Earth's near-surface environments, gold biogeochemical cycling involves gold dissolution and precipitation processes, which are partly attributed to bacteria. These biogeochemical processes as well as abrasion (via physical transport) are known to act upon gold particles, thereby resulting in particle transformation including the development of pure secondary gold and altered morphology, respectively. While previous studies have inferred gold biogeochemical cycling from gold particles obtained from natural environments, little is known about how metal contamination in an environment could impact this cycle. Therefore, this study aims to infer how potentially toxic metal contaminants could affect the structure and chemistry of gold particles and therefore the biogeochemical cycling of gold. In doing so, river sediments and gold particles from the De Kaap Valley, South Africa, were analysed using both microanalytical and molecular techniques. Of the metal contaminants detected in the sediment, mercury can chemically interact with gold particles thereby directly altering particle morphology and "erasing" textural evidence indicative of particle transformation. Other metal contaminants (including mercury) indirectly affect gold cycling by exerting a selective pressure on bacteria living on the surface of gold particles. Particles harbouring gold-tolerant bacteria with diverse metal resistant genes, such as Arthrobacter sp. and Pseudomonas sp., contained nearly two times more secondary gold relative to particles harbouring bacteria with less gold-tolerance. In conclusion, metal contaminants can have a direct or indirect effect on gold biogeochemical cycling in natural environments impacted by anthropogenic activity.

A genetically unique Chinese cattle population shows evidence of common ancestry with wild species when analysed with a reduced ascertainment bias SNP panel.

In Hong Kong, there is a cattle population of ~1,200 individuals of uncertain origin and genetic diversity. This population shows heterogeneous morphology, both in body type and pigmentation. Once used as draught animals by the local farmers, they were abandoned around the 1970s due to changes in the economy, and since then have lived as feral populations. To explore the origins of these cattle, we analysed ~50k genotype data of 21 Hong Kong feral cattle, along with data from 703 individuals of 36 cattle populations of European, African taurine, and Asian origin, the wild x domestic hybrid gayal, plus two wild bovine species, gaur and banteng. To reduce the effect of ascertainment bias ~4k loci that are polymorphic in the two wild species were selected for further analysis. The stringent SNP selection we applied resulted in increased heterozygosity across all populations studies, compared with the full panel of SNP, thus reducing the impact of ascertainment bias and facilitating the comparison of divergent breeds of cattle. Our results showed that Hong Kong feral cattle have relatively high levels of genetic distinctiveness, possibly due to the low level of artificial selection, and a likely common ancestry with wild species. We found signs of a putative taurine introgression, probably dating to the import of north European breeds during the British colonialism of Hong Kong. We showed that Hong Kong feral cattle, are distinct from Bos taurus and Bos indicus breeds. Our results highlight the distinctiveness of Hong Kong feral cattle and stress the conservation value of this indigenous breed that is likely to harbour adaptive genetic variation, which is a fundamental livestock resource in the face of climate change and diversifying market demands.

New insights into mammalian sex chromosome structure and evolution using high-quality sequences from bovine X and Y chromosomes.

BACKGROUND: Mammalian X chromosomes are mainly euchromatic with a similar size and structure among species whereas Y chromosomes are smaller, have undergone substantial evolutionary changes and accumulated male specific genes and genes involved in sex determination. The pseudoautosomal region (PAR) is conserved on the X and Y and pair during meiosis. The structure, evolution and function of mammalian sex chromosomes, particularly the Y chromsome, is still poorly understood because few species have high quality sex chromosome assemblies. RESULTS: Here we report the first bovine sex chromosome assemblies that include the complete PAR spanning 6.84 Mb and three Y chromosome X-degenerate (X-d) regions. The PAR comprises 31 genes, including genes that are missing from the X chromosome in current cattle, sheep and goat reference genomes. Twenty-nine PAR genes are single-copy genes and two are multi-copy gene families, OBP, which has 3 copies and BDA20, which has 4 copies. The Y chromosome X-d1, 2a and 2b regions contain 11, 2 and 2 gametologs, respectively. CONCLUSIONS: The ruminant PAR comprises 31 genes and is similar to the PAR of pig and dog but extends further than those of human and horse. Differences in the pseudoautosomal boundaries are consistent with evolutionary divergence times. A bovidae-specific expansion of members of the lipocalin gene family in the PAR reported here, may affect immune-modulation and anti-inflammatory responses in ruminants. Comparison of the X-d regions of Y chromosomes across species revealed that five of the X-Y gametologs, which are known to be global regulators of gene activity and candidate sexual dimorphism genes, are conserved.

In vivo passage of Salmonella Typhimurium results in minor mutations in the bacterial genome and increases in vitro invasiveness.

Eggs and raw or undercooked egg-containing food items are frequently identified as the bacterial source during epidemiolocal investigation of Salmonella outbreaks. Multi-locus variable number of tandem repeats analysis (MLVA) is a widely used Salmonella typing method enabling the study of diversity within populations of the same serotype. In vivo passage, however, has been linked with changes in MLVA type and more broadly the Salmonella genome. We sought to investigate whether in vivo passage through layer hens had an effect on MLVA type as well as the bacterial genome and whether any mutations affected bacterial virulence. Layer hens were infected with either Salmonella Typhimurium DT9 (03-24-11-11-523) as part of a single infection or were co-infected with an equal amount of Salmonella Mbandaka. Salmonella shedding in both single and co-infected birds was variable over the course of the 16-week experiment. Salmonella Typhimurium and Salmonella Mbandaka were identified in feces of co-infected birds. Salmonella colonies isolated from fecal samples were subtyped using MLVA. A single change in SSTR-6 was observed in Salmonella Typhimurium strains isolated from co-infected birds. Isolates of Salmonella Typhimurium of both the parent (03-24-11-11-523) and modified (03-24-12-11-523) MLVA type were sequenced and compared with the genome of the parent strain. Sequence analysis revealed that in vivo passaging resulted in minor mutation events. Passaged isolates exhibited significantly higher invasiveness in cultured human intestinal epithelial cells than the parent strain. The microevolution observed in this study suggests that changes in MLVA may arise more commonly and may have clinical significance.

Genetic cause for congenital methemoglobinemia in an Australian Pomeranian dog.

Little is known about genetic causes of congenital methemoglobinemia in dogs. Here, we report a CYB5 R3 mutation in a Pomeranian dog with congenital methemoglobinemia. A 6-year-old neutered female Pomeranian dog was investigated for cyanosis noticed during anesthesia for an orthopedic procedure. The history included lifelong mild exercise intolerance and bluish tongue. Methemoglobinemia was diagnosed using co-oximetry. The CYB5 R3 gene was analyzed by comparing the patient's genomic DNA with the reference canine sequence. Mutation functional significance was investigated using snpEff and multispecies protein homology analyses. A homozygous missense single nucleotide CYB5 R3 mutation (ATC ➔ CTC at codon 194) caused a p.Ile194Leu substitution. The pIle194 residue is highly conserved in other mammals, supporting the likely pathogenicity of the substitution. The mutation described here is identical to that associated with familial methemoglobinemia in a family of Japanese Pomeranian dogs. This observation, together with the homozygous mutation found in our case, indicates that the mutant allele may be widespread within the Pomeranian breed internationally.

Chromosome-level assembly of the water buffalo genome surpasses human and goat genomes in sequence contiguity.

Rapid innovation in sequencing technologies and improvement in assembly algorithms have enabled the creation of highly contiguous mammalian genomes. Here we report a chromosome-level assembly of the water buffalo (Bubalus bubalis) genome using single-molecule sequencing and chromatin conformation capture data. PacBio Sequel reads, with a mean length of 11.5 kb, helped to resolve repetitive elements and generate sequence contiguity. All five B. bubalis sub-metacentric chromosomes were correctly scaffolded with centromeres spanned. Although the index animal was partly inbred, 58% of the genome was haplotype-phased by FALCON-Unzip. This new reference genome improves the contig N50 of the previous short-read based buffalo assembly more than a thousand-fold and contains only 383 gaps. It surpasses the human and goat references in sequence contiguity and facilitates the annotation of hard to assemble gene clusters such as the major histocompatibility complex (MHC).

Transcriptional analysis of liver from chickens with fast (meat bird), moderate (F1 layer x meat bird cross) and low (layer bird) growth potential.

BACKGROUND: Divergent selection for meat and egg production in poultry has resulted in strains of birds differing widely in traits related to these products. Modern strains of meat birds can reach live weights of 2 kg in 35 d, while layer strains are now capable of producing more than 300 eggs per annum but grow slowly. In this study, RNA-Seq was used to investigate hepatic gene expression between three groups of birds with large differences in growth potential; meat bird, layer strain as well as an F1 layer x meat bird. The objective was to identify differentially expressed (DE) genes between all three strains to elucidate biological factors underpinning variations in growth performance. RESULTS: RNA-Seq analysis was carried out on total RNA extracted from the liver of meat bird (n = 6), F1 layer x meat bird cross (n = 6) and layer strain (n = 6), males. Differential expression of genes were considered significant at P < 0.05, and a false discovery rate of < 0.05, with any fold change considered. In total, 6278 genes were found to be DE with 5832 DE between meat birds and layers (19%), 2935 DE between meat birds and the cross (9.6%) and 493 DE between the cross and layers (1.6%). Comparisons between the three groups identified 155 significant DE genes. Gene ontology (GO) enrichment and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analysis of the 155 DE genes showed the FoxO signalling pathway was most enriched (P = 0.001), including genes related to cell cycle regulation and insulin signalling. Significant GO terms included 'positive regulation of glucose import' and 'cellular response to oxidative stress', which is also consistent with FoxOs regulation of glucose metabolism. There were high correlations between FoxO pathway genes and bodyweight, as well as genes related to glycolysis and bodyweight. CONCLUSIONS: This study revealed large transcriptome differences between meat and layer birds. There was significant evidence implicating the FoxO signalling pathway (via cell cycle regulation and altered metabolism) as an active driver of growth variations in chicken. Functional analysis of the FoxO genes is required to understand how they regulate growth and egg production.

Evaluation of three cryoprotectants used with bovine milk affected with Mycoplasma bovis in different freezing conditions.

OBJECTIVES: Currently, there is no consensus protocols regarding the combination of glycerol (GLY), gelatin or foetal bovine serum (FBS) with dimethyl sulphoxide (DMSO) as cryoprotectants for Mycoplasma bovis in bovine milk samples. This study aimed to compare different cryopreservation compounds and storage temperatures for M. bovis. RESULTS: There were significant differences in the survival of M. bovis on different media. Differences were also observed between different storage conditions. All additives improved the survival of M. bovis in comparison to control (CON). The combination of GLY and DMSO was shown to be significantly different to CON with 57.1% (95% CI = 21.43-133.34) and 19.1% (95% CI = 11.73-60.27), respectively at week 16, and its use should be encouraged as a cryoprotectant for M. bovis at - 20 and - 80 °C. GEL/DMSO showed the highest survival rate for M. bovis with 57.14% (95% CI = 21.43-133.34) at 4 °C in comparison with CON 14.29% (95% CI = 9.60-50.39). FBS/DMSO showed the highest survival rate for the short-term preservation similarly to other additives. The evaluated cryopreservative compounds would improve survivability of M. bovis in milk for both transport and long-term storage. Hence, it is recommended to use the mentioned methods for routine transportation or storage purposes for suspicious M. bovis milk samples.