Significant Single Nucleotide Polymorphisms Research Articles

Weighted single-step genome-wide association study to reveal new candidate genes for productive traits of Landrace pig in Korea.

The objective of this study was to identify genomic regions and candidate genes associated with productive traits using a total of 37,099 productive records and 6,683 single nucleotide polymorphism (SNP) data obtained from five Great-Grand-Parents (GGP) farms in Landrace. The estimated of heritabilities for days to 105 kg (AGE), average daily gain (ADG), backfat thickness (BF), and eye muscle area (EMA) were 0.49, 0.49, 0.56, and 0.23, respectively. We identified a genetic window that explained 2.05%-2.34% for each trait of the total genetic variance. We observed a clear partitioning of the four traits into two groups, and the most significant genomic region for AGE and ADG were located on the Sus scrofa chromosome (SSC) 1, while BF and EMA were located on SSC 2. We conducted Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), which revealed results in three biological processes, four cellular component, three molecular function, and six KEGG pathway. Significant SNPs can be used as markers for quantitative trait loci (QTL) investigation and genomic selection (GS) for productive traits in Landrace pig.

Combined Genome-Wide Association Study and Transcriptome Analysis Reveal Candidate Genes for Resistance to Rust (Puccinia graminis) in Dactylis glomerata.

Rust disease is a common plant disease that can cause wilting, slow growth of plant leaves, and even affect the growth and development of plants. Orchardgrass (Dactylis glomerata L.) is native to temperate regions of Europe, which has been introduced as a superior forage grass in temperate regions worldwide. Orchardgrass has rich genetic diversity and is widely distributed in the world, which may contain rust resistance genes not found in other crops. Therefore, we collected a total of 333 orchardgrass accessions from different regions around the world. Through a genome-wide association study (GWAS) analysis conducted in four different environments, 91 genes that overlap or are adjacent to significant single nucleotide polymorphisms (SNPs) were identified as potential rust disease resistance genes. Combining transcriptome data from susceptible (PI292589) and resistant (PI251814) accessions, the GWAS candidate gene DG5C04160.1 encoding glutathione S-transferase (GST) was found to be important for orchardgrass rust (Puccinia graminis) resistance. Interestingly, by comparing the number of GST gene family members in seven species, it was found that orchardgrass has the most GST gene family members, containing 119 GST genes. Among them, 23 GST genes showed significant differential expression after inoculation with the rust pathogen in resistant and susceptible accessions; 82% of the genes still showed significantly increased expression 14 days after inoculation in resistant accessions, while the expression level significantly decreased in susceptible accessions. These results indicate that GST genes play an important role in orchardgrass resistance to rust (P. graminis) stress by encoding GST to reduce its oxidative stress response.

Nucleotide excision repair gene polymorphisms and hepatoblastoma susceptibility in Eastern Chinese children: A five-center case-control study.

Nucleotide excision repair (NER) plays a vital role in maintaining genome stability, and the effect of NER gene polymorphisms on hepatoblastoma susceptibility is still under investigation. This study aimed to evaluate the relationship between NER gene polymorphisms and the risk of hepatoblastoma in Eastern Chinese Han children. In this five-center case-control study, we enrolled 966 subjects from East China (193 hepatoblastoma patients and 773 healthy controls). The TaqMan method was used to genotype 19 single nucleotide polymorphisms (SNPs) in NER pathway genes, including ERCC1, XPA, XPC, XPD, XPF, and XPG. Then, multivariate logistic regression analysis was performed, and odds ratios (ORs) and 95% confidence intervals (95% CIs) were utilized to assess the strength of associations. Three SNPs were related to hepatoblastoma risk. XPC rs2229090 and XPD rs3810366 significantly contributed to hepatoblastoma risk according to the dominant model (adjusted OR=1.49, 95% CI=1.07-2.08, P=0.019; adjusted OR=1.66, 95% CI=1.12-2.45, P=0.012, respectively). However, XPD rs238406 conferred a significantly decreased risk of hepatoblastoma under the dominant model (adjusted OR=0.68, 95% CI=0.49-0.95; P=0.024). Stratified analysis demonstrated that these significant associations were more prominent in certain subgroups. Moreover, there was evidence of functional implications of these significant SNPs suggested by online expression quantitative trait loci (eQTLs) and splicing quantitative trait loci (sQTLs) analysis. In summary, NER pathway gene polymorphisms (XPC rs2229090, XPD rs3810366, and XPD rs238406) are significantly associated with hepatoblastoma risk, and further research is required to verify these findings.

Association of PHACTR1 with Coronary Artery Calcium Differs by Sex and Cigarette Smoking.

Background: Coronary artery calcium (CAC) is a marker of subclinical atherosclerosis and is a complex heritable trait with both genetic and environmental risk factors, including sex and smoking. Methods: We performed genome-wide association (GWA) analyses for CAC among all participants and stratified by sex in the COPDGene study (n = 6144 participants of European ancestry and n = 2589 participants of African ancestry) with replication in the Diabetes Heart Study (DHS). We adjusted for age, sex, current smoking status, BMI, diabetes, self-reported high blood pressure, self-reported high cholesterol, and genetic ancestry (as summarized by principal components computed within each racial group). For the significant signals from the GWA analyses, we examined the single nucleotide polymorphism (SNP) by sex interactions, stratified by smoking status (current vs. former), and tested for a SNP by smoking status interaction on CAC. Results: We identified genome-wide significant associations for CAC in the chromosome 9p21 region [CDKN2B-AS1] among all COPDGene participants (p = 7.1 × 10-14) and among males (p = 1.0 × 10-9), but the signal was not genome-wide significant among females (p = 6.4 × 10-6). For the sex stratified GWA analyses among females, the chromosome 6p24 region [PHACTR1] had a genome-wide significant association (p = 4.4 × 10-8) with CAC, but this signal was not genome-wide significant among all COPDGene participants (p = 1.7 × 10-7) or males (p = 0.03). There was a significant interaction for the SNP rs9349379 in PHACTR1 with sex (p = 0.02), but the interaction was not significant for the SNP rs10757272 in CDKN2B-AS1 with sex (p = 0.21). In addition, PHACTR1 had a stronger association with CAC among current smokers (p = 6.2 × 10-7) than former smokers (p = 7.5 × 10-3) and the SNP by smoking status interaction was marginally significant (p = 0.03). CDKN2B-AS1 had a strong association with CAC among both former (p = 7.7 × 10-8) and current smokers (p = 1.7 × 10-7) and the SNP by smoking status interaction was not significant (p = 0.40). Conclusions: Among current and former smokers of European ancestry in the COPDGene study, we identified a genome-wide significant association in the chromosome 6p24 region [PHACTR1] with CAC among females, but not among males. This region had a significant SNP by sex and SNP by smoking interaction on CAC.

First Insight into the Whole Genome Sequencing Whole Variations in Mycobacterium bovis from Cattle in Morocco.

Six cattle heads which tested positive against bovine tuberculosis (bTB) in Morocco were investigated to confirm the disease and to determine the source(s) of infection. Polymerase Chain Reaction (PCR) was directly performed on tissue samples collected from slaughtered animals. All investigated animals tested positive to PCR for the Mycobacterium bovis sub-type. Bacteriological isolation was conducted according to the technique recommended by WOAH for the cultivation of the Mycobacterium tuberculosis Complex (MBTC). Whole genome sequencing (WGS) was carried out on six mycobacterial isolates and the phylogenic tree was constructed. The six Moroccan isolates fit with clades II, III, IV, V and VII and were confirmed to belong to the clonal complexes Eu2, Unknown 2 and 7 as well as to sublineages La1.7.1, La1.2 and La1.8.2. The significant Single Nucleotide Polymorphism (SNPs) ranged from 84 to 117 between the isolates and the reference M. bovis strain and from 17 to 212 between the six isolates. Considering the high resolution of WGS, these results suggests that the source of infection of the bTB could be linked to imported animals as five of the investigated reactor animals were imported a few months prior. WGS can be a useful component to the Moroccan strategy to control bTB.

Relationship between maternal and fetal ERAP-1 gene polymorphism and pre-eclampsia

Objective: To investigate the relationship between the polymorphism of endoplasmic reticulum aminopeptidase 1 (ERAP-1) gene and the occurrence of pre-eclampsia (PE). Methods: A case-control study was conducted in Beijing Obstetrics and Gynecology Hospital from October 2018 to October 2021. A total of 51 PE pregnant women with onset gestational age<34 weeks were selected as the PE group, and 48 normal pregnant women during the same period were selected as the control group. Venous blood samples were collected from the pregnant women before delivery and umbilical cord within 5 minutes after delivery. Single nucleotide polymorphisms (SNP) of ERAP-1 gene in the pregnant women and their fetus were detected by next-generation sequencing. Univariate analysis and multivariate logistic regression analysis were used to analyze all the SNP loci and alleles detected in the two groups, and the significant SNP were screened. Results: (1) A total of 13 target SNP loci of maternal ERAP-1 gene were selected by univariate analysis. Among them, the frequency distribution of genotypes at 96096828, 96121524, 96121715, 96122260 and 96122281 showed statistically significant differences between PE group and control group (all P<0.05). Multivariate logistic regression analysis showed that the risk of PE in pregnant women with TC genotype at locus 96121524 was 2.002 times higher than those with TT genotype (95%CI: 0.687-5.831, P=0.020). (2) A total of 4 target SNP loci of ERAP-1 gene in fetal were selected by univariate analysis, and there was no statistical significance in gene polymorphism of the 4 loci between PE group and control group (all P>0.05). Multivariate logistic regression analysis showed that the risk of PE in fetus with genotype AA at locus 96121406 was 0.236 times that of fetus with genotype GG (95%CI: 0.055-1.025, P=0.016). Conclusion: ERAP-1 gene with TC genotype at 96121524 in the mother and GG genotype at 96121406 in the fetus might be related to the incidence of PE.

Genome-wide association study comparison analysis based on Hanwoo full-sib family.

The improvement of carcass traits is essential for the Hanwoo industry because of the Hanwoo grade determination system, and genome-wide association study (GWAS) analysis is an instrumental tool for identifying the genetic factors that impact these traits. While GWAS analysis utilizing family data offers advantages in minimizing genetic bias, research on family-based GWAS in Hanwoo is currently lacking. This study classified Group A using both parental and offspring genetic information, and Group B based solely on offspring genetic information, to compare GWAS analysis results of Hanwoo carcass traits. A total of 16 significant single nucleotide polymorphism (SNP) markers were identified in Group A, comprising 7 for carcass weight (CWT), 3 for back fat thickness (BFT), and 6 for marbling score (MS). In Group B, 7 significant SNP markers were identified, including 3 for CWT, 1 for eye muscle area, 1 for BFT, and 2 for MS. Functional annotation analysis revealed only one common function related to carcass traits between the groups, while protein-protein interaction analysis indicated more gene interactions in Group A. The reliability of estimated values for common SNP markers identified between the groups was higher in Group A. GWAS analysis utilizing parental genetic information holds greater potential for application, owing to its higher reliability of estimated values and the ability to explore numerous candidate genes.

Interaction between Risk Single-Nucleotide Polymorphisms of Developmental Dyslexia and Parental Education on Reading Ability: Evidence for Differential Susceptibility Theory.

While genetic and environmental factors have been shown as predictors of children's reading ability, the interaction effects of identified genetic risk susceptibility and the specified environment for reading ability have rarely been investigated. The current study assessed potential gene-environment (G×E) interactions on reading ability in 1477 school-aged children. The gene-environment interactions on character recognition were investigated by an exploratory analysis between the risk single-nucleotide polymorphisms (SNPs), which were discovered by previous genome-wide association studies of developmental dyslexia (DD), and parental education (PE). The re-parameterized regression analysis suggested that this G×E interaction conformed to the strong differential susceptibility model. The results showed that rs281238 exhibits a significant interaction with PE on character recognition. Children with the "T" genotype profited from high PE, whereas they performed worse in low PE environments, but "CC" genotype children were not malleable in different PE environments. This study provided initial evidence for how the significant SNPs in developmental dyslexia GWA studies affect children's reading performance by interacting with the environmental factor of parental education.

Linkage Mapping and Genome-Wide Association Study Identified Two Peanut Late Leaf Spot Resistance Loci, PLLSR-1 and PLLSR-2, Using Nested Association Mapping.

Identification of candidate genes and molecular markers for late leaf spot (LLS) disease resistance in peanut (Arachis hypogaea) has been a focus of molecular breeding for the U.S. industry-funded peanut genome project. Efforts have been hindered by limited mapping resolution due to low levels of genetic recombination and marker density available in traditional biparental mapping populations. To address this, a multi-parental nested association mapping population has been genotyped with the peanut 58K single-nucleotide polymorphism (SNP) array and phenotyped for LLS severity in the field for 3 years. Joint linkage-based quantitative trait locus (QTL) mapping identified nine QTLs for LLS resistance with significant phenotypic variance explained up to 47.7%. A genome-wide association study identified 13 SNPs consistently associated with LLS resistance. Two genomic regions harboring the consistent QTLs and SNPs were identified from 1,336 to 1,520 kb (184 kb) on chromosome B02 and from 1,026.9 to 1,793.2 kb (767 kb) on chromosome B03, designated as peanut LLS resistance loci, PLLSR-1 and PLLSR-2, respectively. PLLSR-1 contains 10 nucleotide-binding site leucine-rich repeat disease resistance genes. A nucleotide-binding site leucine-rich repeat disease resistance gene, Arahy.VKVT6A, was also identified on homoeologous chromosome A02. PLLSR-2 contains five significant SNPs associated with five different genes encoding callose synthase, pollen defective in guidance protein, pentatricopeptide repeat, acyl-activating enzyme, and C2 GRAM domains-containing protein. This study highlights the power of multi-parent populations such as nested association mapping for genetic mapping and marker-trait association studies in peanuts. Validation of these two LLS resistance loci will be needed for marker-assisted breeding.

Integrating multilocus genome-wide association studies in chickpea landraces to discern the genetics of drought tolerance

In chickpea breeding, drought is a major concern and a complex trait controlled by several genes. To develop drought-tolerant varieties, it is essential to use the available germplasm and genomic resources. Over the years, the landraces have proven to be a good source for the dissection of genes for different yield and yield-related traits. The present investigation for marker–trait associations (MTAs) and candidate gene identification was conducted by studying 125 chickpea landraces collected from the West Asia and North Africa (WANA) region, along with 4 varieties suitable for irrigated and rainfed environments. This study analyzed 13 physio-morphological traits in 2 consecutive years at two isolated locations (IARI, New Delhi, and Dharwad). A strong correlation coefficient was observed between the trait seed yield (SY) and biological yield (BY) under both conditions. The Drought Susceptibility Index (DSI) ranged from 0.02 to 1.84 and 0.10 to 2.04 at the IARI, New Delhi and Dharwad locations, respectively. The genotypic data of 6,367 single nucleotide polymorphisms (SNPs) distributed across the genome were used for genetic diversity study, population structure, and genome-wide association study (GWAS). The average polymorphic information content (PIC) value observed was 0.25, and the average linkage disequilibrium (LD) decay distance was 152,269 bp across the genome. A total of four subgroups were observed within the population for genotypic data. Fixed and random model Circulating Probability Unification (FarmCPU) was used for the GWAS analysis, which considered both fixed- and random-effect models. A total of 52 significant SNPs were reported in both irrigated and rainfed conditions at low locations; 7 SNPs were associated with more than one trait, which may have pleiotropic effects. Significant SNPs were annotated in the pulse database. The identified genomic region found in or near MTA under rainfed conditions encodes for guard cell hydrogen peroxide-resistant1 (GHR1), late embryogenesis-abundant, E3 ubiquitin-protein ligase, walls are thin1 (WAT1), and beta-galactosidase that are known to be associated with drought tolerance.

Whole-genome resequencing of grass carp (Ctenopharyngodon idella) for genome-wide association study on GCRV resistance

Grass carp (Ctenopharyngodon idella) is a highly productive freshwater fish species that annually contributes significant economic value to the aquaculture industry. However, its production is frequently affected by grass carp haemorrhagic disease (GCHD), which is caused by grass carp reovirus (GCRV). To elucidate the molecular mechanism of GCRV resistance in C. idella, whole-genome resequencing of both susceptible and resistant C. idella was conducted, yielding a total of 3,941,776 high quality single nucleotide polymorphisms (SNPs). A genome-wide association study (GWAS) then identified five significant SNPs, located on chromosomes 2, 3, 5, 6 and 10, associated with the trait of GCRV resistance. An integrative analysis of the GWAS and a previous transcriptome study was performed, and five common genes were identified, including b4galt2 (β-1,4-galactosyltransferase 2), kif2c (Kinesin family member 2C), jmjd8 (JmjC domain containing 8), selenom (Selenoprotein M), and znfx1 (Zinc finger NFX1-type containing 1). This study serves as a valuable basis for implementing marker-assisted selection (MAS) to improve GCRV resistance in C. idella breeding programs, and provides insight into developing methods to control the spread of GCHD.

Combining genome-wide association study and linkage mapping in the genetic dissection of amylose content in maize (Zea mays L.).

Integrated linkage and association analysis revealed genetic basis across multiple environments. The genes Zm00001d003102 and Zm00001d015905 were further verified to influence amylose content using gene-based association study. Maize kernel amylose is an important source of human food and industrial raw material. However, the genetic basis underlying maize amylose content is still obscure. Herein, we used an intermated B73 × Mo17 (IBM) Syn10 doubled haploid population composed of 222 lines and a germplasm set including 305 inbred lines to uncover the genetic control for amylose content under four environments. Linkage mapping detected 16 unique QTL, among which four were individually repeatedly identified across multiple environments. Genome-wide association study revealed 17 significant (P = 2.24E-06) single-nucleotide polymorphisms, of which two (SYN19568 and PZE-105090500) were located in the intervals of the mapped QTL (qAC2 and qAC5-3), respectively. According to the two population co-localized loci, 20 genes were confirmed as the candidate genes for amylose content. Gene-based association analysis indicated that the variants in Zm00001d003102 (Beta-16-galactosyltransferase GALT29A) and Zm00001d015905 (Sugar transporter 4a) affected amylose content across multi-environment. Tissue expression analysis showed that the two genes were specifically highly expressed in the ear and stem, respectively, suggesting that they might participate in sugar transport from source to sink organs. Our study provides valuable genetic information for breeding maize varieties with high amylose.

Genome-wide association study of stem structural characteristics that extracted by a high-throughput phenotypic analysis "LabelmeP rice" in rice.

Stem is important for assimilating transport and plant strength; however, less is known about the genetic basis of its structural characteristics. In this study, a high-throughput method, "LabelmeP rice" was developed to generate 14 traits related to stem regions and vascular bundles, which allows the establishment of a stem cross-section phenotype dataset containing anatomical information of 1738 images from hand-cut transections of stems collected from 387 rice germplasm accessions grown over two successive seasons. Then, the phenotypic diversity of the rice accessions was evaluated. Genome-wide association studies identified 94, 83, and 66 significant single nucleotide polymorphisms (SNPs) for the assayed traits in 2 years and their best linear unbiased estimates, respectively. These SNPs can be integrated into 29 quantitative trait loci (QTL), and 11 of them were common in 2 years, while correlated traits shared 19. In addition, 173 candidate genes were identified, and six located at significant SNPs were repeatedly detected and annotated with a potential function in stem development. By using three introgression lines (chromosome segment substitution lines), four of the 29 QTLs were validated. LOC_Os01g70200, located on the QTL uq1.4, is detected for the area of small vascular bundles (SVB) and the rate of large vascular bundles number to SVB number. Besides, the CRISPR/Cas9 editing approach has elucidated the function of the candidate gene LOC_Os06g46340 in stem development. In conclusion, the results present a time- and cost-effective method that provides convenience for extracting rice stem anatomical traits and the candidate genes/QTL, which would help improve rice.

Genome-wide identification of quantitative trait loci and candidate genes for seven carcass traits in a four-way intercross porcine population

BackgroundCarcass traits are essential economic traits in the commercial pig industry. However, the genetic mechanism of carcass traits is still unclear. In this study, we performed a genome-wide association study (GWAS) based on the specific-locus amplified fragment sequencing (SLAF-seq) to study seven carcass traits on 223 four-way intercross pigs, including dressing percentage (DP), number of ribs (RIB), skin thinkness (ST), carcass straight length (CSL), carcass diagonal length (CDL), loin eye width (LEW), and loin eye thickness (LET).ResultsA total of 227,921 high-quality single nucleotide polymorphisms (SNPs) were detected to perform GWAS. A total of 30 SNPs were identified for seven carcass traits using the mixed linear model (MLM) (p < 1.0 × 10− 5), of which 9 SNPs were located in previously reported quantitative trait loci (QTL) regions. The phenotypic variation explained (PVE) by the significant SNPs was from 2.43 to 16.32%. Furthermore, 11 candidate genes (LYPLAL1, EPC1, MATN2, ZFAT, ZBTB10, ZNF704, INHBA, SMYD3, PAK1, SPTBN2, and ACTN3) were found for carcass traits in pigs.ConclusionsThe GWAS results will improve our understanding of the genetic basis of carcass traits. We hypothesized that the candidate genes associated with these discovered SNPs would offer a biological basis for enhancing the carcass quality of pigs in swine breeding.

Genome-wide association analysis and transgenic characterization for amylose content regulating gene in tuber of Dioscorea zingiberensis

BackgroundAmylose, a prebiotic found in yams is known to be beneficial for the gut microflora and is particularly advantageous for diabetic patients’ diet. However, the genetic machinery underlying amylose production remains elusive. A comprehensive characterization of the genetic basis of amylose content in yam tubers is a prerequisite for accelerating the genetic engineering of yams with respect to amylose content variation.ResultsTo uncover the genetic variants underlying variation in amylose content, we evaluated amylose content in freshly harvested tubers from 150 accessions of Dioscorea zingibensis. With 30,000 high-quality single nucleotide polymorphisms (SNP), we performed a genome-wide association analysis (GWAS). The population structure analysis classified the D. zingiberensis accessions into three groups. A total of 115 significant loci were detected on four chromosomes. Of these, 112 significant SNPs (log10(p) = 5, q-value < 0.004) were clustered in a narrow window on the chromosome 6 (chr6). The peak SNP at the position 75,609,202 on chr6 could explain 63.15% of amylose variation in the population and fell into the first exon of the ADP-glucose pyrophosphorylase (AGPase) small subunit gene, causing a non-synonymous modification of the resulting protein sequence. Allele segregation analysis showed that accessions with the rare G allele had a higher amylose content than those harboring the common A allele. However, AGPase, a key enzyme precursor of amylose biosynthesis, was not expressed differentially between accessions with A and G alleles. Overexpression of the two variants of AGPase in Arabidopsis thaliana resulted in a significantly higher amylose content in lines transformed with the AGPase-G allele.ConclusionsOverall, this study showed that a major genetic variant in AGPase probably enhances the enzyme activity leading to high amylose content in D. zingiberensis tuber. The results provide valuable insights for the development of amylose-enriched genotypes.

Genotyping-by-Sequencing Analysis Reveals Associations between Agronomic and Oil Traits in Gamma Ray-Derived Mutant Rapeseed (Brassica napus L.).

Rapeseed (Brassica napus L.) holds significant commercial value as one of the leading oil crops, with its agronomic features and oil quality being crucial determinants. In this investigation, 73,226 single nucleotide polymorphisms (SNPs) across 95 rapeseed mutant lines induced by gamma rays, alongside the original cultivar ('Tamra'), using genotyping-by-sequencing (GBS) analysis were examined. This study encompassed gene ontology (GO) analysis and a genomewide association study (GWAS), thereby concentrating on agronomic traits (e.g., plant height, ear length, thousand-seed weight, and seed yield) and oil traits (including fatty acid composition and crude fat content). The GO analysis unveiled a multitude of genes with SNP variations associated with cellular processes, intracellular anatomical structures, and organic cyclic compound binding. Through GWAS, we detected 320 significant SNPs linked to both agronomic (104 SNPs) and oil traits (216 SNPs). Notably, two novel candidate genes, Bna.A05p02350D (SFGH) and Bna.C02p22490D (MDN1), are implicated in thousand-seed weight regulation. Additionally, Bna.C03p14350D (EXO70) and Bna.A09p05630D (PI4Kα1) emerged as novel candidate genes associated with erucic acid and crude fat content, respectively. These findings carry implications for identifying superior genotypes for the development of new cultivars. Association studies offer a cost-effective means of screening mutants and selecting elite rapeseed breeding lines, thereby enhancing the commercial viability of this pivotal oil crop.

A high-density linkage map and sex-determination loci in Pacific white shrimp (Litopenaeus vannamei)

BackgroundExpansion of genomic resources for the Pacific white shrimp (Litopenaeus vannamei), such as the construction of dense genetic linkage maps, is crucial for the application of genomic tools in order to improve economically relevant traits. Sexual dimorphism exists in Pacific white shrimp, and the mapping of the sex-determination region in this species may help in future reproductive applications. We have constructed male, female, and sex-averaged high-density genetic maps using a 50 K single-nucleotide polymorphism (SNP) array, followed by a genome-wide association study (GWAS) to identify genomic regions associated with sex in white shrimp.ResultsThe genetic map yielded 15,256 SNPs assigned to 44 linkage groups (LG). The lengths of the male, female, and sex-averaged maps were 5,741.36, 5,461.20 and 5,525.26 cM, respectively. LG18 was found to be the largest for both sexes, whereas LG44 was the shortest for males and LG31 for females. A sex-determining region was found in LG31 with 21 statistically significant SNPs. The most important SNP was previously identified as a sex-linked marker and was able to identify 99% of the males and 88% of the females. Although other significant markers had a lower ability to determine sex, putative genes were intercepted or close to them. The oplophorus-luciferin 2-monooxygenase, serine/arginine repetitive matrix protein and spermine oxidase genes were identified as candidates with possible participation in important processes of sexual differentiation in shrimp.ConclusionsOur results provide novel genomic resources for shrimp, including a high-density linkage map and new insights into the sex-determining region in L. vannamei, which may be usefulfor future genetics and reproduction applications.

Gene-gene interaction analysis for age at onset of bipolar disorder in a Korean population

BackgroundMultiple genes might interact to determine the age at onset of bipolar disorder. We investigated gene-gene interactions related to age at onset of bipolar disorder in the Korean population, using genome-wide association study (GWAS) data. MethodsThe study population consisted of 303 patients with bipolar disorder. First, the top 1000 significant single-nucleotide polymorphisms (SNPs) associated with age at onset of bipolar disorder were selected through single SNP analysis by simple linear regression. Subsequently, the QMDR method was used to find gene-gene interactions. ResultsThe best 10 SNPs from simple regression were located in chromosome 1, 2, 3, 10, 11, 14, 19, and 21. Only five SNPs were found in several genes, such as FOXN3, KIAA1217, OPCML, CAMSAP2, and PTPRS. On QMDR analyses, five pairs of SNPs showed significant interactions with a CVC exceeding 1/5 in a two-locus model. The best interaction was found for the pair of rs60830549 and rs12952733 (CVC = 1/5, P < 1E−07). In three-locus models, four combinations of SNPs showed significant associations with age at onset, with a CVC of >1/5. The best three-locus combination was rs60830549, rs12952733, and rs12952733 (CVC = 2/5, P < 1E−6). The SNPs showing significant interactions were located in the KIAA1217, RBFOX3, SDK2, CYP19A1, NTM, SMYD3, and RBFOX1 genes. ConclusionsOur analysis confirmed genetic interactions influencing the age of onset for bipolar disorder and identified several potential candidate genes. Further exploration of the functions of these promising genes, which may have multiple roles within the neuronal network, is necessary.

Structural variants in linkage disequilibrium with GWAS-significant SNPs

With the recent expansion of structural variant identification in the human genome, understanding the role of these impactful variants in disease architecture is critically important. Currently, a large proportion of genome-wide-significant genome-wide association study (GWAS) single nucleotide polymorphisms (SNPs) are functionally unresolved, raising the possibility that some of these SNPs are associated with disease through linkage disequilibrium with causal structural variants. Hence, understanding the linkage disequilibrium between newly discovered structural variants and statistically significant SNPs may provide a resource for further investigation into disease-associated regions in the genome. Here we present a resource cataloging structural variant-significant SNP pairs in high linkage disequilibrium. The database is composed of (i) SNPs that have exhibited genome-wide significant association with traits, primarily disease phenotypes, (ii) newly released structural variants (SVs), and (iii) linkage disequilibrium values calculated from unphased data. All data files including those detailing SV and GWAS SNP associations and results of GWAS-SNP-SV pairs are available at the SV-SNP LD Database and can be accessed at ‵https://github.com/hliang-SchrodiLab/SV_SNPs. Our analysis results represent a useful fine mapping tool for interrogating SVs in linkage disequilibrium with disease-associated SNPs. We anticipate that this resource may play an important role in subsequent studies which investigate incorporating disease causing SVs into disease risk prediction models.

Genome-wide association studies of root system architecture traits in a broad collection of Brassica genotypes.

The root systems of Brassica species are complex. Eight root system architecture (RSA) traits, including total root length, total root surface area, root average diameter, number of tips, total primary root length, total lateral root length, total tertiary root length, and basal link length, were phenotyped across 379 accessions representing six Brassica species (B. napus, B. juncea, B. carinata, B. oleracea, B. nigra, and B. rapa) using a semi-hydroponic system and image analysis software. The results suggest that, among the assessed species, B. napus and B. oleracea had the most intricate and largest root systems, while B. nigra exhibited the smallest roots. The two species B. juncea and B. carinata shared comparable root system complexity and had root systems with larger root diameters. In addition, 313 of the Brassica accessions were genotyped using a 19K Brassica single nucleotide polymorphism (SNP) array. After filtering by TASSEL 5.0, 6,213 SNP markers, comprising 5,103 markers on the A-genome (covering 302,504 kb) and 1,110 markers on the C-genome (covering 452,764 kb), were selected for genome-wide association studies (GWAS). Two general linear models were tested to identify the genomic regions and SNPs associated with the RSA traits. GWAS identified 79 significant SNP markers associated with the eight RSA traits investigated. These markers were distributed across the 18 chromosomes of B. napus, except for chromosome C06. Sixty-five markers were located on the A-genome, and 14 on the C-genome. Furthermore, the major marker-trait associations (MTAs)/quantitative trait loci (QTLs) associated with root traits were located on chromosomes A02, A03, and A06. Brassica accessions with distinct RSA traits were identified, which could hold functional, adaptive, evolutionary, environmental, pathological, and breeding significance.