Three‐dimensional genetic networks among seed oil‐related traits,metabolites and genes reveal the genetic foundations of oil synthesis in soybean

Although the biochemical and genetic basis of lipid metabolism is clear in Arabidopsis, there is limited information concerning the relevant genes in Glycine max (soybean). To address this issue, we constructed three‐dimensional genetic networks using six seed oil‐related traits, 52 lipid metabolism‐related metabolites and 54 294 SNPs in 286 soybean accessions in total. As a result, 284 and 279 candidate genes were found to be significantly associated with seed oil‐related traits and metabolites by phenotypic and metabolic genome‐wide association studies and multi‐omics analyses, respectively. Using minimax concave penalty (MCP) and smoothly clipped absolute deviation (SCAD) analyses, six seed oil‐related traits were found to be significantly related to 31 metabolites. Among the above candidate genes, 36 genes were found to be associated with oil synthesis (27 genes), amino acid synthesis (four genes) and the tricarboxylic acid (TCA) cycle (five genes), and four genes (GmFATB1a, GmPDAT, GmPLDα1 and GmDAGAT1) are already known to be related to oil synthesis. Using this information, 133 three‐dimensional genetic networks were constructed, 24 of which are known, e.g. pyruvate–GmPDAT–GmFATA2–oil content. Using these networks, GmPDAT, GmAGT and GmACP4 reveal the genetic relationships between pyruvate and the three major nutrients, and GmPDAT, GmZF351 and GmPgs1 reveal the genetic relationships between amino acids and seed oil content. In addition, GmCds1, along with average temperature in July and the rainfall from June to September, influence seed oil content across years. This study provides a new approach for the construction of three‐dimensional genetic networks and reveals new information for soybean seed oil improvement and the identification of gene function.     

Genetic dissection of bull fertility in US Jersey dairy cattle

The service sire has been recognized as an important factor affecting herd fertility in dairy cattle. Recent studies suggest that genetic factors explain part of the difference in fertility among Holstein sires. The main objective of this study was to dissect the genetic architecture of sire fertility in US Jersey cattle. The dataset included 1.5 K Jersey bulls with sire conception rate (SCR) records and 96 K single nucleotide polymorphism (SNP) markers spanning the whole genome. The analysis included whole‐genome scans for both additive and non‐additive effects and subsequent functional enrichment analyses using KEGG Pathway, Gene Ontology (GO) and Medical Subject Headings (MeSH) databases. Ten genomic regions located on eight different chromosomes explained more than 0.5% of the additive genetic variance for SCR. These regions harbor genes, such as PKDREJ, EPB41L2, PDGFD, STX2, SLC25A20 and IP6K1, that are directly implicated in testis development and spermatogenesis, sperm motility and the acrosome reaction. In addition, the genomic scan for non‐additive effects identified two regions on BTA11 and BTA25 with marked recessive effects. These regions harbor three genes—FER1L5, CNNM4 and DNAH3—with known roles in sperm biology. Moreover, the gene‐set analysis revealed terms associated with calcium regulation and signaling, membrane fusion, sperm cell energy metabolism, GTPase activity and MAPK signaling. These gene sets are directly implicated in sperm physiology and male fertility. Overall, this integrative genomic study unravels genetic variants and pathways affecting Jersey bull fertility. These findings may contribute to the development of novel genomic strategies for improving sire fertility in Jersey cattle.     

Discovery of selection-driven genetic differences of Duroc,Landrace, and Yorkshire pig breeds by EigenGWAS and Fst analyses

Pigs are one of the earliest domesticated animals and multiple breeds have been developed to meet the various demands of consumers. EigenGWAS is a novel strategy to identify candidate genes that underlying population genetic differences and to infer candidate regions under selection as well. In this study, EigenGWAS and F_st analyses were performed using the public re-sequencing data of three typical commercial pig breeds, Duroc, Landrace and Yorkshire. The intersection of genome-wide significant SNPs detected by EigenGWAS and top-ranked 1% SNPs of F_st results were treated as signals under selection. Using the data of all three breeds, 3062 signals under selection were detected and the nearby genomic regions within 300 kb upstream and downstream covered 6.54% of whole genome. Pairs of breeds were analysed along with the pathway analysis. The gene function enrichment results indicated that many candidate genes located in the genomic regions of the signals under selection were associated with biological processes related to growth, metabolism, reproduction, sensory perception, etc. Among the candidate genes, the FSHB, AHR, PTHLH, KDR and FST genes were reported to be associated with reproductive performance; the KIT, KITLG, MITF, MC1R and EDNRB genes were previously identified to affect coat colour; the RETREG1, TXNIP, BMP5, PPARD and RBP4 genes were reported to be associated with lipid metabolism and growth traits. The identified genetic differences across the three commercial breeds will advance understanding of the artificial selection history of pigs and the signals under selection will suggest potential uses in pig genomic breeding programmes.     

Search for informative polymorphisms in candidate genes: clock genes and circadian behaviour in blue tits

The identification of functional polymorphisms in genes that underlie behavioural trait variation is a challenging but intriguing task in evolutionary biology. Given the wealth of genomic data and the increasing number of genotype–phenotype association studies in model organisms, one can ask whether and how this information can be used for non-model organisms. Here we describe two strategies to search for likely functional polymorphisms in candidate genes in a bird species that has been intensively studied by behavioural and population ecologists, the blue tit Cyanistes caeruleus. In the first approach we searched for repeating elements in coding regions of the genome using information about repeats in Gallus gallus genes. The rationale is that tandem-repeat elements have a high potential to be polymorphic and functional. The second strategy aimed to replicate reported genotype–phenotype association studies by extrapolating results from model organisms to our study species. Both strategies showed high success rates with respect to finding homologous gene regions and potentially informative genetic variants in the genes AANAT, ADCYAP1, CKIε, CLOCK, CREB1, NPAS2 and PERIOD2. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Detecting selection signatures between Duroc and Duroc synthetic pig populations using high‐density SNP chip

The development of high throughput genotyping techniques has facilitated the identification of selection signatures of pigs. The detection of genomic selection signals in a population subjected to differential selection pressures may provide insights into the genes associated with economically and biologically important traits. To identify genomic regions under selection, we genotyped 488 Duroc (D) pigs and 155 D × Korean native pigs (DKNPs) using the Porcine SNP70K BeadChip. By applying the F_ST and extended haplotype homozygosity (EHH‐Rsb) methods, we detected genes under directional selection associated with growth/stature (DOCK7, PLCB4, HS2ST1, FBP2 and TG), carcass and meat quality (TG, COL14A1, FBXO5, NR3C1, SNX7, ARHGAP26 and DPYD), number of teats (LOC100153159 and LRRC1), pigmentation (MME) and ear morphology (SOX5), which are all mostly near or at fixation. These results could be a basis for investigating the underlying mutations associated with observed phenotypic variation. Validation using genome‐wide association analysis would also facilitate the inclusion of some of these markers in genetic evaluation programs.     

A targeted genotyping approach enhances identification of variants in taste receptor and appetite/reward genes of potential functional importance for obesity‐related porcine traits

Taste receptors (TASRs) and appetite and reward (AR) mechanisms influence eating behaviour, which in turn affects food intake and risk of obesity. In a previous study, we used next generation sequencing to identify potentially functional mutations in TASR and AR genes and found indications for genetic associations between identified variants and growth and fat deposition in a subgroup of animals (n = 38) from the UNIK resource pig population. This population was created for studying obesity and obesity‐related diseases. In the present study we validated results from our previous study by investigating genetic associations between 24 selected single nucleotide variants in TASR and AR gene variants and 35 phenotypes describing obesity and metabolism in the entire UNIK population (n = 564). Fifteen variants showed significant association with specific obesity‐related phenotypes after Bonferroni correction. Six of the 15 genes, namely SIM1, FOS, TAS2R4, TAS2R9, MCHR2 and LEPR, showed good correlation between known biological function and associated phenotype. We verified a genetic association between potentially functional variants in TASR/AR genes and growth/obesity and conclude that the combination of identification of potentially functional variants by next generation sequencing followed by targeted genotyping and association studies is a powerful and cost‐effective approach for increasing the power of genetic association studies.     

A genome‐wide association study of heat stress‐associated SNPs in catfish

Heat tolerance is a complex and economically important trait for catfish genetic breeding programs. With global climate change, it is becoming an increasingly important trait. To better understand the molecular basis of heat stress, a genome‐wide association study (GWAS) was carried out using the 250 K catfish SNP array with interspecific backcross progenies, which derived from crossing female channel catfish with male F1 hybrid catfish (female channel catfish × male blue catfish). Three significant associated SNPs were detected by performing an EMMAX approach for GWAS. The SNP located on linkage group 14 explained 12.1% of phenotypical variation. The other two SNPs, located on linkage group 16, explained 11.3 and 11.5% of phenotypical variation respectively. A total of 14 genes with heat stress related functions were detected within the significant associated regions. Among them, five genes—TRAF2, FBXW5, ANAPC2, UBR1 and KLHL29— have known functions in the protein degradation process through the ubiquitination pathway. Other genes related to heat stress include genes involved in protein biosynthesis (PRPF4 and SYNCRIP), protein folding (DNAJC25), molecule and iron transport (SLC25A46 and CLIC5), cytoskeletal reorganization (COL12A1) and energy metabolism (COX7A2, PLCB1 and PLCB4) processes. The results provide fundamental information about genes and pathways that is useful for further investigation into the molecular mechanisms of heat stress. The associated SNPs could be promising candidates for selecting heat‐tolerant catfish lines after validating their effects on larger and various catfish populations.     

Genome‐wide association analysis identifies the genetic basis of fat deposition in the tails of sheep (Ovis aries)

Fat‐tailed sheep (Ovis aries) can survive in harsh environments and satisfy human's intake of dietary fat. However, the animals require more feed, which increases the cost of farming. Thus, most farmers currently prefer thin‐tailed, short‐tailed or docked sheep. To date, the molecular mechanism of the formation of fat tails in sheep has not been completely elucidated. Here, we conducted a genome‐wide association study using phenotypes and genotypes (the Ovine Infinium HD SNP BeadChip genotype data) of two breeds of contrasting tail types (78 Small‐tailed and 78 Large‐tailed Han sheep breeds) to identify functional genes and variants associated with fat deposition. We identified four significantly (rs416433540, rs409848439, rs408118325 and rs402128848) and three approximately associated autosomal SNPs (rs401248376, rs402445895 and rs416201901). Gene annotation indicated that the surrounding genes (CREB1, STEAP4, CTBP1 and RIP140, also known as NRIP1) function in lipid storage or fat cell regulation. Furthermore, through an X‐chromosome‐wide association analysis, we detected significantly associated SNPs in the OARX: 88–89 Mb region, which could be a strong candidate genomic region for fat deposition in tails of sheep. Our results represent a new genomic resource for sheep genetics and breeding. In addition, the findings provide novel insights into genetic mechanisms of fat deposition in the tail of sheep and other mammals.     

Genome‐wide association study reveals candidate genes associated with body measurement traits in Chinese Wagyu beef cattle

We performed a genome‐wide association study to identify candidate genes for body measurement traits in 463 Wagyu beef cattle typed with the Illumina Bovine HD 770K SNP array. At the genome‐wide level, we detected 18, five and one SNPs associated with hip height, body height and body length respectively. In total, these SNPs are within or near 11 genes, six of which (PENK, XKR4, IMPAD1, PLAG1, CCND2 and SNTG1) have been reported previously and five of which (CSMD3, LAP3, SYN3, FAM19A5 and TIMP3) are novel candidate genes that we found to be associated with body measurement traits. Further exploration of these candidate genes will facilitate genetic improvement in Chinese Wagyu beef cattle.     

Genomic features of the fall armyworm (Spodoptera frugiperda) (J.E. Smith) yield insights into its defense system and flight capability

The fall armyworm (FAW), Spodoptera frugiperda (J.E. Smith) is an agriculturally important insect pest and causes significant economic loss due to its resistance to insecticide and high dispersal ability (i.e. long‐distance flyer). We performed comparative genomics analysis based on the recently released chromosome‐scale genome of the FAW to investigate the genetic mechanisms of FAWs defense system and flight ability. The gene families encoding receptors for bitter or toxic substances and detoxification enzymes, such as cytochrome P450, carboxylesterase and glutathione‐S‐transferase, were expanded in FAW, enabling its ability to detect and detoxify many plant secondary compounds. The chromosomal locations of five detoxification‐related genes (i.e. ABC, CCE, GST, P450, and UGT) were analyzed, which revealed the density of these genes across the FAW genome. Several genes like Mio, Sik2, CRTC, and Lpin were positively selected in FAW, which were associated with energy metabolism, lipid metabolism, and/or muscle development, and thus possibly contributed to the high flight ability of FAW. We also analyzed the homologs to antimicrobial genes reported in the Drosophila genome, and identified 2 attacins, 6 cecropins, 3 gloverins, and 2 moricins in FAW. The findings of this study will broadly expand our knowledge on the biology of the devastating FAW, and contribute to the development of new pest management strategies for the control of major agricultural pests like FAW.     

Association of metabolic gene polymorphisms with alcohol consumption in controls

The objectives were to study the association between metabolic genes involved in alcohol metabolism (CYP2E1 RsaI, CYP2E1 DraI, ADH1C, NQO1) and alcohol consumption in a large sample of healthy controls. Healthy subjects were selected from the International Collaborative Study on Genetic Susceptibility to Environmental Carcinogens (GSEC). Subjects with information on both alcohol consumption and at least one of the studied polymorphisms were included in the analysis (n=2224). Information on the amount of alcohol consumption was available for a subset of subjects (n=844). None of the studied genes was significantly associated with drinking habits. A significant heterogeneity with age was observed when studying the association between CYP2E1 RsaI and alcohol drinking. CYP2E1 RsaI polymorphism was significantly associated with being a never drinker at older ages (odds ratio [OR] 2.4, 95% confidence interval [CI] 1.2–4.8; at ages above 68 years), while the association was reversed at ages below 47 years (OR 0.5, 95% CI 0.2–1.4). For subjects with detailed information on alcohol intake, no association between alcohol quantity and polymorphisms in metabolic genes was observed; subjects carrying the NQO1 polymorphism tended to drink more than subjects carrying the wild-type alleles. Therefore, no significant association between CYP2E1 RsaI, CYP2E1 DraI, ADH1C, NQO1 polymorphisms and alcohol consumption was observed in healthy controls.     

The Role of Obesity‐associated Loci Identified in Genome‐wide Association Studies in the Determination of Pediatric BMI

The prevalence of obesity in children and adults in the United States has increased dramatically over the past decade. Besides environmental factors, genetic factors are known to play an important role in the pathogenesis of obesity. A number of genetic determinants of adult BMI have already been established through genome‐wide association (GWA) studies. In this study, we examined 25 single‐nucleotide polymorphisms (SNPs) corresponding to 13 previously reported genomic loci in 6,078 children with measures of BMI. Fifteen of these SNPs yielded at least nominally significant association to BMI, representing nine different loci including INSIG2, FTO, MC4R, TMEM18, GNPDA2, NEGR1, BDNF, KCTD15, and 1q25. Other loci revealed no evidence for association, namely at MTCH2, SH2B1, 12q13, and 3q27. For the 15 associated variants, the genotype score explained 1.12% of the total variation for BMI z‐score. We conclude that among 13 loci that have been reported to associate with adult BMI, at least nine also contribute to the determination of BMI in childhood as demonstrated by their associations in our pediatric cohort.     

Dissection of genotype-phenotype associations in rice grains using metabolome quantitative trait loci analysis

A comprehensive and large‐scale metabolome quantitative trait loci (mQTL) analysis was performed to investigate the genetic backgrounds associated with metabolic phenotypes in rice grains. The metabolome dataset consisted of 759 metabolite signals obtained from the grains of 85 lines of rice (Oryza sativa, Sasanishiki × Habataki back‐crossed inbred lines). Metabolome analysis was performed using four mass spectrometry pipelines to enhance detection of different classes of metabolites. This mQTL analysis of a wide range of metabolites highlighted an uneven distribution of 802 mQTLs on the rice genome, as well as different modes of metabolic trait (m‐trait) control among various types of metabolites. The levels of most metabolites within rice grains were highly sensitive to environmental factors, but only weakly associated with mQTLs. Coordinated control was observed for several groups of metabolites, such as amino acids linked to the mQTL hotspot on chromosome 3. For flavonoids, m‐trait variation among the experimental lines was tightly governed by genetic factors that alter the glycosylation of flavones. Many loci affecting levels of metabolites were detected by QTL analysis, and plausible gene candidates were evaluated by in silico analysis. Several mQTLs profoundly influenced metabolite levels, providing insight into the control of rice metabolism. The genomic region and genes potentially responsible for the biosynthesis of apigenin‐6,8‐di‐C‐α‐l‐ arabinoside are presented as an example of a critical mQTL identified by the analysis.     

Genetics of fat tissue accumulation in pigs: a comparative approach

Fatness traits are important in pig production since they influence meat quality and fattening efficiency. On the other hand, excessive fat accumulation in humans has become a serious health problem due to worldwide spread of obesity. Since the pig is also considered as an animal model for numerous human diseases, including obesity and metabolic syndrome, comparative genomic studies may bring new insights into genetics of fatness/obesity. Input of genetic factors into phenotypic variability of these traits is rather high and the heritability coefficient (h ²) of these traits oscillates around 0.5. Genome scanning revealed the presence of more than 500 QTLs for fatness in the pig genome. In addition to QTL studies, many candidate gene polymorphisms have been analyzed in terms of their associations with pig fatness, including genes encoding leptin (LEP) and its receptor (LEPR), insulin-like growth factor 2 (IGF-2), fatty acid-binding proteins (FABP3 andFABP4), melanocortin receptor type 4 (MC4R), and theFTO (fat mass and obesity-associated) gene. Among them, a confirmed effect on pig fatness was found for a well-known polymorphism of theIGF-2 gene. In humans the strongest association with predisposition to obesity was shown for polymorphism of theFTO gene, while in pigs such an association seems to be doubtful. The development of functional genomics has revealed a large number of genes whose expression is associated with fat accumulation and lipid metabolism, so far not studied extensively in terms of the association of their polymorphism with pig fatness. Recently, epigenomic mechanisms, mainly RNA interference, have been considered as a potential source of information on genetic input into the fat accumulation process. The rather limited progress in studies focused on the identification of gene polymorphism related with fatness traits shows that their genetic background is highly complex.     

A genomic map of climate adaptation in Mediterranean cattle breeds

Domestic species such as cattle (Bos taurus taurus and B. t. indicus) represent attractive biological models to characterize the genetic basis of short‐term evolutionary response to climate pressure induced by their post‐domestication history. Here, using newly generated dense SNP genotyping data, we assessed the structuring of genetic diversity of 21 autochtonous cattle breeds from the whole Mediterranean basin and performed genome‐wide association analyses with covariables discriminating the different Mediterranean climate subtypes. This provided insights into both the demographic and adaptive histories of Mediterranean cattle. In particular, a detailed functional annotation of genes surrounding variants associated with climate variations highlighted several biological functions involved in Mediterranean climate adaptation such as thermotolerance, UV protection, pathogen resistance or metabolism with strong candidate genes identified (e.g., NDUFB3, FBN1, METTL3, LEF1, ANTXR2 and TCF7). Accordingly, our results suggest that main selective pressures affecting cattle in Mediterranean area may have been related to variation in heat and UV exposure, in food resources availability and in exposure to pathogens, such as anthrax bacteria (Bacillus anthracis). Furthermore, the observed contribution of the three main bovine ancestries (indicine, European and African taurine) in these different populations suggested that adaptation to local climate conditions may have either relied on standing genomic variation of taurine origin, or adaptive introgression from indicine origin, depending on the local breed origins. Taken together, our results highlight the genetic uniqueness of local Mediterranean cattle breeds and strongly support conservation of these populations.     

Metabolism of halophilic archaea

In spite of their common hypersaline environment, halophilic archaea are surprisingly different in their nutritional demands and metabolic pathways. The metabolic diversity of halophilic archaea was investigated at the genomic level through systematic metabolic reconstruction and comparative analysis of four completely sequenced species: Halobacterium salinarum, Haloarcula marismortui, Haloquadratum walsbyi, and the haloalkaliphile Natronomonas pharaonis. The comparative study reveals different sets of enzyme genes amongst halophilic archaea, e.g. in glycerol degradation, pentose metabolism, and folate synthesis. The carefully assessed metabolic data represent a reliable resource for future system biology approaches as it also links to current experimental data on (halo)archaea from the literature. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Screening for signatures of selection of Tianzhu white yak using genome‐wide re‐sequencing

The Tianzhu white yak, a domestic yak indigenous to the Qilian Mountains, migrated inland from the Qinghai‐Tibet Plateau. Specific ecological and long‐term artificial selection influenced the evolution of its pure white coat and physiological characteristics. Therefore, it is not only a natural population that represents a genomic selective region of environmental adaptability but is also an animal model for studying the pigmentation of the yak coat. A total of 24 261 829 variants, including 22 445 252 SNPs, were obtained from 29 yaks by genome‐wide re‐sequencing. According to the results of a selective sweep analysis of Tianzhu white yak in comparison to Tibetan yaks, nine candidate genes under selection in Tianzhu white yak were identified by combining π, Tajima's D, πA/πB and F_ST statistics, with threshold standards of 5%. These genes include PDCD1, NUP210, ABCG8, NEU4, LOC102287650, D2HGDH, COL4A1, RTP5 and HDAC11. Five of the nine genes were classified into 12 molecular signaling pathways, and most of these signaling pathways are involved in environmental information processing, organismal systems and metabolism. A majority of these genes has not been implicated in previous studies of yak coat color and high‐altitude animals. Our findings are helpful not only for explaining the molecular mechanism of yak coat pigmentation but also for exploring the genetic changes in Tianzhu white yak due to environmental adaptation.