The transcriptomic response of bovine uterine tissue is altered in response to sperm from high- and low-fertility bulls†.

Despite stringent quality control checks, some bulls with apparently normal semen quality yield lower than expected pregnancy rates. This study profiled the transcriptome and performed histological analysis of the bovine uterus in response to sperm from high-fertility (HF) and low-fertility (LF) bulls. Postmortem uterine biopsies and uterine explants were collected from heifers 12 h after a fixed time artificial insemination (AI) to a synchronized estrus with frozen-thawed semen from five HF (fertility rate 4.01% ± 0.25) and five LF (fertility rate - 11.29% ± 1.11; mean ± SEM) bulls. Uterine biopsies were also collected from control (CTRL) heifers, which were not inseminated. RNA-sequencing and histological analysis were performed for differential gene expression and neutrophil quantification. In the HF treatment relative to CTRL heifers, there were 376 genes significantly differentially expressed in the endometrium with just one gene differentially expressed in the LF treatment relative to CTRL heifers. Comparing the HF and LF treatments directly, there were 40 significantly differentially expressed genes (P < 0.05). Transcriptomic analysis shows a predominant role for the inflammatory marker Interleukin-1 alpha, which was further confirmed by immunohistochemistry. Quantification of neutrophils in the endometrium showed a significant effect of sperm; however, there was no difference in neutrophil numbers between HF and LF groups. In conclusion, this novel study clearly shows a distinct inflammatory response to sperm in the endometrium and a divergent transcriptomic response to semen from HF and LF bulls.

Immunoproteomic analysis of the serum IgG response to cell wall-associated proteins of Staphylococcus aureus strains belonging to CC97 and CC151.

CC97 and CC151 are two of the most common Staphylococcus aureus lineages associated with bovine intramammary infection. The genotype of the infecting S. aureus strain influences virulence and the progression of intramammary disease. Strains from CC97 and CC151 encode a distinct array of virulence factors. Identification of proteins elaborated in vivo will provide insights into the molecular mechanism of pathogenesis of these lineages, as well as facilitating the development of tailored treatments and pan-lineage vaccines and diagnostics. The repertoire of genes encoding cell wall-anchored (CWA) proteins was identified for S. aureus strains MOK023 (CC97) and MOK124 (CC151); MOK023 encoded more CWA proteins than MOK124. Serum collected during an in vivo challenge trial was used to investigate whether the humoral response to cell wall proteins was strain-specific. Immunoproteomic analysis demonstrated that the humoral response in MOK023-infected cows predominantly targeted high molecular weight proteins while the response in MOK124-infected cows targeted medium or low molecular weight proteins. Antigenic proteins were identified by two-dimensional serum blotting followed by mass spectometry-based identification of immunoreactive spots, with putative antigens subsequently validated. The CWA proteins ClfB, SdrE/Bbp and IsdA were identified as immunogenic regardless of the infecting strain. In addition, a number of putative strain-specific imunogens were identified. The variation in antigens produced by different strains may indicate that these strains have different strategies for exploiting the intramammary niche. Such variation should be considered when developing novel control strategies including vaccines, therapeutics and diagnostics.

Ewe breed differences in the cervical transcriptome at the follicular phase of a synchronised oestrous cycle.

BACKGROUND: Cervical artificial insemination (AI) with frozen-thawed semen results in unacceptably low pregnancy rates internationally. The exception is in Norway, where vaginal deposition of frozen-thawed semen to a natural oestrous routinely yields pregnancy rates in excess of 70%. Previous studies by our group has demonstrated that this is due to differences in cervical sperm transport. However, a potentially important contributory factor is that ewes are inseminated to a natural oestrous in Norway but to a synchronised oestrous across most of the rest of the world. In this study, we interrogated the gene expression of the sheep cervix of four ewe breeds with known differences in pregnancy rates following cervical AI using frozen-thawed semen under the effect of exogenous hormones to synchronise the oestrous cycle. These four ewe breeds (n = 8 to 11 ewes per breed) are from two countries: Ireland (Belclare and Suffolk; medium and low fertility, respectively) and Norway (Norwegian White Sheep (NWS) and Fur; both with high fertility compared to the Irish ewe breeds). RESULTS: RNA extracted from cervical biopsies collected from these breeds was analysed by RNA-sequencing and differential gene expression analysis. Using the low-fertility Suffolk breed as a reference level; 27, 1827 and 2641 genes were differentially expressed in Belclare, Fur and NWS ewes, respectively (P <  0.05 and FC > 1.5). Gene ontology (GO) analysis revealed that Fur and NWS had an up-regulation of enriched pathways involved in muscle contraction and development compared to Suffolk. However, there was a down-regulation of the immune response pathway in NWS compared to Suffolk. In addition, GO analysis showed similar expression patterns involved in muscle contraction, extracellular matrix (ECM) development and cell-cell junction in both Norwegian ewe breeds, which differed to the Irish ewe breeds. CONCLUSIONS: This novel study has identified a number of conserved and breed-specific biological processes under the effect of oestrous synchronisation that may impact cervical sperm transport during the follicular phase of the reproductive cycle.

Cervical immune activation during the luteal phase may compromise subsequent trans-cervical ram sperm transport†.

Worldwide, cervical artificial insemination using frozen-thawed semen yields low pregnancy rates. The only exception to this is in Norway, where vaginal insemination with frozen-thawed semen yields pregnancy rates in excess of 60% and which has been attributed to the specific ewe breed used. Our previous work demonstrated differences in cervical gene expression at the follicular phase of the estrous cycle in ewe breeds with known differences in pregnancy rates. In this study, we characterized the cervical transcriptome of the same ewe breeds [Suffolk, Belclare, Fur, and Norwegian White Sheep (NWS)] during the luteal phase, as an optimal environment at the luteal phase could better prepare the cervix for sperm migration through the cervix at the subsequent follicular phase. High-quality RNA extracted from postmortem cervical tissue was analyzed by RNA sequencing. After stringent filtering, 1051, 1924, and 611 differentially expressed genes (DEGs) were detected in the low-fertility Suffolk breed compared with Belclare, Fur, and NWS, respectively. Gene ontology analysis identified increased humoral adaptive immune response pathways in Suffolk. Increased expression of multiple immune genes supports the presence of an active immune response in the cervix of Suffolk ewes, which differentiates them significantly from the other three ewe breeds. Inflammatory pathways were upregulated in the Suffolk, resulting in higher expression of the potent pro-inflammatory cytokines. Therefore, higher levels of pro-inflammatory cytokines indicate unresolved inflammation in the cervix of the low-fertility Suffolk breed that could contribute to reduced cervical sperm transport in the next follicular phase.

Biochemical and molecular characterization of sialylated cervical mucins in sheep†.

Sialic acid occupies terminal positions on O-glycans of cervical mucins, where they contribute to the increased viscosity of mucin thereby regulating sperm transport. This study characterized the sialylated cervical mucins from follicular phase mucus of six European ewe breeds with known differences in pregnancy rates following cervical artificial insemination (AI) using frozen-thawed semen at both synchronized and natural estrus cycles. These were Suffolk (low fertility) and Belclare (medium fertility) in Ireland, Ile de France and Romanov (both with medium fertility) in France, and Norwegian White Sheep (NWS) and Fur (both with high fertility) in Norway. Expression of mucin and sialic acid related genes was quantified using RNA-sequencing in cervical tissue from Suffolk, Belclare, Fur, and NWS only. Cervical tissue was also assessed for the percentage of cervical epithelial populated by mucin secreting goblet cells in the same four ewe breeds. Biochemical analysis showed that there was an effect of ewe breed on sialic acid species, which was represented by Suffolk having higher levels of Neu5,9Ac2 compared with NWS (P < 0.05). Suffolk ewes had a lower percentage of goblet cells than Fur and NWS (P < 0.05). Gene expression analysis identified higher expression of MUC5AC, MUC5B, ST6GAL1, and ST6GAL2 and lower expression of ST3GAL3, ST3GAL4, and SIGLEC10 in Suffolk compared with high fertility ewe breeds (P < 0.05). Our results indicate that specific alterations in sialylated mucin composition may be related to impaired cervical sperm transport.

Bovine milk somatic cell transcriptomic response to Staphylococcus aureus is dependent on strain genotype.

BACKGROUND: Mastitis is an economically important disease of dairy cows with Staphylococcus aureus a major cause worldwide. Challenge of Holstein-Friesian cows demonstrated that S. aureus strain MOK124, which belongs to Clonal Complex (CC)151, caused clinical mastitis, while strain MOK023, belonging to CC97, caused mild or subclinical mastitis. The aim of this study was to elucidate the molecular mechanisms of the host immune response utilising a transcriptomic approach. Milk somatic cells were collected from cows infected with either S. aureus MOK023 or MOK124 at 0, 24, 48, 72 and 168 h post-infection (hpi) and analysed for differentially expressed (DE) genes in response to each strain. RESULTS: In response to MOK023, 1278, 2278, 1986 and 1750 DE genes were found at 24, 48, 72 and 168 hpi, respectively, while 2293, 1979, 1428 and 1544 DE genes were found in response to MOK124 at those time points. Genes involved in milk production (CSN1, CSN10, CSN1S2, CSN2, a-LACTA and PRLR) were downregulated in response to both strains, with a more pronounced decrease in the MOK124 group. Immune response pathways such as NF-κB and TNF signalling were overrepresented in response to both strains at 24 hpi. These immune pathways continued to be overrepresented in the MOK023 group at 48 and 72 hpi, while the Hippo signalling, extracellular matrix interaction (ECM) and tight junction pathways were overrepresented in the MOK124 group between 48 and 168 hpi. Cellular composition analysis demonstrated that a neutrophil response was predominant in response to MOK124, while M1 macrophages were the main milk cell type post-infection in the MOK023 group. CONCLUSIONS: A switch from immune response pathways to pathways involved in maintaining the integrity of the epithelial cell layer was observed in the MOK124 group from 48 hpi, which coincided with the occurrence of clinical signs in the infected animals. The higher proportion of M1 macrophages in the MOK023 group and lack of substantial neutrophil recruitment in response to MOK023 may indicate immune evasion by this strain. The results of this study highlight that the somatic cell transcriptomic response to S. aureus is dependent on the genotype of the infecting strain.

Conserved and breed-specific differences in the cervical transcriptome of sheep with divergent fertility at the follicular phase of a natural oestrus cycle.

BACKGROUND: The outcome of cervical artificial insemination (AI) with frozen-thawed semen in sheep is limited by the inability of sperm to traverse the cervix of some ewe breeds. Previous research has demonstrated that cervical sperm transport is dependent on ewe breed, as sperm can traverse the cervix in greater numbers in some higher fertility ewe breeds. However, the molecular mechanisms underlying ewe breed differences in sperm transport through the cervix remain unknown. In this study, we aimed to characterise the cervical transcriptome of four European ewe breeds with known differences in pregnancy rates following cervical AI using frozen-thawed semen at the follicular phase of a natural oestrous cycle. Cervical post mortem tissue samples were collected from two Irish ewe breeds (Belclare and Suffolk; medium and low fertility, respectively) and from two Norwegian ewe breeds (Norwegian White Sheep (NWS) and Fur; high fertility compared to both Irish breeds) at the follicular phase of a natural oestrous cycle (n = 8 to 10 ewes per breed). RESULTS: High-quality RNA extracted from biopsies of the mid-region of the cervix was analysed by RNA-sequencing and Gene Ontology (GO). After stringent filtering (P <  0.05 and FC > 1.5), a total of 11, 1539 and 748 differentially expressed genes (DEGs) were identified in Belclare, Fur and NWS compared to the low fertility Suffolk breed, respectively. Gene ontology analysis identified significantly enriched biological processes involved in muscle contraction, extracellular matrix (ECM) development and the immune response. Gene co-expression analysis revealed similar patterns in muscle contraction and ECM development modules in both Norwegian ewe breeds, which differed to the Irish ewe breeds. CONCLUSIONS: These breed-specific biological processes may account for impaired cervical sperm transport through the cervix in sheep during the follicular phase of the reproductive cycle. This novel and comprehensive dataset provides a rich foundation for future targeted initiatives to improve cervical AI in sheep.

Identifying schizophrenia patients who carry pathogenic genetic copy number variants using standard clinical assessment: retrospective cohort study.

BACKGROUND: Copy number variants (CNVs) play a significant role in disease pathogenesis in a small subset of individuals with schizophrenia (~2.5%). Chromosomal microarray testing is a first-tier genetic test for many neurodevelopmental disorders. Similar testing could be useful in schizophrenia. AIMS: To determine whether clinically identifiable phenotypic features could be used to successfully model schizophrenia-associated (SCZ-associated) CNV carrier status in a large schizophrenia cohort. METHOD: Logistic regression and receiver operating characteristic (ROC) curves tested the accuracy of readily identifiable phenotypic features in modelling SCZ-associated CNV status in a discovery data-set of 1215 individuals with psychosis. A replication analysis was undertaken in a second psychosis data-set (n = 479). RESULTS: In the discovery cohort, specific learning disorder (OR = 8.12; 95% CI 1.16-34.88, P = 0.012), developmental delay (OR = 5.19; 95% CI 1.58-14.76, P = 0.003) and comorbid neurodevelopmental disorder (OR = 5.87; 95% CI 1.28-19.69, P = 0.009) were significant independent variables in modelling positive carrier status for a SCZ-associated CNV, with an area under the ROC (AUROC) of 74.2% (95% CI 61.9-86.4%). A model constructed from the discovery cohort including developmental delay and comorbid neurodevelopmental disorder variables resulted in an AUROC of 83% (95% CI 52.0-100.0%) for the replication cohort. CONCLUSIONS: These findings suggest that careful clinical history taking to document specific neurodevelopmental features may be informative in screening for individuals with schizophrenia who are at higher risk of carrying known SCZ-associated CNVs. Identification of genomic disorders in these individuals is likely to have clinical benefits similar to those demonstrated for other neurodevelopmental disorders.

Characterization of the bovine salivary gland transcriptome associated with Mycobacterium avium subsp. paratuberculosis experimental challenge.

BACKGROUND: Mycobacterium avium subsp. paratuberculosis (MAP), the etiologic agent of Johne's disease is spread between cattle via the fecal-oral route, yet the functional changes in the salivary gland associated with infection remain uncharacterized. In this study, we hypothesized that experimental challenge with MAP would induce stable changes in gene expression patterns in the salivary gland that may shed light on the mucosal immune response as well as the regional variation in immune capacity of this extensive gland. Holstein-Friesian cattle were euthanized 33 months' post oral challenge with MAP strain CIT003 and both the parotid and mandibular salivary glands were collected from healthy control (n = 5) and MAP exposed cattle (n = 5) for histopathological and transcriptomic analysis. RESULTS: A total of 205, 21, 61, and 135 genes were significantly differentially expressed between control and MAP exposed cattle in dorsal mandibular (M1), ventral mandibular (M2), dorsal parotid (P1) and ventral parotid salivary glands (P2), respectively. Expression profiles varied between the structurally divergent parotid and mandibular gland sections which was also reflected in the enriched biological pathways identified. Changes in gene expression associated with MAP exposure were detected with significantly elevated expression of BoLA DR-ALPHA, BOLA-DRB3 and complement factors in MAP exposed cattle. In contrast, reduced expression of genes such as polymeric immunoglobin receptor (PIGR), TNFSF13, and the antimicrobial genes lactoferrin (LF) and lactoperoxidase (LPO) was detected in MAP exposed animals. CONCLUSIONS: This first analysis of the transcriptomic profile of salivary glands in cattle adds an important layer to our understanding of salivary gland immune function. Transcriptomic changes associated with MAP exposure have been identified including reduced LF and LPO. These critical antimicrobial and immunoregulatory proteins are known to be secreted into saliva and their downregulation may contribute to disease susceptibility. Future work will focus on the validation of their expression levels in saliva from additional cattle of known infection status as a potential strategy to augment disease diagnosis.

Transcriptome characterization and differentially expressed genes under flooding and drought stress in the biomass grasses Phalaris arundinacea and Dactylis glomerata.

BACKGROUND AND AIMS: Perennial grasses are a global resource as forage, and for alternative uses in bioenergy and as raw materials for the processing industry. Marginal lands can be valuable for perennial biomass grass production, if perennial biomass grasses can cope with adverse abiotic environmental stresses such as drought and waterlogging. METHODS: In this study, two perennial grass species, reed canary grass (Phalaris arundinacea) and cocksfoot (Dactylis glomerata) were subjected to drought and waterlogging stress to study their responses for insights to improving environmental stress tolerance. Physiological responses were recorded, reference transcriptomes established and differential gene expression investigated between control and stress conditions. We applied a robust non-parametric method, RoDEO, based on rank ordering of transcripts to investigate differential gene expression. Furthermore, we extended and validated vRoDEO for comparing samples with varying sequencing depths. KEY RESULTS: This allowed us to identify expressed genes under drought and waterlogging whilst using only a limited number of RNA sequencing experiments. Validating the methodology, several differentially expressed candidate genes involved in the stage 3 step-wise scheme in detoxification and degradation of xenobiotics were recovered, while several novel stress-related genes classified as of unknown function were discovered. CONCLUSIONS: Reed canary grass is a species coping particularly well with flooding conditions, but this study adds novel information on how its transcriptome reacts under drought stress. We built extensive transcriptomes for the two investigated C3 species cocksfoot and reed canary grass under both extremes of water stress to provide a clear comparison amongst the two species to broaden our horizon for comparative studies, but further confirmation of the data would be ideal to obtain a more detailed picture.

Effect of early calf-hood nutrition on the transcriptomic profile of subcutaneous adipose tissue in Holstein-Friesian bulls.

BACKGROUND: Adipose tissue is a major endocrine organ and is thought to play a central role in the metabolic control of reproductive function in cattle. Plane of nutrition during early life has been shown to influence the timing of puberty in both male and female cattle, though the exact biological mechanisms involved are currently unknown. The aim of this study was to investigate the effect of early calf-hood nutrition on the transcriptomic profile of subcutaneous adipose tissue in Holstein-Friesian bulls to identify possible downstream effects on reproductive physiology. RESULTS: Holstein-Friesian bull calves with a mean (±S.D.) age and bodyweight of 19 (±8.2) days and 47.5 (±5.3) kg, respectively, were assigned to either a high (n = 10) or low (n = 10) plane of nutrition. Calves were fed in order to achieve an overall growth rate of 1.08 and 0.57 kg/day for the high and low plane of nutrition treatments, respectively. At 126 days of age, the bulls were euthanized, subcutaneous adipose tissue samples were harvested and RNAseq analysis was performed. There were 674 genes differentially expressed in adipose tissue of calves on the low compared with the high plane of nutrition (P < 0.05; FDR < 0.05; fold change > 2.0). High plane of nutrition positively altered the expression of genes across an array of putative biological processes but the most dominant cellular processes affected were cellular energy production and branched chain amino acid degradation. A high plane of nutrition caused upregulation of genes such as leptin (LEP) and adiponectin (ADIPOQ), which are known to directly affect reproductive function. CONCLUSIONS: These results provide an insight into the effect of augmenting the plane of nutrition of Holstein-Friesian bull calves in the prepubertal period on the transcriptome of adipose tissue.

Genome-wide comparative analysis of atopic dermatitis and psoriasis gives insight into opposing genetic mechanisms.

Atopic dermatitis and psoriasis are the two most common immune-mediated inflammatory disorders affecting the skin. Genome-wide studies demonstrate a high degree of genetic overlap, but these diseases have mutually exclusive clinical phenotypes and opposing immune mechanisms. Despite their prevalence, atopic dermatitis and psoriasis very rarely co-occur within one individual. By utilizing genome-wide association study and ImmunoChip data from >19,000 individuals and methodologies developed from meta-analysis, we have identified opposing risk alleles at shared loci as well as independent disease-specific loci within the epidermal differentiation complex (chromosome 1q21.3), the Th2 locus control region (chromosome 5q31.1), and the major histocompatibility complex (chromosome 6p21-22). We further identified previously unreported pleiotropic alleles with opposing effects on atopic dermatitis and psoriasis risk in PRKRA and ANXA6/TNIP1. In contrast, there was no evidence for shared loci with effects operating in the same direction on both diseases. Our results show that atopic dermatitis and psoriasis have distinct genetic mechanisms with opposing effects in shared pathways influencing epidermal differentiation and immune response. The statistical analysis methods developed in the conduct of this study have produced additional insight from previously published data sets. The approach is likely to be applicable to the investigation of the genetic basis of other complex traits with overlapping and distinct clinical features.

Illumina MiSeq 16S amplicon sequence analysis of bovine respiratory disease associated bacteria in lung and mediastinal lymph node tissue.

BACKGROUND: Bovine respiratory disease (BRD) is caused by growth of single or multiple species of pathogenic bacteria in lung tissue following stress and/or viral infection. Next generation sequencing of 16S ribosomal RNA gene PCR amplicons (NGS 16S amplicon analysis) is a powerful culture-independent open reference method that has recently been used to increase understanding of BRD-associated bacteria in the upper respiratory tract of BRD cattle. However, it has not yet been used to examine the microbiome of the bovine lower respiratory tract. The objective of this study was to use NGS 16S amplicon analysis to identify bacteria in post-mortem lung and lymph node tissue samples harvested from fatal BRD cases and clinically healthy animals. Cranial lobe and corresponding mediastinal lymph node post-mortem tissue samples were collected from calves diagnosed as BRD cases by veterinary laboratory pathologists and from clinically healthy calves. NGS 16S amplicon libraries, targeting the V3-V4 region of the bacterial 16S rRNA gene were prepared and sequenced on an Illumina MiSeq. Quantitative insights into microbial ecology (QIIME) was used to determine operational taxonomic units (OTUs) which corresponded to the 16S rRNA gene sequences. RESULTS: Leptotrichiaceae, Mycoplasma, Pasteurellaceae, and Fusobacterium were the most abundant OTUs identified in the lungs and lymph nodes of the calves which died from BRD. Leptotrichiaceae, Fusobacterium, Mycoplasma, Trueperella and Bacteroides had greater relative abundances in post-mortem lung samples collected from fatal cases of BRD in dairy calves, compared with clinically healthy calves without lung lesions. Leptotrichiaceae, Mycoplasma and Pasteurellaceae showed higher relative abundances in post-mortem lymph node samples collected from fatal cases of BRD in dairy calves, compared with clinically healthy calves without lung lesions. Two Leptotrichiaceae sequence contigs were subsequently assembled from bacterial DNA-enriched shotgun sequences. CONCLUSIONS: The microbiomes of the cranial lung lobe and mediastinal lymph node from calves which died from BRD and from clinically healthy H-F calves have been characterised. Contigs corresponding to the abundant Leptotrichiaceae OTU were sequenced and found not to be identical to any known bacterial genus. This suggests that we have identified a novel bacterial species associated with BRD.

Effect of dietary restriction and subsequent re-alimentation on the transcriptional profile of hepatic tissue in cattle.

BACKGROUND: Compensatory growth (CG) is an accelerated growth phenomenon observed in animals upon re-alimentation following a period of dietary restriction. It is typically utilised in livestock systems to reduce feed costs during periods of reduced feed availability. The biochemical mechanisms controlling this phenomenon, however, are yet to be elucidated. This study aimed to uncover the molecular mechanisms regulating the hepatic expression of CG in cattle, utilising RNAseq. RNAseq was performed on hepatic tissue of bulls following 125 days of dietary restriction (RES) and again following 55 days of subsequent re-alimentation during which the animals exhibited significant CG. The data were compared with those of control animals offered the same diet on an ad libitum basis throughout (ADLIB). Elucidation of the molecular control of CG may yield critical information on genes and pathways which could be targeted as putative molecular biomarkers for the selection of animals with improved CG potential. RESULTS: Following a period of differential feeding, body-weight and liver weight were 161 and 4 kg higher, respectively, for ADLIB compared with RES animals. At this time RNAseq analysis of liver tissue revealed 1352 significantly differentially expressed genes (DEG) between the two treatments. DEGs indicated down-regulation of processes including nutrient transport, cell division and proliferation in RES. In addition, protein synthesis genes were up-regulated in RES following a period of restricted feeding. The subsequent 55 days of ad libitum feeding for both groups resulted in the body-weight difference reduced to 84 kg, with no difference in liver weight between treatment groups. At the end of 55 days of unrestricted feeding, 49 genes were differentially expressed between animals undergoing CG and their continuously fed counterparts. In particular, hepatic expression of cell proliferation and growth genes were greater in animals undergoing CG. CONCLUSIONS: Greater expression of cell cycle and cell proliferation genes during CG was associated with a 100 % recovery of liver weight during re-alimentation. Additionally, an apparent up-regulation in capacity for cellular protein synthesis during restricted feeding may contribute to and sustain CG during re-alimentation. DEGs identified are potential candidate genes for the identification of biomarkers for CG, which may be incorporated into future breeding programmes.

CNV analysis in a large schizophrenia sample implicates deletions at 16p12.1 and SLC1A1 and duplications at 1p36.33 and CGNL1.

Large and rare copy number variants (CNVs) at several loci have been shown to increase risk for schizophrenia. Aiming to discover novel susceptibility CNV loci, we analyzed 6882 cases and 11 255 controls genotyped on Illumina arrays, most of which have not been used for this purpose before. We identified genes enriched for rare exonic CNVs among cases, and then attempted to replicate the findings in additional 14 568 cases and 15 274 controls. In a combined analysis of all samples, 12 distinct loci were enriched among cases with nominal levels of significance (P < 0.05); however, none would survive correction for multiple testing. These loci include recurrent deletions at 16p12.1, a locus previously associated with neurodevelopmental disorders (P = 0.0084 in the discovery sample and P = 0.023 in the replication sample). Other plausible candidates include non-recurrent deletions at the glutamate transporter gene SLC1A1, a CNV locus recently suggested to be involved in schizophrenia through linkage analysis, and duplications at 1p36.33 and CGNL1. A burden analysis of large (>500 kb), rare CNVs showed a 1.2% excess in cases after excluding known schizophrenia-associated loci, suggesting that additional susceptibility loci exist. However, even larger samples are required for their discovery.

An inherited duplication at the gene p21 Protein-Activated Kinase 7 (PAK7) is a risk factor for psychosis.

Identifying rare, highly penetrant risk mutations may be an important step in dissecting the molecular etiology of schizophrenia. We conducted a gene-based analysis of large (>100 kb), rare copy-number variants (CNVs) in the Wellcome Trust Case Control Consortium 2 (WTCCC2) schizophrenia sample of 1564 cases and 1748 controls all from Ireland, and further extended the analysis to include an additional 5196 UK controls. We found association with duplications at chr20p12.2 (P = 0.007) and evidence of replication in large independent European schizophrenia (P = 0.052) and UK bipolar disorder case-control cohorts (P = 0.047). A combined analysis of Irish/UK subjects including additional psychosis cases (schizophrenia and bipolar disorder) identified 22 carriers in 11 707 cases and 10 carriers in 21 204 controls [meta-analysis Cochran-Mantel-Haenszel P-value = 2 × 10(-4); odds ratio (OR) = 11.3, 95% CI = 3.7, ∞]. Nineteen of the 22 cases and 8 of the 10 controls carried duplications starting at 9.68 Mb with similar breakpoints across samples. By haplotype analysis and sequencing, we identified a tandem ~149 kb duplication overlapping the gene p21 Protein-Activated Kinase 7 (PAK7, also called PAK5) which was in linkage disequilibrium with local haplotypes (P = 2.5 × 10(-21)), indicative of a single ancestral duplication event. We confirmed the breakpoints in 8/8 carriers tested and found co-segregation of the duplication with illness in two additional family members of one of the affected probands. We demonstrate that PAK7 is developmentally co-expressed with another known psychosis risk gene (DISC1) suggesting a potential molecular mechanism involving aberrant synapse development and plasticity.

Differentially Expressed Genes in Endometrium and Corpus Luteum of Holstein Cows Selected for High and Low Fertility Are Enriched for Sequence Variants Associated with Fertility.

Despite the importance of fertility in humans and livestock, there has been little success dissecting the genetic basis of fertility. Our hypothesis was that genes differentially expressed in the endometrium and corpus luteum on Day 13 of the estrous cycle between cows with either good or poor genetic merit for fertility would be enriched for genetic variants associated with fertility. We combined a unique genetic model of fertility (cattle that have been selected for high and low fertility and show substantial difference in fertility) with gene expression data from these cattle and genome-wide association study (GWAS) results in ∼20,000 cattle to identify quantitative trait loci (QTL) regions and sequence variants associated with genetic variation in fertility. Two hundred and forty-five QTL regions and 17 sequence variants associated primarily with prostaglandin F2alpha, steroidogenesis, mRNA processing, energy status, and immune-related processes were identified. Ninety-three of the QTL regions were validated by two independent GWAS, with signals for fertility detected primarily on chromosomes 18, 5, 7, 8, and 29. Plausible causative mutations were identified, including one missense variant significantly associated with fertility and predicted to affect the protein function of EIF4EBP3. The results of this study enhance our understanding of 1) the contribution of the endometrium and corpus luteum transcriptome to phenotypic fertility differences and 2) the genetic architecture of fertility in dairy cattle. Including these variants in predictions of genomic breeding values may improve the rate of genetic gain for this critical trait.

Comparative exomics of Phalaris cultivars under salt stress.

BACKGROUND: Reed canary grass (Phalaris arundinacea) is an economically important forage and bioenergy grass of the temperate regions of the world. Despite its economic importance, it is lacking in public genomic data. We explore comparative exomics of the grass cultivars in the context of response to salt exposure. The limited data set poses challenges to the computational pipeline. METHODS: As a prerequisite for the comparative study, we generate the Phalaris reference transcriptome sequence, one of the first steps in addressing the issue of paucity of processed genomic data in this species. In addition, the differential expression (DE) and active-but-stable genes for salt stress conditions were analyzed by a novel method that was experimentally verified on human RNA-seq data. For the comparative exomics, we focus on the DE and stable genic regions, with respect to salt stress, of the genome. RESULTS AND CONCLUSIONS: In our comparative study, we find that phylogeny of the DE and stable genic regions of the Phalaris cultivars are distinct. At the same time we find the phylogeny of the entire expressed reference transcriptome matches the phylogeny of only the stable genes. Thus the behavior of the different cultivars is distinguished by the salt stress response. This is also reflected in the genomic distinctions in the DE genic regions. These observations have important implications in the choice of cultivars, and their breeding, for bio-energy fuels. Further, we identified genes that are representative of DE under salt stress and could provide vital clues in our understanding of the stress handling mechanisms in general.

Identification of mechanosensitive genes during skeletal development: alteration of genes associated with cytoskeletal rearrangement and cell signalling pathways.

BACKGROUND: Mechanical stimulation is necessary for regulating correct formation of the skeleton. Here we test the hypothesis that mechanical stimulation of the embryonic skeletal system impacts expression levels of genes implicated in developmentally important signalling pathways in a genome wide approach. We use a mutant mouse model with altered mechanical stimulation due to the absence of limb skeletal muscle (Splotch-delayed) where muscle-less embryos show specific defects in skeletal elements including delayed ossification, changes in the size and shape of cartilage rudiments and joint fusion. We used Microarray and RNA sequencing analysis tools to identify differentially expressed genes between muscle-less and control embryonic (TS23) humerus tissue. RESULTS: We found that 680 independent genes were down-regulated and 452 genes up-regulated in humeri from muscle-less Spd embryos compared to littermate controls (at least 2-fold; corrected p-value ≤0.05). We analysed the resulting differentially expressed gene sets using Gene Ontology annotations to identify significant enrichment of genes associated with particular biological processes, showing that removal of mechanical stimuli from muscle contractions affected genes associated with development and differentiation, cytoskeletal architecture and cell signalling. Among cell signalling pathways, the most strongly disturbed was Wnt signalling, with 34 genes including 19 pathway target genes affected. Spatial gene expression analysis showed that both a Wnt ligand encoding gene (Wnt4) and a pathway antagonist (Sfrp2) are up-regulated specifically in the developing joint line, while the expression of a Wnt target gene, Cd44, is no longer detectable in muscle-less embryos. The identification of 84 genes associated with the cytoskeleton that are down-regulated in the absence of muscle indicates a number of candidate genes that are both mechanoresponsive and potentially involved in mechanotransduction, converting a mechanical stimulus into a transcriptional response. CONCLUSIONS: This work identifies key developmental regulatory genes impacted by altered mechanical stimulation, sheds light on the molecular mechanisms that interpret mechanical stimulation during skeletal development and provides valuable resources for further investigation of the mechanistic basis of mechanoregulation. In particular it highlights the Wnt signalling pathway as a potential point of integration of mechanical and molecular signalling and cytoskeletal components as mediators of the response.

Delineating the genetic heterogeneity of ALS using targeted high-throughput sequencing.

BACKGROUND: Over 100 genes have been implicated in the aetiology of amyotrophic lateral sclerosis (ALS). A detailed understanding of their independent and cumulative contributions to disease burden may help guide various clinical and research efforts. METHODS: Using targeted high-throughput sequencing, we characterised the variation of 10 Mendelian and 23 low penetrance/tentative ALS genes within a population-based cohort of 444 Irish ALS cases (50 fALS, 394 sALS) and 311 age-matched and geographically matched controls. RESULTS: Known or potential high-penetrance ALS variants were identified within 17.1% of patients (38% of fALS, 14.5% of sALS). 12.8% carried variants of Mendelian disease genes (C9orf72 8.78%; SETX 2.48%; ALS2 1.58%; FUS 0.45%; TARDBP 0.45%; OPTN 0.23%; VCP 0.23%. ANG, SOD1, VAPB 0%), 4.7% carried variants of low penetrance/tentative ALS genes and 9.7% (30% of fALS, 7.1% of sALS) carried previously described ALS variants (C9orf72 8.78%; FUS 0.45%; TARDBP 0.45%). 1.6% of patients carried multiple known/potential disease variants, including all identified carriers of an established ALS variant (p<0.01); TARDBP:c.859G>A(p.[G287S]) (n=2/2 sALS). Comparison of our results with those from studies of other European populations revealed significant differences in the spectrum of disease variation (p=1.7×10(-4)). CONCLUSIONS: Up to 17% of Irish ALS cases may carry high-penetrance variants within the investigated genes. However, the precise nature of genetic susceptibility differs significantly from that reported within other European populations. Certain variants may not cause disease in isolation and concomitant analysis of disease genes may prove highly important.