GWAS reveals determinants of mobilization rate and dynamics of an active endogenous retrovirus of cattle

Five to ten percent of mammalian genomes is occupied by multiple clades of endogenous retroviruses (ERVs), that may count thousands of members. New ERV clades arise by retroviral infection of the germline followed by expansion by reinfection and/or retrotransposition. ERV mobilization is a source of deleterious variation, driving the emergence of ERV silencing mechanisms, leaving "DNA fossils". Here we show that the ERVK[2-1-LTR] clade is still active in the bovine and a source of disease-causing alleles. We develop a method to measure the rate of ERVK[2-1-LTR] mobilization, finding an average of 1 per ~150 sperm cells, with >10-fold difference between animals. We perform a genome-wide association study and identify eight loci affecting ERVK[2-1-LTR] mobilization. We provide evidence that polymorphic ERVK[2-1-LTR] elements in four of these loci cause the association. We generate a catalogue of full length ERVK[2-1-LTR] elements, and show that it comprises 15% of C-type autonomous elements, and 85% of D-type non-autonomous elements lacking functional genes. We show that >25% of the variance of mobilization rate is determined by the number of C-type elements, yet that de novo insertions are dominated by D-type elements. We propose that D-type elements act as parasite-of-parasite gene drives that may contribute to the observed demise of ERV elements.

The rate of de novo structural variation is increased in in vitro-produced offspring and preferentially affects the paternal genome

Assisted reproductive technologies (ARTs), including in vitro maturation and fertilization (IVF), are increasingly used in human and animal reproduction. Whether these technologies directly affect the rate of de novo mutation (DNM), and to what extent, has been a matter of debate. Here we take advantage of domestic cattle, characterized by complex pedigrees that are ideally suited to detect DNMs and by the systematic use of ART, to study the rate of de novo structural variation (dnSV) in this species and how it is impacted by IVF. By exploiting features of associated de novo point mutations (dnPMs) and dnSVs in clustered DNMs, we provide strong evidence that (1) IVF increases the rate of dnSV approximately fivefold, and (2) the corresponding mutations occur during the very early stages of embryonic development (one- and two-cell stage), yet primarily affect the paternal genome.

High-resolution structural variants catalogue in a large-scale whole genome sequenced bovine family cohort data

BACKGROUND

Structural variants (SVs) are chromosomal segments that differ between genomes, such as deletions, duplications, insertions, inversions and translocations. The genomics revolution enabled the discovery of sub-microscopic SVs via array and whole-genome sequencing (WGS) data, paving the way to unravel the functional impact of SVs. Recent human expression QTL mapping studies demonstrated that SVs play a disproportionally large role in altering gene expression, underlining the importance of including SVs in genetic analyses. Therefore, this study aimed to generate and explore a high-quality bovine SV catalogue exploiting a unique cattle family cohort data (total 266 samples, forming 127 trios).

RESULTS

We curated 13,731 SVs segregating in the population, consisting of 12,201 deletions, 1,509 duplications, and 21 multi-allelic CNVs (> 50-bp). Of these, we validated a subset of copy number variants (CNVs) utilising a direct genotyping approach in an independent cohort, indicating that at least 62% of the CNVs are true variants, segregating in the population. Among gene-disrupting SVs, we prioritised two likely high impact duplications, encompassing ORM1 and POPDC3 genes, respectively. Liver expression QTL mapping results revealed that these duplications are likely causing altered gene expression, confirming the functional importance of SVs. Although most of the accurately genotyped CNVs are tagged by single nucleotide polymorphisms (SNPs) ascertained in WGS data, most CNVs were not captured by individual SNPs obtained from a 50K genotyping array.

CONCLUSION

We generated a high-quality SV catalogue exploiting unique whole genome sequenced bovine family cohort data. Two high impact duplications upregulating the ORM1 and POPDC3 are putative candidates for postpartum feed intake and hoof health traits, thus warranting further investigation. Generally, CNVs were in low LD with SNPs on the 50K array. Hence, it remains crucial to incorporate CNVs via means other than tagging SNPs, such as investigation of tagging haplotypes, direct imputation of CNVs, or direct genotyping as done in the current study. The SV catalogue and the custom genotyping array generated in the current study will serve as valuable resources accelerating utilisation of full spectrum of genetic variants in bovine genomes.

Transcriptome analysis reveals tumor microenvironment changes in glioblastoma

A better understanding of transcriptional evolution of IDH-wild-type glioblastoma may be crucial for treatment optimization. Here, we perform RNA sequencing (RNA-seq) (n = 322 test, n = 245 validation) on paired primary-recurrent glioblastoma resections of patients treated with the current standard of care. Transcriptional subtypes form an interconnected continuum in a two-dimensional space. Recurrent tumors show preferential mesenchymal progression. Over time, hallmark glioblastoma genes are not significantly altered. Instead, tumor purity decreases over time and is accompanied by co-increases in neuron and oligodendrocyte marker genes and, independently, tumor-associated macrophages. A decrease is observed in endothelial marker genes. These composition changes are confirmed by single-cell RNA-seq and immunohistochemistry. An extracellular matrix-associated gene set increases at recurrence and bulk, single-cell RNA, and immunohistochemistry indicate it is expressed mainly by pericytes. This signature is associated with significantly worse survival at recurrence. Our data demonstrate that glioblastomas evolve mainly by microenvironment (re-)organization rather than molecular evolution of tumor cells.

Comparing the transcriptome of developing native and iPSC-derived mouse retinae by single cell RNA sequencing

We report the generation and analysis of single-cell RNA-Seq data (> 38,000 cells) from mouse native retinae and induced pluripotent stem cell (iPSC)-derived retinal organoids at four matched stages of development spanning the emergence of the major retinal cell types. We combine information from temporal sampling, visualization of 3D UMAP manifolds, pseudo-time and RNA velocity analyses, to show that iPSC-derived 3D retinal organoids broadly recapitulate the native developmental trajectories. However, we observe relaxation of spatial and temporal transcriptome control, premature emergence and dominance of photoreceptor precursor cells, and susceptibility of dynamically regulated pathways and transcription factors to culture conditions in retinal organoids. We demonstrate that genes causing human retinopathies are enriched in cell-type specifying genes and identify a subset of disease-causing genes with expression profiles that are highly conserved between human retinae and murine retinal organoids. This study provides a resource to the community that will be useful to assess and further improve protocols for ex vivo recapitulation and study of retinal development.

The virulome of Streptomyces scabiei in response to cello-oligosaccharide elicitors

The development of spots or lesions symptomatic of common scab on root and tuber crops is caused by few pathogenic Streptomyces with Streptomyces scabiei 87–22 as the model species. Thaxtomin phytotoxins are the primary virulence determinants, mainly acting by impairing cellulose synthesis, and their production in S. scabiei is in turn boosted by cello-oligosaccharides released from host plants. In this work we aimed to determine which molecules and which biosynthetic gene clusters (BGCs) of the specialized metabolism of S. scabiei 87–22 show a production and/or a transcriptional response to cello-oligosaccharides. Comparative metabolomic analyses revealed that molecules of the virulome of S. scabiei induced by cellobiose and cellotriose include (i) thaxtomin and concanamycin phytotoxins, (ii) desferrioxamines, scabichelin and turgichelin siderophores in order to acquire iron essential for housekeeping functions, (iii) ectoine for protection against osmotic shock once inside the host, and (iv) bottromycin and concanamycin antimicrobials possibly to prevent other microorganisms from colonizing the same niche. Importantly, both cello-oligosaccharides reduced the production of the spore germination inhibitors germicidins thereby giving the ‘green light’ to escape dormancy and trigger the onset of the pathogenic lifestyle. For most metabolites - either with induced or reduced production - cellotriose was revealed to be a slightly stronger elicitor compared to cellobiose, supporting an earlier hypothesis which suggested the trisaccharide was the real trigger for virulence released from the plant cell wall through the action of thaxtomins. Interestingly, except for thaxtomins, none of these BGCs’ expression seems to be under direct control of the cellulose utilization repressor CebR suggesting the existence of a yet unknown mechanism for switching on the virulome. Finally, a transcriptomic analysis revealed nine additional cryptic BGCs that have their expression awakened by cello-oligosaccharides, suggesting that other and yet to be discovered metabolites could be part of the virulome of S. scabiei.

A 12 kb multi-allelic copy number variation encompassing a GC gene enhancer is associated with mastitis resistance in dairy cattle

Clinical mastitis (CM) is an inflammatory disease occurring in the mammary glands of lactating cows. CM is under genetic control, and a prominent CM resistance QTL located on chromosome 6 was reported in various dairy cattle breeds. Nevertheless, the biological mechanism underpinning this QTL has been lacking. Herein, we mapped, fine-mapped, and discovered the putative causal variant underlying this CM resistance QTL in the Dutch dairy cattle population. We identified a ~12 kb multi-allelic copy number variant (CNV), that is in perfect linkage disequilibrium with a lead SNP, as a promising candidate variant. By implementing a fine-mapping and through expression QTL mapping, we showed that the group-specific component gene (GC), a gene encoding a vitamin D binding protein, is an excellent candidate causal gene for the QTL. The multiplicated alleles are associated with increased GC expression and low CM resistance. Ample evidence from functional genomics data supports the presence of an enhancer within this CNV, which would exert cis-regulatory effect on GC. We observed that strong positive selection swept the region near the CNV, and haplotypes associated with the multiplicated allele were strongly selected for. Moreover, the multiplicated allele showed pleiotropic effects for increased milk yield and reduced fertility, hinting that a shared underlying biology for these effects may revolve around the vitamin D pathway. These findings together suggest a putative causal variant of a CM resistance QTL, where a cis-regulatory element located within a CNV can alter gene expression and affect multiple economically important traits.

Single-cell RNA sequencing of human liver reveals hepatic stellate cell heterogeneity

Background & Aims

The multiple vital functions of the human liver are performed by highly specialised parenchymal and non-parenchymal cells organised in complex collaborative sinusoidal units. Although crucial for homeostasis, the cellular make-up of the human liver remains to be fully elucidated. Here, single-cell RNA-sequencing was used to unravel the heterogeneity of human liver cells, in particular of hepatocytes (HEPs) and hepatic stellate cells (HSCs).

Method

The transcriptome of ~25,000 freshly isolated human liver cells was profiled using droplet-based RNA-sequencing. Recently published data sets and RNA in situ hybridisation were integrated to validate and locate newly identified cell populations.

Results

In total, 22 cell populations were annotated that reflected the heterogeneity of human parenchymal and non-parenchymal liver cells. More than 20,000 HEPs were ordered along the portocentral axis to confirm known, and reveal previously undescribed, zonated liver functions. The existence of 2 subpopulations of human HSCs with unique gene expression signatures and distinct intralobular localisation was revealed (i.e. portal and central vein-concentrated GPC3⁺ HSCs and perisinusoidally located DBH⁺ HSCs). In particular, these data suggest that, although both subpopulations collaborate in the production and organisation of extracellular matrix, GPC3⁺ HSCs specifically express genes involved in the metabolism of glycosaminoglycans, whereas DBH⁺ HSCs display a gene signature that is reminiscent of antigen-presenting cells.

Conclusions

This study highlights metabolic zonation as a key determinant of HEP transcriptomic heterogeneity and, for the first time, outlines the existence of heterogeneous HSC subpopulations in the human liver. These findings call for further research on the functional implications of liver cell heterogeneity in health and disease.

Lay summary

This study resolves the cellular landscape of the human liver in an unbiased manner and at high resolution to provide new insights into human liver cell biology. The results highlight the physiological heterogeneity of human hepatic stellate cells.

•

A cell atlas from the near-native transcriptome of >25,000 human liver cells is presented.

•

Hepatocytes were ordered along the portocentral axis to reveal previously undescribed gene expression patterns and zonated liver functions.

•

Two subpopulations of human hepatic stellate cells (HSCs) are reported, characterised by different spatial distribution in the native tissue.

•

Characteristic gene signatures of HSC subpopulations are suggestive of far-reaching functional differences.

SNP-based quantitative deconvolution of biological mixtures: application to the detection of cows with subclinical mastitis by whole-genome sequencing of tank milk

Biological products of importance in food (e.g., milk) and medical (e.g., donor blood-derived products) sciences often correspond to mixtures of samples contributed by multiple individuals. Identifying which individuals contributed to the mixture and in what proportions may be of interest in several circumstances. We herein present a method that allows to do this by shallow whole-genome sequencing of the DNA in mixed samples from hundreds of donors. We show the efficacy of the approach for the detection of cows with subclinical mastitis by analysis of farms' tank mixtures containing milk from as many as 500 cows.

Detection and identification of transgenic events by next generation sequencing combined with enrichment technologies

Next generation sequencing (NGS) is a promising tool for analysing the quality and safety of food and feed products. The detection and identification of genetically modified organisms (GMOs) is complex, as the diversity of transgenic events and types of structural elements introduced in plants continue to increase. In this paper, we show how a strategy that combines enrichment technologies with NGS can be used to detect a large panel of structural elements and partially or completely reconstruct the new sequence inserted into the plant genome in a single analysis, even at low GMO percentages. The strategy of enriching sequences of interest makes the approach applicable even to mixed products, which was not possible before due to insufficient coverage of the different genomes present. This approach is also the first step towards a more complete characterisation of agrifood products in a single analysis.

Analysis of Genes Associated With Monogenic Primary Immunodeficiency Identifies Rare Variants in XIAP in Patients With Crohn's Disease

BACKGROUND & AIMS

A few rare monogenic primary immunodeficiencies (PIDs) are characterized by chronic intestinal inflammation that resembles Crohn's disease (CD). We investigated whether 23 genes associated with 10 of these monogenic disorders contain common, low-frequency, or rare variants that increase risk for CD.

METHODS

Common and low frequency variants in 1 Mb loci centered on the candidate genes were analyzed using meta-data corresponding to genotypes of approximately 17,000 patients with CD or without CD (controls) in Europe. The contribution of rare variants was assessed by high-throughput sequencing of 4750 individuals, including 660 early-onset and/or familial cases among the 2390 patients with CD. Variants were expressed from vectors in SW480 or HeLa cells and functions of their products were analyzed in immunofluorescence, luciferase, immunoprecipitation, and immunoblot assays.

RESULTS

We reproduced the association of the interleukin 10 locus with CD (P = .007), although none of the significantly associated variants modified the coding sequence of interleukin 10. We found XIAP to be significantly enriched for rare coding mutations in patients with CD vs controls (P = .02). We identified 4 previously unreported missense variants associated with CD. Variants in XIAP cause the PID X-linked lymphoproliferative disease type 2, yet none of the carriers of these variants had all the clinical features of X-linked lymphoproliferative disease type 2. Identified XIAP variants S123N, R233Q, and P257A were associated with an impaired activation of NOD2 signaling after muramyl dipeptide stimulation.

CONCLUSIONS

In a systematic analysis of variants in 23 PID-associated genes, we confirmed the association of variants in XIAP with CD. Further screenings for CD-associated variants and analyses of their functions could increase our understanding of the relationship between PID-associated genes and CD pathogenesis.

Isolation, Characterization, and Antibacterial Activity of Hard-to-Culture Actinobacteria from Cave Moonmilk Deposits

Cave moonmilk deposits host an abundant and diverse actinobacterial population that has a great potential for producing novel natural bioactive compounds. In our previous attempt to isolate culturable moonmilk-dwelling Actinobacteria, only Streptomyces species were recovered, whereas a metagenetic study of the same deposits revealed a complex actinobacterial community including 46 actinobacterial genera in addition to streptomycetes. In this work, we applied the rehydration-centrifugation method to lessen the occurrence of filamentous species and tested a series of strategies to achieve the isolation of hard-to-culture and rare Actinobacteria from the moonmilk deposits of the cave “Grotte des Collemboles”. From the “tips and tricks” that were tested, separate autoclaving of the components of the International Streptomyces Project (ISP) medium number 5 (ISP5) medium, prolonged incubation time, and dilution of the moonmilk suspension were found to most effectively improve colony forming units. Taxonomic analyses of the 40 isolates revealed new representatives of the Agromyces, Amycolatopsis, Kocuria, Micrococcus, Micromonospora, Nocardia, and Rhodococcus species, as well as additional new streptomycetes. The applied methodologies allowed the isolation of strains associated with both the least and most abundant moonmilk-dwelling actinobacterial operational taxonomic units. Finally, bioactivity screenings revealed that some isolates displayed high antibacterial activities, and genome mining uncovered a strong potential for the production of natural compounds.

Transcriptome-wide analysis of natural antisense transcripts shows their potential role in breast cancer

Non-coding RNAs (ncRNA) represent 1/5 of the mammalian transcript number, and 90% of the genome length is transcribed. Many ncRNAs play a role in cancer. Among them, non-coding natural antisense transcripts (ncNAT) are RNA sequences that are complementary and overlapping to those of either protein-coding (PCT) or non-coding transcripts. Several ncNATs were described as regulating protein coding gene expression on the same loci, and they are expected to act more frequently in cis compared to other ncRNAs that commonly function in trans. In this work, 22 breast cancers expressing estrogen receptors and their paired adjacent non-malignant tissues were analyzed by strand-specific RNA sequencing. To highlight ncNATs potentially playing a role in protein coding gene regulations that occur in breast cancer, three different data analysis methods were used: differential expression analysis of ncNATs between tumor and non-malignant tissues, differential correlation analysis of paired ncNAT/PCT between tumor and non-malignant tissues, and ncNAT/PCT read count ratio variation between tumor and non-malignant tissues. Each of these methods yielded lists of ncNAT/PCT pairs that were enriched in survival-associated genes. This work highlights ncNAT lists that display potential to affect the expression of protein-coding genes involved in breast cancer pathology.

Population-based cohort study of variation in the use of emergency cholecystectomy for benign gallbladder diseases

Background

The aims of this prospective population-based cohort study were to identify the patient and hospital characteristics associated with emergency cholecystectomy, and the influences of these in determining variations between hospitals.

Methods

Data were collected for consecutive patients undergoing cholecystectomy in acute UK and Irish hospitals between 1 March and 1 May 2014. Potential explanatory variables influencing the performance of emergency cholecystectomy were analysed by means of multilevel, multivariable logistic regression modelling using a two-level hierarchical structure with patients (level 1) nested within hospitals (level 2).

Results

Data were collected on 4744 cholecystectomies from 165 hospitals. Increasing age, lower ASA fitness grade, biliary colic, the need for further imaging (magnetic retrograde cholangiopancreatography), endoscopic interventions (endoscopic retrograde cholangiopancreatography) and admission to a non-biliary centre significantly reduced the likelihood of an emergency cholecystectomy being performed. The multilevel model was used to calculate the probability of receiving an emergency cholecystectomy for a woman aged 40 years or over with an ASA grade of I or II and a BMI of at least 25·0 kg/m2, who presented with acute cholecystitis with an ultrasound scan showing a thick-walled gallbladder and a normal common bile duct. The mean predicted probability of receiving an emergency cholecystectomy was 0·52 (95 per cent c.i. 0·45 to 0·57). The predicted probabilities ranged from 0·02 to 0·95 across the 165 hospitals, demonstrating significant variation between hospitals.

Conclusion

Patients with similar characteristics presenting to different hospitals with acute gallbladder pathology do not receive comparable care.

NGS-based reverse genetic screen for common embryonic lethal mutations compromising fertility in livestock

We herein report the result of a large-scale, next generation sequencing (NGS)-based screen for embryonic lethal (EL) mutations in Belgian beef and New Zealand dairy cattle. We estimated by simulation that cattle might carry, on average, ∼0.5 recessive EL mutations. We mined exome sequence data from >600 animals, and identified 1377 stop-gain, 3139 frame-shift, 1341 splice-site, 22,939 disruptive missense, 62,399 benign missense, and 92,163 synonymous variants. We show that cattle have a comparable load of loss-of-function (LoF) variants (defined as stop-gain, frame-shift, or splice-site variants) as humans despite having a more variable exome. We genotyped >40,000 animals for up to 296 LoF and 3483 disruptive missense, breed-specific variants. We identified candidate EL mutations based on the observation of a significant depletion in homozygotes. We estimated the proportion of EL mutations at 15% of tested LoF and 6% of tested disruptive missense variants. We confirmed the EL nature of nine candidate variants by genotyping 200 carrier × carrier trios, and demonstrating the absence of homozygous offspring. The nine identified EL mutations segregate at frequencies ranging from 1.2% to 6.6% in the studied populations and collectively account for the mortality of ∼0.6% of conceptuses. We show that EL mutations preferentially affect gene products fulfilling basic cellular functions. The resulting information will be useful to avoid at-risk matings, thereby improving fertility.

A stop-gain in the laminin, alpha 3 gene causes recessive junctional epidermolysis bullosa in Belgian Blue cattle

Four newborn purebred Belgian Blue calves presenting a severe form of epidermolysis bullosa were recently referred to our heredo-surveillance platform. SNP array genotyping followed by autozygosity mapping located the causative gene in a 8.3-Mb interval on bovine chromosome 24. Combining information from (i) whole-genome sequencing of an affected calf, (ii) transcriptomic data from a panel of tissues and (iii) a list of functionally ranked positional candidates pinpointed a private G to A nucleotide substitution in the LAMA3 gene that creates a premature stop codon (p.Arg2609*) in exon 60, truncating 22% of the corresponding protein. The LAMA3 gene encodes the alpha 3 subunit of the heterotrimeric laminin-332, a key constituent of the lamina lucida that is part of the skin basement membrane connecting epidermis and dermis layers. Homozygous loss-of-function mutations in this gene are known to cause severe junctional epidermolysis bullosa in human, mice, horse, sheep and dog. Overall, our data strongly support the causality of the identified gene and mutation.

Genome-wide next-generation DNA and RNA sequencing reveals a mutation that perturbs splicing of the phosphatidylinositol glycan anchor biosynthesis class H gene (PIGH) and causes arthrogryposis in Belgian Blue cattle

BACKGROUND

Cattle populations are characterized by regular outburst of genetic defects as a result of the extensive use of elite sires. The causative genes and mutations can nowadays be rapidly identified by means of genome-wide association studies combined with next generation DNA sequencing, provided that the causative mutations are conventional loss-of-function variants. We show in this work how the combined use of next generation DNA and RNA sequencing allows for the rapid identification of otherwise difficult to identify splice-site variants.

RESULTS

We report the use of haplotype-based association mapping to identify a locus on bovine chromosome 10 that underlies autosomal recessive arthrogryposis in Belgian Blue Cattle. We identify 31 candidate mutations by resequencing the genome of four cases and 15 controls at ~10-fold depth. By analyzing RNA-Seq data from a carrier fetus, we observe skipping of the second exon of the PIGH gene, which we confirm by RT-PCR to be fully penetrant in tissues from affected calves. We identify - amongst the 31 candidate variants - a C-to-G transversion in the first intron of the PIGH gene (c211-10C > G) that is predicted to affect its acceptor splice-site. The resulting PIGH protein is likely to be non-functional as it lacks essential domains, and hence to cause arthrogryposis.

CONCLUSIONS

This work illustrates how the growing arsenal of genome exploration tools continues to accelerate the identification of an even broader range of disease causing mutations, therefore improving the management and control of genetic defects in livestock.

Early-onset type 2 diabetes impairs skeletal acquisition in the male TALLYHO/JngJ mouse

Type 2 diabetes (T2D) incidence in adolescents is rising and may interfere with peak bone mass acquisition. We tested the effects of early-onset T2D on bone mass, microarchitecture, and strength in the TALLYHO/JngJ mouse, which develops T2D by 8 weeks of age. We assessed metabolism and skeletal acquisition in male TALLYHO/JngJ and SWR/J controls (n = 8-10/group) from 4 weeks to 8 and 17 weeks of age. Tallyho mice were obese; had an approximately 2-fold higher leptin and percentage body fat; and had lower bone mineral density vs SWR at all time points (P < .03 for all). Tallyho had severe deficits in distal femur trabecular bone volume fraction (-54%), trabecular number (-27%), and connectivity density (-82%) (P < .01 for all). Bone formation was higher in Tallyho mice at 8 weeks but lower by 17 weeks of age vs SWR despite similar numbers of osteoblasts. Bone marrow adiposity was 7- to 50-fold higher in Tallyho vs SWR. In vitro, primary bone marrow stromal cell differentiation into osteoblast and adipocyte lineages was similar in SWR and Tallyho, suggesting skeletal deficits were not due to intrinsic defects in Tallyho bone-forming cells. These data suggest the Tallyho mouse might be a useful model to study the skeletal effects of adolescent T2D.

A 660-Kb deletion with antagonistic effects on fertility and milk production segregates at high frequency in Nordic Red cattle: additional evidence for the common occurrence of balancing selection in livestock

In dairy cattle, the widespread use of artificial insemination has resulted in increased selection intensity, which has led to spectacular increase in productivity. However, cow fertility has concomitantly severely declined. It is generally assumed that this reduction is primarily due to the negative energy balance of high-producing cows at the peak of lactation. We herein describe the fine-mapping of a major fertility QTL in Nordic Red cattle, and identify a 660-kb deletion encompassing four genes as the causative variant. We show that the deletion is a recessive embryonically lethal mutation. This probably results from the loss of RNASEH2B, which is known to cause embryonic death in mice. Despite its dramatic effect on fertility, 13%, 23% and 32% of the animals carry the deletion in Danish, Swedish and Finnish Red Cattle, respectively. To explain this, we searched for favorable effects on other traits and found that the deletion has strong positive effects on milk yield. This study demonstrates that embryonic lethal mutations account for a non-negligible fraction of the decline in fertility of domestic cattle, and that associated positive effects on milk yield may account for part of the negative genetic correlation. Our study adds to the evidence that structural variants contribute to animal phenotypic variation, and that balancing selection might be more common in livestock species than previously appreciated.

We report the identification of a large deletion encompassing four genes and the demonstration of its negative effect on fertility in Nordic Red dairy cattle. We show that this deletion is recessively lethal (homozygous embryos die) and therefore, when carrier cows are mated to carrier bulls, there is a high risk of embryonic mortality. As a result, chances of insemination failure are higher for such matings. Surprisingly, despite its negative effect, the deletion is frequent in Nordic Red cattle. We show that this high frequency may be a consequence of the fact that the deletion is associated with increased milk production and therefore selected for. Due to increased levels of inbreeding resulting from the widespread use of artificial insemination, such recessive lethal alleles may account for a non-negligible fraction of the reduction in fertility observed in cattle.

A genome scan conducted in a multigenerational pedigree with convergent strabismus supports a complex genetic determinism

A genome-wide linkage scan was conducted in a Northern-European multigenerational pedigree with nine of 40 related members affected with concomitant strabismus. Twenty-seven members of the pedigree including all affected individuals were genotyped using a SNP array interrogating > 300,000 common SNPs. We conducted parametric and non-parametric linkage analyses assuming segregation of an autosomal dominant mutation, yet allowing for incomplete penetrance and phenocopies. We detected two chromosome regions with near-suggestive evidence for linkage, respectively on chromosomes 8 and 18. The chromosome 8 linkage implied a penetrance of 0.80 and a rate of phenocopy of 0.11, while the chromosome 18 linkage implied a penetrance of 0.64 and a rate of phenocopy of 0. Our analysis excludes a simple genetic determinism of strabismus in this pedigree.

A deletion in the bovine FANCI gene compromises fertility by causing fetal death and brachyspina

Fertility is one of the most important traits in dairy cattle, and has been steadily declining over the last decades. We herein use state-of-the-art genomic tools, including high-throughput SNP genotyping and next-generation sequencing, to identify a 3.3 Kb deletion in the FANCI gene causing the brachyspina syndrome (BS), a rare recessive genetic defect in Holstein dairy cattle. We determine that despite the very low incidence of BS (<1/100,000), carrier frequency is as high as 7.4% in the Holstein breed. We demonstrate that this apparent discrepancy is likely due to the fact that a large proportion of homozygous mutant calves die during pregnancy. We postulate that several other embryonic lethals may segregate in livestock and significantly compromise fertility, and propose a genotype-driven screening strategy to detect the corresponding deleterious mutations.

Serial translocation by means of circular intermediates underlies colour sidedness in cattle

Colour sidedness is a dominantly inherited phenotype of cattle characterized by the polarization of pigmented sectors on the flanks, snout and ear tips. It is also referred to as 'lineback' or 'witrik' (which means white back), as colour-sided animals typically display a white band along their spine. Colour sidedness is documented at least since the Middle Ages and is presently segregating in several cattle breeds around the globe, including in Belgian blue and brown Swiss. Here we report that colour sidedness is determined by a first allele on chromosome 29 (Cs(29)), which results from the translocation of a 492-kilobase chromosome 6 segment encompassing KIT to chromosome 29, and a second allele on chromosome 6 (Cs(6)), derived from the first by repatriation of fused 575-kilobase chromosome 6 and 29 sequences to the KIT locus. We provide evidence that both translocation events involved circular intermediates. This is the first example, to our knowledge, of a phenotype determined by homologous yet non-syntenic alleles that result from a novel copy-number-variant-generating mechanism.

Variants modulating the expression of a chromosome domain encompassing PLAG1 influence bovine stature

We report mapping of a quantitative trait locus (QTL) with a major effect on bovine stature to a ∼780-kb interval using a Hidden Markov Model-based approach that simultaneously exploits linkage and linkage disequilibrium. We re-sequenced the interval in six sires with known QTL genotype and identified 13 clustered candidate quantitative trait nucleotides (QTNs) out of >9,572 discovered variants. We eliminated five candidate QTNs by studying the phenotypic effect of a recombinant haplotype identified in a breed diversity panel. We show that the QTL influences fetal expression of seven of the nine genes mapping to the ∼780-kb interval. We further show that two of the eight candidate QTNs, mapping to the PLAG1-CHCHD7 intergenic region, influence bidirectional promoter strength and affect binding of nuclear factors. By performing expression QTL analyses, we identified a splice site variant in CHCHD7 and exploited this naturally occurring null allele to exclude CHCHD7 as single causative gene.

Genetic and functional confirmation of the causality of the DGAT1 K232A quantitative trait nucleotide in affecting milk yield and composition

We recently used a positional cloning approach to identify a nonconservative lysine to alanine substitution (K232A) in the bovine DGAT1 gene that was proposed to be the causative quantitative trait nucleotide underlying a quantitative trait locus (QTL) affecting milk fat composition, previously mapped to the centromeric end of bovine chromosome 14. We herein generate genetic and functional data that confirm the causality of the DGAT1 K232A mutation. We have constructed a high-density single-nucleotide polymorphism map of the 3.8-centimorgan BULGE30-BULGE9 interval containing the QTL and show that the association with milk fat percentage maximizes at the DGAT1 gene. We provide evidence that the K allele has undergone a selective sweep. By using a baculovirus expression system, we have expressed both DGAT1 alleles in Sf9 cells and show that the K allele, causing an increase in milk fat percentage in the live animal, is characterized by a higher Vmax in producing triglycerides than the A allele.

Molecular dissection of a quantitative trait locus: a phenylalanine-to-tyrosine substitution in the transmembrane domain of the bovine growth hormone receptor is associated with a major effect on milk yield and composition

We herein report on our efforts to improve the mapping resolution of a QTL with major effect on milk yield and composition that was previously mapped to bovine chromosome 20. By using a denser chromosome 20 marker map and by exploiting linkage disequilibrium using two distinct approaches, we provide strong evidence that a chromosome segment including the gene coding for the growth hormone receptor accounts for at least part of the chromosome 20 QTL effect. By sequencing individuals with known QTL genotype, we identify an F to Y substitution in the transmembrane domain of the growth hormone receptor gene that is associated with a strong effect on milk yield and composition in the general population.

Simultaneous mining of linkage and linkage disequilibrium to fine map quantitative trait loci in outbred half-sib pedigrees: revisiting the location of a quantitative trait locus with major effect on milk production on bovine chromosome 14

A maximum-likelihood QTL mapping method that simultaneously exploits linkage and linkage disequilibrium and that is applicable in outbred half-sib pedigrees is described. The method is applied to fine map a QTL with major effect on milk fat content in a 3-cM marker interval on proximal BTA14. This proximal location is confirmed by applying a haplotype-based association method referred to as recombinant ancestral haplotype analysis. The origin of the discrepancy between the QTL position derived in this work and that of a previous analysis is examined and shown to be due to the existence of distinct marker haplotypes associated with QTL alleles having large substitution effects.

Positional candidate cloning of a QTL in dairy cattle: identification of a missense mutation in the bovine DGAT1 gene with major effect on milk yield and composition

We recently mapped a quantitative trait locus (QTL) with a major effect on milk composition--particularly fat content--to the centromeric end of bovine chromosome 14. We subsequently exploited linkage disequilibrium to refine the map position of this QTL to a 3-cM chromosome interval bounded by microsatellite markers BULGE13 and BULGE09. We herein report the positional candidate cloning of this QTL, involving (1) the construction of a BAC contig spanning the corresponding marker interval, (2) the demonstration that a very strong candidate gene, acylCoA:diacylglycerol acyltransferase (DGAT1), maps to that contig, and (3) the identification of a nonconservative K232A substitution in the DGAT1 gene with a major effect on milk fat content and other milk characteristics.

Human-ovine comparative sequencing of a 250-kb imprinted domain encompassing the callipyge (clpg) locus and identification of six imprinted transcripts: DLK1, DAT, GTL2, PEG11, antiPEG11, and MEG8

Two ovine BAC clones and a connecting long-range PCR product, jointly spanning approximately 250 kb and representing most of the MULGE5-OY3 marker interval known to contain the clpg locus, were completely sequenced. The resulting genomic sequence was aligned with its human ortholog and extensively annotated. Six transcripts, four of which were novel, were predicted to originate from within the analyzed region and their existence confirmed experimentally: DLK1, DAT, GTL2, PEG11, antiPEG11, and MEG8. RT-PCR experiments performed on a range of tissues sampled from an 8-wk-old animal demonstrated the preferential expression of all six transcripts in skeletal muscle, which suggests that they are under control of common regulatory elements. The six transcripts were also shown to be subject to parental imprinting: DLK1, DAT, and PEG11 were shown to be paternally expressed and GTL2, antiPEG11, and MEG8 to be maternally expressed.

The callipyge mutation enhances the expression of coregulated imprinted genes in cis without affecting their imprinting status

The callipyge (CLPG) phenotype (from kappa(alpha)lambda(iota), "beautiful," and pi(iota)gamma(epsilon), "buttocks") described in sheep is an inherited muscular hypertrophy that is subject to an unusual parent-of-origin effect referred to as polar overdominance: only heterozygous individuals having inherited the CLPG mutation from their sire exhibit the muscular hypertrophy. The callipyge (clpg) locus was mapped to a chromosome segment of approximately 400 kb (refs. 2-4), which was shown to contain four genes (DLK1, GTL2, PEG11 and MEG8) that are preferentially expressed in skeletal muscle and subject to parental imprinting in this tissue. Here we describe the effect of the CLPG mutation on the expression of these four genes, and demonstrate that callipyge individuals have a unique expression profile that may account for the observed polar overdominance.

Convenient genotyping of six myostatin mutations causing double-muscling in cattle using a multiplex oligonucleotide ligation assay

We herein describe a procedure that allows for simultaneous genotyping of six loss-of-function mutations in the bovine myostatin gene associated with the double-muscling phenotype. The proposed method relies on a multiplex oligonucleotide ligation assay and detection of the fluorescently labelled products using automatic sequencers.

Extensive genome-wide linkage disequilibrium in cattle

A genome-wide linkage disequilibrium (LD) map was generated using microsatellite genotypes (284 autosomal microsatellite loci) of 581 gametes sampled from the dutch black-and-white dairy cattle population. LD was measured between all marker pairs, both syntenic and nonsyntenic. Analysis of syntenic pairs revealed surprisingly high levels of LD that, although more pronounced for closely linked marker pairs, extended over several tens of centimorgan. In addition, significant gametic associations were also shown to be very common between nonsyntenic loci. Simulations using the known genealogies of the studied sample indicate that random drift alone is likely to account for most of the observed disequilibrium. No clear evidence was obtained for a direct effect of selection ("Bulmer effect"). The observation of long range disequilibrium between syntenic loci using low-density marker maps indicates that LD mapping has the potential to be very effective in livestock populations. The frequent occurrence of gametic associations between nonsyntenic loci, however, encourages the combined use of linkage and linkage disequilibrium methods to avoid false positive results when mapping genes in livestock.

Single photon emission tomography with 99m Tc-HMPAO in Arab patients with depression

BACKGROUND

This study investigates the rate of cerebral blood flow (rCBF) in Arab patients wth depression.

METHODS

Forty-four patients with DSM-III-R major depressive disorders were studied at rest using single photon emission computerized tomography (SPECT) with 99m Tc-HMPAO in comparison with 20 normal controls. All patients were assessed using the Hamilton Rating Scale for Depression (HRSD).

RESULTS

The depressed group showed greater rCBF in left and right posterior frontal and parietal cortical regions than normal controls. Within the depressed group, patients with the least severe illness (HRSD < 20) had significantly lower rCBF than normal controls, whilst those with moderately severe (HRSD 20-29) and severe (HRSD > 30) had significantly greater rCBF in most cortical regions than normal controls. Symptom scores, derived from the HRSD were predicted by rCBF principally increased rCBF in the left frontal cortex.

CONCLUSIONS

These results suggest a generalized cerebral activation principally in the frontal cortex which is in contrast to the results of most previous studies but more in line with the results of studies of induced affect and some studies of depression subsyndromes.

Fine-mapping of quantitative trait loci by identity by descent in outbred populations: application to milk production in dairy cattle

We previously mapped a quantitative trait locus (QTL) affecting milk production to bovine chromosome 14. To refine the map position of this QTL, we have increased the density of the genetic map of BTA14q11-16 by addition of nine microsatellites and three single nucleotide polymorphisms. Fine-mapping of the QTL was accomplished by a two-tiered approach. In the first phase, we identified seven sires heterozygous "Qq" for the QTL by marker-assisted segregation analysis in a Holstein-Friesian pedigree comprising 1,158 individuals. In a second phase, we genotyped the seven selected sires for the newly developed high-density marker map and searched for a shared haplotype flanking an hypothetical, identical-by-descent QTL allele with large substitution effect. The seven chromosomes increasing milk fat percentage were indeed shown to carry a common chromosome segment with an estimated size of 5 cM predicted to contain the studied QTL. The same haplotype was shown to be associated with increased fat percentage in the general population as well, providing additional support in favor of the location of the QTL within the corresponding interval.

The role of pterins and related factors in the biology of early postpartum depression

Plasma tryptophan and other putative amino acids, cortisol, folate and vitamin B12 and urinary biopterin (B) and neopterins (N) were measured in three groups of women: 62 women in the early postpartum period, 23 pregnant and 38 non-gravid controls. Sixty-two postpartum women were screened for depression by the Edinburgh postnatal depression scale (EPDS) on day 7 after delivery. Postpartum women had significantly lower tryptophan, vitamin B12 and significantly greater levels of cortisol, folate, neopterins and biopterins than controls. Comparisons between women who were classified on the EPDS as cases and non-cases revealed only a statistically significant difference for lower N:B (P<0.01) and lower folate (P<0.01) ratio in cases than non-cases. Multiple regression analysis showed a significant contribution for low tryptophan to increased EPDS which also showed significant correlations with low methionine, low tyrosine, low N:B ratio and high vitamin B12.

Hormonal aspects of postpartum depression

Plasma cortisol, prolactin, oestrogen, progesterone, thyroxine, thyrotrophin (TSH) were collected from 23 pregnant, 70 postpartum women at 7 days postpartum, and 38 non-gravid controls. Sixty two postpartum women were screened for depression by the Edinburgh Postnatal Depression Scale (EPDS) on day 7 after delivery and 34 of them were assessed by the Present State Examination (PSE) at 8 +/- 2 weeks after delivery. Postpartum women had a significantly greater level of cortisol, prolactin, thyroxine and oestrogen than non-puerperal women. Postpartum women with current depression (EPDS > or = 11) had significantly lower plasma prolactin levels than those without depression and those who developed depression within 6-10 weeks after delivery (PSE level > or = 5) had significantly lower plasma prolactin and significantly greater progesterone levels than those who were not depressed. There were significant correlations between age and plasma cortisol and prolactin levels. Higher thyroxine levels predicted greater severity of concurrent symptoms of depression (total EPDS score) whilst higher progesterone and lower prolactin levels predicted the occurrence of depression (total PSE score) 6-10 weeks after delivery. Women who breastfed had significantly lower EPDS and total PSE scores and higher plasma prolactin levels than those who did not breastfed their infants whilst women who had previous episodes of depression had significantly greater EPDS and PSE scores, lower prolactin and higher TSH levels than those who had not suffered from previous episodes of depression.

A QTL with major effect on milk yield and composition maps to bovine chromosome 14

A whole genome scan was undertaken in a granddaughter design comprising 1158 progeny-tested bulls in order to map QTL influencing milk yield and composition. In this paper we report the identification of a locus on the centromeric end of bovine Chromosome (Chr) 14, with major effect on fat and protein percentage as well as milk yield. The genuine nature of this QTL was verified using the grand2-daughter design, that is, by tracing the segregating QTL alleles from heterozygous grandsires to their maternal grandsons and confirming the predicted QTL allele substitution effect.

A QTL affecting milk yield and composition maps to bovine chromosome 20: a confirmation

As part of a whole genome scan undertaken to detect quantitative trait loci (QTL) affecting milk yield and composition, we have genotyped a granddaughter design comprising 1152 sons for six microsatellite markers spanning bovine chromosome 20. An analysis performed across families provided strong evidence (experiment-wise P-values < 0.01) for the presence of a QTL affecting primarily protein percentage towards the telomeric end of the chromosome. A founder sire, shown in a previous study to segregate for a similar QTL in the corresponding chromosome region, was characterized by 29 and 57 sons and maternal grandsons, respectively, in the present design. Sorting corresponding sons and grandsons by paternal or grandpaternal allele provided significant evidence for the segregation of a QTL on chromosome 20. Altogether these results confirm the location of a QTL affecting milk production on bovine chromosome 20.

Anorexia nervosa in an Arab culture

The available literature indicates that anorexia nervosa is rare in Arab culture. We report 5 cases of anorexia nervosa: 3 female and 2 male patients. Two are UAE nationals, 2 are Omanis, and 1 is Sudanese. The occurrence of these cases is discussed in the context of sociocultural changes and the increasing westernization of the UAE society.

Quantitative trait loci analysis for five milk production traits on chromosome six in the Dutch Holstein-Friesian population

Twenty Dutch Holstein-Friesian families, with a total of 715 sires, were evaluated in a granddaughter experiment design for marker-QTL associations. Five traits-milk, fat and protein yield and fat and protein percent-were analyzed. Across-family analysis was undertaken using multimarker regression principles. One and two QTL models were fitted. Critical values for the test statistic were calculated empirically by permuting the data. Individual trait distributions of permuted test statistics differed and, thus distributions, had to be calculated for each trait. Experimentwise critical values, which account for evaluating marker-QTL associations on all 29 autosomal-bovine chromosomes and for five traits, were calculated. A QTL for protein percent was identified in one, and two QTL models and was significant at the 1 and 2% level, respectively. Extending the multimarker regression approach to an analysis including two QTL was limited by families not being informative at all markers, which resulted in singularity. Below average heterozygosity for the first and last marker lowered information content for the first and last marker bracket. Highly informative markers at the ends of the mapped chromosome would overcome the decrease in information content in the first and last marker bracket and singularity for the two QTL model.

Efficacy and safety of risperidone in psychotic patients: an open study

This open prospective study was undertaken to determine the efficacy and safety of a fixed dose (6 mg) of risperidone in psychotic patients. Hospital in-patients who fulfilled DSM-111-R criteria for schizophrenia, schizoaffective and bipolar disorders were eligible for entry into the study (n = 15). Patients who were on other antipsychotics had a washout period of 1 week before they were started on the drug. A fixed dose of risperidone was administered (6 mg). The Brief Psychiatric Rating Scale (BPRS), Negative Symptom Rating Scale (NSRS) and Abnormal Involuntary Movement Scale were used to measure psychopathology and extrapyramidal side-effects. Five patients dropped out of the study. Two patients became very agitated and potentially aggressive, one patient became very restless and did not respond to benzodiazepines, and one dropped out because of restlessness that did not respond to clonazepam. Of the 10 patients who completed the study, 50 per cent reduction on BPRS and NSRS was achieved by five and six patients respectively. There was a marginally significant trend towards a greater reduction in the magnitude of negative symptoms. Four patients required treatment with anticholinergic drugs. Risperidone was effective in resistent psychotic patients, but agitated and impulsive psychotic patients with positive symptoms may not be best candidates for treatment with risperidone. On average, negative symptoms respond better than positive symptoms.

Predictors of employment status of treated patients with DSM-III-R diagnosis. Can logistic regression model find a solution?

To investigate the predictors of employment status of patients with DSM-III-R diagnosis, 55 patients were selected by a simple random technique from the main psychiatric clinic in Al Ain, United Arab Emirates. Structured and formal assessments were carried out to extract the potential predictors of outcome of schizophrenia. Logistic regression model revealed that being married, absence of schizoid personality, free or with minimum symptoms of the illness, later age of onset, and higher educational attainment were the most significant predictors of employment outcome. The implications of the results of this study are discussed in the text.