A genome-wide association study on growth traits in orangespotted grouper(Epinephelus coioides) 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 which110 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, col22 a1,creb5, dnd1, dzank1, dnai1, npy2 r, fat3, lrrk2, lrp5, map3 k9, and so on. Among these candidate genes, npy2 r(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.     

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.

     

Genomewide association study for economic traits in the large yellow croaker with different numbers of extreme phenotypes

A traditional genomewide association study (GWAS) detects genotype–phenotype associations by the vast number of genotyped individuals. This method requires large-scale samples and considerable sequencing costs. Extreme phenotypic sampling proposes make GWAS more cost-efficient and are applied more widely. With extreme phenotypic sampling, we performed a GWAS for n-3 highly unsaturated fatty acids (HUFA) and eviscerated weight (EW) traits in the large yellow croaker population. Of the 32,249 and 29,748 detected SNPs for the two traits, three candidate regions were found in each trait. Three candidate regions associated with HUFA were known near genes on chromosomes 4 and 11, and three candidate regions were on chromosome 6, and 15 for the EW trait. By combing through our GWAS results and the biological functional analysis of the genes, we suggest that the FABP, DGAT, ATP8B1, FAF2 and CERS2 genes,  as well as the IGF2, BORA, CYP1A1, GRTP1 and HOX genes are promising candidate genes for n-3 HUFA and EW, respectively, in the large yellow croaker. Moreover, compared with the different numbers of the extreme phenotypic sampling, we conclude that 60% of the extreme phenotypic subsample can obtain a similar result as GWAS with whole phenotypes. Thus, extreme phenotypic sampling could save 40% of the cost for genotyping and DNA extraction without loss of the candidate regions and functional genes. Our study may provide a basis for further genomic breeding and a reference for others who want to perform GWAS with extreme phenotypes.     

Evaluation of European Schizophrenia GWAS Loci in Asian Populations via Comprehensive Meta-Analyses

Schizophrenia is a severe and highly heritable neuropsychiatric disorder. Recent genetic analyses including genome-wide association studies (GWAS) have implicated multiple genome-wide significant variants for schizophrenia among European populations. However, many of these risk variants were not largely validated in other populations of different ancestry such as Asians. To validate whether these European GWAS significant loci are associated with schizophrenia in Asian populations, we conducted a systematic literature search and meta-analyses on 19 single nucleotide polymorphisms (SNPs) in Asian populations by combining all available case-control and family-based samples, including up to 30,000 individuals. We employed classical fixed (or random) effects inverse variance weighted methods to calculate summary odds ratios (ORs) and 95 % confidence intervals (CIs). Among the 19 GWAS loci, we replicated the risk associations of nine markers (e.g., SNPs at VRK2, ITIH3/4, NDST3, NOTCH4) surpassing significance level (two-tailed P?<?0.05), and three additional SNPs in MIR137 and ZNF804A also showed trend associations (one-tailed P?<?0.05). These risk associations are in the same directions of allelic effects between Asian replication samples and initial European GWAS findings, and the successful replications of these GWAS loci in a different ethnic group provide stronger evidence for their clinical associations with schizophrenia. Further studies, focusing on the molecular mechanisms of these GWAS significant loci, will become increasingly important for understanding of the pathogenesis to schizophrenia.     

SNP markers associated with body size and pelt length in American mink (<Emphasis Type="Italic">Neovison vison</Emphasis>)

Background

Identification of genes underlying production traits is a key aim of the mink research community. Recent availability of genomic tools have opened the possibility for faster genetic progress in mink breeding. Availability of mink genome assembly allows genome-wide association studies in mink.

Results

In this study, we used genotyping-by-sequencing to obtain single nucleotide polymorphism (SNP) genotypes of 2496 mink. After multiple rounds of filtering, we retained 28,336 high quality SNPs and 2352 individuals for a genome-wide association study (GWAS). We performed the first GWAS for body weight, behavior, along with 10 traits related to fur quality in mink.

Conclusions

Combining association results with existing functional information of genes and mammalian phenotype databases, we proposed WWC3, MAP2K4, SLC7A1 and USP22 as candidate genes for body weight and pelt length in mink.

     

Impact of variants on type-2 diabetes risk genes identified through genomewide association studies in polycystic ovary syndrome: a case–control study

Polycystic ovary syndrome (PCOS) is a common endocrine disorder in females, and is associated with altered metabolic processes in particular insulin resistance and diabetes mellitus. PCOS shares with type-2 diabetes (T2D) a number of features, including beta cell dysfunction, impaired glucose tolerance and dyslipidaemia. Recently, genomewide association studies (GWAS) have reported a number of genes with reproducible associations and susceptibilities to T2D. To address this, we examined the association between the T2D GWAS candidate genes (CDKAL1, CDKN2B, COL8A1, HHEX, IGF2BP2, KCNJ1, KCNQ1 and SLC30A8) and PCOS in Saudi women. A case–control study, includes 162 cases and 162 controls was enrolled. Genotyping was carried out by the allelic discrimination method. Our results showed that the variants including rs792837 of COL8A1, rs61873498 of KCNQ1 and rs13266634 of SLC30A8 genes to be significantly more frequent in PCOS patients than in controls. Our results suggest that COL8A1, KCNQ1 and SLC30A8, which are previously identified through GWAS as T2D-associated genes, are associated with PCOS.     

Genome-wide association study of body weight in Wenshang Barred chicken based on the SLAF-seq technology

Chicken body weight (BW) is an economically important trait, and many studies have been conducted on genetic selection for BW. However, previous studies have detected functional chromosome mutations or regions using gene chips. The present study used the specific-locus amplified fragment sequencing (SLAF-seq) technology to perform a genome-wide association study (GWAS) on purebred Wengshang Barred chicken. A total of 1,286,715 single-nucleotide polymorphisms (SNPs) were detected, and 175,211 SNPs were selected as candidate SNPs for genome-wide association analysis using TASSEL general linear models. Six SNP markers reached genome-wide significance. Of these, rs732048524, rs735522839, rs738991545, and rs15837818 were significantly associated with body weight at 28 days (BW28), while rs314086457 and rs315694878 were significantly associated with BW120. These SNPs are close to seven genes (PRSS23, ME3, FAM181B, NABP1, SDPR, TSSK6L2, and RBBP8). Moreover, 24 BW-associated SNPs reached “suggestive” genome-wide significance. Of these, 6, 13, 1, and 4 SNPs were associated with BW28, BW56, BW80, and BW120, respectively. These results would enrich the studies on BW and promote the use of Chinese chicken, especially the Wenshang Barred chicken.     

Discovery and characterization of two new stem rust resistance genes in <Emphasis Type="Italic">Aegilops sharonensis</Emphasis>

Key message

We identified two novel wheat stem rust resistance genes, Sr-1644-1Sh and Sr-1644-5Sh in Aegilops sharonensis that are effective against widely virulent African races of the wheat stem rust pathogen.

Abstract

Stem rust is one of the most important diseases of wheat in the world. When single stem rust resistance (Sr) genes are deployed in wheat, they are often rapidly overcome by the pathogen. To this end, we initiated a search for novel sources of resistance in diverse wheat relatives and identified the wild goatgrass species Aegilops sharonesis (Sharon goatgrass) as a rich reservoir of resistance to wheat stem rust. The objectives of this study were to discover and map novel Sr genes in Ae. sharonensis and to explore the possibility of identifying new Sr genes by genome-wide association study (GWAS). We developed two biparental populations between resistant and susceptible accessions of Ae. sharonensis and performed QTL and linkage analysis. In an F₆ recombinant inbred line and an F₂ population, two genes were identified that mapped to the short arm of chromosome 1S^sh, designated as Sr-1644-1Sh, and the long arm of chromosome 5S^sh, designated as Sr-1644-5Sh. The gene Sr-1644-1Sh confers a high level of resistance to race TTKSK (a member of the Ug99 race group), while the gene Sr-1644-5Sh conditions strong resistance to TRTTF, another widely virulent race found in Yemen. Additionally, GWAS was conducted on 125 diverse Ae. sharonensis accessions for stem rust resistance. The gene Sr-1644-1Sh was detected by GWAS, while Sr-1644-5Sh was not detected, indicating that the effectiveness of GWAS might be affected by marker density, population structure, low allele frequency and other factors.

     

Genome-wide association study for soybean cyst nematode resistance in Chinese elite soybean cultivars

Soybean cyst nematode (SCN) (Heterodera glycines Ichinohe) is a highly recalcitrant endoparasite of soybean roots, causing more yield loss than any other pest. To identify quantitative trait loci (QTL) controlling resistance to SCN (HG type 2.5.7, race 1), a genome-wide association study (GWAS) was performed. The association panel, consisting of 120 Chinese soybean cultivars, was genotyped with 7189 single nucleotide polymorphism (SNPs). A total of 6204 SNPs with minor allele frequency >0.05 were used to estimate linkage disequilibrium (LD) and population structure. The mean level of LD measured by r ² declined very rapidly to half its maximum value (0.51) at 220 kb. The overall population structure was approximately coincident with geographic origin. The GWAS results identified 13 SNPs in 7 different genomic regions significantly associated with SCN resistance. Of these, three SNPs were localized in previously mapped QTL intervals, including rhg1 and Rhg4. The GWAS results also detected 10 SNPs in 5 different genomic regions associated with SCN resistance. The identified loci explained an average of 95.5% of the phenotypic variance. The proportion of phenotypic variance was due to additive genetic variance of the validated SNPs. The present study identified multiple new loci and refined chromosomal regions of known loci associated with SCN resistance. The loci and trait-associated SNPs identified in this study can be used for developing soybean cultivars with durable resistance against SCN.     

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%.     

Genome-wide association study of heading and flowering dates and construction of its prediction equation in Chinese common wheat

Heading date is one of the most important traits in wheat breeding as it affects adaptation and yield potential. A genome-wide association study (GWAS) using the 90 K iSelect SNP genotyping assay indicated that a total of 306 loci were significantly associated with heading and flowering dates in 13 environments in Chinese common wheat from the Yellow and Huai wheat region. Of these, 105 loci were significantly correlated with both heading and flowering dates and were found in clusters on chromosomes 2, 5, 6, and 7. Based on differences in distribution of the vernalization and photoperiod genes among chromosomes, arms, or block regions, 13 novel, environmentally stable genetic loci were associated with heading and flowering dates, including RAC875_c41145_189 on 1DS, RAC875_c50422_299 on 2BL, and RAC875_c48703_148 on 2DS, that accounted for more than 20% phenotypic variance explained (PVE) of the heading/flowering date in at least four environments. GWAS and t test of a combination of SNPs and vernalization and photoperiod alleles indicated that the Vrn-B1, Vrn-D1, and Ppd-D1 genes significantly affect heading and flowering dates in Chinese common wheat. Based on the association of heading and flowering dates with the vernalization and photoperiod alleles at seven loci and three significant SNPs, optimal linear regression equations were established, which show that of the seven loci, the Ppd-D1 gene plays the most important role in modulating heading and flowering dates in Chinese wheat, followed by Vrn-B1 and Vrn-D1. Additionally, three novel genetic loci (RAC875_c41145_189, Excalibur_c60164_137, and RAC875_c50422_299) also show important effect on heading and flowering dates. Therefore, Ppd-D1, Vrn-B1, Vrn-D1, and the novel genetic loci should be further investigated in terms of improving heading and flowering dates in Chinese wheat. Further quantitative analysis of an F₁₀ recombinant inbred lines population identified a major QTL that controls heading and flowering dates within the Ppd-D1 locus with PVEs of 28.4% and 34.0%, respectively; this QTL was also significantly associated with spike length, peduncle length, fertile spikelets number, cold resistance, and tiller number.     

Single nucleotide polymorphisms in <Emphasis Type="Italic">TaER</Emphasis> genes and their association with carbon isotope discrimination in wheat genotypes under drought

Candidate gene association studies implicate the detection of contributing single nucleotide polymorphism (SNP) for the target traits and have been recommended as a promising technique to anatomize the complex characters in plants. The ERECTA gene in plants controls different physiological functions. In this study, we identified SNPs in 1.1 kb partial sequences of TaER-1 and TaER-2 of wheat (Triticum aestivum L.). Thirty-nine SNPs were identified in the coding regions of TaER-1 gene in 33 wheat genotypes, of which 20 SNPs caused non-synonymous mutations while 19 SNPs produced synonymous mutations; 31 SNPs were located in the coding regions of TaER-2 gene in 26 genotypes, of which 18 SNPs caused non-synonymous mutations and 13 SNPs caused synonymous mutations. In addition, 32 SNPs in TaER-1 and 9 SNPs in TaER-2 were also identified in the non-coding regions. Moreover, the significant genetic associations of SNPs of TaER-1 and TaER-2 genes with carbon isotope discrimination, stomatal conductance, photosynthetic rate, transpiration rate, intrinsic water use efficiency (iWUE), leaf length, leaf width, stomatal density, epidermal cell density, and stomatal index were noted in wheat genotypes. This study confirms the importance of TaER-1 and TaER-2 genes which could improve iWUE of wheat by regulating leaf gas exchange and leaf structural traits. These identified SNPs may play a critical role in molecular breeding by means of marker-assisted selection.     

Association between expression levels and growth trait-related SNPs located in promoters of the <Emphasis Type="Italic">MC4R</Emphasis> and <Emphasis Type="Italic">MSTN</Emphasis> genes in <Emphasis Type="Italic">Spinibarbus hollandi</Emphasis>

Melanocortin 4 receptor: (MC4R) and Myostatin (MSTN) are two important growth trait-related genes in animals. In this study, we showed that two SNPs, MC4R-719A>G and MSTN-519C>T, found in the promoters of the MC4R and MSTN genes, respectively, are both associated with growth traits in Spinibarbus hollandi. Furthermore, we observed that there were significant associations between the expression levels of the MC4R and MSTN genes and these two growth trait-related SNPs. The expression level of MC4R gene in brain was lower in GG genotype fish with extremely high growth performance than that in AA genotype fish with extremely low growth performance. Expression level of the MSTN gene in muscle was lower in TT genotype fish with extremely high growth performance than that in CC and CT genotype fish with lower growth performance. The results indicated that these SNPs located in the promoters of MC4R and MSTN are associated with growth-related traits through modification of gene expression levels. The MSTN and MC4R SNPs may have useful application in effective marker-assisted selection aimed to increase output in S. hollandi.     

MicroRNA Biogenesis Pathway Gene Polymorphisms Are Associated with Breast Cancer Risk

MicroRNAs (miRNAs) play an important role as epigenetic regulators in cancer initiation and progression. One of the mechanisms of miRNA dysregulation is altered functioning of proteins involved in miRNA processing machinery. It has been suggested that single nucleotide polymorphisms (SNPs) within miRNA gene regions, miRNA target genes, and miRNA machinery genes may affect the miRNAs regulation. We selected 25 SNPs in the key genes of miRNA biosynthesis, including DROSHA/RNASEN, DGCR8, DICER1, XPO5, RAN, PIWIL1/HIWI, AGO1/EIF2C1, AGO2, GEMIN4, GEMIN3/DDX20, and DDX5, and investigated the association between these SNPs and the risk of breast cancer. The total number of breast cancer cases and cancer-free controls enrolled in the investigation were 778 (417 breast cancer patients and 361 healthy women). We found that rs11060845 and rs10773771 in the PIWIL1 gene, rs3809142/RAN, rs10719/DROSHA, rs1640299/DGCR8, rs563002/DDX20, rs595055/AGO1, and rs2740348/GEMIN4 were associated with breast cancer risk in Russians.     

Evaluation of antifungal,phosphate solubilisation,and siderophore and chitinase release activities of endophytic fungi from <Emphasis Type="Italic">Pistacia vera</Emphasis>

Fungal endophytes use different strategies to protect host plants from abiotic and biotic stress. In this study, we isolated endophytic fungi from Pistacia vera and characterised their antifungal activity against Aspergillus flavus, Rhizoctonia solani and Sclerotinia sclerotiorum, and their release of some factors that can alter plant growth capability. Trichoderma harzianum TH 5-1-2, T. harzianum TH 10-2-2 and T. atroviride TA 2-2-1 exhibited the highest growth inhibition percentages in dual culture assays against A. flavus, R. solani and S. sclerotiorum, respectively. Among the fungal endophyte cultures, ethyl acetate extracts of T. harzianum TH 10-2-2, T. harzianum TH 5-1-2 and T. atroviride TA 2-2-1 exhibited the highest growth inhibition of S. sclerotiorum, R. solani and A. flavus, respectively. Phosphate solubilisation was induced by Byssochlamys nivea BN 1-1-1 in culture. Large amounts of siderophore production were observed with Quambalaria cyanescens QC 11-3-2 and Epicoccum nigrum EN1, but Trichoderma spp. also produced siderophore in lower amounts. Trichoderma harzianum TH 5-1-2 produced the highest chitinase activity (2.92 U/mL). In general, among the endophytes isolated, Trichoderma spp. appear to have the most promise for promoting healthy growth of P. vera.