Revisiting the association of sudden infant death syndrome (SIDS) with polymorphisms of NHE3 and IL13

OBJECTIVES

Disturbances of the central nervous system and immune system are thought to play a role in sudden infant death syndrome (SIDS). Dysregulated expression of sodium (Na+)/hydrogen (H+) exchanger 3 (NHE3) in the brainstem and of interleukin 13 (IL13) in the lungs has been observed in SIDS. An association of single-nucleotide polymorphisms (SNPs) in NHE3 and IL13 with SIDS has been proposed, but controversial results were reported. Therefore, there is a need to revisit the association of SNPs in NHE3 and IL13 with SIDS.

METHODS

Genotyping of rs71597645 (G1131A) and rs2247114 (C2405T) in NHE3 and rs20541 (+ 4464A/G) in IL13 was performed in 201 SIDS cases and 338 controls. A meta-analysis was performed after merging our data with previously published data (all from European populations).

RESULTS

Polymorphisms rs2247114 (NHE3) and rs20541 (IL13) were significantly associated with SIDS overall and in multiple subgroups, but no association was found for rs71597645 (NHE3). After combining our data with previously published data, a fixed-effect meta-analysis showed that rs2247114 in NHE3 retained a significant association with SIDS under a recessive model (OR 2.78, 95%CI 1.53 to 5.06; p = 0.0008).

CONCLUSION

Our findings suggest an association of NHE3 variant rs2247114 (C2405T), though not rs71597645 (NHE3), with SIDS. A potential role of rs20541 (IL13) still has to be elucidated. Especially NHE3 seems to be an interesting topic for future SIDS research.

A mutation in LPAR2 activates the miR-939-5p-LPAR2-PI3K/AKT axis to regulate the proliferation and apoptosis of granulosa cells in sheep

Lysophosphatidic acid receptor-2 (LPAR2) is a G protein-coupled receptor, which is involved in various physiological processes such as cell development, proliferation, and apoptosis, and is thought to play an important role in follicular development and reproduction. There is evidence that miRNA recognition elements (MRE) in the gene 3'UTR often contain single nucleotide polymorphisms (SNPs) that can alter the binding affinity of the target miRNA, leading to dysregulation of gene expression. In this study, we detected a SNP in LPAR2 3 'UTR (rs410670692, c.*701C > T) in 384 small-tailed Han sheep using Sequenom MassARRAY®SNP genotyping. Association analysis showed that the SNP was significantly associated with litter size. Then, the effect of LPAR2 rs410670692 mutation on gene expression in sheep hosts was studied by molecular biotechnology. The results showed that the expression of LPAR2 in the TT genotype was significantly higher than that in the CC genotype, which confirmed the existence of rs410670692, a functional SNP, in LPAR2 3'UTR. We then used bioinformatics methods and double luciferase reporter gene assay to predict and confirm LPAR2 SNP rs410670692 as the direct targeting regulatory element of miR-939-5p. Cell transfection experiments further found that SNP rs410670692 down-regulated the mRNA and protein levels of LPAR2 by influencing the binding of miR-939-5p. To understand the function and mechanism of miR-939-5p in sheep granulosa cells (GCs), we conducted cell proliferation and apoptosis experiments which showed inhibited GCs proliferation along with promoted GCs apoptosis upon overexpression of miR-939-5p. Moreover, overexpression of miR-939-5p promotes apoptosis of granulosa cells by blocking the LPAR2-dependent PI3K/Akt signaling pathway. In conclusion, these results indicate that the SNP rs410670692 of LPAR2 is related to the litter size of small-tailed cold sheep, and miR-939-5p can act as a regulatory element binding to the C mutation of rs410670692 to regulate the expression of LPAR2, affect the development of GCs, and thus indirectly affect the litter size of sheep. These studies provide evidence for the involvement of LPAR2 polymorphism in sheep reproduction and are expected to provide new insights into the molecular genetic mechanisms of litter size traits in sheep.

The potential up-regulation risk of 3' UTR SNP (rs10787760 G > A) for the VAX1 gene is associated with NSCLP in the northwest Chinese population

AIMS

To investigate the association between single nucleotide polymorphisms (SNPs) in 3'UTR region of VAX1, SYT14 and PAX7 genes and the risk of non-syndromic cleft palate (NSCLP) in a northwest Chinese population.

MAIN METHODS

A case-control study was conducted in 406 normal controls and 399 NSCLP patients. Using iMLDRTM genotyping technology, eight SNPs of three genes ((rs10787760, rs7086344 at VAX1), (rs1010113, rs851114, and rs485874 at PAX7), and (rs61820397, rs4609425, rs12133399 at SYT14)) were genotyped to investigate the differences in alleles and genotype distribution frequencies between NSCLP patients and healthy controls. RNA Folding Form software was used to predict RNA secondary structure and expression vectors were constructed to explore the function of the relevant SNP. The effect of SNP polymorphism of gene transcription and translation was assessed using qPCR and Western blot analysis.

KEY FINDINGS

Among the eight SNPs of three genes, rs10787760 of VAX1 gene was found to be associated with an increased risk of NSCLP (OR = 1.341,CI = 1.004-1.790) and the GA genotype of rs10787760 increased the risk of cleft lip and/or palate (CL/P) about 1.42 times (p < 0.05), and carrying the A allele might increase the risk of NSCL/P in male (OR = 1.356, 95 % CI = 1.010-1.823). But there was no association observed with cleft palate only (CPO). Cell function experiments revealed that the G to A mutation in rs10787760 up-regulated GFP-VAX1 transcriptional level by 2.39 and 3.13 times in two cell lines respectively, and enhance the protein expression of the VAX1 gene further. RNA secondary structure study showed that the rs10787760 (G > A) had two different secondary structures in 3'UTR region.

SIGNIFICANCE

The rs10787760 variant in the 3'UTR region of VAX1 gene is associated with CL/P in northwest Chinese population. We hypothesize that the machanism of it might be caused by the RNA differenct fold in the 3'UTR region caused by the polymorphism of the gene.

LEVEL OF EVIDENCE

Original Reports.

Extended application of PGT-M strategies for small pathogenic CNVs

PURPOSE

The preimplantation genetic testing for aneuploidy (PGT-A) platform is not currently available for small copy-number variants (CNVs), especially those < 1 Mb. Through strategies used in PGT for monogenic disease (PGT-M), this study intended to perform PGT for families with small pathogenic CNVs.

METHODS

Couples who carried small pathogenic CNVs and underwent PGT at the Reproductive and Genetic Hospital of CITIC-Xiangya (Hunan, China) between November 2019 and April 2023 were included in this study. Haplotype analysis was performed through two platforms (targeted sequencing and whole-genome arrays) to identify the unaffected embryos, which were subjected to transplantation. Prenatal diagnosis using amniotic fluid was performed during 18-20 weeks of pregnancy.

RESULTS

PGT was successfully performed for 20 small CNVs (15 microdeletions and 5 microduplications) in 20 families. These CNVs distributed on chromosomes 1, 2, 6, 7, 13, 15, 16, and X with sizes ranging from 57 to 2120 kb. Three haplotyping-based PGT-M strategies were applied. A total of 89 embryos were identified in 25 PGT cycles for the 20 families. The diagnostic yield was 98.9% (88/89). Nineteen transfers were performed for 17 women, resulting in a 78.9% (15/19) clinical pregnancy rate after each transplantation. Of the nine women who had healthy babies, eight accepted prenatal diagnosis and the results showed no related pathogenic CNVs.

CONCLUSION

Our results show that the extended haplotyping-based PGT-M strategy application for small pathogenic CNVs compensated for the insufficient resolution of PGT-A. These three PGT-M strategies could be applied to couples with small pathogenic CNVs.

Genomic sequencing revealed recombination event between clade 1 and clade 2 occurs in circulating varicella-zoster virus in China

Varicella-zoster virus (VZV), a member of the Alphaherpesvirinae subfamily, causes varicella in primary infections and establishing a latent stage in sensory ganglia. Upon reactivation, VZV causes herpes zoster with severe neuralgia, especially in elderly patients. The mutation rate for VZV is comparatively lower than the other members of other alpha herpesviruses. Due to geographic isolation, different genotypes of VZV are circulating on separate continents. Here, we successfully isolated a VZV from the vesicular fluid of a youth zoster patient. Based on the single-nucleotide polymorphism profiles of different open reading frames that define the genotype, this newly isolated VZV primarily represents genotype clade 2 but also has characteristics of genotype clade 1. The next-generation sequencing provided a nearly full-length sequence, and further phylogenetic analysis revealed that this VZV isolate is distinct from clades 1 and 2. The Recombination Detection Program indicates that a possible recombinant event may occur between the VZV isolate and clade 1. In summary, we found that there is a circulating VZV isolate in China that may represent a recombinant between clade 1 and clade 2, providing new concerns that need to be considered in the future VZV vaccination program.

Metabolic syndrome and risk of colorectal cancer: A Mendelian randomization study

Background

Observational studies have previously demonstrated a significant relationship among both metabolic syndrome (Mets) and colorectal cancer (CRC). Whether there is a causal link remains controversial.

Objective

To clarify whether Mets and their components have a causal effect on colorectal cancer, we have carried out a bidirectional Mendelian randomization analysis (MR).

Methods

This study started from genome-wide association data for Mets and its 5 components (hypertension, waist circumference, fasting blood glucose, serum triglycerides, and serum high-density lipoprotein cholesterol) and colorectal cancer. Mendelian randomization (MR) techniques were used in the study to examine their associations.

Results

After Benjamini-Hochberg multiple corrections, genetically predicted significant causal link exists between WC (waist circumference) and CRC. The OR was 1.35 (95 % CI: 1.08-1.69; p = 0.0096). Other Mets components (HBP, FBG, TG, HDL), on the other hand, found no evidence of a genetic link between CRC and Mets. In addition, MR results showed that CRC was not causally related to either Mets or the components. We get the same result in the validated dataset.

Conclusion

According to the bidirectional MR investigation shows a significant causal relationship among obesity and CRC in the Mets component but no causal relationship in the opposite direction.

Mutations in Coding and Non-Coding Regions in Varicella-Zoster Virus Causing Fatal Hemorrhagic Fever Without Rash in an Immunocompetent Patient: Case Report

INTRODUCTION

We report the case of a fatal hemorrhagic varicella primary infection in an immunocompetent man and whole-genome characterization of the virus for the investigation of biomarkers of virulence.

CASE

A 38-year-old patient born in Nigeria presented to the emergency department with abdominal pain and subsequently developed fatal hemorrhagic disease without skin rash. Extensive laboratory tests including serology and PCR for arenaviruses, bunyaviruses and ebolaviruses were negative. Varicella-zoster virus (VZV) PCR of sera, liver and spleen tissue samples from autopsy revealed the presence of VZV DNA. Primary infection by varicella-zoster virus with hemorrhagic manifestations was diagnosed after virological testing. The VZV genome was sequenced using a mWGS approach. Bioinformatic analysis showed 53 mutations across the genome, 33 of them producing non-synonymous variants affecting up to 14 genes. Some of them, such as ORF11 and ORF 62, encoded for essential functions related to skin or neurotropism. To our knowledge, the mutations reported here have never been described in a VZV causing such a devastating outcome.

DISCUSSION

In immunocompetent patients, viral factors should be considered in patients with uncommon symptoms or severe diseases. Some relevant mutations revealed by using whole genome sequencing (WGS) directly from clinical samples may be involved in this case and deserves further investigation.

CONCLUSION

Differential diagnosis of varicella-zoster virus in immunocompetent adults should be considered among patients with suspected VHF, even if the expected vesicular rash is not present at admission and does not arise thereafter. Whole genome sequencing of strains causing uncommon symptoms and/or mortality is needed for epidemiological surveillance and further characterization of putative markers of virulence. Additionally, this report highlights the recommendation for a VZV vaccination policy in non-immunized migrants from developing countries.

A homozygous missense mutation in the fibroblast growth factor 5 gene is associated with the long-hair trait in Angora rabbits

BACKGROUND

Rabbits are well-domesticated animals. As a crucial economic animal, rabbit has been successfully bred into wool-use, meat-use and fur-use breeds. Hair length is one of the most economically important traits affecting profitability in wool rabbits. In this study, to identify selection signatures with the long-hair trait, whole-genomic resequencing of long-haired rabbits (Angora rabbits) and short-haired rabbits (Rex and New Zealand rabbits) was performed.

RESULTS

By genome-wide selective sweeping analysis based on population comparison, we identified a total of 5.85 Mb regions (containing 174 candidate genes) with strong selection signals. Six of these genes (Dusp1, Ihh, Fam134a, Map3k1, Spata16, and Fgf5) were enriched in the MAPK signalling and Hedgehog signalling pathways, both of which are closely associated with hair growth regulation. Among these genes, Fgf5 encodes the FGF5 protein, which is a well-established regulator of hair growth. There was a nonsynonymous nucleotide substitution (T19234C) in the Fgf5 gene. At this locus, the C allele was present in all of the tested Angora rabbits, while the T allele was dominant in New Zealand and Rex rabbits. We further confirmed that the C allele was conserved in Angora rabbits by screening an additional 135 rabbits. Moreover, the results of functional predictions and co-immunoprecipitation revealed that the T19234C mutation impaired the binding capacity of FGF5 to its receptor FGFR1.

CONCLUSIONS

We discovered that the homozygous missense mutation T19234C within Fgf5 might contribute to the long-hair trait of Angora rabbits by reducing its receptor binding capacity. This finding will provide new insights into the genetic basis underlying the genetic improvement of Angora rabbits and benefit the improvement of rabbit breeding in the future.

Introduction of Genetic Mutations Into Mice by Base Editor and Target-AID

Generating genetically modified animal models that precisely recapitulate disease characteristics forms an integral and indispensable tool to understanding disease pathophysiology. Recently, important advances in genome editing technologies have enabled us to efficiently create sophisticated animal models in short periods of time. Base editing is a modified CRISPR/Cas system that induces base substitution at targeted genomic regions. Here I describe a basic protocol to introduce disease-relevant pathogenic mutations into mice utilizing two representative base editing tools, Base Editor and Target-AID.

Is the Association of the Rare rs35667974 IFIH1 Gene Polymorphism With Autoimmune Diseases a Case of RNA Epigenetics?

Interferon induced with helicase C domain-containing protein 1 (IFIH1) gene encodes a cytoplasmic RNA helicase otherwise known as melanoma differentiation-associated 5 (MDA5), a RIG-1-like RNA helicase that recognizes viral RNA and is involved in innate immunity through recognition of viral RNA. Upon binding to double-stranded (ds) RNA, MDA5 forms a filamentous assembly along the length of dsRNA and utilizes molecular signatures to discriminate self, versus non-self on the basis of dsRNA length and methylation. Its missense variant rs35667974 is protective for type 1 diabetes, psoriasis, and psoriatic arthritis, but is also found to be associated with an increased risk for ankylosing spondylitis, Crohn's disease, and ulcerative colitis. To gain insight into the complex role of this variant we performed a structural analysis of MDA5 in complex with dsRNA using molecular dynamics simulations. Our data suggest that while the Ile923Val mutation of the rs35667974 variant does not affect binding to native dsRNA significantly, it displays a destabilizing effect in the presence of 2'-O uridine methylation. Thus, the presence of 2'-O-methylation at the dsRNA introduces a sensing signature that leads to selective reduction of the overall MDA catalytic activity. This study represents an evaluation of the role of the shared rs35667974 variant of autoimmune locus IFIH1, reported to lead to selectively reduced catalytic activity of the modified MDA5 phenotype and, as a consequence, reduced negative feedback on cytokine and chemokine signaling and selectively protection against autoimmunity.

Detection of Allele-Specific 3D Chromatin Interactions Using High-Resolution In-Nucleus 4C-seq

Chromosome conformation capture techniques are a set of methods used to determine 3D genome organization through the capture and identification of physical contacts between pairs of genomic loci. Among them, 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing) allows for the identification and quantification of the sequences interacting with a preselected locus of interest. 4C-seq has been widely used in the literature, mainly to study chromatin loops between enhancers and promoters or between CTCF binding sites and to identify chromatin domain boundaries. As 3D-contacts may be established in an allele-specific manner, we describe an up-to-date allele-specific 4C-seq protocol, starting from the selection of allele-specific viewpoints to Illumina sequencing. This protocol has mainly been optimized for cultured mammalian cells, but can be adapted for other cell types with relatively minor changes in initial steps.

Association of PARK-2 Non-synonyms Polymorphisms and Their In Silico Validation Among North Indian Colorectal Cancer Patients

PURPOSE

PARK2 is a potential tumour suppressor gene and its genetic alterations (regionic loss) are common across many human cancers. The association of PARK2 germline variations (SNPs) with Parkinson's has been shown, but their association in development and progression of cancer remains elusive. The aim of this study was to identify association of PARK2 polymorphisms (rs1801474, rs1801334) with colorectal cancer in a case control study design.

METHODS

This case control study included a total of 650 genetically unrelated subjects comprising 300 colorectal cancer cases and 350 healthy controls belonging to North Indian. Both SNPs were analyzed using the PCR-RFLP assay. Statistical analysis for describing risk and association was performed using SPSS-17 software. Structural deviations due to non- synonymous substitutions (S167N and D394N) were analyzed using MD simulations.

RESULTS

The genotype distributions of both the SNPs were in Hardy-Weinberg equilibrium. For both the polymorphisms, the allelic model showed statistically significant risk with OR ~ 1.3. Many of the associations remained significant even after Bonferroni correction (P < 0.00125). The result suggested that both S167N and D394N were deviated from wild type and structures and were stable after 5 ns. The average value of RMSD for backbone atoms was calculated from 5 to 10 ns molecular dynamics simulation data.

CONCLUSION

In conclusion, our study revealed a significant association of PARK2 SNPs with colorectal cancer as well as their relations with other clinical parameters highlighting their contribution towards colorectal cancer susceptibility in North Indian population.

Splicing and editing of ionotropic glutamate receptors: a comprehensive analysis based on human RNA-Seq data

Ionotropic glutamate receptors (iGluRs) play key roles for signaling in the central nervous system. Alternative splicing and RNA editing are well-known mechanisms to increase iGluR diversity and to provide context-dependent regulation. Earlier work on isoform identification has focused on the analysis of cloned transcripts, mostly from rodents. We here set out to obtain a systematic overview of iGluR splicing and editing in human brain based on RNA-Seq data. Using data from two large-scale transcriptome studies, we established a workflow for the de novo identification and quantification of alternative splice and editing events. We detected all canonical iGluR splice junctions, assessed the abundance of alternative events described in the literature, and identified new splice events in AMPA, kainate, delta, and NMDA receptor subunits. Notable events include an abundant transcript encoding the GluA4 amino-terminal domain, GluA4-ATD, a novel C-terminal GluD1 (delta receptor 1) isoform, GluD1-b, and potentially new GluK4 and GluN2C isoforms. C-terminal GluN1 splicing may be controlled by inclusion of a cassette exon, which shows preference for one of the two acceptor sites in the last exon. Moreover, we identified alternative untranslated regions (UTRs) and species-specific differences in splicing. In contrast, editing in exonic iGluR regions appears to be mostly limited to ten previously described sites, two of which result in silent amino acid changes. Coupling of proximal editing/editing and editing/splice events occurs to variable degree. Overall, this analysis provides the first inventory of alternative splicing and editing in human brain iGluRs and provides the impetus for further transcriptome-based and functional investigations.

Developmental validation of the MGIEasy Signature Identification Library Prep Kit, an all-in-one multiplex system for forensic applications

Analyzing genetic markers in nuclear and mitochondrial genomes is helpful in various forensic applications, such as individual identifications and kinship analyses. However, most commercial kits detect these markers separately, which is time-consuming, laborious, and more error-prone (mislabelling, contamination, ...). The MGIEasy Signature Identification Library Prep Kit (hereinafter "MGIEasy identification system"; MGI Tech, Shenzhen, China) has been designed to provide a simple, fast, and robust way to detect appropriate markers in one multiplex PCR reaction: 52 autosomal STRs, 27 X-chromosomal STRs, 48 Y-chromosomal STRs, 145 identity-informative SNPs, 53 ancestry-informative SNPs, 29 phenotype-informative SNPs, and the hypervariable regions of mitochondrial DNA (mtDNA). Here, we validated the performance of MGIEasy identification system following the guidelines of the Scientific Working Group on DNA Analysis Methods (SWGDAM), assessing species specificity, sensitivity, mixture identification, stability under non-optimal conditions (degraded samples, inhibitor contamination, and various substrates), repeatability, and concordance. Libraries prepared using MGIEasy identification system were sequenced on a MGISEQ-2000 instrument (MGI Tech). MGIEasy-derived STR, SNP, and mtDNA genotypes were highly concordant with CE-based STR genotypes (99.79%), MiSeq FGx-based SNP genotypes (99.78%), and Sanger-based mtDNA genotypes (100%), respectively. This system was strongly human-specific, resistant to four common PCR inhibitors, and reliably amplified both low quantities of DNA (as low as 0.125 ng) and degraded DNA (~ 150 nt). Most of the unique alleles from the minor contributor were detected in 1:10 male-female and male-male mixtures; some minor Y-STR alleles were even detected in 1:1000 male-female mixtures. MGIEasy also successfully directly amplified markers from blood stains on FTA cards, filter papers, and swabs. Thus, our results demonstrated that MGIEasy identification system was suitable for use in forensic analyses due to its robust and reliable performance on samples of varying quality and quantity.

Ancestry-specific profiles of genetic determinants of severe hypertriglyceridemia

BACKGROUND

Susceptibility to severe hypertriglyceridemia (HTG), defined as plasma triglyceride (TG) levels ≥10 mmol/L (880 mg/dL), is conferred by both heterozygous rare variants in five genes involved in TG metabolism and numerous common single-nucleotide polymorphisms (SNPs) associated with TG levels.

OBJECTIVE

To date, these genetic susceptibility factors have been comprehensively assessed primarily in severe HTG patients of European ancestry. Here, we expand our analysis to HTG patients of East Asian and Hispanic ancestry.

METHODS

The genomic DNA of 336, 63 and 199 severe HTG patients of European, East Asian and Hispanic ancestry, respectively, was evaluated using a targeted next-generation sequencing panel to screen for: 1) rare variants in LPL, APOA5, APOC2, GPIHBP1 and LMF1; 2) common, small-to-moderate effect SNPs, quantified using a polygenic score; and 3) common, large-effect polymorphisms, APOA5 p.G185C and p.S19W.

RESULTS

While the proportion of individuals with high polygenic scores was similar, frequency of rare variant carriers varied across ancestries. Compared with ancestry-matched controls, Hispanic patients were the most likely to have a rare variant (OR = 5.02; 95% CI 3.07-8.21; p < 0.001), while European patients were the least likely (OR = 2.56; 95% CI 1.58-4.13; p < 0.001). The APOA5 p.G185C polymorphism, exclusive to East Asians, was significantly enriched in patients compared with controls (OR = 10.1; 95% CI 5.6-18.3; p < 0.001), showing the highest enrichment among the measured genetic factors.

CONCLUSION

While TG-associated rare variants and common SNPs are both found in statistical excess in severe HTG patients of different ancestral backgrounds, the overall genetic profiles of each ancestry group were distinct.

COVID-19 susceptibility: potential of ACE2 polymorphisms

Background

Angiotensin-converting enzyme 2 (ACE2) is a metallopeptidase that primarily functions as a negative regulator of renin angiotensin system (RAS) by converting angiotensin II (Ang II) to angiotensin 1-7. Contrary to this, another RAS component, angiotensin-converting enzyme (ACE) catalyzes synthesis of Ang II from angiotensin I (Ang I) that functions as active compound in blood pressure regulation. This indicates importance of ACE/ACE2 level in regulating blood pressure by targeting Ang II. An outbreak of severe acute respiratory syndrome (SARS) highlighted the additional role of ACE2 as a receptor for SARS coronavirus (SARS-CoV) infection.

Main body of the abstract

ACE2 is a functional receptor for SARS-CoV and SARS-CoV-2. Activation of spike (S)-protein by either type II transmembrane serine proteases (TTSPs) or cathepsin-mediated cleavage initiates receptor binding and viral entry. In addition to TTSPs, ACE2 can also be trimmed by ADAM 17 (a disintegrin and metalloproteinase 17) that competes for the same receptor. Cleavage by TTSPs activates ACE2 receptor for binding, whereas ADAM17 releases extracellular fragment called soluble ACE2 (sACE2). Structural studies of both ACE2 and S-protein have found critical sites involved in binding mechanism. In addition to studies on structural motifs, few single-nucleotide polymorphism (SNPs) studies have been done to find an association between genetic variants and SARS susceptibility. Though no association was found in those reports, but seeing the non-reproducibility of SNP studies among different ethnic groups, screening of ACE2 SNPs in individual population can be undertaken.

Short conclusion

Thus, screening for novel SNPs focussing on recently identified critical regions of ACE2 can be targeted to monitor susceptibility towards coronavirus disease 2019 (COVID-19).

A novel forensic panel of 186-plex SNPs and 123-plex STR loci based on massively parallel sequencing

The MiSeq® FGX Forensic system and the HID-Ion AmpliSeq Panel were previously developed for massively parallel sequencing (MPS) for forensic casework. Among the three major sequencing platforms, BGISEQ-500^TM, which is based on multiple PCRs, is still lacking in forensics. Here, a novel forensic panel was constructed to detect 186 single-nucleotide polymorphisms (SNPs) and 123 short tandem repeats (STRs) with MPS technology on the BGISEQ-500™ platform. First, the library preparation, sequencing process, and data analysis were performed, focusing on the average depth of coverage and heterozygote balance. We calculated the allelic frequencies and forensic parameters of STR and SNP loci in 73 unrelated Chinese Han individuals. In addition, performance was evaluated with accuracy, uniformity, sensitivity, PCR inhibitor, repeatability and reproducibility, mixtures, degraded samples, case-type samples, and pedigree analyses. The results showed that 100% accurate and concordant genotypes can be obtained, and the loci with an abundance in the interquartile range accounted for 92.90% of the total, suggesting reliable uniformity in this panel. We obtained a locus detection rate that was higher than 98.78% from 78 pg of input DNA, and the optimal amount was 1.25-10 ng. The maximum concentrations of hematin and humic acid were 200 and 100 μM, respectively (the ratios of detected loci were 96.52% and 92.41%), in this panel. As a mixture, compared with those of SNPs, minor-contributor alleles of STRs could be detected at higher levels. For the degraded sample, the ratio of detected loci was 98.41%, and most profiles from case-type samples were not significantly different in abundance in our studies. As a whole, this panel showed high-performance, reliable, robust, repeatable, and reproducible results, which are sufficient for paternity testing, individual identification, and use for potentially degraded samples in forensic science.

Interaction of two functional genetic variants LOXL1 rs1048661 and VEGFA rs3025039 on the risk of age-related macular degeneration in Chinese women

Background

Cumulative evidence indicates that LOXL1 and VEGF-a play important roles in extracellular matrix formation and angiogenesis, respectively. The disorder of extracellular matrix and angiogenesis are the key factors of pathogenesis of age-related macular degeneration (AMD). We hypothesized that rs1048661 (T>G) in the LOXL1 gene and rs3025039 (C>T) in the VEGFA gene might be associated with risk of AMD.

Methods

A total of 533 unrelated Chinese subjects, 286 cases (247 with early AMD and 39 with late neovascular AMD) and 247 controls, were included in the study. The gene sequences of LOXL1 rs1048661 and VEGFA rs3025039 were amplified by polymerase chain reaction and genotyped. Interaction between rs1048661 and rs3025039 on AMD risk was also assessed.

Results

LOXL1 rs1048661 but not VEGFA rs3025039 was associated with a significantly increased risk of AMD. The adjusted odds ratio was 1.6 (95% CI, 1.1-2.5) for rs1048661 TT + GT genotype compared with GG homozygotes in the dominant model analysis. Moreover, there was a significant gene-gene interaction between these two polymorphic loci. In VEGFA rs3025039 CC + CT genotype which indicated sufficient expression of VEGF-a, LOXL1 rs1048661 had odds ratios of 1.7 (95% CI, 1.1-2.7) for early AMD and 3.6 (95% CI, 1.1-12.3) for late neovascular AMD in the dominant model analysis. However, LOXL1 rs1048661 did not confer the risk of AMD in subjects harboring VEGFA rs3025039 TT genotype which indicated decreased expression of VEGF-a.

Conclusions

Our findings suggest that LOXL1 rs1048661 (T>G) may be involved in the risk of AMD. In addition, LOXL1 rs1048661 and VEGFA rs3025039 interacted to confer the development of AMD, especially for late-stage neovascular AMD. Our data need to be further validated.

Genetic associations and serum paraoxonase levels with atherosclerosis in western Iranian patients

The oxidative modification of low-density lipoprotein (LDL) in the arterial wall plays a pivotal role in the initiation and progression of atherosclerosis which is a complex and progressive disorder. Paraoxonase1 (PON1), which is required for lipid metabolism, is believed to protect LDL from oxidation. The relationship between PON1 gene Leusin55Methionin (L55M) and Glutamine192Arginine (Q192R) polymorphisms in western Iranians with atherosclerosis and its association with enzyme activity and oxidized low-density lipoprotein (oxLDL) were examined in the present study. In this study, blood specimens were collected from 145 healthy individuals and 154 patients with atherosclerosis proven by angiography referred to Shahid Madani Hospital, Khorramabad, Iran. Genomic deoxy ribonucleic acid (DNA) was extracted from whole blood. For all the subjects, restriction fragment length polymorphism-polymerase chain reaction (RFLP-PCR) was carried out for the detection of L55M and Q192R polymorphisms. PON1 enzyme activity and the level of oxLDL were also evaluated. There was a 3.114-fold increase in the risk of developing atherosclerosis in the subjects presenting the PON1L55M, MM genotype compared to those with the LL genotype (OR 3.114; 95% CI 1.412-6.870). PON1Q192R polymorphism in the PON1 gene was not associated with atherosclerosis. Patients with atherosclerosis had significantly higher oxLDL and reduced PON1 enzyme activity (P < 0.05) compared to the controls. There was no association between the type of genotype, enzyme activity, and oxLDL level. It has been concluded that PON1L55M polymorphism and MM genotype are associated with an increased risk of coronary artery disease (CAD) in Iranian patients with atherosclerosis. We did not find any relationship between PON1Q192R polymorphism and atherosclerosis.

Age-dependent association of polymorphisms in the promoter and 5'-untranslated region of the norepinephrine transporter gene with generalized anxiety disorder

BACKGROUND

Norepinephrine transporter (NET), which regulates synaptic norepinephrine for noradrenergic signaling, is involved in the pathogenesis of anxiety, while expression of the NET gene differs at different ages. Here, we examine whether genetic variants in the NET gene are associated, in an age-specific manner, with increased risk of generalized anxiety disorder (GAD), one of the most disabling anxiety disorders.

METHODS

Three common single-nucleotide polymorphisms (SNPs) in the promoter (rs168924: A/G; rs2242446: T/C) and 5'-untranslated region (5'-UTR) (rs2397771: G/C) of the NET gene were genotyped in 2,317 Han-Chinese participants (791 GAD patients and 1,526 controls; age: 20-65). Potential confounding factors, such as gender, stress levels and psychiatric comorbidities, were included as covariates.

RESULTS

An interaction between age and NET genotypes and haplotypes was found for the risk of GAD. In the younger participants, rs168924 minor allele G homozygotes had the lowest incidence of GAD; however, older subjects displayed an inverse pattern, with homozygous G/G carriers presenting the highest prevalence of GAD. Additionally, younger individuals carrying 2 copies of the GGT haplotype composed of rs2397771-rs168924-rs2242446 had the lowest rate of GAD. However, those with 2 copies of the same haplotype exhibited the highest risk of GAD in the older groups.

LIMITATIONS

Only 3 common SNPs in the promoter and 5'-UTR of the NET gene were analyzed.

CONCLUSIONS

Our findings are the first to demonstrate that potentially functional SNPs in the NET promoter and 5'-UTR are associated with an increased risk of GAD, and that such associations are determined in an age-specific way.

Lack of association between interleukin-22 gene polymorphisms and cancer risk: a case-control study and a meta-analysis

BACKGROUND

Interleukin-22 (IL22) has been implicated in inflammation and tumorigenesis. The association between IL22 gene polymorphisms and cancer risk has been widely explored. However, the limited sample sizes of previous studies may produce inadequate statistical power and conflicting results, which calls for further investigations. In this study, we recruited a total of 1490 cancer patients (480 liver cancer patients, 550 lung cancer patients, and 460 gastric cancer patients) and 800 normal controls to explore the associations between IL22 gene polymorphisms (rs1179251, rs2227485, rs2227511, and rs2227473) and cancer risk.

METHOD

The genotyping was performed with polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) and Sanger sequencing.

RESULTS

Our results showed that none of the four IL22 gene polymorphisms was associated with the risk of liver, lung or gastric cancer in Hubei Han Chinese population. To improve the statistical strength, a meta-analysis was further conducted. The results further confirmed our present findings and showed that rs1179251, rs2227485, and rs2227473 were not associated with cancer risk in total or stratified analysis.

CONCLUSION

Consequently, the rs1179251, rs2227485, rs2227511, and rs2227473 polymorphisms may not be associated with cancer risk. However, further investigations using larger samples in different ethnic populations are required.

Genetic association of FKBP5 with PTSD in US service members deployed to Iraq and Afghanistan

Post-traumatic stress disorder (PTSD) is a debilitating mental disorder with a prevalence of more than 7% in the US population and 12% in the military. An interaction of childhood trauma with FKBP5 (a glucocorticoid-regulated immunophilin) has been reported to be associated with PTSD in the general population. However, there are few reports on the association of FKBP5 with PTSD, particularly in important high-risk population such as the military. Here, we examined the association between four single-nucleotide polymorphisms (SNPs; rs3800373, rs9296158, rs1360780, rs9470080) covering the FKBP5 gene and probable PTSD in US service members deployed to Iraq and Afghanistan, a high-risk military population (n = 3890) (Hines et al., 2014). We found that probable PTSD subjects were significantly more likely to carry the A-allele of rs3800373, G-allele of rs9296158, C-allele of rs1360780, and C-allele of rs9470080. Furthermore, the four SNPs were in one block of strong pairwise linkage disequilibrium (r = 0.91-0.96). Within the block there were two major haplotypes of CATT and AGCC (rs3800373-rs9296158-rs1360780-rs9470080) that account for 99% of haplotype diversity. The distribution of the AGCC haplotype was significantly higher in probable PTSD subjects compared to non-PTSD (p<.05). The diplotype-based analysis indicated that the AGCC carriers tended to be probable PTSD. In this study, we demonstrated the association between FKBP5 and probable PTSD in US service members deployed to Iraq and Afghanistan, indicating that FKBP5 might be a risk factor for PTSD.

The polygenic nature of mild-to-moderate hypertriglyceridemia

BACKGROUND

Patients with mild-to-moderate hypertriglyceridemia (HTG) are thought to share specific genetic susceptibility factors that are also present in patients with severe HTG, but no data have been reported on this issue.

OBJECTIVE

The objective of this study was to characterize genetic profiles of patients with mild-to-moderate HTG and compare them to patients with severe HTG.

METHODS

DNA from patients with mild-to-moderate HTG was sequenced using our targeted sequencing panel, "LipidSeq". For each patient, we assessed 1) rare variants disrupting five TG metabolism genes and 2) the accumulation of 16 common single-nucleotide polymorphisms (SNPs) using a polygenic risk score. The genetic profiles for these patients were then compared with normolipidemic controls from the 1000 Genomes Project and with patients with severe HTG.

RESULTS

Across 134 patients with mild-to-moderate HTG, 9.0% carried heterozygous rare variants and 26.9% had an excess accumulation of common SNPs. Patients with mild-to-moderate HTG were 2.38 times (95% CI [1.13-4.99]; P = .021) more likely to carry a rare variant and 3.26 times (95% CI [2.02-5.26]; P < .0001) more likely to have an extreme polygenic risk score compared with the 1000 Genomes Project. In addition, patients with severe HTG were 1.86 times (95% CI [0.98-3.51]; P = .032) more likely to carry a rare variant and 1.63 times (95% CI [1.07-2.48]; P = .013) more likely to have an extreme polygenic risk score than patients with mild-to-moderate HTG.

CONCLUSIONS

We report an increased prevalence of genetic determinants in patients with an increased severity of the HTG phenotype when considering either rare variants disrupting TG metabolism genes or an excess accumulation of common SNPs. As well, the findings confirm that the most prevalent genetic contributor to HTG, regardless of severity, is polygenic SNP accumulation.

Relationship between HLA-DPA1 genetic polymorphism and anembryonic pregnancy

BACKGROUND

Human leukocyte antigen (HLA)-DP is an HLA class II molecule. Overexpression of HLA class II molecules in placental trophoblast cells may induce pregnancy loss. However, the association between HLA-DP and pregnancy loss remains unclear. HLA-DPA1 is an HLA-DP peptide chain. The objective of this study was to assess the association between HLA-DPA1 genetic polymorphism and anembryonic pregnancy, a type of early pregnancy loss, in the Chinese population.

METHODS

A case-control study was designed to compare the frequencies of HLA-DPA1 gene polymorphisms in an anembryonic pregnancy group and a control group. Sixty-eight cases and 122 controls were recruited. Statistical analysis was performed to assess the correlation between single-nucleotide polymorphisms (SNPs) and anembryonic pregnancy susceptibility. MassARRAY high-throughput DNA analysis was used to analyze 19 HLA-DPA1 SNPs. To explore how HLA-DPA1 polymorphism could affect anembryonic pregnancy, HLA-DPA1 serum levels were analyzed by ELISA.

RESULTS

Homozygous typing of rs1431403 (CC and TT) significantly increased the risk of anembryonic pregnancy in the case group (ORCC  = 3.13, 95% CI: 1.50-6.53; ORTT  = 2.96, 95% CI: 1.31-6.66; ORCC+TT  = 3.06, 95% CI: 1.62-5.78). In samples with high HLA-DPA1 levels (≥1,500 pg/ml), the homozygous rs1431403 genotypes (nCC  = 21, 43.8%; nTT  = 20, 57.1%) were observed more frequently than were heterozygous genotypes.

CONCLUSION

HLA-DPA1 rs1431403 may be a risk factor for anembryonic pregnancy in the Chinese population. Homozygous rs1431403 genotypes (CC and TT) may increase the risk of anembryonic pregnancy by aberrantly increasing the HLA-DPA1 levels.

Genetic variant rs16944 in IL1B gene is a risk factor for early-onset sepsis susceptibility and outcome in preterm infants

OBJECTIVE

Interleukin-1-B (IL1B) is a proinflammatory cytokine that plays an important role in sepsis. The aim of this study was to evaluate the relationships between IL1B-511G/A polymorphism and susceptibility and outcome of early-onset sepsis (EOS) in preterm infants.

METHODS

DNA was extracted from the buccal swabs of 471 (285 with EOS and 186 control) preterm infants. Genotypes of rs16944 polymorphism were determined with real-time PCR method.

RESULTS

We found statistically significant higher frequency of IL1B-511AA genotype in EOS group than in control group (p = 0.012). Also, IL1B-511AA genotype is statistically significantly more frequent in patients with lethal EOS outcome (p = 0.011).

CONCLUSION

Genotype IL1B-511AA was associated with susceptibility to EOS and it is a significant predictor of lethal outcome in preterm infants with EOS.

Genetic polymorphisms in part of intron 7 and exon 8 of HSP90AA1 gene and its association with heat tolerance traits in two exotic layer chicken strains

Chickens, like other vertebrates, react to stress conditions through their cultured cells by expressing heat shock proteins (HSPs). Genetic association of single-nucleotide polymorphisms (SNPs) in HSPs with desirable traits will reveal their importance as potential genetic markers. Blood samples were collected from 50 birds per strain for DNA extraction, polymerase chain reaction amplification and sequencing of the HSP90AA1 gene. SNPs were detected using Codoncode Aligner. Association of each SNP with heat tolerance traits was analysed using generalized linear model procedure of SAS. A total of seven SNPs were detected, four SNPs; A7T, A160T, T223A and C134T were detected in part of intron 7 to exon 8 of HSP90AA1 gene of BRD while three A160T, T223A and C134T were detected in HYL. SNPC134T, a synonymous variant, was detected in exon 8. Only SNPA7T was in Hardy-Weinberg equilibrium (X² = 0.03) but had no association with the traits measured. Polymorphic information content calculated showed SNPA160T to be moderately polymorphic; other SNPs were lowly polymorphic. Heterozygosity for SNPs-A160T and T223A of BRD showed moderate genetic variation while the other SNPs and those in HYL recorded low genetic variation. The study concluded that the SNPs detected were majorly lowly polymorphic and also the SNP locus A7T in intron 7 of HSP90AA1 of BRD had no genetic association with heat tolerance traits.

Genomic organization of the molt-inhibiting hormone gene in the red swamp crayfish Procambarus clarkii and characterization of single-nucleotide polymorphisms associated with growth

Molt-inhibiting hormone (MIH), a neuropeptide synthesized in the eyestalk in crustaceans, is mainly responsible for the molting by negatively controlling the ecdysteroids secretion. Although there are several reports of the isolation and protein sequencing of MIH in the red swamp crayfish, little is known about the nucleotide sequence and gene organization of this neuropeptide, even less about the association of MIH polymorphisms and growth traits. Here, a 1237 bp full-length MIH cDNA was obtained from the crayfish eyestalk, which encodes a putative protein of 106 amino acids, with a 191 bp 5'-UTR and a 728 bp 3'-UTR. The MIH genomic DNA sequence is 4205 bp in length, which includes three exons interrupted by two introns, and a 929 bp 5'-flanking region. Potential transcription initiation site and transcription factor binding sites were identified in the 5'-flanking region, implying a potential role in transcriptional regulation. Seventeen SNPs in the 5'-flanking region and 3'-UTR were identified, and the associations between these SNPs and growth traits were evaluated with a two-stage design. A SNPs g. -12C > G that showed a significant association with body weight was identified. Individuals with GG genotype had a significantly higher body weight than those with CC genotype (43.98 ± 9.82 g vs. 34.27 ± 6.87 g; P ﹤ 0.001), indicating a beneficial effect of the G allele on the growth of red swamp crayfish. The obtained MIH gene, as well as the identified SNPs, may serve as targets for molecular marker-aided selection in growth improvement of the red swamp crayfish in future studies.

Correlation of EYS polymorphisms with lumbar disc herniation risk among Han Chinese population

BACKGROUND

Lumbar disc herniation (LDH) is a common spinal disease in clinical practice. Once lumbar disc herniation occurs, it seriously reduces patient's quality of life. The EYS (eyes shut homolog) was discovered in recent years and it may be related to lumbar disc herniation. So we conducted a case-control study to explore the relationship between EYS polymorphism and lumbar disc herniation risk.

METHODS

We selected 5 single-nucleotide polymorphisms (SNPs) of EYS gene in a case-control study with 508 cases and 508 healthy controls to evaluate the relatedness by using genetic model, haplotype, and stratification analysis.

RESULTS

We found that the minor alleles of rs62413038 (OR = 1.21, 95%CI: 1.01-1.43, p = .036) and rs9450607 (OR = 1.26, 95% CI: 1.05-1.53, p = .016) were associated with an increased risk of lumbar disc herniation in the allelic model analysis. In the genotypic model analysis, rs62413038 displayed a significantly increased risk of lumbar disc herniation in log-additive models (OR = 1.20, 95% CI: 1.01-1.43, p = .039). While the rs9450607 was also obviously associated with an increased lumbar disc herniation risk in recessive (OR = 1.98, 95% CI: 1.24-3.13, p = .004) and log-additive models (OR = 1.27, 95% CI: 1.05-1.55, p = .014). In addition, in the haplotype analyses of the SNPs, we found that the "CGGA" haplotype of rs1482456, rs9342097, rs9450607, and rs7757884 was associated with lumbar disc herniation. (OR = 0.52, 95% CI: 0.30-0.89, p = .017).

CONCLUSION

These results suggest that EYS polymorphism may be associated with lumbar disc herniation among Han Chinese population. It also opens up a new exploration direction for the etiology of lumbar disc herniation.

Functional characterization of SLC26A3 c.392C>G (p.P131R) mutation in intestinal barrier function using CRISPR/CAS9-created cell models

Background

Congenital chloride diarrhea (CCD) in a newborn is a rare autosomal recessive disorder with life-threatening complications, requiring early diagnostics and treatment to prevent severe dehydration and infant mortality. SLC26A3 rs386833481 (c.392C>G; p.P131R) gene polymorphism is an important genetic determinant of CCD. Here, we report the influence of the non-synonymous SLC26A3 variant rs386833481 gene polymorphism on the function of the epithelial barrier and the potential mechanisms of these effects.

Results

We found that P131R-SLC26A3 increased dysfunction of the epithelial barrier compared with wild type SLC26A3 in human colonic Caco-2 and mouse colonic CMT-93 cells. When P131R-SLC26A3 was subsequently reverted to wild type, the epithelial barrier function was restored similar to wild type cells. Further study demonstrated that variant P131R-SLC26A3 disrupts function of epithelial barrier through two distinct molecular mechanisms: (a) decreasing SLC26A3 expression through a ubiquitination pathway and (b) disrupting a key interaction with its partner ZO-1/CFTR, thereby increasing the epithelial permeability.

Conclusion

Our study provides an important insight of SLC26A3 SNPs in the regulation of the epithelial permeability and indicates that SLC26A3 rs386833481 is likely a causative mutation in the dysfunction of epithelial barrier of CCD, and correction of this SNP or increasing SLC26A3 function could be therapeutically beneficial for chronic diarrhea diseases.

Comparative transcriptomics analysis uncovers alternative splicing events and molecular markers in cabbage (Brassica oleracea L.)

Alternative splicing (AS) events were identified and verified in cabbage by comparative transcriptome analysis. The corresponding markers were developed and the germplasm resources were identified. Alternative splicing (AS) is a central regulatory mechanism that greatly contributes to plant gene expression and transcriptome diversity. A large body of evidence has shown that AS complexity is relevant for plant development, evolution, complexity, and adaptation. Both insertion/deletion (InDel) and single nucleotide polymorphism (SNP) are typically co-dominant inheritance markers and have abundant polymorphisms. These have been widely used for marker-assisted selection, genetic mapping, and germplasm identification in plants. However, little is known about the molecular mechanisms underlying AS events and the development of markers including SNP and InDel from the cabbage transcriptome. In this study, three cabbage transcriptome datasets were collected and aligned to the cabbage reference genome to analyze AS events and marker development. 31,524 AS events were identified from three cabbage genotypes, accounting for 20.8% of the total cabbage genes. Alternative 3' splice site donor (A3SS) was the most frequent type of the four main AS events in cabbage. 70,475 InDels and 706,269 SNPs were identified with average frequencies of 1 InDel/6.9 kb and 1 SNP/0.7 kb, respectively. 71,942 potential SSRs were identified in 53,129 assembled unigenes with a density of 1 SSR/6.8 kb. The ratio of SNPs with synonymous/non-synonymous mutations was 1:0.65. 142 InDels and 36 SNPs were randomly selected and validated via Sanger sequencing and polymorphism was found among 66.2% of the InDels and 78.6% of the SNPs. Furthermore, 35 informative InDel markers were successfully used for genetic diversity analysis on 36 cabbage accessions. These results facilitate understanding of the molecular regulation mechanism underlying AS events in cabbage. They also provide molecular marker resource data for genetic mapping construction and germplasm identification, and facilitate the genetic improvement of cabbage via breeding.

A schizophrenia associated CMYA5 allele displays differential binding with desmin

CMYA5 is a candidate gene for schizophrenia because of the genetic association of variant rs10043986 (C > T) to this severe mental disorder. Studies of CMYA5 and its gene product, myospryn, in the brain and neuronal cells have not been previously reported. The SNP rs10043986 changes the 4,063rd amino acid from Pro to Leu, which is likely to alter protein function. To understand its potential role in the brain, we examined the neuronal expression of myospryn and its binding partner, desmin, an intermediate filament (IF) protein, and investigated how the two alleles of myospryn affect its binding to desmin. Myospryn and desmin are shown to be expressed in the brain and myospryn is shown to localize to the cytoplasm and nucleus of myoblast, neuroblastoma, and glioblastoma cell lines. Peripherin and vimentin, known brain IF proteins, have high protein similarity to desmin but were found not to interact with myospryn using yeast two-hybrid (Y2H). Using a quantitative Y2H assay and surface plasmon resonance, the T allele (Leu) of rs10043986 was found to have stronger binding to desmin than the C allele (Pro). Based on findings described in this report, we hypothesize that the interaction between myospryn to IF provides structural support and efficient rearrangement of the cytoskeleton network during early neuritogenesis.

PEX10, SIRPA-SIRPG, and SOX5 gene polymorphisms are strongly associated with nonobstructive azoospermia susceptibility

PURPOSE

Male infertility is a multifactorial syndrome encompassing a wide variety of disorders. A previous Chinese genome-wide single-nucleotide polymorphism (SNP) association studies have identified four SNPs (rs12097821 in PRMT6 gene, rs2477686 in PEX10 gene, rs6080550 in SIRPA-SIRPG, and rs10842262 in SOX5 gene) as being significantly associated with risk factors for nonobstructive azoospermia (NOA). However, the results were not fully repeated in later studies, which calls for further investigations.

METHODS

We here performed a case-control study in a central Chinese population to explore the association between the four SNPs and male infertility, which included 631 infertile men (NOA and oligozoospermia) and 720 healthy fertile men. The genotyping was performed using the polymerase chain reaction-restriction fragment length polymorphism and confirmed by sequencing.

RESULTS

The results showed that rs12097821 and rs10842262 were strongly associated with the risk of NOA but not total male infertility or oligozoospermia, while rs2477686 and rs6080550 were not associated with the risk of total male infertility, NOA, or oligozoospermia. To improve the statistical strength, a meta-analysis was conducted. The results suggested that rs2477686, rs6080550, and rs10842262 were significantly associated with male infertility, especially with NOA, while rs12097821 was only found to be associated with total male infertility.

CONCLUSIONS

Collectively, the rs2477686, rs6080550, and rs10842262 may indeed be the genetic risk factors for NOA, which requires further investigation using larger independent sets of samples in different ethnic populations.

Mutational analysis of theFAM175A gene in patients with premature ovarian insufficiency

RESEARCH QUESTION

The family with sequence similarity 175 member A gene (FAM175A; also known as ABRAXAS1, CCDC98 and ABRA1), a member of the DNA repair family, contributes to the BRCA1 (BRCA1 DNA repair associated)-dependent DNA damage response and is associated with age at natural menopause. However, it remains poorly understood whether sequence variants in FAM175A are causative for premature ovarian insufficiency (POI). The aim of this study was to investigate whether mutations in the gene FAM175A were present in patients with POI.

DESIGN

A total of 400 women with idiopathic POI and 498 control women with regular menstruation (306 age-matched women and 192 women over 40 years old) were recruited. After Sanger sequencing of FAM175A, functional experiments were carried out to explore the deleterious effects of the identified variation. DNA damage was subsequently induced by mitomycin C (MMC), and DNA repair capacity and G2-M checkpoint activation were evaluated by examining the phosphorylation level of H2AX (H2A histone family, member X) and the percentage of mitotic cells, respectively.

RESULTS

One rare single-nucleotide polymorphism, rs755187051 in gene FAM175A, c.C727G (p.L243V), was identified in two patients but absent in the 498 controls. The functional experiments demonstrated that overexpression of variant p.L243V in HeLa cells resulted in a similar sensitivity to MMC-induced damage compared with cells transfected with wild-type FAM175A. Moreover, after treatment with MMC, there were no differences in DNA repair capacity and G2-M checkpoint activation between the mutant and wild-type genes.

CONCLUSION

Our results suggest that the p.L243V variant of FAM175A may not be causative for POI. The contribution of FAM175A to POI needs further exploration.

Association between expression levels and growth trait-related SNPs located in promoters of the MC4R and MSTN genes in Spinibarbus hollandi

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.

Evaluating the breast cancer predisposition role of rare variants in genes associated with low-penetrance breast cancer risk SNPs

Background

Genome-wide association studies (GWASs) have identified numerous single-nucleotide polymorphisms (SNPs) associated with small increases in breast cancer risk. Studies to date suggest that some SNPs alter the expression of the associated genes, which potentially mediates risk modification. On this basis, we hypothesised that some of these genes may be enriched for rare coding variants associated with a higher breast cancer risk.

Methods

The coding regions and exon-intron boundaries of 56 genes that have either been proposed by GWASs to be the regulatory targets of the SNPs and/or located < 500 kb from the risk SNPs were sequenced in index cases from 1043 familial breast cancer families that previously had negative test results for BRCA1 and BRCA2 mutations and 944 population-matched cancer-free control participants from an Australian population. Rare (minor allele frequency ≤ 0.001 in the Exome Aggregation Consortium and Exome Variant Server databases) loss-of-function (LoF) and missense variants were studied.

Results

LoF variants were rare in both the cases and control participants across all the candidate genes, with only 38 different LoF variants observed in a total of 39 carriers. For the majority of genes (n = 36), no LoF variants were detected in either the case or control cohorts. No individual gene showed a significant excess of LoF or missense variants in the cases compared with control participants. Among all candidate genes as a group, the total number of carriers with LoF variants was higher in the cases than in the control participants (26 cases and 13 control participants), as was the total number of carriers with missense variants (406 versus 353), but neither reached statistical significance (p = 0.077 and p = 0.512, respectively). The genes contributing most of the excess of LoF variants in the cases included TET2, NRIP1, RAD51B and SNX32 (12 cases versus 2 control participants), whereas ZNF283 and CASP8 contributed largely to the excess of missense variants (25 cases versus 8 control participants).

Conclusions

Our data suggest that rare LoF and missense variants in genes associated with low-penetrance breast cancer risk SNPs may contribute some additional risk, but as a group these genes are unlikely to be major contributors to breast cancer heritability.

Electronic supplementary material

The online version of this article (doi:10.1186/s13058-017-0929-z) contains supplementary material, which is available to authorized users.

Accelerating public sector rice breeding with high-density KASP markers derived from whole genome sequencing of indica rice

Few public sector rice breeders have the capacity to use NGS-derived markers in their breeding programmes despite rapidly expanding repositories of rice genome sequence data. They rely on > 18,000 mapped microsatellites (SSRs) for marker-assisted selection (MAS) using gel analysis. Lack of knowledge about target SNP and InDel variant loci has hampered the uptake by many breeders of Kompetitive allele-specific PCR (KASP), a proprietary technology of LGC genomics that can distinguish alleles at variant loci. KASP is a cost-effective single-step genotyping technology, cheaper than SSRs and more flexible than genotyping by sequencing (GBS) or array-based genotyping when used in selection programmes. Before this study, there were 2015 rice KASP marker loci in the public domain, mainly identified by array-based screening, leaving large proportions of the rice genome with no KASP coverage. Here we have addressed the urgent need for a wide choice of appropriate rice KASP assays and demonstrated that NGS can detect many more KASP to give full genome coverage. Through re-sequencing of nine indica rice breeding lines or released varieties, this study has identified 2.5 million variant sites. Stringent filtering of variants generated 1.3 million potential KASP assay designs, including 92,500 potential functional markers. This strategy delivers a 650-fold increase in potential selectable KASP markers at a density of 3.1 per 1 kb in the indica crosses analysed and 377,178 polymorphic KASP design sites on average per cross. This knowledge is available to breeders and has been utilised to improve the efficiency of public sector breeding in Nepal, enabling identification of polymorphic KASP at any region or quantitative trait loci in relevant crosses. Validation of 39 new KASP was carried out by genotyping progeny from a range of crosses to show that they detected segregating alleles. The new KASP have replaced SSRs to aid trait selection during marker-assisted backcrossing in these crosses, where target traits include rice blast and BLB resistance loci. Furthermore, we provide the software for plant breeders to generate KASP designs from their own datasets.

Distinguishing the disease-associated SNPs based on composition frequency analysis

Single-nucleotide polymorphism (SNP) is a basical variation in genome. When SNPs occur at the binding sites of microRNA, they can influence the binding efficiency, cause a fluctuation of the mRNA in vivo, and thus arouse posttranscriptional level abnormality. Therefore, SNP has a strong correlation with diseases. Although enormous SNPs have been experimentally identified, only a tiny proportion of them are truly disease-associated SNPs (dSNPs) that relate to microRNA modification and then are involved in disease causing process. Therefore, it is important to distinguish dSNPs from the usual SNPs. Analysis here shows that composition is different between sequence segments centered by dSNP and SNP. Inspired by the composition, transition and distribution features which are meaningful and effective in characterizing proteins' sequence information, we improved and applied it to represent the frequency and physicochemical properties of a gene sequence. Binary encoding scheme was also used for further labelling four nucleic acids (A, T, C, and G). First, clustering analysis was performed to gain reasonable negative samples. Then, optimization tests were implemented on different ratios of positive vs negative samples and different feature subsets retrieved by evaluation method of F score. The optimal model constructed by random forest achieves an accuracy of more than 90% on the testing data set. Moreover, the promising results of the external validation also demonstrate the practical applicability of our method. Finally, principal component analysis on the features indicates that all features in our method gain the gross contribution to the prediction model.

A common variant of MAF/c-MAF, transcriptional factor gene in the kidney, is associated with gout susceptibility

Gout is a multifactorial disease characterized by acute inflammatory arthritis, and it is caused as a consequence of hyperuricemia. A recent meta-analysis of genome-wide association studies has newly identified the relationship between serum uric acid (SUA) levels and rs889472, a single nucleotide polymorphism of musculoaponeurotic fibrosarcoma oncogene (MAF/c-MAF). However, it remained unclear whether rs889472 is associated with gout susceptibility. In the present study, we investigate the association between c-MAF rs889472 and gout in Japanese male population. We genotyped 625 male patients who were clinically diagnosed as gout and 1221 male control subjects without hyperuricemia or a history of gout by TaqMan method. As a result, the major allele (C), which reportedly increases SUA levels, had a higher frequency in the gout cases (58.8%) than in the controls (55.0%). A logistic regression analysis showed a significant association between rs889472 and gout (p = 0.029, odds ratio = 1.17; 95% confidence interval 1.02-1.34). C-MAF is reported as a pivotal transcriptional factor in the development and differentiation of renal proximal tubular cells. Because urate is mainly regulated in renal proximal tubular cells, c-MAF may have an important role in urate regulation in the kidney and influence not only SUA but also gout susceptibility. Our finding shows that rs889472 of c-MAF is associated with gout susceptibility.

The haplotype of UBE2L3 gene is associated with Hashimoto's thyroiditis in a Chinese Han population

BACKGROUND

The ubiquitin conjugating enzyme E2L3 (UBE2L3) gene is associated with susceptibility to many autoimmune diseases. The aim of this study was to investigate the association between UBE2L3 gene and autoimmune thyroid diseases (AITDs) and their clinical phenotypes.

METHODS

We genotyped five single-nucleotide polymorphisms (SNPs) rs131654, rs5754217, rs2298428, rs140489 and rs5998672 of UBE2L3 gene in case groups including 1028 patients with AITDs [676 cases of Graves' disease (GD) and 352 cases of Hashimoto's thyroiditis (HT)] and control group including 897 healthy individuals. The genotyping was performed with the method of polymerase chain reaction-ligase detection reaction (PCR-LDR).

RESULTS

The frequencies of allele and genotype of five SNPs in gene UBE2L3 showed no statistically significant difference between case groups and control group, respectively. Moreover, no significant differences in frequencies of allele and genotype of five SNPs of the gene were found between clinical subphenotypes of AITDs and control group. Such subphenotypes included GD, HT, and thyroid associated ophthalmopathy (TAO). The negative results were also found in the frequency of other haplotypes of the gene except the haplotype of TCGGC, which was significantly higher in HT group than in control group (P = 0.031, OR = 1.441).

CONCLUSIONS

The present findings indicate that TCGGC haplotype is associated with an increased risk of HT and UBE2L3 gene is likely to be a susceptibility factor to HT in a Chinese Han population.

An integrated analysis tool for analyzing hybridization intensities and genotypes using new-generation population-optimized human arrays

BACKGROUND

Affymetrix Axiom single nucleotide polymorphism (SNP) arrays provide a cost-effective, high-density, and high-throughput genotyping solution for population-optimized analyses. However, no public software is available for the integrated genomic analysis of hybridization intensities and genotypes for this new-generation population-optimized genotyping platform.

RESULTS

A set of statistical methods was developed for an integrated analysis of allele frequency (AF), allelic imbalance (AI), loss of heterozygosity (LOH), long contiguous stretch of homozygosity (LCSH), and copy number variation or alteration (CNV/CNA) on the basis of SNP probe hybridization intensities and genotypes. This study analyzed 3,236 samples that were genotyped using different SNP platforms. The proposed AF adjustment method considerably increased the accuracy of AF estimation. The proposed quick circular binary segmentation algorithm for segmenting copy number reduced the computation time of the original segmentation method by 30-67 %. The proposed CNV/CNA detection, which integrates AI and LOH/LCSH detection, had a promising true positive rate and well-controlled false positive rate in simulation studies. Moreover, our real-time quantitative polymerase chain reaction experiments successfully validated the CNVs/CNAs that were identified in the Axiom data analyses using the proposed methods; some of the validated CNVs/CNAs were not detected in the Affymetrix Array 6.0 data analysis using the Affymetrix Genotyping Console. All the analysis functions are packaged into the ALICE (AF/LOH/LCSH/AI/CNV/CNA Enterprise) software.

CONCLUSIONS

ALICE and the used genomic reference databases, which can be downloaded from http://hcyang.stat.sinica.edu.tw/software/ALICE.html , are useful resources for analyzing genomic data from the Axiom and other SNP arrays.

Single-Nucleotide Polymorphism to Associate Cancer Risk

Genetic heterogeneity explains variation in predisposition for cancer. Whole-genome analysis allows risk to be quantified, giving better targeted screening and quantification of the personalized risk posed by environmental factors. Array-based approaches to whole-genome analysis are rapidly being overtaken by next-generation sequencing (NGS). In this review the different platforms currently available for NGS are compared and the opportunities and risks of this approach are discussed: including the informatics packages required and the ethical issues. Methods applicable to the personal genome machine (PGM) are given as an example of workflows.

Targeting the Plasmodium vivax equilibrative nucleoside transporter 1 (PvENT1) for antimalarial drug development

Infection with Plasmodium falciparum and vivax cause most cases of malaria. Emerging resistance to current antimalarial medications makes new drug development imperative. Ideally a new antimalarial drug should treat both falciparum and vivax malaria. Because malaria parasites are purine auxotrophic, they rely on purines imported from the host erythrocyte via Equilibrative Nucleoside Transporters (ENTs). Thus, the purine import transporters represent a potential target for antimalarial drug development. For falciparum parasites the primary purine transporter is the P. falciparum Equilibrative Nucleoside Transporter Type 1 (PfENT1). Recently we identified potent PfENT1 inhibitors with nanomolar IC₅₀ values using a robust, yeast-based high throughput screening assay. In the current work we characterized the Plasmodium vivax ENT1 (PvENT1) homologue and its sensitivity to the PfENT1 inhibitors. We expressed a yeast codon-optimized PvENT1 gene in Saccharomyces cerevisiae. PvENT1-expressing yeast imported both purines ([³H]adenosine) and pyrimidines ([³H]uridine), whereas wild type (fui1Δ) yeast did not. Based on radiolabel substrate uptake inhibition experiments, inosine had the lowest IC₅₀ (3.8 μM), compared to guanosine (14.9 μM) and adenosine (142 μM). For pyrimidines, thymidine had an IC₅₀ of 183 μM (vs. cytidine and uridine; mM range). IC₅₀ values were higher for nucleobases compared to the corresponding nucleosides; hypoxanthine had a 25-fold higher IC₅₀ than inosine. The archetypal human ENT1 inhibitor 4-nitrobenzylthioinosine (NBMPR) had no effect on PvENT1, whereas dipyridamole inhibited PvENT1, albeit with a 40 μM IC₅₀, a 1000-fold less sensitive than human ENT1 (hENT1). The PfENT1 inhibitors blocked transport activity of PvENT1 and the five known naturally occurring non-synonymous single nucleotide polymorphisms (SNPs) with similar IC₅₀ values. Thus, the PfENT1 inhibitors also target PvENT1. This implies that development of novel antimalarial drugs that target both falciparum and vivax ENT1 may be feasible.

•

PvENT1 can be functionally expressed in Saccharomyces cerevisiae.

•

PvENT1 transports purine and pyrimidine nucleosides and nucleobases but does not transport nucleotides.

•

PvENT1 is inhibited by recently described PfENT1 inhibitors with similar potency.

•

Identified PvENT1 non-synonymous SNPs do not change PfENT1 inhibitor potency.

•

Plasmodium ENTs may be feasible target for development of novel antimalarial drugs.

Association of collagen I, IX and vitamin D receptor gene polymorphisms with radiological severity of intervertebral disc degeneration in Southern European Ancestor

PURPOSE

Several genomic loci have been previously found to be associated with intervertebral disc degeneration, so far. Data are mostly derived from northern European countries whereas data derived from Southern European Ancestor are limited. This study aimed to evaluate the association between radiological disease severity of lumbar disc degeneration and certain genetic loci in a sample of participants from Southern Europe.

METHODS

Seventy-five patients with mild to severe lumbar disc degeneration and 25 healthy controls were enrolled into the study. In each subject, each lumbar intervertebral disc was separately examined to obtain a total radiological score for disease severity. In addition, single-nucleotide polymorphisms of predefined genetic samples were analyzed in all participants: COL1A1 Sp1, COL9a2 Trp2, COL9a3 Trp3, and VDR TaqI.

RESULTS

Degeneration scores were significantly worse in cases with COL1A1 Sp1, COL9a3 Trp3, and VDR TaqI mutations; however, COL9a2 Trp2 mutation was not associated with a difference in the severity of disc degeneration. In addition, subjects with mutation in more than one gene sample (n = 20) had significantly worse degeneration scores than the remaining study participants (n = 80) (17.70 ± 2.72 vs. 21.81 ± 1.81, p < 0.001).

CONCLUSION

Single-nucleotide polymorphisms occurring in COL1A1, COL9a3 and VDR genes seem to be associated with the development of lumbar disc degeneration in this cohort, possibly with even more pronounced association when multiple mutations are present in the same individual. By further prospective twin studies in associated genes and analyses of their relationship with environmental factors in an internationally sampled large cohort will make a more clear-minded conclusion about their association with disc degeneration, which would yield better appreciation and clinical planning of some predisposed people for these pathologies.

Metallothionein 2A core promoter region genetic polymorphism and its impact on the risk, tumor behavior, and recurrences of sinonasal inverted papilloma (Schneiderian papilloma)

Inverted papillomas are a unique group of locally aggressive benign epithelial neoplasms in the nasal cavity and paranasal sinuses arising from the Schneiderian mucosa. Metallothioneins are sulfhydryl-rich heavy metal-binding proteins required for metal toxicity protection and regulation of biological mechanisms including proliferation and invasion. The goal of this study was to identify three SNPs at loci −5 A/G (rs28366003) and −209 A/G (rs1610216) in the core promoter region and at locus +838 C/G (rs10636) in 3′UTR region of the MT2A gene with IP risk and with tumor invasiveness according to Krouse staging. Genotyping was performed using the PCR restriction fragment length polymorphism technique in 130 genetically unrelated IP individuals, and 418 randomly selected healthy volunteers. The presence of the rs28366003 SNP was significantly related to the risk of IP within the present population-based case-control study. Compared to homozygous common allele carriers, heterozygosity and homozygosity for the G variant had a significantly increased risk of IP (adjusted odds ratio [OR] = 7.71, 95 % confidence interval [CI]: 4.01–14.91, p_dominant < 0.001). Moreover, risk allele carriers demonstrated higher Krouse stage (pT1 vs. pT2-4) (OR = 19.32; 95 % CI, 2.30–173.53; p < 0.0001), diffuse tumor growth (OR = 4.58; 95 % CI, 1.70–12.11; p = 0.0008), bone destruction (OR = 4.13; 95 % CI, 1.50–11.60; p = 0.003), and higher incidence of tumor recurrences (OR = 5.11; 95 % CI, 1.68–15.20; p = 0.001). The findings suggest that MT2A gene variation rs28366003 may be implicated in the etiology of sinonasal inverted papilloma in a Polish population.

The effect of metallothionein 2A core promoter region single-nucleotide polymorphism on accumulation of toxic metals in sinonasal inverted papilloma tissues

Metallothioneins (MTs) are intracellular thiol-rich heavy metal-binding proteins which join trace metal ions protecting cells against heavy metal toxicity and regulate metal distribution and donation to various enzymes and transcription factors. The goal of this study was to identify the -5 A/G (rs28366003) single-nucleotide polymorphism (SNP) in the core promoter region of the MT2A gene, and to investigate its effect on allele-specific gene expression and Cd, Zn, Cu and Ni content in sinonasal inverted papilloma tissue (IP), with non-cancerous sinonasal mucosa (NCM) as a control. The MT2A promoter region -5 A/G SNP was identified by restriction fragment length polymorphism using 117 IP and 132 NCM. MT2A gene analysis was performed by quantitative real-time PCR. Metal levels were analyzed by flame atomic absorption spectrometry. The frequency of A allele carriage was 99.2% and 100% in IP and NCM, respectively. The G allele carriage was detected in 23.9% of IP and in 12.1% of the NCM samples. As a result, a significant association of -5 A/G SNP in MT2A gene with mRNA expression in both groups was determined. A significant association was identified between the -5 A/G SNP in the MT2A gene with mRNA expression in both groups. A highly significant association was detected between the rs28366003 genotype and Cd and Zn content in IP. Furthermore, significant differences were identified between A/A and A/G genotype with regard to the type of metal contaminant. The Spearman rank correlation results showed the MT2A gene expression and both Cd and Cu levels were negatively correlated. The results obtained in this study suggest that the -5 A/G SNP in the MT2A gene may have an effect on allele-specific gene expression and toxic metal accumulation in sinonasal inverted papilloma.

IL-18 polymorphisms in hepatitis B virus related liver disease

Interleukine-18 (IL-18) was originally called interferon (INF-γ) inducing factor and plays a critical dual role in Th1 polarization and viral clearance. We aimed to explore whether single-nucleotide promoter polymorphisms (SNPs) are associated with the outcome of hepatitis B virus (HBV) infection. 271 HBV infected patients were recruited in this study out of these 109 were spontaneously recovered and 162 were diagnosed to be having persistent HBV infection which includes 48 chronic hepatitis, 84 liver cirrhosis, 30 HCC cases and were compared with 280 healthy controls. IL-18 promoter genotyping was performed with sequence-specific primers. The results demonstrated the significant involvement of genotype AA at position -607 in healthy controls (38.6%) when compared to cases (26.0%) (OR=0.54 (0.385-0.797)) and also associated with spontaneous clearance (37.6%) compared to persistent HBV infections (17.9%) (OR=2.76 (1.582-4.832)). Whereas, genotype CC at position -607 in cases (18.0%) when compared to healthy controls (6.7%) (OR=3.03 (1.734-5.303)) also associated with persistent HBV infections (24.1%) compared to spontaneous clearance (9.2%) (OR=0.31 (0.151-0.67)). And genotype GC at position -137 in cases (49.5%) compared to healthy controls (38.5%) (OR=1.55 (1.11-2.18)). Whereas, genotype GG at position -137 in healthy controls (56.8%) compared to cases (45.4%) (OR=0.63 (0.451-0.885)). No significant difference at position -137 was observed between spontaneous clearance and persistent HBV infections. These polymorphisms of the IL-18 gene promoter region at position -607 and -137 could be associated with different outcomes of HBV infection. The people with allele A at position -607 may be protected against HBV infection; moreover AA genotype is associated with spontaneous clearance.

Scalable privacy-preserving data sharing methodology for genome-wide association studies

The protection of privacy of individual-level information in genome-wide association study (GWAS) databases has been a major concern of researchers following the publication of "an attack" on GWAS data by Homer et al. (2008). Traditional statistical methods for confidentiality and privacy protection of statistical databases do not scale well to deal with GWAS data, especially in terms of guarantees regarding protection from linkage to external information. The more recent concept of differential privacy, introduced by the cryptographic community, is an approach that provides a rigorous definition of privacy with meaningful privacy guarantees in the presence of arbitrary external information, although the guarantees may come at a serious price in terms of data utility. Building on such notions, Uhler et al. (2013) proposed new methods to release aggregate GWAS data without compromising an individual's privacy. We extend the methods developed in Uhler et al. (2013) for releasing differentially-private χ(2)-statistics by allowing for arbitrary number of cases and controls, and for releasing differentially-private allelic test statistics. We also provide a new interpretation by assuming the controls' data are known, which is a realistic assumption because some GWAS use publicly available data as controls. We assess the performance of the proposed methods through a risk-utility analysis on a real data set consisting of DNA samples collected by the Wellcome Trust Case Control Consortium and compare the methods with the differentially-private release mechanism proposed by Johnson and Shmatikov (2013).

Applications and challenges of next-generation sequencing in Brassica species

Next-generation sequencing (NGS) produces numerous (often millions) short DNA sequence reads, typically varying between 25 and 400 bp in length, at a relatively low cost and in a short time. This revolutionary technology is being increasingly applied in whole-genome, transcriptome, epigenome and small RNA sequencing, molecular marker and gene discovery, comparative and evolutionary genomics, and association studies. The Brassica genus comprises some of the most agro-economically important crops, providing abundant vegetables, condiments, fodder, oil and medicinal products. Many Brassica species have undergone the process of polyploidization, which makes their genomes exceptionally complex and can create difficulties in genomics research. NGS injects new vigor into Brassica research, yet also faces specific challenges in the analysis of complex crop genomes and traits. In this article, we review the advantages and limitations of different NGS technologies and their applications and challenges, using Brassica as an advanced model system for agronomically important, polyploid crops. Specifically, we focus on the use of NGS for genome resequencing, transcriptome sequencing, development of single-nucleotide polymorphism markers, and identification of novel microRNAs and their targets. We present trends and advances in NGS technology in relation to Brassica crop improvement, with wide application for sophisticated genomics research into agronomically important polyploid crops.

Identification and characterization of a novel variant of the human P2X(7) receptor resulting in gain of function

The P2X(7) receptor exhibits significant allelic polymorphism in humans, with both loss and gain of function variants potentially impacting on a variety of infectious and inflammatory disorders. At least five loss-of-function polymorphisms (G150R, R307Q, T357S, E496A, and I568N) and two gain-of-function polymorphisms (H155Y and Q460R) have been identified and characterized to date. In this study, we used RT-PCR cloning to isolate and characterize P2X(7) cDNA clones from human PBMCs and THP-1 cells. A previously unreported variant with substitutions of V80M and A166G was identified. When expressed in HEK293 cells, this variant exhibited heightened sensitivity to the P2X(7) agonist (BzATP) relative to the most frequent allele, as shown by pore formation measured by fluorescent dye uptake into cells. Mutational analyses showed that A166G alteration was critical for the gain-of-function change, while V80M was not. Full-length variants with multiple previously identified nonsynonymous SNPs (H155Y, H270R, A348T, and E496A) were also identified. Distinct functional phenotypes of the P2X(7) variants or mutants constructed with multiple polymorphisms were observed. Gain-of-function variations (A166G or H155Y) could not rescue the loss-of-function E496A polymorphism. Synergistic effects of the gain-of-function variations were also observed. We also identified the A348T alteration as a weak gain-of-function variant. Thus, these results identify the new gain-of-function variant A166G and demonstrate that multiple-gene polymorphisms contribute to functional phenotypes of the human P2X(7) receptor. Furthermore, the results demonstrate that the C-terminal of the cysteine-rich domain 1 of P2X(7) is critical for regulation of P2X(7)-mediated pore formation.