Context The alpaca roan pattern is characterised by white and coloured fibre interspersed together, with a distinctive lighter body and darker extremities, and commonly is believed to be inherited in an autosomal dominant manner. It is of interest to the alpaca fibre industry as it causes ‘contamination’ of coloured fibre with white fibres, but cannot be detected in white or light fawn animals. Other livestock species, such as horses, cattle, goats, and pigs, exhibit comparable phenotypes, which are associated with candidate variant(s) in either KIT or KITLG. Aims To identify a region or regions of the genome that is/are causative of the roan pattern in alpacas. Methods We conducted a genome-wide association study (GWAS) by using 13 roan and 14 non-roan alpacas sampled from the USA, Australia, and New Zealand. Regions of genome-wide significance were examined for variants that correlated with the roan phenotype. Key results A novel candidate single-nucleotype polymorphism (SNP; Super-Scaffold_15:39 742 851T > A), located 272 kb upstream of KITLG, was identified in 1 of 12 regions with genome-wide significant association (P ≤ 5 × 10−8). We identified the candidate SNP-containing region (Super-Scaffold_15:39 742 096–39 887 419) to be a 145 kb copy number variant (CNV) that is likely to be a tandem duplication. All 13 roan alpacas had one or two copies of the roan-associated T allele and all except three non-roans had zero copies. Furthermore, we determined the Mendelian inheritance of copy number haplotypes and their allelic composition in a roan and a non-roan family. Conclusions Our data support the hypothesised autosomal incomplete dominant mode of inheritance of the roan pattern in alpacas and suggests that the effect of the T allele CNV version is likely to be suppressed when in cis with the A allele CNV version. However, additional verification is required to validate the finding and determine the functional effect. Implications Identification of the cause, or a marker for roan pattern will allow alpaca breeders to select for or against the roan pattern, even when the phenotype is hidden, and therefore increase production output and profitability.

Human cytomegalovirus (HCMV) is carried lifelong by similar to 80 % of adults worldwide, generating distinct disease syndromes in transplant recipients, people with HIV (PWH) and neonates. Amino acids 15-23 encoded by the HCMV gene UL40 match positions 3-11 of HLA-A and HLA-C, and constitute a "signal peptide" able to stabilise cell surface HLA-E as a restriction element and a ligand of NKG2A and NKG2C. We present next generation sequencing of UL40 amplified from 15 Australian renal transplant recipients (RTR), six healthy adults and four neonates, and 21 Indonesian PWH. We found no groupwise associations between the presence of multiple sequences and HCMV burden (highest in PWH) or HCMV-associated symptoms in neonates. Homology between UL40 and corresponding HLA-C and HLA-A peptides in 11 RTR revealed perfect matches with HLA-C in three individuals, all carrying HCMV encoding only VMAPRTLIL - a peptide previously associated with viremia. However indices of the burden of HCMV did not segregate in our cohort.

Context Currently, there is a growing consumer demand for more ecologically sustainable practices in the textile industry. Fabric dyeing is highly pollutive, and one way to avoid dyeing is to use naturally coloured fibres. Alpacas exhibit a wide range of fleece colours, thereby making them a good source of fibre for sustainable textile production. Our understanding of the colour genetics of alpacas is improving, but there is still no explanation for all the variation seen in alpaca coat colours. Aims To identify a region or regions in the alpaca genome that contribute to differences in pigment intensity. Methods Colorimetric analysis using L*a*b* colourspace of fibre from white and black alpacas, that had been genotyped using the Neogen Australasia alpaca coat colour test, was used to classify the samples into intense and dilute groups for each colour. Illumina short-read genotyping by sequencing of the DNA from these alpacas was used to identify single nucleotide polymorphisms that were subsequently used in a case–control genome-wide association study (GWAS) comparing the extreme dilute and intensely pigmented animals. Key results Intense eumelanic fibre is darker (P = 0.0003), less red (P = 0.004), and more blue (P = 0.001) than is dilute eumelanic fibre. Intense pheomelanic fibre is darker (P = 1 × 10−7), more red (P = 3 × 10−20), and more yellow (P = 2 × 10−6) than is dilute pheomelanic fibre. The GWAS showed six regions of genome-wide significance. After manual inspection of these six regions, the best candidate region was upstream of KITLG, a gene previously associated with pigment intensity in dogs. Conclusions In combination with ASIP genotype, a regulatory mutation in a region upstream of KITLG in alpacas potentially has a major effect on pigment intensity in the species. Implications Successful identification of a marker for pigment intensity will allow breeders to select more precisely for breeding and production animals that will assist them in supplying the desired fibre colours to the textile industry.

We propose a general deep variational model (reduced version, full version as well as the extension) via a comprehensive fusion approach in this paper. It is able to realize various image tasks in a completely unsupervised way without learning from samples. Technically, it can properly incorporate the CNN based deep image prior (DIP) architecture into the classic variational image processing models. The minimization problem solving strategy is transformed from iteratively minimizing the sub-problem for each variable to automatically minimizing the loss function by learning the generator network parameters. The proposed deep variational (DV) model contributes to the high order image edition and applications such as image restoration, inpainting, decomposition and texture segmentation. Experiments conducted have demonstrated significant advantages of the proposed deep variational model in comparison with several powerful techniques including variational methods and deep learning approaches.

Human cytomegalovirus (HCMV) is a beta-herpesvirus carried by ~80% of adults worldwide. Acute infections are often asymptomatic in healthy individuals but generate diverse syndromes in neonates, renal transplant recipients (RTR), and people with HIV (PWH). The HCMV gene UL111a encodes a homolog of human interleukin-10 (IL-10) that interacts with the human IL-10 receptor. Deep sequencing technologies were used to sequence UL111a directly from 59 clinical samples from Indonesian PWH and Australian RTR, healthy adults, and neonates. Overall, 93% of samples contained more than one variant of HCMV, as defined by at least one nonsynonymous variation. Carriage of these variants differed between neonates and adults, Australians and Indonesians, and between saliva and blood leukocytes. The variant alleles of N41D and S71Y occurred together in Australian RTR and were associated with higher T-cell responses to HCMV pp65. The variant P122S was associated with lower levels of antibodies reactive with a lysate of HCMV-infected fibroblasts. L174F was associated with increased levels of antibodies reactive with HCMV lysate, immediate-early 1 (IE-1), and glycoprotein B (gB) in Australian RTR and Indonesians PWH, suggesting a higher viral burden. We conclude that variants of UL111a are common in all populations and may influence systemic responses to HCMV.

Around 80% of adults worldwide carry human cytomegaloviris (HCMV). The HCMV gene UL18 is a homolog of HLA class I genes and encodes a protein with high affinity for the NK and T-cell cytotoxicity inhibitor LIR-1. UL18 was deep sequenced from blood, saliva or urine from Indonesian people with HIV (PWH) (n = 28), Australian renal transplant recipients (RTR) (n = 21), healthy adults (n = 7) and neonates (n = 4). 95% of samples contained more than one variant of HCMV UL18, as defined by carriage of nonsynonymous variations. When aligned with immunological markers of the host’s burden of HCMV, the S318N variation associated with high levels of antibody reactive with HCMV lysate in PWH over 12 months on antiretroviral therapy. The A107T variation associated with HCMV antibody levels and inflammatory biomarkers in PWH at early timepoints. Variants D32G, D248N, V250A and E252D aligned with elevated HCMV antibody levels in RTR, while M191K, E196Q and F165L were associated with HCMV-reactive T-cells and proportions of Vδ2− γδ T-cells—populations linked with high burdens of HCMV. We conclude that UL18 is a highly variable gene, where variation may alter the persistent burden of HCMV and/or the host response to that burden.

Chronic psychosocial stress is increasingly being recognised as a risk factor for sporadic Alzheimer's disease (AD). The hypothalamic-pituitary-adrenal axis (HPA axis) is the major stress response pathway in the body and tightly regulates the production of cortisol, a glucocorticoid hormone. Dysregulation of the HPA axis and increased levels of cortisol are commonly found in AD patients and make a major contribution to the disease process. The underlying mechanisms remain poorly understood. In addition, within the general population there are interindividual differences in sensitivities to glucocorticoid and stress responses, which are thought to be due to a combination of genetic and environmental factors. These differences could ultimately impact an individuals' risk of AD. The purpose of this review is first to summarise the literature describing environmental and genetic factors that can impact an individual's HPA axis reactivity and function and ultimately AD risk. Secondly, we propose a mechanism by which genetic factors that influence HPA axis reactivity may also impact inflammation, a key driver of neurodegeneration. We hypothesize that these factors can mediate glucocorticoid priming of the immune cells of the brain, microglia, to become pro-inflammatory and promote a neurotoxic environment resulting in neurodegeneration. Understanding the underlying molecular mechanisms and identifying these genetic factors has implications for evaluating stress-related risk/progression to neurodegeneration, informing the success of interventions based on stress management and potential risks associated with the common use of glucocorticoids.

The subfamily Acastinae contains a diverse group of barnacles that are obligate symbionts of sponges and alcyonacean and antipatharian corals. Integrating morphological and genetic (COI) data to compare against known species, this paper reports on nine species of sponge-inhabiting barnacles of the subfamily Acastinae, including three undescribed species (Acasta caveata sp. nov., Euacasta acutaflava sp. nov., and E. excoriatrix sp. nov.) and three species previously not recorded in Australian waters (A. sandwichi, Pectinoacasta cancellorum, and P. sculpturata). The new species are distinguished from similar species by a suite of morphological characters as well as genetic distances. A lectotype for Pectinoacasta cancellorum is designated. Sponge hosts were identified for all specimens where possible and are represented by 19 species from eight families and five orders.

Human cytomegalovirus (HCMV) is a beta-herpesvirus carried by ∼80% of the world's population. Acute infections are asymptomatic in healthy individuals but generate diverse syndromes in neonates, solid organ transplant recipients, and HIV-infected individuals. The HCMV gene US28 encodes a homolog of a human chemokine receptor that is able to bind several chemokines and HIV gp120. Deep sequencing technologies were used to sequence US28 directly from 60 clinical samples from Indonesian HIV patients and Australian renal transplant recipients, healthy adults, and neonates. Molecular modeling approaches were used to predict whether nine nonsynonymous mutations in US28 may alter protein binding to a panel of six chemokines and two variants of HIV gp120. Ninety-two percent of samples contained more than one variant of HCMV, as defined by at least one nonsynonymous mutation. Carriage of these variants differed between neonates and adults, Australian and Indonesian samples, and saliva samples and blood leukocytes. Two nonsynonymous mutations (N170D and R267K) were associated with increased levels of immediate early protein 1 (IE-1) and glycoprotein B (gB) HCMV-reactive antibodies, suggesting a higher viral burden. Seven of the nine mutations were predicted to alter binding of at least one ligand. Overall, HCMV variants are common in all populations and have the potential to affect US28 interactions with human chemokines and/or gp120 and alter responses to the virus. The findings relied on deep sequencing technologies applied directly to clinical samples, so the variants exist in vivo. IMPORTANCE Human cytomegalovirus (HCMV) is a common viral pathogen of solid organ transplant recipients, neonates, and HIV-infected individuals. HCMV encodes homologs of several host genes with the potential to influence viral persistence and/or pathogenesis. Here, we present deep sequencing of an HCMV chemokine receptor homolog, US28, acquired directly from clinical specimens. Carriage of these variants differed between patient groups and was associated with different levels of circulating HCMV-reactive antibodies. These features are consistent with a role for US28 in HCMV persistence and pathogenesis. This was supported by in silico analyses of the variant sequences demonstrating altered ligand-binding profiles. The data delineate a novel approach to understanding the pathogenesis of HCMV and may impact the development of an effective vaccine.

Chestnut coat color in horses is determined by a missense mutation within the MC1R gene. However, the intensity of the chestnut color can vary widely within individuals possessing this genotype. Here, we investigated this variation using standardized photographs of 96 horses. Each horse was ranked lightest to darkest within the cohort for phenotype by 3 blinded observers. A genome-wide association study utilizing the relative shade ranking as the phenotype and using 268 487 single-nucleotide polymorphisms (SNPs) genotyped using the Affymetrix Equine 670k array identified a single significantly associated region on chromosome 3 (P = 2.934 × 10-8). Analysis of whole-genome sequences for horses spanning the diverse range of chestnut color identified candidate SNPs within the coding sequence of the only gene in the region: SALL1. The function of SALL1 is largely unknown, though it is predicted to interact with the Hermansky-Pudlak Syndrome type 1 (HPS1) protein, which causes partial albinism in humans. However, with only one study suggesting a circumstantial influence of the SALL1 protein on pigmentation, additional work is needed to confirm this new coat color locus in larger populations and investigate the function of this protein for impacts on equine health.

Context Breech strike is a serious disease for wool sheep. Skin wrinkle and dags are known predisposing factors for breech strike; however, a large part of the variation among sheep is unknown. Aims We studied the natural diversity and difference in microbial populations in the skin around the breech area in Merino sheep genetically resistant and susceptible to breech strike, by using 16S rRNA gene sequence analysis. Methods The sheep were from the breech strike flocks at the Mount Barker research station in Western Australia and from the CSIRO research station near Armidale in New South Wales. Skin samples were collected from the breech of all 2013-born progeny in both flocks before they were struck. Yearling ewes and rams were then naturally exposed to challenge by Lucilia cuprina blowflies. Breeding values for breech strike were estimated and used with phenotypic data to identify breech strike-resistant and -susceptible sheep. Skin samples of 78 unstruck and 73 struck sheep were selected, their microbiomes were analysed using 16S rRNA meta-barcoding, and operational taxonomic unit counts were analysed. Results The diversity analyses showed that the two flocks in the different environments had different microbiome profiles, but no difference was found between sexes or between breech strike-resistant and -susceptible sheep in either flock. Conclusions The results indicated that microbial differences on the skin of sheep are not associated with differences in susceptibility to breech strike. Implications Microbial differences do not offer opportunities to manage breech strike in Merino sheep.

In order to characterize the phenotype and to examine the effects of sun exposure on the color and structure of eumelanin (EM) and pheomelanin (PM) in alpaca fibers, we applied Soluene-350 solubilization, alkaline hydrogen peroxide oxidation (AHPO) and hydroiodic acid (HI) hydrolysis to the base and tip fibers of 20 true-black (TB) and 20 warm-black (WB) alpacas. We analyzed absorbances at 500 nm (A500) and 650 nm (A650), Free and Total pyrrole-2,3,5-tricarboxylic acid (PTCA), 2,3,4,5-tetracarboxylic acid (PTeCA) as degradative products from EM, and 4-amino-3-hydroxyphenylalanine (4-AHP), 3-amino-4-hydroxyphenylalanine (3-AHP) and thiazole-2,4,5-tricarboxylic acid (TTCA) as degradative products from PM. We found that the ratio of PTeCA/Total PTCA increased significantly from the base to the tip in both colors of alpaca fibers, while the ratios of A650/A500 and 4-AHP/3-AHP decreased significantly. These results show that structures made of both EM and PM in alpaca fibers are modified significantly by sun exposure inducing color change. This study indicates that the ratios of A650/A500, PTeCA/Total PTCA and 4-AHP/3-AHP are highly sensitive markers of color change and photodegradation of EM and PM, respectively.

Sponges represent one of the most species-rich hosts for commensal barnacles yet host utilisation and diversity have not been thoroughly examined. This study investigated the diversity and phylogenetic relationships of sponge-inhabiting barnacles within a single, targeted host group, primarily from Western Australian waters. Specimens of the sponge order Dictyoceratida were surveyed and a total of 64 host morphospecies, representing four families, were identified as barnacle hosts during the study. Utilising molecular (COI, 12S) and morphological methods 42 molecular operational taxonomic units (MOTUs) of barnacles, representing Acasta, Archiacasta, Euacasta and Neoacasta were identified. Comparing inter- and intra-MOTU genetic distances showed a barcode gap between 2.5% and 5% for COI, but between 1% and 1.5% in the 12S dataset, thus demonstrating COI as a more reliable barcoding region. These sponge-inhabiting barnacles were demonstrated to show high levels of host specificity with the majority being found in a single sponge species (74%), a single genus (83%) or a single host family (93%). Phylogenetic relationships among the barnacles were reconstructed using mitochondrial (12S, COI) and nuclear (H3, 28S) markers. None of the barnacle genera were recovered as monophyletic. Euacasta was paraphyletic in relation to the remaining Acastinae genera, which were polyphyletic. Six well-supported clades of molecular operational taxonomic units, herein considered to represent species complexes, were recovered, but relationships between them were not well supported. These complexes showed differing patterns of host usage, though most were phylogenetically conserved with sister lineages typically occupying related hosts within the same genus or family of sponge. The results show that host specialists are predominant, and the dynamics of host usage have played a significant role in the evolutionary history of the Acastinae.

Additional file 5 FR-AgENCODE genes and transcripts (structure, expression, positional and coding classes).• bos_taurus_cuff_tpm0.1_2sample_complete.gff• bos_taurus_cuff_tpm0.1_2sample_trid_4posclasses_3codingclasses_booleans.tsv• bos_taurus.frag.gnid.posclasslist.codclasslist.tsv• bos_taurus.fraggn.tpm.tsv• capra_hircus_cuff_tpm0.1_2sample_complete.gff• capra_hircus_cuff_tpm0.1_2sample_trid_4posclasses_3codingclasses_booleans.tsv• capra_hircus.frag.gnid.posclasslist.codclasslist.tsv• capra_hircus.fraggn.tpm.tsv• gallus_gallus_cuff_tpm0.1_2sample_complete.gff• gallus_gallus_cuff_tpm0.1_2sample_trid_4posclasses_3codingclasses_booleans.tsv• gallus_gallus.frag.gnid.posclasslist.codclasslist.tsv• gallus_gallus.fraggn.tpm.tsv• sus_scrofa_cuff_tpm0.1_2sample_complete.gff• sus_scrofa_cuff_tpm0.1_2sample_trid_4posclasses_3codingclasses_booleans.tsv• sus_scrofa.frag.gnid.posclasslist.codclasslist.tsv• sus_scrofa.fraggn.tpm.tsv

Human cytomegalovirus (HCMV) infections are common following renal transplantation and may have long-lasting effects. HCMV can be measured directly by viral DNA or indirectly via host immune responses. HCMV-encoded microRNA (miRNA) may alter the pathobiology of HCMV infections and contribute to the progression of HCMV disease. HCMV-encoded miRNAs can be detected in blood but have not been sought in saliva. We investigated saliva samples from 32 renal transplant recipients (RTR) and 12 seropositive healthy controls for whom immunological data was available. Five HCMV-encoded miRNAs (miR-UL112-5p, miR-US5-2-3p, miR-UL36, miR-US25-2-3p and miR-UL22A) were sought using primer probe assays. HCMV miRNA species were detected in saliva from 15 RTR and 3 healthy controls, with miR-US5-2-3p most commonly detected. The presence of HCMV miRNAs associated with increased T-cell responses to HCMV IE-1 in RTR, suggesting a link with frequent reactivations of HCMV.

BACKGROUND:The Djallonke sheep is well adapted to harsh environmental conditions, and is relatively resistant to Haemonchosis and resilient to animal trypanosomiasis. The larger Sahelian sheep, which cohabit the same region, is less well adapted to these disease challenges. Haemonchosis and Trypanosomiasis collectively cost the worldwide animal industry billions of dollars in production losses annually.RESULTS:Here, we separately sequenced and then pooled according to breed the genomes from five unrelated individuals from each of the Djallonke and Sahelian sheep breeds (sourced from Ghana), at greater than 22-fold combined coverage for each breed. A total of approximately 404 million (97%) and 343 million (97%) sequence reads from the Djallonke and Sahelian breeds respectively, were successfully mapped to the sheep reference genome Oar v3.1. We identified approximately 11.1 million and 10.9 million single nucleotide polymorphisms (SNPs) in the Djallonke and Sahelian breeds, with approximately 15 and 16% respectively of these not previously reported in sheep. Multiple regions of reduced heterozygosity were also found; 70 co-localised within genomic regions harbouring genes that mediate disease resistance, immune response and adaptation in sheep or cattle. Thirty- three of the regions of reduced heterozygosity co-localised with previously reported genes for resistance to haemonchosis and trypanosomiasis.CONCLUSIONS:Our analyses suggest that these regions of reduced heterozygosity may be signatures of selection for these economically important diseases.

Abstract Improving the functional annotation of animal genomes is a key challenge in bridging the gap between genotype and phenotype, thus enabling predictive biology. Regarding livestock production, major outcomes are expected from a better understanding of the genetic architecture underlying quantitative traits. As part of the Functional Annotation of ANimal Genomes action (FAANG: www.faang.org), the FR-AgENCODE project generated omics data to improve the reference annotation of the cattle, pig, goat and chicken genome. High-throughput molecular assays have been performed on tissues/cells relevant to immune and metabolic traits. From two males and two females per species (pig, cattle, goat, chicken), strand-oriented RNA-seq gene expression and ATAC-seq chromatin accessibility assays were performed on liver and two PBMC-sorted T-cell types (CD4+ and CD8+). Chromosome Conformation Capture (in situ Hi-C) was also carried out on liver samples. About 4,000 samples have been collected at the INRA biorepository and registered at the EBI BioSamples registry. More than 80% of the planned experiments could be completed, generating ~11.5 billions of sequencing reads over the 3 assays. While most (50–80%) RNA-seq reads mapped to annotated exons, thousands of novel transcripts were found, with ~60K mRNAs and ~22K lncRNAs in cattle. Differentially expressed genes between cell types were enriched for immunity- or metabolism-related terms, and differentially accessible chromatin regions were identified as potential regulatory sites. Interestingly, correlations between gene expression and promoter accessibility across samples were skewed towards both positive and negative values, suggesting distinct regulatory mechanisms of gene expression. These patterns have been further investigated using human data from the Epigenome Roadmap Mapping Consortium. Altogether, this study illustrates the interest of a coordinated effort to tackle the genome-to-phenome challenge and provides a useful resource to the community. Availability: www.fragencode.org.

Additional file 2 Reference genes and transcripts (structure, expression) of the 4 species. Archive content: • bos_taurus.gtf • bos_taurus.refgn.tpm.tsv • capra_hircus.gtf • capra_hircus.refgn.tpm.tsv • gallus_gallus.gtf • gallus_gallus.refgn.tpm.tsv • sus_scrofa.gtf •sus_scrofa.refgn.tpm.tsv

The alpaca classic grey phenotype is of particular interest to the industry. Until now, there were only indirect data suggesting that the KIT gene was involved in the classic grey phenotype. All exons of KIT in three black and three classic silvergrey alpacas were sequenced. Five non-synonymous SNPs were observed. There was only one SNP found that was present only in the silvergrey alpacas, and this was also the only SNP predicted to be damaging. This variant results in a change of a glycine (Gly) to an arginine (Arg) at amino acid position 126 (c.376G>A), occurring in the second Ig-like domain of the extracellular domain of KIT. Basic protein modelling predicted that this variant is likely destabilising. Therefore, an additional 488 alpacas were genotyped for this SNP using the tetra-primer amplification refractory mutation system PCR (Tetra-primer ARMS-PCR). All classic grey alpacas were observed to be heterozygous, and 99.3% of non-grey dark base colour alpacas were found to be homozygous for the wildtype allele in this position. These results confirm that the classic grey phenotype in alpacas is the result of a c.376G>A (p.Gly126Arg) SNP in exon 3 of KIT. These data also support the hypothesis that the grey phenotype is autosomal dominant and that the mutation is most likely homozygous lethal.