Putative SNPs Research Articles

Functional Polymorphism in the MSI1 Gene Promoter Confers a Decreased Risk of Lung Cancer in Chinese by Reducing MSI1 Expression.

Background: Musashi1 (MSI1) is a characteristic stem cell marker that regulates the balance between cell self-renewal and differentiation. Evidence has identified MSI1 as a pivotal oncogenic regulator in diverse malignancies. However, little evidence uncovers the role of genetic variations of MSI1 gene in cancer etiology.Objective: The aim of this study was to investigate the association between genetic variants in the MSI1 gene and lung cancer risk.Methods: Based on a two-stage retrospective study with a total of 1559 patients with lung cancer and 1667 healthy controls, we evaluated the relevance between three putative functional SNPs in the MSI1 promoter (i.e., -2696T>C[rs7959801], -2297T>C[rs3742038] and -1081C>T[rs34570155]) and lung cancer risk.Results: We found that the SNP rs7959801T>C was significantly associated with lung cancer susceptibility. Compared to those with rs7959801TT wild-genotype, individuals with CT/CC variant genotypes exerted consistently beneficial roles in lung cancer risk in the discovery set (adjusted odd ratios [OR] = 0.67; 95% confidence interval [CI] = 0.57-0.80), and in the validation set (OR=0.69; 95%CI=0.54-0.88). Functional assays indicated that the allele transformation from T to C in rs7959801 of MSI1 gene arrestingly decreased its transcription activity in vitro. Furthermore, the expression levels of MSI1 were significantly lower in the patients with CT/CC variants than in those who were with TT genotype.Conclusion: Our findings suggested that the rs7959801T>C polymorphism in the MSI1 promoter conferred a decreased risk to lung cancer by reducing the expression of MSI1 and it may be a promising indicator for lung cancer predisposition.

TIA: algorithms for development of identity-linked SNP islands for analysis by massively parallel DNA sequencing

BackgroundSingle nucleotide polymorphisms (SNPs) located within the human genome have been shown to have utility as markers of identity in the differentiation of DNA from individual contributors. Massively parallel DNA sequencing (MPS) technologies and human genome SNP databases allow for the design of suites of identity-linked target regions, amenable to sequencing in a multiplexed and massively parallel manner. Therefore, tools are needed for leveraging the genotypic information found within SNP databases for the discovery of genomic targets that can be evaluated on MPS platforms.ResultsThe SNP island target identification algorithm (TIA) was developed as a user-tunable system to leverage SNP information within databases. Using data within the 1000 Genomes Project SNP database, human genome regions were identified that contain globally ubiquitous identity-linked SNPs and that were responsive to targeted resequencing on MPS platforms. Algorithmic filters were used to exclude target regions that did not conform to user-tunable SNP island target characteristics. To validate the accuracy of TIA for discovering these identity-linked SNP islands within the human genome, SNP island target regions were amplified from 70 contributor genomic DNA samples using the polymerase chain reaction. Multiplexed amplicons were sequenced using the Illumina MiSeq platform, and the resulting sequences were analyzed for SNP variations. 166 putative identity-linked SNPs were targeted in the identified genomic regions. Of the 309 SNPs that provided discerning power across individual SNP profiles, 74 previously undefined SNPs were identified during evaluation of targets from individual genomes. Overall, DNA samples of 70 individuals were uniquely identified using a subset of the suite of identity-linked SNP islands.ConclusionsTIA offers a tunable genome search tool for the discovery of targeted genomic regions that are scalable in the population frequency and numbers of SNPs contained within the SNP island regions. It also allows the definition of sequence length and sequence variability of the target region as well as the less variable flanking regions for tailoring to MPS platforms. As shown in this study, TIA can be used to discover identity-linked SNP islands within the human genome, useful for differentiating individuals by targeted resequencing on MPS technologies.

RNAseq analysis of bronchial epithelial cells to identify COPD-associated genes and SNPs

BackgroundThere is a need for more powerful methods to identify low-effect SNPs that contribute to hereditary COPD pathogenesis. We hypothesized that SNPs contributing to COPD risk through cis-regulatory effects are enriched in genes comprised by bronchial epithelial cell (BEC) expression patterns associated with COPD.MethodsTo test this hypothesis, normal BEC specimens were obtained by bronchoscopy from 60 subjects: 30 subjects with COPD defined by spirometry (FEV1/FVC < 0.7, FEV1% < 80%), and 30 non-COPD controls. Targeted next generation sequencing was used to measure total and allele-specific expression of 35 genes in genome maintenance (GM) genes pathways linked to COPD pathogenesis, including seven TP53 and CEBP transcription factor family members. Shrinkage linear discriminant analysis (SLDA) was used to identify COPD-classification models. COPD GWAS were queried for putative cis-regulatory SNPs in the targeted genes.ResultsOn a network basis, TP53 and CEBP transcription factor pathway gene pair network connections, including key DNA repair gene ERCC5, were significantly different in COPD subjects (e.g., Wilcoxon rank sum test for closeness, p-value = 5.0E-11). ERCC5 SNP rs4150275 association with chronic bronchitis was identified in a set of Lung Health Study (LHS) COPD GWAS SNPs restricted to those in putative regulatory regions within the targeted genes, and this association was validated in the COPDgene non-hispanic white (NHW) GWAS. ERCC5 SNP rs4150275 is linked (D’ = 1) to ERCC5 SNP rs17655 which displayed differential allelic expression (DAE) in BEC and is an expression quantitative trait locus (eQTL) in lung tissue (p = 3.2E-7). SNPs in linkage (D’ = 1) with rs17655 were predicted to alter miRNA binding (rs873601). A classifier model that comprised gene features CAT, CEBPG, GPX1, KEAP1, TP73, and XPA had pooled 10-fold cross-validation receiver operator characteristic area under the curve of 75.4% (95% CI: 66.3%–89.3%). The prevalence of DAE was higher than expected (p = 0.0023) in the classifier genes.ConclusionsGM genes comprised by COPD-associated BEC expression patterns were enriched for SNPs with cis-regulatory function, including a putative cis-rSNP in ERCC5 that was associated with COPD risk. These findings support additional total and allele-specific expression analysis of gene pathways with high prior likelihood for involvement in COPD pathogenesis.

Association of μ-Calpain and Calpastatin Polymorphisms with Meat Tenderness in a Brahman-Angus Population.

Autogenous proteolytic enzymes of the calpain family are implicated in myofibrillar protein degradation. As a result, the μ-calpain gene and its specific inhibitor, calpastatin, have been repeatedly investigated for their association with meat quality traits in cattle; however, no functional mutation has been identified for these two genes. The objectives of this study were: (1) to assess breed composition effect on tenderness; (2) to perform a linkage disequilibrium (LD) analysis in μ-calpain and calpastatin genes as well as an association analyses with tenderness; and (3) to analyze putative functional SNPs inside the significant LD block for an effect on tenderness. Tenderness measurements and genotypes for 16 SNPs in μ-calpain gene and 28 SNPs in calpastatin gene from 673 steers were analyzed. A bioinformatic analysis identified “putative functional SNPs” inside the associated LD block – polymorphisms able to produce a physical and/or chemical change in the DNA, mRNA, or translated protein in silico. Breed composition had a significant (P < 0.0001) effect on tenderness where animals with more than 80% Angus composition had the most tender meat. One 11-kb LD-block and three LD-blocks of 37, 17, and 14 kb in length were identified in the μ-calpain and calpastatin genes, respectively. Out of these, the LD-block 3 in calpastatin, tagged by SNPs located at 7-98566391 and 7-98581038, had a significant effect on tenderness with the TG-CG diplotype being approximately 1 kg more tender than the toughest diplotype, TG-CG. A total of 768 SNPs in the LD-block 3 of calpastatin were included in the bioinformatic analysis, and 28 markers were selected as putative functional SNPs inside the LD-block 3 of calpastatin; however, none of them were polymorphic in this population. Out of 15 initial polymorphisms segregating inside the LD-block 3 of calpastatin in this population, markers ARSUSMARC116, Cast5, rs730723459, and rs210861835 were found to be significantly associated with tenderness.

Inferring causal relationships between phenotypes using summary statistics from genome-wide association studies.

Genome-wide association studies (GWAS) have successfully identified numerous genetic variants associated with diverse complex phenotypes and diseases, and provided tremendous opportunities for further analyses using summary association statistics. Recently, Pickrell et al. developed a robust method for causal inference using independent putative causal SNPs. However, this method may fail to infer the causal relationship between two phenotypes when only a limited number of independent putative causal SNPs identified. Here, we extended Pickrell's method to make it more applicable for the general situations. We extended the causal inference method by replacing the putative causal SNPs with the lead SNPs (the set of the most significant SNPs in each independent locus) and tested the performance of our extended method using both simulation and empirical data. Simulations suggested that when the same number of genetic variants is used, our extended method had similar distribution of test statistic under the null model as well as comparable power under the causal model compared with the original method by Pickrell et al. But in practice, our extended method would generally be more powerful because the number of independent lead SNPs was often larger than the number of independent putative causal SNPs. And including more SNPs, on the other hand, would not cause more false positives. By applying our extended method to summary statistics from GWAS for blood metabolites and femoral neck bone mineral density (FN-BMD), we successfully identified ten blood metabolites that may causally influence FN-BMD. We extended a causal inference method for inferring putative causal relationship between two phenotypes using summary statistics from GWAS, and identified a number of potential causal metabolites for FN-BMD, which may provide novel insights into the pathophysiological mechanisms underlying osteoporosis.

DNA methylation rather than single nucleotide polymorphisms regulates the production of an aberrant splice variant of IL6R in mastitic cows.

Interleukin-6 receptor-alpha (IL6R) interacts with IL6 and forms a ligand-receptor complex, which can stimulate various cellular responses, such as cell proliferation, cell differentiation, and activation of inflammatory processes. Both genetic mutation and epigenetic modification regulate gene transcription. We identified a novel splice variant of bovine IL6R, designated as IL6R-TV, which is characterized by the skipping of exon 2 of the NCBI-referenced IL6R gene (IL6R-reference). The expression levels of IL6R-TV and IL6R-reference transcripts were lower in normal mammary gland tissues. These transcripts play a potential role during inflammatory infection. We also detected two putative functional SNPs (g.19711 T > C and g.19731 G > C) located within the upstream 100bp of exon 2. These SNPs formed two haplotypes (T-G and C-C). Two mutant pSPL3 exon-trapping plasmids (pSPL3-T-G and pSPL3-C-C) were transferred into the bovine mammary epithelial cells (MAC-T) and human embryonic kidney 293 T cells (HEK293T) to investigate the relationship between the two SNPs and the aberrant splicing of IL6R. DNA methylation levels of the alternatively spliced exon in normal and mastitis-infected mammary gland tissues were quantified through nested bisulfate sequencing PCR (BSP) and cloning sequencing. We found that DNA methylation regulated IL6R transcription. The DNA methylation level was high in mastitis-infected mammary gland tissues and stimulated IL6R expression, thereby promoting the inclusion of the alternatively spliced exon. The upregulated expression of the two transcripts was due to DNA methylation modification rather than genetic mutations.

Identification and in silico analysis of functional SNPs of human TAGAP protein: A comprehensive study.

Genetic polymorphisms in TAGAP gene have been associated with many diseases including rheumatoid arthritis, multiple sclerosis and other autoimmune disorders. Identifying functional SNPs in such disease associated genes is an uphill task hence before planning larger population study, it is better to scrutinize putative functional SNPs. In this study we used various computational approaches to identify nsSNPs which are deleterious to the structure and/or function of TAGAP protein that might be causing these diseases. Computational analysis was performed by five different in silico tools including SIFT, PROVEAN, PolyPhen-2, PhD-SNP and SNPs&GO. The study concludes that mutations of Glycine → Glutamic Acid at position 120, Glycine → Tryptophan at position 141 and Valine → Methionine at position 151 are major mutations in native TAGAP protein which might contribute to its malfunction and ultimately causing disease. The study also proposed 3D structures of native TAGAP protein and its three mutants. Future studies should consider these nsSNPs as main target mutations in various diseases involving TAGAP malfunction. This is the first comprehensive study, where TAGAP gene variants were analyzed using in silico tools hence will be of great help while considering large scale studies and also in developing precision medicines for cure of diseases related to these polymorphisms. Furthermore, animal models of various autoimmune diseases and having these mutations might be of help in exploring their precise roles.

Association between 3'UTR polymorphisms in genes ACVR2A, AGTR1 and RGS2 and preeclampsia.

Preeclampsia (PE) is a pregnancy specific disease with several risk factors such as genetic polymorphisms, environmental and social factors participating in its development. The aim of this study was to investigate whether distribution of three putative regulatory SNPs rs13430086, rs5186, rs4606 in 3'UTR of genes ACVR2A, AGTR1 and RGS2, respectively, that have been associated with hypertension and regulation of trophoblast invasion differ between women with PE and control group. The associations of rs13430086, rs5186 and rs4606 with preeclampsia were tested in two groups - the group of 50 women with PE and the control group of 42 healthy pregnant women at term. DNA was isolated from blood samples and the determination of genotypes was performed using Real-Time PCR. Power analysis for the size of the cohort was performed and the results were analyzed using Fisher exact test. The AA genotype of ACVR2A rs13430086 was significantly associated with higher risk to preeclampsia compared with TT genotype (p = 0.026, OR: 5.39, 95%CI: 1.21-31.54). Results showed no association between genotypes and preeclampsia for polymorphisms rs5186, rs4606. Further studies are important in order to better understand the role of ACVR2A in the pathogenesis of PE.

Development of nuclear SNP markers for the timber tracking of the African tree species Sapelli, Entandrophragma cylindricum

We describe the development of new nuclear SNP markers for the genetic timber tracking of the geographical origin of Sapelli, Entandrophragma cylindricum (Meliaceae). Restriction associated DNA sequencing (RADseq) of two reference individuals yielded 1131 putative SNPs. Among those, 131 were selected to design four MassARRAY multiplexes and screened at 178 individuals. Seventy-two loci were selected for further use in genetic tracking.

The developmental transcriptome landscape of receptive endometrium during embryo implantation in dairy goats

Under natural conditions, some embryos cannot implant successfully because of the dysfunction of receptive endometrium (RE). Thus, it is imperative for us to study the molecular mechanisms involved in the formation of the RE from pre-receptive endometrium (PE). In this study, the endometrium from gestational day 5 (D5, PE) and gestational day 15 (D15, RE) dairy goats were selected to systematically analyze the transcriptome using strand-specific Ribo-Zero RNA-Seq, >120 million high-quality paired-end reads were generated and 47,616 transcripts were identified in the endometrium of dairy goats. A total of 810 mRNAs were differentially expressed genes (DEGs) between the RE and PE meeting the criteria of P-values<0.05. Bioinformatics analysis of the DEGs revealed that a number of biological processes and pathways were potentially involved in the establishment of the RE, notably energy metabolism and amino acid metabolism. Furthermore, we speculated that CXCL14, IGFBP3, and LGALS15 potentially participated in the development of endometrium. What's more, putative SNPs, InDels and AS events were identified and analyzed in the endometrium. In a word, this resulting view of the transcriptome greatly enhances the comprehensive transcript catalog and uncovers the global trends in gene expression during the formation of receptive endometrium in dairy goats.

Identification of SNPs associated with muscle yield and quality traits using allelic-imbalance analyses of pooled RNA-Seq samples in rainbow trout

BackgroundCoding/functional SNPs change the biological function of a gene and, therefore, could serve as “large-effect” genetic markers. In this study, we used two bioinformatics pipelines, GATK and SAMtools, for discovering coding/functional SNPs with allelic-imbalances associated with total body weight, muscle yield, muscle fat content, shear force, and whiteness. Phenotypic data were collected for approximately 500 fish, representing 98 families (5 fish/family), from a growth-selected line, and the muscle transcriptome was sequenced from 22 families with divergent phenotypes (4 low- versus 4 high-ranked families per trait).ResultsGATK detected 59,112 putative SNPs; of these SNPs, 4798 showed allelic imbalances (>2.0 as an amplification and <0.5 as loss of heterozygosity). SAMtools detected 87,066 putative SNPs; and of them, 4962 had allelic imbalances between the low- and high-ranked families. Only 1829 SNPs with allelic imbalances were common between the two datasets, indicating significant differences in algorithms. The two datasets contained 7930 non-redundant SNPs of which 4439 mapped to 1498 protein-coding genes (with 6.4% non-synonymous SNPs) and 684 mapped to 295 lncRNAs. Validation of a subset of 92 SNPs revealed 1) 86.7–93.8% success rate in calling polymorphic SNPs and 2) 95.4% consistent matching between DNA and cDNA genotypes indicating a high rate of identifying SNPs with allelic imbalances. In addition, 4.64% SNPs revealed random monoallelic expression. Genome distribution of the SNPs with allelic imbalances exhibited high density for all five traits in several chromosomes, especially chromosome 9, 20 and 28. Most of the SNP-harboring genes were assigned to important growth-related metabolic pathways.ConclusionThese results demonstrate utility of RNA-Seq in assessing phenotype-associated allelic imbalances in pooled RNA-Seq samples. The SNPs identified in this study were included in a new SNP-Chip design (available from Affymetrix) for genomic and genetic analyses in rainbow trout.

Genomic footprints of adaptation in a cooperatively breeding tropical bird across a vegetation gradient.

Identifying the genetic basis of phenotypic variation and its relationship with the environment is key to understanding how local adaptations evolve. Such patterns are especially interesting among populations distributed across habitat gradients, where genetic structure can be driven by isolation by distance (IBD) and/or isolation by environment (IBE). Here, we used variation in ~1,600 high-quality SNPs derived from paired-end sequencing of double-digest restriction site-associated DNA (ddRAD-Seq) to test hypotheses related to IBD and IBE in the Yucatan jay (Cyanocorax yucatanicus), a tropical bird endemic to the Yucatán Peninsula. This peninsula is characterized by a precipitation and vegetation gradient-from dry to evergreen tropical forests-that is associated with morphological variation in this species. We found a moderate level of nucleotide diversity (π=.008) and little evidence for genetic differentiation among vegetation types. Analyses of neutral and putatively adaptive SNPs (identified by complementary genome-scan approaches) indicate that IBD is the most reliable explanation to account for frequency distribution of the former, while IBE has to be invoked to explain those of the later. These results suggest that selective factors acting along a vegetation gradient can promote local adaptation in the presence of gene flow in a vagile, nonmigratory and geographically restricted species. The putative candidate SNPs identified here are located within or linked to a variety of genes that represent ideal targets for future genomic surveys.

Development and characterization of SNP derived from spinyhead croaker (Collichthys lucidus) by RNA-seq

A novel set of SNP markers were detected from RNA-seq data of spinyhead croaker (Collichthys lucidus). Totally 49,617,970 clean reads and 78,671 putative SNPs were obtained. Among them, 40 SNPs were randomly selected to perform validation in 30 wild individuals of C. lucidus captured from 3 geographic stocks in the coast of China. Eventually, 29 SNP loci were successfully validated and showed bi-allelic polymorphism in genotyping by Sequenom MassARRAY. These novel SNPs provided useful genomic resources for C. lucidus, especially in population genetics, molecular assisted selection and conservation studies.

Deciphering the distance to antibiotic resistance for the pneumococcus using genome sequencing data

Advances in genome sequencing technologies and genome-wide association studies (GWAS) have provided unprecedented insights into the molecular basis of microbial phenotypes and enabled the identification of the underlying genetic variants in real populations. However, utilization of genome sequencing in clinical phenotyping of bacteria is challenging due to the lack of reliable and accurate approaches. Here, we report a method for predicting microbial resistance patterns using genome sequencing data. We analyzed whole genome sequences of 1,680 Streptococcus pneumoniae isolates from four independent populations using GWAS and identified probable hotspots of genetic variation which correlate with phenotypes of resistance to essential classes of antibiotics. With the premise that accumulation of putative resistance-conferring SNPs, potentially in combination with specific resistance genes, precedes full resistance, we retrogressively surveyed the hotspot loci and quantified the number of SNPs and/or genes, which if accumulated would confer full resistance to an otherwise susceptible strain. We name this approach the ‘distance to resistance’. It can be used to identify the creep towards complete antibiotics resistance in bacteria using genome sequencing. This approach serves as a basis for the development of future sequencing-based methods for predicting resistance profiles of bacterial strains in hospital microbiology and public health settings.

High-Density SNP Genotyping Array for Hexaploid Wheat and Its Relatives.

A lack of genetic diversity between wheat breeding lines has been recognized as a significant block to future yield increases. Wheat breeding and prebreeding strategies are increasingly using material from wheat ancestors or wild relatives to reintroduce diversity. Where molecular markers are polymorphic between the host and introgressed material, they may be used to track the size and location of the introgressed material through generations of backcrossing. To generate markers for this purpose, sequence capture targeted resequencing was carried out for a range of wheat varieties, wheat relatives, and wheat progenitors. From these sequences, putative SNPs were identified and used to generate the Axiom® Wheat HD array. A selection of varieties representing a selection of elite wheat breeding material, progenitor species, and wild relatives were used to validate the array. The procedures used are described here in detail.

Genome-Wide SNPs Identification and Determination of Proteins Associated with Stress Response in Sorghum (&amp;lt;i&amp;gt;Sorghum bicolor&amp;lt;/i&amp;gt; L. Monech) Accessions

Current efforts in sorghum breeding programs are exploiting genotyping-by-sequencing (GBS) data to provide full-genome scans for desired traits. The aim of this study was to utilize GBS approach for the identification of genomic regions associated with stress response in sorghum (Sorghum bicolor L. Monech) accessions. DNA samples of twenty sorghum accessions, having different response to drought, were used to prepare GBS libraries for sequencing. SNPs were called using the TASSELGBS pipeline and the tags that present at least 10 times in the dataset were considered and aligned to the reference genome of Sorghum bicolor. The identified SNPs were all compared with the published sorghum transcript related to stress response gene activity. Overall; 94.40% tags were aligned and 69,736 putative SNPs positions were identified. Blast search revealed homology to annotated heat and drought–tolerance associated genes which code for ATPases, Peroxidase, Hydrophobic protein LTI6A, Aquaporin SIP2-1, Aconitate hydratase and phosphatidylinositol-4-phosphate-5-kinase. The phylogeny of the 20 accessions was constructed using the generated SNPs data. Phylogenetic analysis data showed that the phenotypically tolerant line (El9) makes a separate cluster and the same for the accessions HSD8653 and HSD5612 near to the cluster that includes most accessions with known post-flowering drought tolerance (HSD7410, HDS10033, HSD8552, GESHEISH and HSD8849). Post-flowering drought sensitive accessions (Tabat, Wadahmed, HSD6468 and HSD6478) formed a separate cluster while the sensitive accession HSD9959 and the tolerant accessions HSD8511 and HSD9566 were distributed between the two clusters. Thus, cluster analysis confirmed the variation among accessions in post-flowering drought tolerance. With further validation, these markers may be used for marker assisted selection for breeding new sorghum genotypes with stress adaptation.

Generalization and fine mapping of European ancestry-based central adiposity variants in African ancestry populations.

Background/ObjectivesCentral adiposity measures such as waist circumference (WC) and waist-to-hip ratio (WHR) are associated with cardiometabolic disorders independently of BMI and are gaining clinically utility. Several studies report genetic variants associated with central adiposity, but most utilize only European ancestry populations. Understanding whether the genetic associations discovered among mainly European descendants are shared with African ancestry populations will help elucidate the biological underpinnings of abdominal fat deposition.Subjects/MethodsTo identify the underlying functional genetic determinants of body fat distribution, we conducted an array-wide association meta-analysis among persons of African ancestry across seven studies/consortia participating in the Population Architecture using Genomics and Epidemiology (PAGE) consortium. We used the Metabochip array, designed for fine mapping cardiovascular associated loci, to explore novel array-wide associations with WC and WHR among 15 945 African descendants using all and sex-stratified groups. We further interrogated 17 known WHR regions for African ancestry-specific variants.ResultsOf the 17 WHR loci, eight SNPs located in four loci were replicated in the sex-combined or sex-stratified meta-analyses. Two of these eight independently associated with WHR after conditioning on the known variant in European descendants (rs12096179 in TBX15-WARS2 and rs2059092 in ADAMTS9). In the fine mapping assessment, the putative functional region was reduced across all four loci but to varying degrees (average 40% drop in number of putative SNPs and 20% drop in genomic region). Similar to previous studies, the significant SNPs in the female stratified analysis were stronger than the significant SNPs from the sex-combined analysis. No novel associations were detected in the array-wide analyses.ConclusionsOf 17 previously identified loci, four loci replicated in the African ancestry populations of this study. Utilizing different linkage disequilibrium patterns observed between European and African ancestries, we narrowed the suggestive region containing causative variants for all four loci.

Large-Scale Genotyping-by-Sequencing Indicates High Levels of Gene Flow in the Deep-Sea Octocoral Swiftia simplex (Nutting 1909) on the West Coast of the United States.

Deep-sea corals are a critical component of habitat in the deep-sea, existing as regional hotspots for biodiversity, and are associated with increased assemblages of fish, including commercially important species. Because sampling these species is so difficult, little is known about the connectivity and life history of deep-sea octocoral populations. This study evaluates the genetic connectivity among 23 individuals of the deep-sea octocoral Swiftia simplex collected from Eastern Pacific waters along the west coast of the United States. We utilized high-throughput restriction-site associated DNA (RAD)-tag sequencing to develop the first molecular genetic resource for the deep-sea octocoral, Swiftia simplex. Using this technique we discovered thousands of putative genome-wide SNPs in this species, and after quality control, successfully genotyped 1,145 SNPs across individuals sampled from California to Washington. These SNPs were used to assess putative population structure across the region. A STRUCTURE analysis as well as a principal coordinates analysis both failed to detect any population differentiation across all geographic areas in these collections. Additionally, after assigning individuals to putative population groups geographically, no significant FST values could be detected (FST for the full data set 0.0056), and no significant isolation by distance could be detected (p = 0.999). Taken together, these results indicate a high degree of connectivity and potential panmixia in S. simplex along this portion of the continental shelf.

Development of SNP markers and validation 31 SNPs in Coilia nasus

Estuarine tapertail anchovy (Coilia nasus) is a widely distributed and commercially important aquaculture species. However, excessive fishing and changes in aquatic ecology have almost caused extinction of the species in the middle reaches of the Yangtze River. Therefore, it is very urgent to perform genetic research on C. nasus to protect the wild population. In this study, we sequenced the liver transcriptome of C. nasus using Illumina Hiseq 2500 platform. We obtained 65,129 unigenes with an average length of 606 bp and a total of 62,947 putative SNPs were detected. Among them, 31 SNPs were selected for the validation experiments. All of the 31 loci were found to be polymorphic and showed bi-allelic in 58 individuals of C. nasus. These SNP markers should not only be useful for population conservation, but also for construction of genetic linkage map and economic performance improvement of C. nasus.

Exploring genomic databases for in silico discovery of Pht1 genes in high syntenic close related grass species with focus in sugarcane (Saccharum spp.)

The genomic sequences available at several public databases and the local alignments allow a fast and accurate search for nucleotide homology among different species, revealing the structure, composition and function of nucleotides and amino acids. The Phosphate Transporter 1 gene family (Pht1) has crucial role in several metabolic processes from plants and has been identified in different crops. Here, the CDS sequence of 28 Pht1 genes described for Oryza sativa, Zea mays and Arabidopsis thaliana were used to identify nine homologous sequences in Sorghum bicolor and Setaria italica and the S. bicolor sequences were used as reference to guide the identification of those homologous sequences in Saccharum spp. Using five different libraries of paired-end Illumina reads, nine sugarcane genes (ScPht1) were assembled and the rate of gene recover, the average of gene length, the sequencing depth and the GC content were estimated in 93.72%, 2053bp, 16.08X and 58.70%, respectively. Putative SNPs and microsatellites loci were found for those ScPht1 genes and three of them seem to be under positive selection pressure. A phylogenetic tree corroborates with orthology relationship among Saccharum spp., O. sativa and Z. mays genes, emphasizing the levels of synteny described for grass species. This study discuss the potential use of homology in genomic studies as one of the most useful tools for gene discovery and annotation, illustrating how useful is the exploitation of public databases for gene identification and characterization, allowing further gene manipulation in crop improvements.