AmelHap: Leveraging drone whole-genome sequence data to create a honey bee HapMap.

Honey bee, Apis mellifera, drones are typically haploid, developing from an unfertilized egg, inheriting only their queen's alleles and none from the many drones she mated with. Thus the ordered combination or 'phase' of alleles is known, making drones a valuable haplotype resource. We collated whole-genome sequence data for 1,407 drones, including 45 newly sequenced Scottish drones, collectively representing 19 countries, 8 subspecies and various hybrids. Following alignment to Amel_HAv3.1, variant calling and quality filtering, we retained 17.4 M high quality variants across 1,328 samples with a genotyping rate of 98.7%. We demonstrate the utility of this haplotype resource, AmelHap, for genotype imputation, returning >95% concordance when up to 61% of data is missing in haploids and up to 12% of data is missing in diploids. AmelHap will serve as a useful resource for the community for imputation from low-depth sequencing or SNP chip data, accurate phasing of diploids for association studies, and as a comprehensive reference panel for population genetic and evolutionary analyses.

Longitudinal transcriptome analysis of cattle infected with Theileria parva.

The apicomplexan cattle parasite Theileria parva is a major barrier to improving the livelihoods of smallholder farmers in Africa, killing over one million cattle on the continent each year. Although exotic breeds not native to Africa are highly susceptible to the disease, previous studies have illustrated that such breeds often show innate tolerance to infection by the parasite. The mechanisms underlying this tolerance remain largely unclear. To better understand the host response to T. parva infection we characterised the transcriptional response over 15 days in tolerant and susceptible cattle (n = 29) naturally exposed to the parasite. We identify key genes and pathways activated in response to infection as well as, importantly, several genes differentially expressed between the animals that ultimately survived or succumbed to infection. These include genes linked to key cell proliferation and infection pathways. Furthermore, we identify response expression quantitative trait loci containing genetic variants whose impact on the expression level of nearby genes changes in response to the infection. These therefore provide an indication of the genetic basis of differential host responses. Together these results provide a comprehensive analysis of the host transcriptional response to this under-studied pathogen, providing clues as to the mechanisms underlying natural tolerance to the disease.

Optical mapping compendium of structural variants across global cattle breeds.

Structural variants (SV) have been linked to important bovine disease phenotypes, but due to the difficulty of their accurate detection with standard sequencing approaches, their role in shaping important traits across cattle breeds is largely unexplored. Optical mapping is an alternative approach for mapping SVs that has been shown to have higher sensitivity than DNA sequencing approaches. The aim of this project was to use optical mapping to develop a high-quality database of structural variation across cattle breeds from different geographical regions, to enable further study of SVs in cattle. To do this we generated 100X Bionano optical mapping data for 18 cattle of nine different ancestries, three continents and both cattle sub-species. In total we identified 13,457 SVs, of which 1,200 putatively overlap coding regions. This resource provides a high-quality set of optical mapping-based SV calls that can be used across studies, from validating DNA sequencing-based SV calls to prioritising candidate functional variants in genetic association studies and expanding our understanding of the role of SVs in cattle evolution.

A cattle graph genome incorporating global breed diversity.

Despite only 8% of cattle being found in Europe, European breeds dominate current genetic resources. This adversely impacts cattle research in other important global cattle breeds, especially those from Africa for which genomic resources are particularly limited, despite their disproportionate importance to the continent's economies. To mitigate this issue, we have generated assemblies of African breeds, which have been integrated with genomic data for 294 diverse cattle into a graph genome that incorporates global cattle diversity. We illustrate how this more representative reference assembly contains an extra 116.1 Mb (4.2%) of sequence absent from the current Hereford sequence and consequently inaccessible to current studies. We further demonstrate how using this graph genome increases read mapping rates, reduces allelic biases and improves the agreement of structural variant calling with independent optical mapping data. Consequently, we present an improved, more representative, reference assembly that will improve global cattle research.

Using regulatory variants to detect gene-gene interactions identifies networks of genes linked to cell immortalisation.

The extent to which the impact of regulatory genetic variants may depend on other factors, such as the expression levels of upstream transcription factors, remains poorly understood. Here we report a framework in which regulatory variants are first aggregated into sets, and using these as estimates of the total cis-genetic effects on a gene we model their non-additive interactions with the expression of other genes in the genome. Using 1220 lymphoblastoid cell lines across platforms and independent datasets we identify 74 genes where the impact of their regulatory variant-set is linked to the expression levels of networks of distal genes. We show that these networks are predominantly associated with tumourigenesis pathways, through which immortalised cells are able to rapidly proliferate. We consequently present an approach to define gene interaction networks underlying important cellular pathways such as cell immortalisation.

Correction: Correlation between Cu ion migration behaviour and deNOx activity in Cu-SSZ-13 for the standard NH3-SCR reaction.

Correction for 'Correlation between Cu ion migration behaviour and deNOx activity in Cu-SSZ-13 for the standard NH3-SCR reaction' by A. M. Beale et al., Chem. Commun., 2016, 52, 6170-6173.

Genome-wide scans between two honeybee populations reveal putative signatures of human-mediated selection.

Human-mediated selection has left signatures in the genomes of many domesticated animals, including the European dark honeybee, Apis mellifera mellifera, which has been selected by apiculturists for centuries. Using whole-genome sequence information, we investigated selection signatures in spatially separated honeybee subpopulations (Switzerland, n = 39 and France, n = 17). Three different test statistics were calculated in windows of 2 kb (fixation index, cross-population extended haplotype homozygosity and cross-population composite likelihood ratio) and combined into a recently developed composite selection score. Applying a stringent false discovery rate of 0.01, we identified six significant selective sweeps distributed across five chromosomes covering eight genes. These genes are associated with multiple molecular and biological functions, including regulation of transcription, receptor binding and signal transduction. Of particular interest is a selection signature on chromosome 1, which corresponds to the WNT4 gene, the family of which is conserved across the animal kingdom with a variety of functions. In Drosophila melanogaster, WNT4 alleles have been associated with differential wing, cross vein and abdominal phenotypes. Defining phenotypic characteristics of different Apis mellifera ssp., which are typically used as selection criteria, include colour and wing venation pattern. This signal is therefore likely to be a good candidate for human mediated-selection arising from different applied breeding practices in the two managed populations.

Pitfalls in metal-organic framework crystallography: towards more accurate crystal structures.

Single crystal X-ray diffraction (SC-XRD) is the principal method for determining the crystal structures of metal-organic frameworks (MOFs). This tutorial deals with the handling of MOF crystals and analysis of crystallographic data obtained from single-crystal X-ray diffraction, focusing on two features that are particularly important in MOF crystallography and have a large impact on the quality and reliability of the final crystal structures: (1) the treatment of pore-occupying entities (both in the physical crystals and in the crystallographic model) and (2) crystallographic twinning. Proper handling of samples and data will reduce the need for using solvent masking software (e.g. SQUEEZE) to obtain acceptable crystal structures. If SC-XRD is to retain its position as the definitive method of MOF structure determination, these issues must be addressed when a new MOF structure is determined and reported. The issues addressed in this review is also valid for other porous, crystalline solids such as porous organic cages, metal-organic polyhedra, covalent organic frameworks and zeotype materials.

Two new Cu(ii) and La(iii) 2D coordination polymers, synthesis and in situ structural analysis by X-ray diffraction.

Two new coordination polymers were synthesized solvothermally using 4,4'-dimethoxy-3,3'-biphenyldicarboxylic acid (H2dmbpdc), and di- and trivalent metal salts (Cu(NO3)2·2.5H2O and La(NO3)3·6H2O). Their structures were determined by single-crystal X-ray diffraction analysis, and their thermal stability was evaluated by thermogravimetric analysis. The copper compound Cu(dmbpdc)(DMF; N,N-dimethylformamide), CPO-71-Cu, is based on the well known copper acetate paddlewheel secondary building unit. The asymmetric unit comprises one copper cation with one DMF molecule and one linker molecule coordinated. The lanthanum compound La2(dmbpdc)3(DMF)(H2O)3, CPO-72-La, is formed from a dimer of nine-coordinate, edge sharing lanthanum cations. To this dimer, three water molecules and one DMF molecule are coordinated in an ordered fashion. In addition, the asymmetric unit contains three crystallographically unique linker molecules. Both CPO-71-Cu and CPO-72-La form two-dimensional layered structures, and topological analyses reveal sql topologies with point symbol 4(4)·6(2) and vertex symbol 4·4·4·4·6(2)·6(2). The thermal behavior of CPO-71-Cu was investigated in an in situ structural analysis by variable temperature powder- and single-crystal X-ray diffraction.

Correlation between Cu ion migration behaviour and deNOx activity in Cu-SSZ-13 for the standard NH3-SCR reaction.

Here we present the results of a synchrotron-based in situ, time-resolved PXRD study during activation of two Cu-SSZ-13 catalysts under O2/He and one during standard NH3-SCR reaction conditions to obtain insight into the behaviour of Cu ions. The results obtained indicate that deNOx activity is inexorably linked with occupancy of the zeolite 6r.

Insights into the genetic history of Green-legged Partridgelike fowl: mtDNA and genome-wide SNP analysis.

The Green-legged Partridgelike (GP) fowl, an old native Polish breed, is characterised by reseda green-coloured shanks rather than yellow, white, slate or black commonly observed across most domestic breeds of chicken. Here, we investigate the origin, genetic relationships and structure of the GP fowl using mtDNA D-loop sequencing and genome-wide SNP analysis. Genome-wide association analysis between breeds enables us to verify the genetic control of the reseda green shank phenotype, a defining trait for the breed. Two mtDNA D-loop haplogroups and three autosomal genetic backgrounds are revealed. Significant associations of SNPs on chromosomes GGA24 and GGAZ indicate that the reseda green leg phenotype is associated with recessive alleles linked to the W and Id loci. Our results provide new insights into the genetic history of European chicken, indicating an admixed origin of East European traditional breeds of chicken on the continent, as supported by the presence of the reseda green phenotype and the knowledge that the GP fowl as a breed was developed before the advent of commercial stocks.

Analysis of genome-wide structure, diversity and fine mapping of Mendelian traits in traditional and village chickens.

Extensive phenotypic variation is a common feature among village chickens found throughout much of the developing world, and in traditional chicken breeds that have been artificially selected for traits such as plumage variety. We present here an assessment of traditional and village chicken populations, for fine mapping of Mendelian traits using genome-wide single-nucleotide polymorphism (SNP) genotyping while providing information on their genetic structure and diversity. Bayesian clustering analysis reveals two main genetic backgrounds in traditional breeds, Kenyan, Ethiopian and Chilean village chickens. Analysis of linkage disequilibrium (LD) reveals useful LD (r(2) ≥ 0.3) in both traditional and village chickens at pairwise marker distances of ~10 Kb; while haplotype block analysis indicates a median block size of 11-12 Kb. Association mapping yielded refined mapping intervals for duplex comb (Gga 2:38.55-38.89 Mb) and rose comb (Gga 7:18.41-22.09 Mb) phenotypes in traditional breeds. Combined mapping information from traditional breeds and Chilean village chicken allows the oocyan phenotype to be fine mapped to two small regions (Gga 1:67.25-67.28 Mb, Gga 1:67.28-67.32 Mb) totalling ~75 Kb. Mapping the unmapped earlobe pigmentation phenotype supports previous findings that the trait is sex-linked and polygenic. A critical assessment of the number of SNPs required to map simple traits indicate that between 90 and 110K SNPs are required for full genome-wide analysis of haplotype block structure/ancestry, and for association mapping in both traditional and village chickens. Our results demonstrate the importance and uniqueness of phenotypic diversity and genetic structure of traditional chicken breeds for fine-scale mapping of Mendelian traits in the species, with village chicken populations providing further opportunities to enhance mapping resolutions.