Integrative Transcriptome,Genome and Quantitative Trait Loci Resources Identify Single Nucleotide Polymorphisms in Candidate Genes for Growth Traits in Turbot

Growth traits represent a main goal in aquaculture breeding programs and may be related to adaptive variation in wild fisheries. Integrating quantitative trait loci (QTL) mapping and next generation sequencing can greatly help to identify variation in candidate genes, which can result in marker-assisted selection and better genetic structure information. Turbot is a commercially important flatfish in Europe and China, with available genomic information on QTLs and genome mapping. Muscle and liver RNA-seq from 18 individuals was carried out to obtain gene sequences and markers functionally related to growth, resulting in a total of 20,447 genes and 85,344 single nucleotide polymorphisms (SNPs). Many growth-related genes and SNPs were identified and placed in the turbot genome and genetic map to explore their co-localization with growth-QTL markers. Forty-five SNPs on growth-related genes were selected based on QTL co-localization and relevant function for growth traits. Forty-three SNPs were technically feasible and validated in a wild Atlantic population, where 91% were polymorphic. The integration of functional and structural genomic resources in turbot provides a practical approach for QTL mining in this species. Validated SNPs represent a useful set of growth-related gene markers for future association, functional and population studies in this flatfish species.     

Novel Single Nucleotide Polymorphisms of the Insulin-Like Growth Factor-I Gene and Their Associations with Growth Traits in Common Carp (Cyprinus carpio L.)

Insulin-like growth factor-I (IGF-I) plays an important role in the growth and development of vertebrates. To study polymorphisms of IGF-I, we screened a total of 4555 bp of genomic sequences in four exons and partial introns for the discovery of single nucleotide polymorphism (SNP) in common carp (Cyprinus carpio). Three SNPs (g.3759T>G, g.7627T>A and g.7722T>C) in intron 2 and a nonsynonymous SNP (g.7892C>T) in exon 3 were identified in a pilot population including random parents and their progenies. 289 progenies were further genotyped for studying possible associations between genotypes or combined genotypes and growth traits. The results showed that the locus g.7627T>A was significantly associated with body weight and body length, and fish with genotype AA had a mean body weight 5.9% higher than those with genotype TT. No significant associations were observed between genotypes of other loci and growth traits. However, when both g.7627T>A and g.7722T>C were considered, the combined genotype TT/TT was extremely associated with the lowest values of body length and body weight and the highest K value in comparison with other diplotypes (p < 0.01). These results suggest that genotype AA at g.7627T>A and its combined genotypes with alleles from another locus have positive effects on growth traits, which would be a candidate molecular marker for further studies in marker-assisted selection in common carp.     

Lack of Association between the TSPAN18 Gene and Schizophrenia Based on New Data from Han Chinese and a Meta-Analysis

Tetraspanin-18 (TSPAN18) potentially plays a role in the calcium signaling that is associated with dopamine-induced cortical neuron apoptosis and is considered to be an important mechanism in the pathogenesis of schizophrenia (SCZ). Furthermore, a genome-wide association study (GWAS) identified TSPAN18 as a possible susceptibility gene for SCZ. To validate these findings and reveal the effects of different inheritance models, seven single nucleotide polymorphisms (SNPs) of the TSPAN18 gene were analyzed in 443 patients with SCZ and 628 controls of Han Chinese descent via the SNPscan method. Single SNP, genotype, and association analyses with different models (i.e., additive, dominant, and recessive models) were performed, and the published datasets (2062 cases and 2053 controls) were combined with our results to determine the inheritance effects of the SNPs on SCZ. We observed genotypes and allele distributions of TSPAN18 gene did not show any significant associations in the Han Chinese population based on our experimental and meta-analytical results. Our findings indicate that the TSPAN18 gene is unlikely to be a major susceptibility gene for schizophrenia in Han Chinese.     

Single Nucleotide Polymorphisms in Growth Hormone Gene and Their Association with Growth Traits in Siniperca chuatsi (Basilewsky)

Growth hormone (GH) has been considered as a candidate gene for growth traits in fish. In this study, polymorphisms of the GH gene were evaluated for associations with growth traits in 282 Siniperca chuatsi individuals. Using directly sequencing, four single nucleotide polymorphisms (SNPs) were identified in GH gene, with two mutations in intron 4 (g.4940A>C, g.4948A>T), one mutation in exon 5 (g.5045T>C) and one in intron 5 (g.5234T>G). Notably, three of them were significantly associated with growth performance, particularly for g.4940A>C which was highly correlated with all the four growth traits. In conclusion, our results demonstrated that these SNPs in GH gene could influence growth performance of S.chuatsi and could be used for marker-assisted selection (MAS) in this species.     

Polymorphisms of Leptin-b Gene Associated with Growth Traits in Orange-Spotted Grouper (Epinephelus coioides)

In mammals, leptin has been demonstrated to perform important roles in many physiological activities and to influence development, growth, metabolism and reproduction. However, in fish, its function is still unclear. Duplicate leptin genes, leptin-a and leptin-b, have been identified in the orange-spotted grouper. In the present study, the polymorphisms in the leptin-b gene of the orange-spotted grouper were detected, and the relation between these polymorphisms and 12 growth traits were analyzed. Six polymorphisms (including 3 single nucleotide polymorphisms (c.14G>A, c.93A>G, c.149G>A) in exon 1, 2 SNPs (c.181A>G, c.193G>A) in intron 1, and 1 SNP (c.360C>T) in exon 2) were identified and genotyped from 200 different individuals. The results revealed that the SNP c.149G>A was significantly associated with growth traits, that the heterozygous mutation genotype GA having negative effects on growth traits. However, the other five SNPs (c.14G>A, c.93A>G, c.181A>G, c.193G>A, c.360C>T) did not show significant associations with all the growth traits. Compared with our findings in leptin-a gene, the results suggested that the leptin-a hormone has more important physiological effects in fish bodies than the leptin-b type. Moreover, leptin genes were supposed to be one class of major candidate genes of regulating growth traits in the orange-spotted grouper.     

Association of Allelic Variation in PtoXET16A with Growth and Wood Properties in Populus tomentosa

Xyloglucan endo-transglycosylases (XETs) modify the xyloglucan-cellulose framework of plant cell walls and, thus, affect cell wall expansion and strength. Dissecting the mechanism by which natural variation in XETs affects wood properties can inform breeding efforts to improve wood quality and yield traits. To this end, we isolated a full-length PtoXET16A cDNA clone from Populus tomentosa. Real-time PCR analysis showed that PtoXET16A was maximally expressed in the root, followed by phloem, cambium, and developing xylem, suggesting that PtoXET16A plays important roles in the development of vascular tissues. Nucleotide diversity and linkage disequilibrium analysis revealed that PtoXET16A has high single nucleotide polymorphism (SNP) diversity (π = 0.01266 and θ_w = 0.01392) and low linkage disequilibrium (r² ≥ 0.1, within 900 bp). SNP- and haplotype-based association analyses of 426 individuals from a natural population indicated that nine SNPs (including two non-synonymous markers and one splicing variant) (p ≤ 0.05, false discovery rate Q ≤ 0.01), and nine haplotypes (p ≤ 0.05) were significantly associated with growth and wood properties, each explaining from 3.40%–10.95% of phenotypic variance. This work shows that examination of allelic variation and linkage disequilibrium by a candidate-gene-based approach can help to decipher the genetic basis of wood formation. Moreover, the SNP markers identified in this study can potentially be applied for marker-assisted selection to improve growth and wood-property traits in Populus.     

Genome-Wide Association Study in Bread Wheat Identifies Genomic Regions Associated with Grain Yield and Quality under Contrasting Water Availability

The objectives of this study were to identify genetic loci in the bread wheat genome that would influence yield stability and quality under water stress, and to identify accessions that can be recommended for cultivation in dry and hot regions. We performed a genome-wide association study (GWAS) using a panel of 232 wheat accessions spanning diverse ecogeographic regions. Plants were evaluated in the Israeli Northern Negev, under two environments: water-limited (D; 250 mm) and well-watered (W; 450 mm) conditions; they were genotyped with ~71,500 SNPs derived from exome capture sequencing. Of the 14 phenotypic traits evaluated, 12 had significantly lower values under D compared to W conditions, while the values for two traits were higher under D. High heritability (H² = 0.5–0.9) was observed for grain yield, spike weight, number of grains per spike, peduncle length, and plant height. Days to heading and grain yield could be partitioned based on accession origins. GWAS identified 154 marker-trait associations (MTAs) for yield and quality-related traits, 82 under D and 72 under W, and identified potential candidate genes. We identified 24 accessions showing high and/or stable yields under D conditions that can be recommended for cultivation in regions under the threat of global climate change.     

Deciphering the Genetic Basis of Root and Biomass Traits in Rapeseed (Brassica napus L.) through the Integration of GWAS and RNA-Seq under Nitrogen Stress

An excellent root system is responsible for crops with high nitrogen-use efficiency (NUE). The current study evaluated the natural variations in 13 root- and biomass-related traits under a low nitrogen (LN) treatment in a rapeseed association panel. The studied traits exhibited significant phenotypic differences with heritabilities ranging from 0.53 to 0.66, and most of the traits showed significant correlations with each other. The genome-wide association study (GWAS) found 51 significant and 30 suggestive trait–SNP associations that integrated into 14 valid quantitative trait loci (QTL) clusters and explained 5.7–21.2% phenotypic variance. In addition, RNA sequencing was performed at two time points to examine the differential expression of genes (DEGs) between high and low NUE lines. In total, 245, 540, and 399 DEGs were identified as LN stress-specific, high nitrogen (HN) condition-specific, and HNLN common DEGs, respectively. An integrated analysis of GWAS, weighted gene co-expression network, and DEGs revealed 16 genes involved in rapeseed root development under LN stress. Previous studies have reported that the homologs of seven out of sixteen potential genes control root growth and NUE. These findings revealed the genetic basis underlying nitrogen stress and provided worthwhile SNPs/genes information for the genetic improvement of NUE in rapeseed.     

Polymorphisms in Myostatin Gene and Associations with Growth Traits in the Common Carp (Cyprinus carpio L.)

Myostatin (MSTN) is a member of the transforming growth factor-β superfamily that negatively regulates skeletal muscle development and growth. In the present study, partial genomic fragments of MSTN were screened for single nucleotide polymorphisms (SNPs) in selected common carp individuals from wild populations, and two SNPs in intron 2 (c.371 + 749A > G, c.371 + 781T > C) and two synonymous SNPs in exon 3 (c.42A > G, c.72C > T) were identified. Genotyping by direct sequencing of polymerase chain reaction (PCR) products for these four SNPs were performed in 162 individuals from a commercial hatchery population. Association analysis showed that two SNPs in exon 3 were significantly associated with body weight (BW) and condition factor (K), and haplotype analyses revealed that haplotype H7H8 showed better growth performance. Our results demonstrated that some of the SNPs in MSTN may have positive effects on growth traits and suggested that MSTN could be a candidate gene for growth and marker-assisted selection in common carp.     

Identification of Loci Governing Agronomic Traits and Mutation Hotspots via a GBS-Based Genome-Wide Association Study in a Soybean Mutant Diversity Pool

In this study, we performed a genotyping-by-sequencing analysis and a genome-wide association study of a soybean mutant diversity pool previously constructed by gamma irradiation. A GWAS was conducted to detect significant associations between 37,249 SNPs, 11 agronomic traits, and 6 phytochemical traits. In the merged data set, 66 SNPs on 13 chromosomes were highly associated (FDR p < 0.05) with the following 4 agronomic traits: days of flowering (33 SNPs), flower color (16 SNPs), node number (6 SNPs), and seed coat color (11 SNPs). These results are consistent with the findings of earlier studies on other genetic features (e.g., natural accessions and recombinant inbred lines). Therefore, our observations suggest that the genomic changes in the mutants generated by gamma irradiation occurred at the same loci as the mutations in the natural soybean population. These findings are indicative of the existence of mutation hotspots, or the acceleration of genome evolution in response to high doses of radiation. Moreover, this study demonstrated that the integration of GBS and GWAS to investigate a mutant population derived from gamma irradiation is suitable for dissecting the molecular basis of complex traits in soybeans.     

Identification of Genetic Variations and Candidate Genes Responsible for Stalk Sugar Content and Agronomic Traits in Fresh Corn via GWAS across Multiple Environments

The stem and leaves of fresh corn plants can be used as green silage or can be converted to biofuels, and the stalk sugar content and yield directly determine the application value of fresh corn. To identify the genetic variations and candidate genes responsible for the related traits in fresh corn, the genome-wide scan and genome-wide association analysis (GWAS) were performed. A total of 32 selective regions containing 172 genes were detected between sweet and waxy corns. Using the stalk sugar content and seven other agronomic traits measured in four seasons over two years, the GWAS identified ninety-two significant single nucleotide polymorphisms (SNPs). Most importantly, seven SNPs associated with the stalk sugar content were detected across multiple environments, which could explain 13.68–17.82% of the phenotypic variation. Accessions differing in genotype for certain significant SNPs showed significant variation in the stalk sugar content and other agronomic traits, and the expression levels of six important candidate genes were significantly different between two materials with different stalk sugar content. The genetic variations and candidate genes provide valuable resources for future studies of the molecular mechanism of the stalk sugar content and establish the foundation for molecular marker-assisted breeding of fresh corn.     

Associations of MTHFR C677T and MTRR A66G Gene Polymorphisms with Metabolic Syndrome: A Case-Control Study in Northern China

Prior evidence indicates that homocysteine plays a role in the development of metabolic syndrome (MetS). Methylenetetrahydrofolate reductase (MTHFR) C677T and methionine synthase reductase (MTRR) A66G polymorphisms are common genetic determinants of homocysteine levels. To investigate the associations of the MTHFR C677T and MTRR A66G polymorphisms with MetS, 692 Chinese Han subjects with MetS and 878 controls were recruited. The component traits of MetS and the MTHFR C677T and MTRR A66G genotypes were determined. A significant association was observed between the MTHFR 677T allele and increased risk of MetS, high fasting blood glucose, high waist circumference, and increasing number of MetS components. The MTRR A66G polymorphism was associated with an increased risk of MetS when combined with the MTHFR 677TT genotype, although there was no association found between MetS and MTRR A66G alone. Furthermore, the MTRR 66GG genotype was associated with high fasting blood glucose and triglycerides. Our data suggest that the MTHFR 677T allele may contribute to an increased risk of MetS in the northern Chinese Han population. The MTRR A66G polymorphism is not associated with MetS. However, it may exacerbate the effect of the MTHFR C677T variant alone. Further large prospective population-based studies are required to confirm our findings.     

Risk-Association of Five SNPs in TOX3/LOC643714 with Breast Cancer in Southern China

The specific mechanism by which low-risk genetic variants confer breast cancer risk is currently unclear, with contradictory evidence on the role of single nucleotide polymorphisms (SNPs) in TOX3/LOC643714 as a breast cancer susceptibility locus. Investigations of this locus using a Chinese population may indicate whether the findings initially identified in a European population are generalizable to other populations, and may provide new insight into the role of genetic variants in the etiology of breast cancer. In this case-control study, 623 Chinese female breast cancer patients and 620 cancer-free controls were recruited to investigate the role of five SNPs in TOX3/LOC643714 (rs8051542, rs12443621, rs3803662, rs4784227, and rs3112612); Linkage disequilibrium (LD) pattern analysis was performed. Additionally, we evaluated how these common SNPs influence the risk of specific types of breast cancer, as defined by estrogen receptor (ER) status, progesterone receptor (PR) status and human epidermal growth factor receptor 2 (HER2) status. Significant associations with breast cancer risk were observed for rs4784227 and rs8051542 with odds ratios (OR) of 1.31 ((95% confidence intervals (CI), 1.10–1.57)) and 1.26 (95% CI, 1.02–1.56), respectively, per T allele. The T-rs8051542 allele was significantly associated with ER-positive and HER2-negative carriers. No significant association existed between rs12443621, rs3803662, and rs3112612 polymorphisms and risk of breast cancer. Our results support the hypothesis that the applicability of a common susceptibility locus must be confirmed among genetically different populations, which may together explain an appreciable fraction of the genetic etiology of breast cancer.     

The Integration of Genome-Wide Association Study and Homology Analysis to Explore the Genomic Regions and Candidate Genes for Panicle-Related Traits in Foxtail Millet

Panicle traits are important factors affecting yield, and their improvement has long been a critical goal in foxtail millet breeding. In order to understand the genetic basis of panicle formation, a large-scale genome-wide association study (GWAS) was performed in this study for six panicle-related traits based on 706,646 high-polymorphism SNP loci in 407 accessions. As a result, 87 quantitative trait loci (QTL) regions with a physical distance of less than 100 kb were detected to be associated with these traits in three environments. Among them, 27 core regions were stably detected in at least two environments. Based on rice–foxtail millet homologous comparison, expression, and haplotype analysis, 27 high-confidence candidate genes in the QTL regions, such as Si3g11200 (OsDER1), Si1g27910 (OsMADS6), Si7g27560 (GS5), etc., affected panicle-related traits by involving multiple plant growth regulator pathways, a photoperiod response, as well as panicle and grain development. Most of these genes showed multiple effects on different panicle-related traits, such as Si3g11200 affecting all six traits. In summary, this study clarified a strategy based on the integration of GWAS, a homologous comparison, and haplotype analysis to discover the genomic regions and candidate genes for important traits in foxtail millet. The detected QTL regions and candidate genes could be further used for gene clone and marker-assisted selection in foxtail millet breeding.     

Genetic Predisposition to the Mortality in Septic Shock Patients: From GWAS to the Identification of a Regulatory Variant Modulating the Activity of a CISH Enhancer

The high mortality rate in septic shock patients is likely due to environmental and genetic factors, which influence the host response to infection. Two genome-wide association studies (GWAS) on 832 septic shock patients were performed. We used integrative bioinformatic approaches to annotate and prioritize the sepsis-associated single nucleotide polymorphisms (SNPs). An association of 139 SNPs with death based on a false discovery rate of 5% was detected. The most significant SNPs were within the CISH gene involved in cytokine regulation. Among the 139 SNPs associated with death and the 1311 SNPs in strong linkage disequilibrium with them, we investigated 1439 SNPs within non-coding regions to identify regulatory variants. The highest integrative weighted score (IW-score) was obtained for rs143356980, indicating that this SNP is a robust regulatory candidate. The rs143356980 region is located in a non-coding region close to the CISH gene. A CRISPR-Cas9-mediated deletion of this region and specific luciferase assays in K562 cells showed that rs143356980 modulates the enhancer activity in K562 cells. These analyses allowed us to identify several genes associated with death in patients with septic shock. They suggest that genetic variations in key genes, such as CISH, perturb relevant pathways, increasing the risk of death in sepsis patients.     

APOE Polymorphisms Contribute to Reduced Atorvastatin Response in Chilean Amerindian Subjects

Genetic factors can determine the high variability observed in response to lipid-lowering therapy with statins. Nonetheless, the frequency of single nucleotide polymorphisms (SNPs) and their impact can vary due to ethnicity. Because the Chilean population carries a strong Amerindian background, the objective of this study was to evaluate the influence of apolipoprotein E (APOE) variants (rs429358, rs7412) and the 1959C>T SNP (rs5925) in the low-density lipoprotein receptor (LDLR) in response to atorvastatin treatment in hypercholesterolemic individuals. A hundred and thirty nine subjects undergoing statin therapy were included. Identification of Amerindian mtDNA haplogroups was determined by polymerase chain reaction (PCR) and PCR followed by restriction fragment length polymorphism (RFLP), respectively. SNPs were determined by PCR-RFLP. Out of the 139 individuals studied, 84.4% had an Amerindian background, according to mtDNA analysis. In relation to APOE variants, carriers of the E3/4 genotype presented lower cholesterol reduction compared to genotype E3/3 (LDL-C: −18% vs. −29%, p ˂ 0.001). On the other hand, the LDLR rs5925 SNP was not related to atorvastatin response (p = 0.5760). Our results suggest that APOE SNPs are potential predictors to atorvastatin therapy in Amerindian Chilean subjects.     

Cis-eQTL Analysis and Functional Validation of Candidate Genes for Carcass Yield Traits in Beef Cattle

Carcass yield traits are of considerable economic importance for farm animals, which act as a major contributor to the world’s food supply. Genome-wide association studies (GWASs) have identified many genetic variants associated with carcass yield traits in beef cattle. However, their functions are not effectively illustrated. In this study, we performed an integrative analysis of gene-based GWAS with expression quantitative trait locus (eQTL) analysis to detect candidate genes for carcass yield traits and validate their effects on bovine skeletal muscle satellite cells (BSCs). The gene-based GWAS and cis-eQTL analysis revealed 1780 GWAS and 1538 cis-expression genes. Among them, we identified 153 shared genes that may play important roles in carcass yield traits. Notably, the identified cis-eQTLs of PON3 and PRIM2 were significantly (p < 0.001) enriched in previous GWAS loci for carcass traits. Furthermore, overexpression of PON3 and PRIM2 promoted the BSCs’ proliferation, increased the expression of MYOD and downregulated the expression of MYOG, which indicated that these genes may inhibit myogenic differentiation. In contrast, PON3 and PRIM2 were significantly downregulated during the differentiation of BSCs. These findings suggested that PON3 and PRIM2 may promote the proliferation of BSCs and inhibit them in the pre-differentiation stage. Our results further contribute to the understanding of the molecular mechanisms of carcass yield traits in beef cattle.     

LILRA3 Is Associated with Benign Prostatic Hyperplasia Risk in a Chinese Population

A recent prostate cancer (PCa) genome-wide association study (GWAS) identified rs103294, a single nucleotide polymorphism (SNP) located on LILRA3, a key component in the regulation of inflammatory inhibition, to be significantly associated with PCa risk in a Chinese population. Because inflammation may be a common etiological risk factor between PCa and benign prostatic hyperplasia (BPH), the current study was conducted to investigate the association of rs103294 with BPH risk. rs103294 was genotyped in a Chinese population of 426 BPH cases and 1,008 controls from Xinhua Hospital in Shanghai, China. Association between rs103294, BPH risk and clinicopathological traits were tested with adjustment for age. rs103294 was significantly associated with BPH risk with a p-value of 0.0067. Individuals with risk allele “C” had increased risk for BPH (OR = 1.34, 95% CI: 1.09–1.66). Stratified analysis revealed a stronger association risk for younger patients who are below 72 years old (OR = 1.51, 95% CI: 1.06–2.16). Our study represents the first effort to demonstrate that LILRA3 gene is significantly associated with BPH risk in a Chinese population. Our results support a common role of inflammation in the development of PCa and BPH. Additional studies are needed to further evaluate our results.     

A Simple Strategy for Development of Single Nucleotide Polymorphisms from Non-Model Species and Its Application in Panax

Single nucleotide polymorphisms (SNPs) are widely employed in the studies of population genetics, molecular breeding and conservation genetics. In this study, we explored a simple route to develop SNPs from non-model species based on screening the library of single copy nuclear genes (SCNGs). Through application of this strategy in Panax, we identified 160 and 171 SNPs from P. quinquefolium and P. ginseng, respectively. Our results demonstrated that both P. ginseng and P. quinquefolium possessed a high level of nucleotide diversity. The number of haplotype per locus ranged from 1 to 12 for P. ginseng and from 1 to 9 for P. quinquefolium, respectively. The nucleotide diversity of total sites (π_T) varied between 0.000 and 0.023 for P. ginseng and 0.000 and 0.035 for P. quinquefolium, respectively. These findings suggested that this approach is well suited for SNP discovery in non-model organisms and is easily employed in standard genetics laboratory studies.