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.

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.

Evaluating the contribution of rare variants to type 2 diabetes and related traits using pedigrees.

A major challenge in evaluating the contribution of rare variants to complex disease is identifying enough copies of the rare alleles to permit informative statistical analysis. To investigate the contribution of rare variants to the risk of type 2 diabetes (T2D) and related traits, we performed deep whole-genome analysis of 1,034 members of 20 large Mexican-American families with high prevalence of T2D. If rare variants of large effect accounted for much of the diabetes risk in these families, our experiment was powered to detect association. Using gene expression data on 21,677 transcripts for 643 pedigree members, we identified evidence for large-effect rare-variant cis-expression quantitative trait loci that could not be detected in population studies, validating our approach. However, we did not identify any rare variants of large effect associated with T2D, or the related traits of fasting glucose and insulin, suggesting that large-effect rare variants account for only a modest fraction of the genetic risk of these traits in this sample of families. Reliable identification of large-effect rare variants will require larger samples of extended pedigrees or different study designs that further enrich for such variants.

Genome sequencing identifies major causes of severe intellectual disability.

Severe intellectual disability (ID) occurs in 0.5% of newborns and is thought to be largely genetic in origin. The extensive genetic heterogeneity of this disorder requires a genome-wide detection of all types of genetic variation. Microarray studies and, more recently, exome sequencing have demonstrated the importance of de novo copy number variations (CNVs) and single-nucleotide variations (SNVs) in ID, but the majority of cases remain undiagnosed. Here we applied whole-genome sequencing to 50 patients with severe ID and their unaffected parents. All patients included had not received a molecular diagnosis after extensive genetic prescreening, including microarray-based CNV studies and exome sequencing. Notwithstanding this prescreening, 84 de novo SNVs affecting the coding region were identified, which showed a statistically significant enrichment of loss-of-function mutations as well as an enrichment for genes previously implicated in ID-related disorders. In addition, we identified eight de novo CNVs, including single-exon and intra-exonic deletions, as well as interchromosomal duplications. These CNVs affected known ID genes more frequently than expected. On the basis of diagnostic interpretation of all de novo variants, a conclusive genetic diagnosis was reached in 20 patients. Together with one compound heterozygous CNV causing disease in a recessive mode, this results in a diagnostic yield of 42% in this extensively studied cohort, and 62% as a cumulative estimate in an unselected cohort. These results suggest that de novo SNVs and CNVs affecting the coding region are a major cause of severe ID. Genome sequencing can be applied as a single genetic test to reliably identify and characterize the comprehensive spectrum of genetic variation, providing a genetic diagnosis in the majority of patients with severe ID.

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.

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.

Detection and phasing of single base de novo mutations in biopsies from human in vitro fertilized embryos by advanced whole-genome sequencing.

Currently, the methods available for preimplantation genetic diagnosis (PGD) of in vitro fertilized (IVF) embryos do not detect de novo single-nucleotide and short indel mutations, which have been shown to cause a large fraction of genetic diseases. Detection of all these types of mutations requires whole-genome sequencing (WGS). In this study, advanced massively parallel WGS was performed on three 5- to 10-cell biopsies from two blastocyst-stage embryos. Both parents and paternal grandparents were also analyzed to allow for accurate measurements of false-positive and false-negative error rates. Overall, >95% of each genome was called. In the embryos, experimentally derived haplotypes and barcoded read data were used to detect and phase up to 82% of de novo single base mutations with a false-positive rate of about one error per Gb, resulting in fewer than 10 such errors per embryo. This represents a ∼ 100-fold lower error rate than previously published from 10 cells, and it is the first demonstration that advanced WGS can be used to accurately identify these de novo mutations in spite of the thousands of false-positive errors introduced by the extensive DNA amplification required for deep sequencing. Using haplotype information, we also demonstrate how small de novo deletions could be detected. These results suggest that phased WGS using barcoded DNA could be used in the future as part of the PGD process to maximize comprehensiveness in detecting disease-causing mutations and to reduce the incidence of genetic diseases.

Whole Genome Sequencing Increases Molecular Diagnostic Yield Compared with Current Diagnostic Testing for Inherited Retinal Disease.

PURPOSE: To compare the efficacy of whole genome sequencing (WGS) with targeted next-generation sequencing (NGS) in the diagnosis of inherited retinal disease (IRD). DESIGN: Case series. PARTICIPANTS: A total of 562 patients diagnosed with IRD. METHODS: We performed a direct comparative analysis of current molecular diagnostics with WGS. We retrospectively reviewed the findings from a diagnostic NGS DNA test for 562 patients with IRD. A subset of 46 of 562 patients (encompassing potential clinical outcomes of diagnostic analysis) also underwent WGS, and we compared mutation detection rates and molecular diagnostic yields. In addition, we compared the sensitivity and specificity of the 2 techniques to identify known single nucleotide variants (SNVs) using 6 control samples with publically available genotype data. MAIN OUTCOME MEASURES: Diagnostic yield of genomic testing. RESULTS: Across known disease-causing genes, targeted NGS and WGS achieved similar levels of sensitivity and specificity for SNV detection. However, WGS also identified 14 clinically relevant genetic variants through WGS that had not been identified by NGS diagnostic testing for the 46 individuals with IRD. These variants included large deletions and variants in noncoding regions of the genome. Identification of these variants confirmed a molecular diagnosis of IRD for 11 of the 33 individuals referred for WGS who had not obtained a molecular diagnosis through targeted NGS testing. Weighted estimates, accounting for population structure, suggest that WGS methods could result in an overall 29% (95% confidence interval, 15-45) uplift in diagnostic yield. CONCLUSIONS: We show that WGS methods can detect disease-causing genetic variants missed by current NGS diagnostic methodologies for IRD and thereby demonstrate the clinical utility and additional value of WGS.

Developmental timing of mutations revealed by whole-genome sequencing of twins with acute lymphoblastic leukemia.

Acute lymphoblastic leukemia (ALL) is the major pediatric cancer. At diagnosis, the developmental timing of mutations contributing critically to clonal diversification and selection can be buried in the leukemia's covert natural history. Concordance of ALL in monozygotic, monochorionic twins is a consequence of intraplacental spread of an initiated preleukemic clone. Studying monozygotic twins with ALL provides a unique means of uncovering the timeline of mutations contributing to clonal evolution, pre- and postnatally. We sequenced the whole genomes of leukemic cells from two twin pairs with ALL to comprehensively characterize acquired somatic mutations in ALL, elucidating the developmental timing of all genetic lesions. Shared, prenatal, coding-region single-nucleotide variants were limited to the putative initiating lesions. All other nonsynonymous single-nucleotide variants were distinct between tumors and, therefore, secondary and postnatal. These changes occurred in a background of noncoding mutational changes that were almost entirely discordant in twin pairs and likely passenger mutations acquired during leukemic cell proliferation.

Whole genome sequences are required to fully resolve the linkage disequilibrium structure of human populations.

BACKGROUND: An understanding of linkage disequilibrium (LD) structures in the human genome underpins much of medical genetics and provides a basis for disease gene mapping and investigating biological mechanisms such as recombination and selection. Whole genome sequencing (WGS) provides the opportunity to determine LD structures at maximal resolution. RESULTS: We compare LD maps constructed from WGS data with LD maps produced from the array-based HapMap dataset, for representative European and African populations. WGS provides up to 5.7-fold greater SNP density than array-based data and achieves much greater resolution of LD structure, allowing for identification of up to 2.8-fold more regions of intense recombination. The absence of ascertainment bias in variant genotyping improves the population representativeness of the WGS maps, and highlights the extent of uncaptured variation using array genotyping methodologies. The complete capture of LD patterns using WGS allows for higher genome-wide association study (GWAS) power compared to array-based GWAS, with WGS also allowing for the analysis of rare variation. The impact of marker ascertainment issues in arrays has been greatest for Sub-Saharan African populations where larger sample sizes and substantially higher marker densities are required to fully resolve the LD structure. CONCLUSIONS: WGS provides the best possible resource for LD mapping due to the maximal marker density and lack of ascertainment bias. WGS LD maps provide a rich resource for medical and population genetics studies. The increasing availability of WGS data for large populations will allow for improved research utilising LD, such as GWAS and recombination biology studies.

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.

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.

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.

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.

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.

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.

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.

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.