Variants In Protein-coding Genes Research Articles

Analysis of genetic variants in protein-coding genes of Aoluguya reindeer based on the whole-genome data.

Analysis of genetic variants in protein-coding genes of Aoluguya reindeer based on the whole-genome data.

Whole Genome Analysis of SNV and Indel Polymorphism in Common Marmosets (Callithrix jacchus).

The common marmoset (Callithrix jacchus) is one of the most widely used nonhuman primate models of human disease. Owing to limitations in sequencing technology, early genome assemblies of this species using short-read sequencing suffered from gaps. In addition, the genetic diversity of the species has not yet been adequately explored. Using long-read genome sequencing and expert annotation, we generated a high-quality genome resource creating a 2.898 Gb marmoset genome in which most of the euchromatin portion is assembled contiguously (contig N50 = 25.23 Mbp, scaffold N50 = 98.2 Mbp). We then performed whole genome sequencing on 84 marmosets sampling the genetic diversity from several marmoset research centers. We identified a total of 19.1 million single nucleotide variants (SNVs), of which 11.9 million can be reliably mapped to orthologous locations in the human genome. We also observed 2.8 million small insertion/deletion variants. This dataset includes an average of 5.4 million SNVs per marmoset individual and a total of 74,088 missense variants in protein-coding genes. Of the 4956 variants orthologous to human ClinVar SNVs (present in the same annotated gene and with the same functional consequence in marmoset and human), 27 have a clinical significance of pathogenic and/or likely pathogenic. This important marmoset genomic resource will help guide genetic analyses of natural variation, the discovery of spontaneous functional variation relevant to human disease models, and the development of genetically engineered marmoset disease models.

Mitochondrial DNA heteroplasmy analysis in keratoconus patients from China.

Background: Mitochondrial DNA (mtDNA) variants have been implicated in keratoconus (KC). The present study aimed to characterize the mtDNA heteroplasmy profile in KC and explore the association of mitochondrial heteroplasmic levels with KC. Methods: Mitochondrial sequencing of peripheral blood samples and corneal tomography were conducted in 300 KC cases and 300 matched controls. The number of heteroplasmic and homoplasmic variants was calculated across the mitochondrial genome. Spearman's correlation was used to analyze the correlation between the number of heteroplasmic variants and age. The association of mtDNA heteroplasmic level with KC was analyzed by logistic regression analysis. Moreover, the relationship between mitochondrial heteroplasmic levels and clinical parameters was determined by linear regression analysis. Results: The distribution of mtDNA heteroplasmic variants showed the highest number of heteroplasmic variants in the non-coding region, while the COX3 gene exhibited the highest number in protein-coding genes. Comparisons of the number of heteroplasmic and homoplasmic non-synonymous variants in protein-coding genes revealed no significant differences between KC cases and controls (all p > 0.05). In addition, the number of heteroplasmic variants was positively associated with age in all subjects (r = 0.085, p = 0.037). The logistic regression analyses indicated that the heteroplasmic levels of m.16180_16181delAA was associated with KC (p < 0.005). Linear regression analyses demonstrated that the heteroplasmic levels of m.16180_16181delAA and m.302A>C were not correlated with thinnest corneal thickness (TCT), steep keratometry (Ks), and flat keratometry (Kf) (all p > 0.05) in KC cases and controls separately. Conclusion: The current study characterized the mtDNA heteroplasmy profile in KC, and revealed that the heteroplasmic levels of m.16180_16181delAA were associated with KC.

Whole-Exome sequencing analysis identified TMSB10/TRABD2A locus to be associated with carfilzomib-related cardiotoxicity among patients with multiple myeloma.

Proteasome inhibitor Carfilzomib (CFZ) is effective in treating patients with refractory or relapsed multiple myeloma (MM) but has been associated with cardiovascular adverse events (CVAE) such as hypertension, cardiomyopathy, and heart failure. This study aimed to investigate the contribution of germline genetic variants in protein-coding genes in CFZ-CVAE among MM patients using whole-exome sequencing (WES) analysis. Exome-wide single-variant association analysis, gene-based analysis, and rare variant analyses were performed on 603,920 variants in 247 patients with MM who have been treated with CFZ and enrolled in the Oncology Research Information Exchange Network (ORIEN) at the Moffitt Cancer Center. Separate analyses were performed in European Americans and African Americans followed by a trans-ethnic meta-analysis. The most significant variant in the exome-wide single variant analysis was a missense variant rs7148 in the thymosin beta-10/TraB Domain Containing 2A (TMSB10/TRABD2A) locus. The effect allele of rs7148 was associated with a higher risk of CVAE [odds ratio (OR) = 9.3 with a 95% confidence interval of 3.9-22.3, p = 5.42*10-7]. MM patients with rs7148 AG or AA genotype had a higher risk of CVAE (50%) than those with GG genotype (10%). rs7148 is an expression quantitative trait locus (eQTL) for TRABD2A and TMSB10. The gene-based analysis also showed TRABD2A as the most significant gene associated with CFZ-CVAE (p = 1.06*10-6). We identified a missense SNP rs7148 in the TMSB10/TRABD2A as associated with CFZ-CVAE in MM patients. More investigation is needed to understand the underlying mechanisms of these associations.

Aberrant phase separation and nucleolar dysfunction in rare genetic diseases

Thousands of genetic variants in protein-coding genes have been linked to disease. However, the functional impact of most variants is unknown as they occur within intrinsically disordered protein regions that have poorly defined functions1–3. Intrinsically disordered regions can mediate phase separation and the formation of biomolecular condensates, such as the nucleolus4,5. This suggests that mutations in disordered proteins may alter condensate properties and function6–8. Here we show that a subset of disease-associated variants in disordered regions alter phase separation, cause mispartitioning into the nucleolus and disrupt nucleolar function. We discover de novo frameshift variants in HMGB1 that cause brachyphalangy, polydactyly and tibial aplasia syndrome, a rare complex malformation syndrome. The frameshifts replace the intrinsically disordered acidic tail of HMGB1 with an arginine-rich basic tail. The mutant tail alters HMGB1 phase separation, enhances its partitioning into the nucleolus and causes nucleolar dysfunction. We built a catalogue of more than 200,000 variants in disordered carboxy-terminal tails and identified more than 600 frameshifts that create arginine-rich basic tails in transcription factors and other proteins. For 12 out of the 13 disease-associated variants tested, the mutation enhanced partitioning into the nucleolus, and several variants altered rRNA biogenesis. These data identify the cause of a rare complex syndrome and suggest that a large number of genetic variants may dysregulate nucleoli and other biomolecular condensates in humans.

Towards a Multi-Omics of Male Infertility.

Infertility is a common problem affecting one in six couples and in 30% of infertile couples, the male factor is a major cause. A large number of genes are involved in spermatogenesis and a significant proportion of male infertility phenotypes are of genetic origin. Studies on infertility have so far primarily focused on chromosomal abnormalities and sequence variants in protein-coding genes and have identified a large number of disease-associated genes. However, it has been shown that a multitude of factors across various omics levels also contribute to infertility phenotypes. The complexity of male infertility has led to the understanding that an integrated, multi-omics analysis may be optimal for unravelling this disease. While there is a vast array of different factors across omics levels associated with infertility, the present review focuses on known factors from the genomics, epigenomics, transcriptomics, proteomics, metabolomics, glycomics, lipidomics, miRNomics, and integrated omics levels. These include: repeat expansions in AR, POLG, ATXN1, DMPK, and SHBG, multiple SNPs, copy number variants in the AZF region, disregulated miRNAs, altered H3K9 methylation, differential MTHFR, MEG3, PEG1, and LIT1 methylation, altered protamine ratios and protein hypo/hyperphosphorylation. This integrative review presents a step towards a multi-omics approach to understanding the complex etiology of male infertility. Currently only a few genetic factors, namely chromosomal abnormalities and Y chromosome microdeletions, are routinely tested in infertile men undergoing intracytoplasmic sperm injection. A multi-omics approach to understanding infertility phenotypes may yield a more holistic view of the disease and contribute to the development of improved screening methods and treatment options. Therefore, beside discovering as of yet unknown genetic causes of infertility, integrating multiple fields of study could yield valuable contributions to the understanding of disease development. Future multi-omics studies will enable to synthesise fragmented information and facilitate biomarker discovery and treatments in male infertility.

Combining genetic constraint with predictions of alternative splicing to prioritize deleterious splicing in rare disease studies

BackgroundDespite numerous molecular and computational advances, roughly half of patients with a rare disease remain undiagnosed after exome or genome sequencing. A particularly challenging barrier to diagnosis is identifying variants that cause deleterious alternative splicing at intronic or exonic loci outside of canonical donor or acceptor splice sites.ResultsSeveral existing tools predict the likelihood that a genetic variant causes alternative splicing. We sought to extend such methods by developing a new metric that aids in discerning whether a genetic variant leads to deleterious alternative splicing. Our metric combines genetic variation in the Genome Aggregate Database with alternative splicing predictions from SpliceAI to compare observed and expected levels of splice-altering genetic variation. We infer genic regions with significantly less splice-altering variation than expected to be constrained. The resulting model of regional splicing constraint captures differential splicing constraint across gene and exon categories, and the most constrained genic regions are enriched for pathogenic splice-altering variants. Building from this model, we developed ConSpliceML. This ensemble machine learning approach combines regional splicing constraint with multiple per-nucleotide alternative splicing scores to guide the prediction of deleterious splicing variants in protein-coding genes. ConSpliceML more accurately distinguishes deleterious and benign splicing variants than state-of-the-art splicing prediction methods, especially in “cryptic” splicing regions beyond canonical donor or acceptor splice sites.ConclusionIntegrating a model of genetic constraint with annotations from existing alternative splicing tools allows ConSpliceML to prioritize potentially deleterious splice-altering variants in studies of rare human diseases.

The relevance of mitochondrial DNA variants fluctuation during reprogramming and neuronal differentiation of human iPSCs.

SummaryThe generation of inducible pluripotent stem cells (iPSCs) is a revolutionary technique allowing production of pluripotent patient-specific cell lines used for disease modeling, drug screening, and cell therapy. Integrity of nuclear DNA (nDNA) is mandatory to allow iPSCs utilization, while quality control of mitochondrial DNA (mtDNA) is rarely included in the iPSCs validation process. In this study, we performed mtDNA deep sequencing during the transition from parental fibroblasts to reprogrammed iPSC and to differentiated neuronal precursor cells (NPCs) obtained from controls and patients affected by mitochondrial disorders. At each step, mtDNA variants, including those potentially pathogenic, fluctuate between emerging and disappearing, and some having functional implications. We strongly recommend including mtDNA analysis as an unavoidable assay to obtain fully certified usable iPSCs and NPCs.

Characterizing the molecular and immune landscape of canine bladder cancer.

Transitional cell carcinoma (TCC), also known as urothelial carcinoma, is the most common bladder cancer in humans and dogs. Approximately one-quarter of human TCCs are muscle-invasive and associated with a high risk of death from metastasis. Canine TCC (cTCC) tumours are typically high-grade and muscle-invasive. Shared similarities in risk factors, histopathology, and clinical presentation suggest that cTCC may serve as a model for the assessment of novel therapeutics that may inform therapies for human muscle-invasive TCC. The goal of this study was to characterize cTCC at the molecular level to identify drivers of oncogenesis and druggable targets. We performed whole exome sequencing (WES) of 11 cTCC tumours and three matched normal samples, identifying 583 variants in protein-coding genes. The most common variant was a V-to-E missense mutation in BRAF, identified in 4 out of 11 samples (36%) via WES. Sanger sequencing identified BRAF variants in 8 out of the same 11 cTCC samples, as well as in 22 out of 32 formalin-fixed paraffin embedded (FFPE) cTCC samples, suggesting an overall prevalence of 70%. RNA-Seq was performed to compare the gene expression profiles of cTCC tumours to normal bladder tissue. cTCC tumours exhibited up-regulation of genes involved in the cell cycle, DNA repair, and antiviral immunity. We also analysed the immune landscape of cTCC using immune gene signatures and immunohistochemical analysis. A subset of tumours had characteristics of a hot tumour microenvironment and exhibited high expression of signatures associated with complete response to PD-1/PD-L1 blockade in human bladder cancer.

Spectrum of splicing variants in disease genes and the ability of RNA analysis to reduce uncertainty in clinical interpretation

SummaryThe complexities of gene expression pose challenges for the clinical interpretation of splicing variants. To better understand splicing variants and their contribution to hereditary disease, we evaluated their prevalence, clinical classifications, and associations with diseases, inheritance, and functional characteristics in a 689,321-person clinical cohort and two large public datasets. In the clinical cohort, splicing variants represented 13% of all variants classified as pathogenic (P), likely pathogenic (LP), or variants of uncertain significance (VUSs). Most splicing variants were outside essential splice sites and were classified as VUSs. Among all individuals tested, 5.4% had a splicing VUS. If RNA analysis were to contribute supporting evidence to variant interpretation, we estimated that splicing VUSs would be reclassified in 1.7% of individuals in our cohort. This would result in a clinically significant result (i.e., P/LP) in 0.1% of individuals overall because most reclassifications would change VUSs to likely benign. In ClinVar, splicing VUSs were 4.8% of reported variants and could benefit from RNA analysis. In the Genome Aggregation Database (gnomAD), splicing variants comprised 9.4% of variants in protein-coding genes; most were rare, precluding unambiguous classification as benign. Splicing variants were depleted in genes associated with dominant inheritance and haploinsufficiency, although some genes had rare variants at essential splice sites or had common splicing variants that were most likely compatible with normal gene function. Overall, we describe the contribution of splicing variants to hereditary disease, the potential utility of RNA analysis for reclassifying splicing VUSs, and how natural variation may confound clinical interpretation of splicing variants.

Brain organoids with Neanderthal genes

Human Evolution The genomes of Neanderthals and modern humans are overall very similar. To understand the impact of genetic variants that are specific to modern humans, Trujillo et al. performed a genome-wide analysis to identify 61 coding variants in protein-coding genes. Identifying the gene encoding the RNA-binding protein NOVA1 as a top candidate for functional analyses, they introduced the archaic gene variant into human pluripotent stem cells and generated brain organoids. These organoids showed alterations in gene expression and splicing as well as morphology and synaptogenesis, suggesting that this method could be used to explore other genetic changes that underlie the phenotypic traits separating our species from extinct relatives. Science , this issue p. [eaax2537][1] [1]: /lookup/doi/10.1126/science.aax2537

Novel insights into non-HLA alloimmunity in kidney transplantation.

SummaryRecognition of non‐self structures on donor cells represents the main immunological barrier in solid organ transplantation. The human leukocyte antigens (HLA) are considered the most important non‐self (allo)antigens in transplantation. Long‐term graft attrition is mainly caused by the formation of alloreactive antibodies that are directed against non‐self structures (i.e., epitopes) on cell surface proteins. Recently published data provided evidence for a similar importance of non‐HLA mismatches between donors and recipients in acute rejection as well as long‐term kidney allograft survival. These data suggest a broader concept of immunological non‐self that goes beyond HLA incompatibility and expands the current concept of polymorphic non‐self epitopes on cell surface molecules from HLA to non‐HLA targets. Amino acid substitutions caused by single nucleotide variants in protein‐coding genes or complete loss of gene expression represent the basis for polymorphic residues in both HLA and non‐HLA molecules. To better understand these novel insights in non‐HLA alloimmunity, we will first review basic principles of the alloimmune response with a focus on the HLA epitope concept in donor‐specific antibody formation before discussing key publications on non‐HLA antibodies.

The Mir-99b/Let-7e/Mir-125a cluster Is Associated with Platelet Count and Mean Platelet Volume and Regulates Megakaryocyte Proplatelet Formation

Platelet count (PC) and mean platelet volume (MPV) are critical measurements when considering the role of platelets in hemostasis and thrombosis, abnormalities in PC and MPV have been shown to be predictors of ischemic cardiovascular disease. Ischemic vascular disease has a known genetic component, and we and others have shown high heritability for both PC and MPV. Many inherited, and most immune, thrombocytopenias show an inverse relationship between PC and MPV, and such a relationship has also been reported in healthy subjects. Although genome-wide association studies have identified numerous single nucleotide variants in protein-coding genes associated with PC and MPV (Geiger Nat Genet 2011), non-coding genes have had rather little consideration. The goal of our studies was to identify microRNAs (miRs) regulating PC and MPV and to begin to understand the molecular mechanisms of this inverse platelet number-volume relationship. Initially, we queried our Platelet RNA and Expression-1 data set and observed that PC and MPV are inversely correlated in these 154 healthy subjects (r=-0.45, p=3.98E-09). We next determined which of the 734 platelet-expressed miRs were significantly (p&lt;0.05) correlated with both PC and MPV. Of these miRs, only miR-125a-5p, miR-99b-5p, let-7e-5p and miR-150-5p had an inverse relationship with PC and MPV. The Table shows the correlations between these miRs and complete blood count data. Subsequent studies focused on miR-125a-5p, miR-99b-5p, let-7e-5p because they (1) were associated with only platelet parameters and not WBC or Hb, (2) were more highly expressed than miR-150-5p and (3) showed the same directionality (e.g., positive r value for PC). We were also intrigued by the fact that these 3 miRs exist as a genomic cluster within 727 bp on chr19q13.41. The expression of each of the 3 miRs in this cluster was strongly correlated with each other in platelets (r=0.93, p=1.02E-68 for miR-125a and let-7e; r=0.91, p=8.53E-61 for miR-99b and let-7e; r=0.89, p=1.19E-52 for miR-125a and miR-99b), consistent with a single, polycistronic transcript containing all genes. To begin to address the potential functionality of this miR cluster, we overexpressed and knocked down their expression in megakaryocytes (MKs), since MK proplatelet formation (PPF) is a necessary step in the generation of peripheral blood platelets. We utilized human umbilical cord blood derived CD34-positive hematopoietic stem and progenitor cells differentiated along the MK lineage with thrombopoietin. Quantification of miR levels with qPCR demonstrated all 3 miRs in the cluster progressively increased from culture day 6 to day 9 to day 13 (the time during which MKs mature in culture). Compared to a scrambled control, when miR activity was knocked down with locked nucleic acids (LNAs), we observed a 61% reduction in PPF for miR-125a-5p (p&lt;0.05), a 54% reduction for let-7e-5p (p&lt;0.05) and a ~28% reduction for miR-99b-5p (p=0.33). PPF was quantified blinded as to treatment group. When all 3 miRs in the cluster were inhibited simultaneously by LNAs, we observed a 74% reduction in PPF (p&lt;0.05) (Figure 1). LNA inhibition of two miRs not associated with PC or MPV did not alter MK PFF. To date, we have only tested the effect of overexpression of one miR, but observed that lentiviral overexpression of miR-125a-5p induced an increase in MK PPF (p&lt;0.05). In summary, we have identified a genomic cluster of miRs (miR-99b/let-7e/miR-125a) that regulates PPF in cultured human MKs and whose expression is positively correlated with human PC and negatively correlated with MPV. We speculate that platelet mass homeostasis, which is determined by PC and MPV, may be controlled in part by this novel genomic locus. Disclosures No relevant conflicts of interest to declare.

Detecting the Population Structure and Scanning for Signatures of Selection in Horses (Equus caballus) From Whole-Genome Sequencing Data.

Animal domestication gives rise to gradual changes at the genomic level through selection in populations. Selective sweeps have been traced in the genomes of many animal species, including humans, cattle, and dogs. However, little is known regarding positional candidate genes and genomic regions that exhibit signatures of selection in domestic horses. In addition, an understanding of the genetic processes underlying horse domestication, especially the origin of Chinese native populations, is still lacking. In our study, we generated whole genome sequences from 4 Chinese native horses and combined them with 48 publicly available full genome sequences, from which 15 341 213 high-quality unique single-nucleotide polymorphism variants were identified. Kazakh and Lichuan horses are 2 typical Asian native breeds that were formed in Kazakh or Northwest China and South China, respectively. We detected 1390 loss-of-function (LoF) variants in protein-coding genes, and gene ontology (GO) enrichment analysis revealed that some LoF-affected genes were overrepresented in GO terms related to the immune response. Bayesian clustering, distance analysis, and principal component analysis demonstrated that the population structure of these breeds largely reflected weak geographic patterns. Kazakh and Lichuan horses were assigned to the same lineage with other Asian native breeds, in agreement with previous studies on the genetic origin of Chinese domestic horses. We applied the composite likelihood ratio method to scan for genomic regions showing signals of recent selection in the horse genome. A total of 1052 genomic windows of 10 kB, corresponding to 933 distinct core regions, significantly exceeded neutral simulations. The GO enrichment analysis revealed that the genes under selective sweeps were overrepresented with GO terms, including “negative regulation of canonical Wnt signaling pathway,” “muscle contraction,” and “axon guidance.” Frequent exercise training in domestic horses may have resulted in changes in the expression of genes related to metabolism, muscle structure, and the nervous system.

Using ALoFT to determine the impact of putative loss-of-function variants in protein-coding genes

Variants predicted to result in the loss of function of human genes have attracted interest because of their clinical impact and surprising prevalence in healthy individuals. Here, we present ALoFT (annotation of loss-of-function transcripts), a method to annotate and predict the disease-causing potential of loss-of-function variants. Using data from Mendelian disease-gene discovery projects, we show that ALoFT can distinguish between loss-of-function variants that are deleterious as heterozygotes and those causing disease only in the homozygous state. Investigation of variants discovered in healthy populations suggests that each individual carries at least two heterozygous premature stop alleles that could potentially lead to disease if present as homozygotes. When applied to de novo putative loss-of-function variants in autism-affected families, ALoFT distinguishes between deleterious variants in patients and benign variants in unaffected siblings. Finally, analysis of somatic variants in >6500 cancer exomes shows that putative loss-of-function variants predicted to be deleterious by ALoFT are enriched in known driver genes.

A GJA9 frameshift variant is associated with polyneuropathy in Leonberger dogs

BackgroundMany inherited polyneuropathies (PN) observed in dogs have clinical similarities to the genetically heterogeneous group of Charcot-Marie-Tooth (CMT) peripheral neuropathies in humans. The canine disorders collectively show a variable expression of progressive clinical signs and ages of onset, and different breed prevalences. Previously in the Leonberger breed, a variant highly associated with a juvenile-onset PN was identified in the canine orthologue of a CMT-associated gene. As this deletion in ARHGEF10 (termed LPN1) does not explain all cases, PN in this breed may encompass variants in several genes with similar clinical and histopathological features.ResultsA genome-wide comparison of 173 k SNP genotypes of 176 cases, excluding dogs homozygous for the ARHGEF10 variant, and 138 controls, was carried out to detect further PN-associated variants. A single suggestive significant association signal on CFA15 was found. The genome of a PN-affected Leonberger homozygous for the associated haplotype was sequenced and variants in the 7.7 Mb sized critical interval were identified. These variants were filtered against a database of variants observed in 202 genomes of various dog breeds and 3 wolves, and 6 private variants in protein-coding genes, all in complete linkage disequilibrium, plus 92 non-coding variants were revealed. Five of the coding variants were predicted to have low or moderate effect on the encoded protein, whereas a 2 bp deletion in GJA9 results in a frameshift of high impact. GJA9 encodes connexin 59, a connexin gap junction family protein, and belongs to a group of CMT-associated genes that have emerged as important components of peripheral myelinated nerve fibers. The association between the GJA9 variant and PN was confirmed in an independent cohort of 296 cases and 312 controls. Population studies showed a dominant mode of inheritance, an average age of onset of approximately 6 years, and incomplete penetrance.ConclusionsThis GJA9 variant represents a highly probable candidate variant for another form of PN in Leonberger dogs, which we have designated LPN2, and a new candidate gene for CMT disease. To date, approximately every third PN-diagnosed Leonberger dog can be explained by the ARHGEF10 or GJA9 variants, and we assume that additional genetic heterogeneity in this condition exists in the breed.

Analysis of protein-coding genetic variation in 60,706 humans.

SummaryLarge-scale reference data sets of human genetic variation are critical for the medical and functional interpretation of DNA sequence changes. We describe the aggregation and analysis of high-quality exome (protein-coding region) sequence data for 60,706 individuals of diverse ethnicities generated as part of the Exome Aggregation Consortium (ExAC). This catalogue of human genetic diversity contains an average of one variant every eight bases of the exome, and provides direct evidence for the presence of widespread mutational recurrence. We have used this catalogue to calculate objective metrics of pathogenicity for sequence variants, and to identify genes subject to strong selection against various classes of mutation; identifying 3,230 genes with near-complete depletion of truncating variants with 72% having no currently established human disease phenotype. Finally, we demonstrate that these data can be used for the efficient filtering of candidate disease-causing variants, and for the discovery of human “knockout” variants in protein-coding genes.

Deep sequencing of Danish Holstein dairy cattle for variant detection and insight into potential loss-of-function variants in protein coding genes.

BackgroundOver the last few years, continuous development of high-throughput sequencing platforms and sequence analysis tools has facilitated reliable identification and characterization of genetic variants in many cattle breeds. Deep sequencing of entire genomes within a cattle breed that has not been thoroughly investigated would be imagined to discover functional variants that are underlying phenotypic differences. Here, we sequenced to a high coverage the Danish Holstein cattle breed to detect and characterize single nucleotide polymorphisms (SNPs), insertion/deletions (Indels), and loss-of-function (LoF) variants in protein-coding genes in order to provide a comprehensive resource for subsequent detection of causal variants for recessive traits.ResultsWe sequenced four genetically unrelated Danish Holstein cows with a mean coverage of 27X using an Illumina Hiseq 2000. Multi-sample SNP calling identified 10,796,794 SNPs and 1,295,036 indels whereof 482,835 (4.5 %) SNPs and 231,359 (17.9 %) indels were novel. A comparison between sequencing-derived SNPs and genotyping from the BovineHD BeadChip revealed a concordance rate of 99.6–99.8 % for homozygous SNPs and 93.3–96.5 % for heterozygous SNPs. Annotation of the SNPs discovered 74,886 SNPs and 1937 indels affecting coding sequences with 2145 being LoF mutations. The frequency of LoF variants differed greatly across the genome, a hot spot with a strikingly high density was observed in a 6 Mb region on BTA18. LoF affected genes were enriched for functional categories related to olfactory reception and underrepresented for genes related to key cellular constituents and cellular and biological process regulation. Filtering using sequence derived genotype data for 288 Holstein animals from the 1000 bull genomes project removing variants containing homozygous individuals retained 345 of the LoF variants as putatively deleterious. A substantial number of the putative deleterious LoF variants had a minor allele frequency >0.05 in the 1000 bull genomes data set.ConclusionsDeep sequencing of Danish Holstein genomes enabled us to identify 12.1 million variants. An investigation into LoF variants discovered a set of variants predicted to disrupt protein-coding genes. This catalog of variants will be a resource for future studies to understand variation underlying important phenotypes, particularly recessively inherited lethal phenotypes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2249-y) contains supplementary material, which is available to authorized users.

Exome Capture with Heterologous Enrichment in Pig (Sus scrofa).

The discovery of new protein-coding DNA variants related to carcass traits is very important for the Italian pig industry, which requires heavy pigs with higher thickness of subcutaneous fat for Protected Designation of Origin (PDO) productions. Exome capture techniques offer the opportunity to focus on the regions of DNA potentially related to the gene and protein expression. In this research a human commercial target enrichment kit was used to evaluate its performances for pig exome capture and for the identification of DNA variants suitable for comparative analysis. Two pools of 30 pigs each, crosses of Italian Duroc X Large White (DU) and Commercial hybrid X Large White (HY), were used and NGS libraries were prepared with the SureSelectXT Target Enrichment System for Illumina Paired-End Sequencing Library (Agilent). A total of 140.2 M and 162.5 M of raw reads were generated for DU and HY, respectively. Average coverage of all the exonic regions for Sus scrofa (ENSEMBL Sus_scrofa.Sscrofa10.2.73.gtf) was 89.33X for DU and 97.56X for HY; and 35% of aligned bases uniquely mapped to off-target regions. Comparison of sequencing data with the Sscrofa10.2 reference genome, after applying hard filtering criteria, revealed a total of 232,530 single nucleotide variants (SNVs) of which 20.6% mapped in exonic regions and 49.5% within intronic regions. The comparison of allele frequencies of 213 randomly selected SNVs from exome sequencing and the same SNVs analyzed with a Sequenom MassARRAY® system confirms that this “human-on-pig” approach offers new potentiality for the identification of DNA variants in protein-coding genes.

Many Inflammatory Bowel Disease Risk Loci Include Regions That Regulate Gene Expression in Immune Cells and the Intestinal Epithelium

The contribution of genetic factors to the pathogenesis of inflammatory bowel disease (IBD) has been established by twin, targeted-sequencing, and genome-wide association studies. These studies identified many risk loci, and research is underway to identify causal variants. These studies have focused mainly on protein-coding genes. We investigated other functional elements in the human genome, such as regulatory regions. Using acetylated histone 3 lysine 27 chromatin immunoprecipitation and sequencing, we identified tens of thousands of potential regulatory regions that are active in intestinal epithelium (primary intestinal crypts and cultured organoids) isolated from resected material and from biopsies collected during ileo-colonoscopies and immune cells (monocytes, macrophages, CD34(+), CD4(+), and CD8(+)). We correlated these regions with susceptibility loci for IBD. We have generated acetylated histone 3 lysine 27 profiles from primary intestinal epithelium and cultured organoids, which we have made publically available. We found that 45 of 163 single nucleotide polymorphisms (SNPs) associated with IBD overlap specifically with active regulatory elements. In addition, by taking strong linkage disequilibrium into account, another 47 IBD-associated SNPs colocalized with active regulatory elements through other SNPs in their vicinity. Altogether, 92 of 163 IBD-associated SNPs correlated with distinct active regulatory elements-a frequency 2.5- to 3.5-fold greater than that expected from random sampling. The variations in these SNPs often create or disrupt known binding motifs; they might affect the binding of transcriptional regulators to alter expression of regulated genes. In addition to variants in protein coding genes, variants in noncoding DNA regulatory regions that are active in intestinal epithelium and immune cells are potentially involved in the pathogenesis of IBD.