Targeted imputation of sequence variants and gene expression profiling identifies twelve candidate genes associated with lactation volume,composition and calving interval in dairy cattle

Dairy cattle are an interesting model for gaining insights into the genes responsible for the large variation between and within mammalian species in the protein and fat content of their milk and their milk volume. Large numbers of phenotypes for these traits are available, as well as full genome sequence of key founders of modern dairy cattle populations. In twenty target QTL regions affecting milk production traits, we imputed full genome sequence variant genotypes into a population of 16,721 Holstein and Jersey cattle with excellent phenotypes. Association testing was used to identify variants within each target region, and gene expression data were used to identify possible gene candidates. There was statistical support for imputed sequence variants in or close to BTRC, MGST1, SLC37A1, STAT5A, STAT5B, PAEP, VDR, CSF2RB, MUC1, NCF4, and GHDC associated with milk production, and EPGN for calving interval. Of these candidates, analysis of RNA-Seq data demonstrated that PAEP, VDR, SLC37A1, GHDC, MUC1, CSF2RB, and STAT5A were highly differentially expressed in mammary gland compared to 15 other tissues. For nine of the other target regions, the most significant variants were in non-coding DNA. Genomic predictions in a third dairy breed (Australian Reds) using sequence variants in only these candidate genes were for some traits more accurate than genomic predictions from 632,003 common SNP on the Bovine HD array. The genes identified in this study are interesting candidates for improving milk production in cattle and could be investigated for novel biological mechanisms driving lactation traits in other mammals.     

κ-casein gene (<Emphasis Type="Italic">CSN3</Emphasis>) allelic polymorphism in Russian cattle breeds and its information value as a genetic marker

The frequencies of the κ-casein gene (CSN3) alleles and genotypes have been determined in five Russian cattle breeds (Bestuzhev, Kalmyk, Russian Black Pied, Yaroslavl, and Yakut breeds) by means of PCR-RFLP analysis using two independent restriction nucleases (HinfI and TaqI) and by allele-specific PCR. Typing alleles A and B of CSN3 is of practical importance, because allele B is correlated with commercially valuable parameters of milk productivity (protein content and milk yield) and improves the cheese yielding capacity. The frequencies of the B allele of CSN3 in the breeds studied vary from 0.16 to 0.50; and those of the AB and BB genotypes, from 0.27 to 0.60 and from 0.02 to 0.23, respectively. The Yaroslavl breed had the highest frequencies of CSN3 allele B and genotype BB (0.50 and 0.23, respectively). The frequencies of the B allele and BB genotype in other breeds studied varied from 0.25 to 0.32 and from 0.03 to 0.09, respectively. In none of the breeds studied have the observed and expected heterozygosities been found to differ from each other significantly. However, the observed genotype distributions significantly differ from the expected one in some herds (in most such cases, an excess of heterozygotes is observed). Two herds of the Yaroslavl breed dramatically differ from each other in the heterozygosity level: a deficit (D = ?0.14) and an excess (D = 0.20) of heterozygotes have been observed at the Mikhailovskoe and Gorshikha farms, respectively. In general, however, the heterozygosity of the Yaroslavl breed corresponds to the expected level (D = 0.04). Analysis of breeds for homogeneity with the use of Kulback’s test has shown that all cattle breeds studied are heterogeneous, the CSN3 diversity within breeds being higher than that among different breeds, which is confirmed by low F _st values (0.0025–0.0431). Thus, a DNA marker based on CSN3 gene polymorphism is extremely important for breeding practice as a marker of milk quality; however, it is inapplicable to marking differences between breeds or phylogenetic relationships between cattle breeds because of the high diversity with respect to this locus within breeds.     

Genome-wide association study of <Emphasis Type="Italic">Mycobacterium avium</Emphasis> subspecies <Emphasis Type="Italic">Paratuberculosis</Emphasis> infection in Chinese Holstein

Background

Paratuberculosis is a contagious, chronic and enteric disease in ruminants, which is caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection, resulting in enormous economic losses worldwide. There is currently no effective cure for MAP infection or a vaccine, it is thus important to explore the genetic variants that contribute to host susceptibility to infection by MAP, which may provide a better understanding of the mechanisms of paratuberculosis and benefit animal genetic improvement. Herein we performed a genome-wide association study (GWAS) to identify genomic regions and candidate genes associated with susceptibility to MAP infection in dairy cattle.

Results

Using Illumina Bovine 50?K (54,609 SNPs) and GeneSeek HD (138,893 SNPs) chips, two analytical approaches were performed, GRAMMAR-GC and ROADTRIPS in 937 Chinese Holstein cows, among which individuals genotyped by the 50?K chip were imputed to HD SNPs with Beagle software. Consequently, 15 and 11 significant SNPs (P?<?5?×?10^??5) were identified with GRAMMAR-GC and ROADTDRIPS, respectively. A total of 10 functional genes were in proximity to (i.e., within 1?Mb) these SNPs, including IL4, IL5, IL13, IRF1, MyD88, PACSIN1, DEF6, TDP2, ZAP70 and CSF2. Functional enrichment analysis showed that these genes were involved in immune related pathways, such as interleukin, T cell receptor signaling pathways and inflammatory bowel disease (IBD), implying their potential associations with susceptibility to MAP infection. In addition, by examining the publicly available cattle QTLdb, a previous QTL for MAP was found to be overlapped with one of regions detected currently at 32.5?Mb on BTA23, where the TDP2 gene was anchored.

Conclusions

In conclusion, we identified 26 SNPs located on 15 chromosomes in the Chinese Holstein population using two GWAS strategies with high density SNPs. Integrated analysis of GWAS, biological functions and the reported QTL information helps to detect positional candidate genes and the identification of regions associated with susceptibility to MAP traits in dairy cattle.

     

Novel SNPs of the bovine <Emphasis Type="Italic">NUCB2</Emphasis> gene and their association with growth traits in three native Chinese cattle breeds

In this study, polymorphism in the exon 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 and 11 of bovine NUCB2 gene was detected by PCR-SSCP and DNA sequencing methods in 686 individuals from three Chinese cattle breeds. Two haplotypes (M and N), three observed genotypes (MM, MN and NN) and two SNPs (NC_007313: g. 27451G>A, NC_007313: g. 27472T>C) were detected. The frequencies of haplotypes M and N in inland Chinese three breeds were 0.531–0.721 and 0.279–0.469 respectively. The studied showed that Nanyang, Jiaxian Red and Qinchuan cattle populations were in Hardy–Weinberg equilibrium at SNPs locus of NUCB2 gene (P > 0.05). Polymorphism of the NUCB2 gene was shown to be associated with growth traits in Qingchuan and Nanyang cattle breed. The linkage of two mutant sites in the bovine NUCB2 gene had significant effects on body length, body weight, heart girth, and average daily gain at 24 months (P < 0.05). Results of this study suggested that the NUCB2-gene-specific SNP may be a useful marker for growth traits in future marker-assisted selection programmes in inland Chinese cattle.     

The effects of polymorphisms in <Emphasis Type="Italic">growth hormone</Emphasis> and <Emphasis Type="Italic">growth hormone receptor</Emphasis> genes on production and reproduction traits in Aberdeen-Angus cattle (<Emphasis Type="Italic">Bos taurus</Emphasis> L., 1758)

The study was aimed to analyze the relation between individual genotypes and allelic variants of SNPs g.2141C>G of growth hormone gene, g.914T>A and g.257A>G of growth hormone receptor gene with growth and reproduction traits and to evaluate the populationgenetic structure in Aberdeen-Angus cattle (Bos taurus L., 1758) sample of Eastern Ukraine according SNPs studied. Allele C of SNP g.2141C>G has a positive correlation with birth weight, body stature, bigger rump, udder and total exterior evaluation score, shorter calving interval and better calve birth weight and negative correlation with calve average daily gain. Allele T of SNP g.914T>A has positive correlation with the muscle and udder size; live weight in each age, average daily gain, weight and average daily gain of calves born conform to the principle AA>TT≥TA. SNP g.257A>G showed a positive correlation for G allele with muscle size. The population is in equilibrium for SNPs g.2141C>G and g.257A>G, and in disequilibrium for SNP g.914T>A. The analysis showed no linkage disequilibrium between SNPs g.914T>A and g.257A>G. Inbreeding coefficient F_ST in Aberdeen-Angus group studied was 16.1%.     

Rank-based genome-wide analysis reveals the association of Ryanodine receptor-2 gene variants with childhood asthma among human populations

Background

The standard approach to determine unique or shared genetic factors across populations is to identify risk alleles in one population and investigate replication in others. However, since populations differ in DNA sequence information, allele frequencies, effect sizes, and linkage disequilibrium patterns, SNP association using a uniform stringent threshold on p values may not be reproducible across populations. Here, we developed rank-based methods to investigate shared or population-specific loci and pathways for childhood asthma across individuals of diverse ancestry. We performed genome-wide association studies on 859,790 SNPs genotyped in 527 affected offspring trios of European, African, and Hispanic ancestry using publically available asthma database in the Genotypes and Phenotypes database.

Results

Rank-based analyses showed that there are shared genetic factors for asthma across populations, more at the gene and pathway levels than at the SNP level. Although the top 1,000 SNPs were not shared, 11 genes (RYR2, PDE4D, CSMD1, CDH13, ROBO2, RBFOX1, PTPRD, NPAS3, PDE1C, SEMA5A, and CTNNA2) mapped by these SNPs were shared across populations. Ryanodine receptor 2 (RYR2, a statin response-related gene) showed the strongest association in European (p value?=?2.55?×?10^?7) and was replicated in African (2.57?×?10^?4) and Hispanic (1.18?×?10^?3) Americans. Imputation analyses based on the 1000 Genomes Project uncovered additional RYR2 variants associated with asthma. Network and functional ontology analyses revealed that RYR2 is an integral part of dermatological or allergic disorder biological networks, specifically in the functional classes involving inflammatory, eosinophilic, and respiratory diseases.

Conclusion

Our rank-based genome-wide analysis revealed for the first time an association of RYR2 variants with asthma and replicated previously discovered PDE4D asthma gene across human populations. The replication of top-ranked asthma genes across populations suggests that such loci are less likely to be false positives and could indicate true associations. Variants that are associated with asthma across populations could be used to identify individuals who are at high risk for asthma regardless of genetic ancestry.

     

Sex influences eQTL effects of SLE and Sjögren’s syndrome-associated genetic polymorphisms

Background

Systemic lupus erythematosus (SLE) and primary Sjögren’s syndrome (pSS) are autoimmune disorders characterized by autoantibodies, dysregulated B cells, and notably high female-to-male incidence ratios. Genome-wide association studies have identified several susceptibility SNPs for both diseases. Many SNPs in the genome are expression quantitative trait loci (eQTLs), with context-dependent effects. Assuming that sex is a biological context, we investigated whether SLE/pSS SNPs act as eQTLs in B cells and used a disease-targeted approach to understand if they display sex-specific effects.

Methods

We used genome-wide genotype and gene expression data from primary B cells from 125 males and 162 females. The MatrixEQTL R package was used to identify eQTLs within a genomic window of 2 Mb centered on each of 22 established SLE and/or pSS susceptibility SNPs. To find sex-specific eQTLs, we used a linear model with a SNP * sex interaction term.

Results

We found ten SNPs affecting the expression of 16 different genes (FDR <?0.05). rs7574865-INPP1, rs7574865-MYO1B, rs4938573-CD3D, rs11755393-SNRPC, and rs4963128-PHRF1 were novel observations for the immune compartment and B cells. By analyzing the SNP * sex interaction terms, we identified six genes with differentially regulated expression in females compared to males, depending on the genotype of SLE/pSS-associated SNPs: SLC39A8 (BANK1 locus), CD74 (TNIP1 locus), PXK, CTSB (BLK/FAM167A locus), ARCN1 (CXCR5 locus), and DHX9 (NCF2 locus).

Conclusions

We identified several unknown sex-specific eQTL effects of SLE/pSS-associated genetic polymorphisms and provide novel insight into how gene-sex interactions may contribute to the sex bias in systemic autoimmune diseases.

     

Sequence polymorphisms in <Emphasis Type="Italic">Zmisa2</Emphasis> gene are significantly associated with starch pasting and gelatinization properties in maize (<Emphasis Type="Italic">Zea mays</Emphasis> L.)

Researches on the genetic basis of starch pasting and gelatinization properties will provide foundation of maize improvement for quality, feed and industrial applications. Maize gene Zmisa2, encoding an isoamylase-type starch debranching enzyme, plays important roles in starch biosynthesis. In this study, the genomic sequences of the gene Zmisa2 in 72 elite maize inbred lines were obtained, and the nucleotide polymorphisms and haplotype diversity were detected. In addition, seven pasting and four gelatinization properties of maize were measured for the tested inbred lines using rapid visco analyzer and differential scanning calorimeter, respectively. A total of 99 sequence variants, including 91 SNPs and 8 indels, were identified at the promoter and coding regions of this gene. Although the frequency of polymorphism in promoter region is much higher than that of coding region, the SNPs in the coding region of maize gene Zmisa2 classified this gene into 21 haplotypes, which encode 11 different ISA2 proteins. Furthermore, the association of the variants of Zmisa2 gene with maize starch pasting and gelatinization properties was estimated, and the results revealed that seven SNPs in coding region, including four nonsynonymous sites, were significantly associated with phenotypic variations of pasting time and enthalpy of transition (ΔH). These results suggested that the polymorphism in maize Zmisa2 locus could be used in molecular marker-assisted selection for improvement of quality in maize breeding programs.     

Genetic risk factors for myocardial infarction more clearly manifest for early age of first onset

Epidemiological genetics established that heritability in determining the risk of myocardial infarction (MI) is substantially greater when MI occurs early in life. However, the genetic architecture of early-onset and late-onset MI was not compared. We analyzed genotype frequencies of SNPs in/near 20 genes whose protein products are involved in the pathogenesis of atherosclerosis in two groups of Russian patients with MI: the first group included patients with age of first MI onset <60 years (N?=?230) and the second group with onset ≥60 years (N?=?174). The control group of corresponding ethnicity consisted of 193 unrelated volunteers without cardiovascular diseases (93 individuals were over 60 years). We found that in the group of patients with age of onset <60 years, SNPs FGB rs1800788*T, TGFB1 rs1982073*T/T, ENOS rs2070744*C and CRP rs1130864*T/T were associated with risk of MI, whereas in patients with age of onset ≥60 years, only TGFB1 rs1982073*T/T was associated with risk of MI. Using APSampler software, we found composite markers associated with MI only in patients with early onset: FGB rs1800788*T?+?TGFB1 rs1982073*T; FGB rs1800788*T?+?LPL rs328*C?+?IL4 rs2243250*C; FGB rs1800788*T?+?ENOS rs2070744*C (Fisher p values of 1.4?×?10^?6 to 2.2?×?10^?5; the permutation p values of 1.1?×?10^?5 to 3.0?×?10^?4; ORs?=?2.67–2.54). Alleles included in the combinations were associated with MI less significantly and with lower ORs than the combinations themselves. The result showed a substantially greater contribution of the genetic component in the development of MI if it occurs early in life, and demonstrated the usefulness of genetic testing for young people.     

Results of <Emphasis Type="Italic">Camelus dromedarius</Emphasis> and <Emphasis Type="Italic">Camelus bactrianus</Emphasis> Genotyping by Alpha-S<Subscript>1</Subscript>-Casein,Kappa-Casein Loci,and DNA Fingerprinting

Genotyping of Kazakh camels Camelus dromedarius (milk breed) (n = 18) and Camelus bactrianus (meat breed) (n = 18) by alpha-S₁-casein (αs₁-CN) and kappa-casein (κ-CN) loci was conducted using the PCR–RFLP analysis method. A new pair of primers was suggested for the amplification of the CSN3 gene fragment with subsequent cleavage of the reaction products by AluI restriction endonuclease in order to identify the gene genetic variants. DNA polymorphism was detected only for the kappa-casein locus; no genetic polymorphism for alpha-S₁-casein gene was found in the studied populations. Analysis of the results of DNA fingerprinting demonstrated that the band sharing (BS) coefficient between the groups was low enough (0.13), and the genetic distance (D) between Dromedary and Bactrian breeds was 0.305. The results of genotyping of Bactrian and Dromedary Kazakh camel breeds by alpha-S₁-casein, kappa-casein loci, and DNA fingerprinting indicate that the Dromedary breed female camels are more polymorphic as compared with Bactrian.     

Polymorphism of the genes for antioxidant defense enzymes and their association with the development of chronic obstructive pulmonary disease in the population of Bashkortostan

In this study, frequencies of the polymorphic variants of the genes encoding antioxidant enzymes, GSTM1, GSTT1, GSTP1, CAT, GPX1, NQO1, SOD1, and SOD3 were examined in three ethnic groups of healthy subjects from the Republic of Bashkortostan (Russians, Tatars, and Bashkirs). An association of these markers with the development of chronic obstructive pulmonary disease (COPD) was tested. Interethnic differences relative to the distribution of the polymorphic variants of the GSTP1 locus Ile105Val and the NQO1 locus 609C/T were revealed. Relative to the genotype distribution at the Ile105Val locus of the GSTP1 gene, ethnic group of Bashkirs was found to be statistically significantly different from Tatars (χ² = 8.819; d.f. = 2; P = 0.012). Relative to the genotype frequency distribution pattern at the NQO1 locus 609C/T, the group of Bashkirs differed from Russians (χ² = 8.913; d.f. = 2; P = 0.012). An association of genotype Val/Val of the GSTP1 Ile105Val locus with the risk of COPD in Russians (χ² = 5.25; P = 0.022; P _cor = 0.044; OR = 4.09), and of the GSTP1 haplotype *D in Tatars, was demonstrated (χ² = 11.575; P = 0.0014; P _cor = 0.0042; OR = 3.178). Genotype TT of the CAT ?262C/T locus marked resistance to the COPD development in Russians (χ² = 6.82; P = 0.0098; P _cor = 0.0196; OR = 0.31; 95%CI, 0.119 to 0.77). The risk for COPD in the ethnic group of Tatars was associated with the CAT haplotype (?262)C/(1167)T (χ² = 6.038; P = 0.0147; P _cor = 0.044; OR = 1.71). Analysis of the NQO1 haplotypes at the 465C/T and 609C/T loci showed that haplotype 465C/609T was associated with COPD in Russians (χ² = 4.571; P = 0.0328; P _cor = 0.01; OR = 1.799). It was demonstrated that Gly allele of the Arg213Gly polymorphic locus of the SOD3 gene marked the risk for COPD in the ethnic group of Tatars (OR = 2.23; 95%CI, 1.22 to 4.1). Thus, GSTP1, CAT, NQO1, and SOD3 polymorphisms play an important role in the development of COPD among the population of Bashkortostan.     

Promoters pro-SmAMP1 and pro-SmAMP2 from Wild Plant <Emphasis Type="Italic">Stellaria media</Emphasis> for the Biotechnology of Dicotyledons

Comparative results of the studied effectiveness of two new promoters, pro-SmAMP1 and pro- SmAMP2, from chickweed (Stellaria media L.) in various types of cultivated plants with transient expression and in stable transformants are given. The effectiveness of the promoters was evaluated through the expression of the reporter uidA gene by measuring the activity of its GUS protein product. It was found that the deletion variant (442 bp) of the pro-SmAMP1 promoter was significantly stronger in plants of Nicotiana benthamiana (Domin) with transient expression than the deletion variant (455 bp) of the pro-SmAMP2 promoter. The effectiveness of these short deletion variants of both promoters under transient expression in the plants of rapeseed (Brassica napus L.) and sugar beet (Beta vulgaris L.) was comparable with that of the viral CaMV35S promoter. The functionality of the pro-SmAMP2 promoter in the calluses of common flax plants (Linum usitatissimum L.) was shown. In the homozygous lines of transgenic tobacco plants (Nicotiana tabacum L.), all deletion variants of the pro-SmAMP1 promoter and the shortest version of pro-SmAMP2 were twice as strong as the CaMV35S viral promoter. The effectiveness of short variants of both promoters from the chickweed in controlling the gene encoding neomycin phosphotransferase II in the transgenic plants of tobacco and arabidopsis (Arabidopsis thaliana L.) growing on media supplemented with recommended concentrations of kanamycin are not inferior to the duplicated 2хCaMV35S viral promoter. The obtained experimental data show that short deletion variants of pro-SmAMP1 (442 bp) and pro-SmAMP2 (455 bp) promoters may be recommended as strong constitutive promoters for use in the biotechnology of crop plants.     

A survey of endogenous retrovirus (ERV) sequences in the vicinity of multiple sclerosis (MS)-associated single nucleotide polymorphisms (SNPs)

Although multiple sclerosis (MS) is one of the most common central nervous system diseases in young adults, little is known about its etiology. Several human endogenous retroviruses (ERVs) are considered to play a role in MS. We are interested in which ERVs can be identified in the vicinity of MS associated genetic marker to find potential initiators of MS. We analysed the chromosomal regions surrounding 58 single nucleotide polymorphisms (SNPs) that are associated with MS identified in one of the last major genome wide association studies. We scanned these regions for putative endogenous retrovirus sequences with large open reading frames (ORFs). We observed that more retrovirus-related putative ORFs exist in the relatively close vicinity of SNP marker indices in multiple sclerosis compared to control SNPs. We found very high homologies to HERV-K, HCML-ARV, XMRV, Galidia ERV, HERV-H/env62 and XMRV-like mouse endogenous retrovirus mERV-XL. The associated genes (CYP27B1, CD6, CD58, MPV17L2, IL12RB1, CXCR5, PTGER4, TAGAP, TYK2, ICAM3, CD86, GALC, GPR65 as well as the HLA DRB1*1501) are mainly involved in the immune system, but also in vitamin D regulation. The most frequently detected ERV sequences are related to the multiple sclerosis-associated retrovirus, the human immunodeficiency virus 1, HERV-K, and the Simian foamy virus. Our data shows that there is a relation between MS associated SNPs and the number of retroviral elements compared to control. Our data identifies new ERV sequences that have not been associated with MS, so far.     

Characteristics of a transverse RF discharge in Xe/Cl<Subscript>2</Subscript> mixtures

Results are presented from the study of the electrical and optical characteristics of a transverse RF discharge in Xe/Cl₂ mixtures at pressures of p≤400 Pa. The working mixture was excited by a modulated RF discharge (f=1.76 MHz) with a transverse electrode configuration (L≤17 cm). The emission spectrum in the spectral range of 210–600 nm and the waveforms of the discharge current, discharge voltage, and plasma emission intensity were investigated. The UV emission power from the discharge was studied as a function of the pressure and composition of a Xe/Cl₂ mixture. It is shown that a discharge in a xenon-chlorine mixture acts as planar excimer-halogen lamp operating in the spectral range of 220–450 nm, which contains a system of overlapping XeCl(D, B-X; B, C-A) and Cl₂(D′-A′) bands. Transverse RF discharges in Xe/Cl₂ mixtures can be used to create a wideband lamp with two 50-cm² planar apertures and the low circulation rate of the working mixture.     

A genome-wide association study on growth traits in orange-spotted grouper (<Emphasis Type="Italic">Epinephelus coioides</Emphasis>) with RAD-seq genotyping

The orange-spotted grouper, Epinephelus coioides, is one of the most popular fish in China and Southeast Asian countries because of its important economic value. However, molecular mechanism underlying the growth of orange-spotted grouper has never been fully understood. Herein, we performed a genome-wide association study (GWAS) on a natural population of 198 individuals aiming to screen the whole genome of orange-spotted grouper for identification of growth-related loci by restrictionsite associated DNA sequencing. In this research, 261,366 single nucleotide polymorphisms (SNPs) were developed, in which 110 SNPs were identified to be correlated with growth and 20 SNPs were further confirmed to be associated with both body weight and total length. From these identified SNPs, we annotated a total of 34 genes, including adgrb2, csnkza1, cers5, col22a1, creb5, dnd1, dzank1, dnai1, npy2r, fat3, lrrk2, lrp5, map3k9, and so on. Among these candidate genes, npy2r (neuropeptide Y receptor Y2) was reported to play a critical role in growth of the orange-spotted grouper. In addition, population structure, principal component analysis, kinship matrix and linkage disequilibrium were examined to verify the accuracy and reliability of our GWAS results. Our data will also provide a valuable genetic resource for further marker-assisted selection program to improve growth quality in groupers.