Identification of genome-wide single nucleotide polymorphisms in allopolyploid crop Brassica napus

Background

Single nucleotide polymorphisms (SNPs) are the most common type of genetic variation. Identification of large numbers of SNPs is helpful for genetic diversity analysis, map-based cloning, genome-wide association analyses and marker-assisted breeding. Recently, identifying genome-wide SNPs in allopolyploid Brassica napus (rapeseed, canola) by resequencing many accessions has become feasible, due to the availability of reference genomes of Brassica rapa (2n = AA) and Brassica oleracea (2n = CC), which are the progenitor species of B. napus (2n = AACC). Although many SNPs in B. napus have been released, the objective in the present study was to produce a larger, more informative set of SNPs for large-scale and efficient genotypic screening. Hence, short-read genome sequencing was conducted on ten elite B. napus accessions for SNP discovery. A subset of these SNPs was randomly selected for sequence validation and for genotyping efficiency testing using the Illumina GoldenGate assay.

Results

A total of 892,536 bi-allelic SNPs were discovered throughout the B. napus genome. A total of 36,458 putative amino acid variants were located in 13,552 protein-coding genes, which were predicted to have enriched binding and catalytic activity as a result. Using the GoldenGate genotyping platform, 94 of 96 SNPs sampled could effectively distinguish genotypes of 130 lines from two mapping populations, with an average call rate of 92%.

Conclusions

Despite the polyploid nature of B. napus, nearly 900,000 simple SNPs were identified by whole genome resequencing. These SNPs were predicted to be effective in high-throughput genotyping assays (51% polymorphic SNPs, 92% average call rate using the GoldenGate assay, leading to an estimated >450 000 useful SNPs). Hence, the development of a much larger genotyping array of informative SNPs is feasible. SNPs identified in this study to cause non-synonymous amino acid substitutions can also be utilized to directly identify causal genes in association studies.     

Endocytosis of the non-catalytic ADAM23: Recycling and long half-life properties

A Disintegrin And Metalloprotease 23 (ADAM23) is a member of the ADAMs family of transmembrane proteins, mostly expressed in nervous system, and involved in traffic and stabilization of Kv1-potassium channels, synaptic transmission, neurite outgrowth, neuronal morphology and cell adhesion. Also, ADAM23 has been linked to human pathological conditions, such as epilepsy, cancer metastasis and cardiomyopathy. ADAM23 functionality depends on the molecule presence at the cell surface and along the secretory pathway, as expected for a cell surface receptor. Because endocytosis is an important functional regulatory mechanism of plasma membrane receptors and no information is available about the traffic or turnover of non-catalytic ADAMs, we investigated ADAM23 internalization, recycling and half-life properties. Here, we show that ADAM23 undergoes constitutive internalization from the plasma membrane, a process that depends on lipid raft integrity, and is redistributed to intracellular vesicles, especially early and recycling endosomes. Furthermore, we observed that ADAM23 is recycled from intracellular compartments back to the plasma membrane and thus has longer half-life and higher cell surface stability compared with other ADAMs. Our findings suggest that regulation of ADAM23 endocytosis/stability could be exploited therapeutically in diseases in which ADAM23 is directly involved, such as epilepsy, cancer progression and cardiac hypertrophy.     

Accurate measurement of gene copy number for human alpha-defensin DEFA1A3

Background

Multi-allelic copy number variants include examples of extensive variation between individuals in the copy number of important genes, most notably genes involved in immune function. The definition of this variation, and analysis of its impact on function, has been hampered by the technical difficulty of large-scale but accurate typing of genomic copy number. The copy-variable alpha-defensin locus DEFA1A3 on human chromosome 8 commonly varies between 4 and 10 copies per diploid genome, and presents considerable challenges for accurate high-throughput typing.

Results

In this study, we developed two paralogue ratio tests and three allelic ratio measurements that, in combination, provide an accurate and scalable method for measurement of DEFA1A3 gene number. We combined information from different measurements in a maximum-likelihood framework which suggests that most samples can be assigned to an integer copy number with high confidence, and applied it to typing 589 unrelated European DNA samples. Typing the members of three-generation pedigrees provided further reassurance that correct integer copy numbers had been assigned. Our results have allowed us to discover that the SNP rs4300027 is strongly associated with DEFA1A3 gene copy number in European samples.

Conclusions

We have developed an accurate and robust method for measurement of DEFA1A3 copy number. Interrogation of rs4300027 and associated SNPs in Genome-Wide Association Study SNP data provides no evidence that alpha-defensin copy number is a strong risk factor for phenotypes such as Crohn’s disease, type I diabetes, HIV progression and multiple sclerosis.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-719) contains supplementary material, which is available to authorized users.     

Identification of a gene associated with avian migratory behaviour

Bird migration is one of the most spectacular and best-studied phenomena in behavioural biology. Yet, while the patterns of variation in migratory behaviour and its ecological causes have been intensively studied, its genetic, physiological and neurological control remains poorly understood. The lack of knowledge of the molecular basis of migration is currently not only limiting our insight into the proximate control of migration, but also into its evolution. We investigated polymorphisms in the exons of six candidate genes for key behavioural traits potentially linked to migration, which had previously been identified in several bird species, and eight control loci in 14 populations of blackcaps (Sylvia atricapilla), representing the whole range of geographical variation in migration patterns found in this species, with the aim of identifying genes controlling variation in migration. We found a consistent association between a microsatellite polymorphism and migratory behaviour only at one candidate locus: the ADCYAP1 gene. This polymorphism explained about 2.6 per cent of the variation in migratory tendency among populations, and 2.7-3.5% of variation in migratory restlessness among individuals within two independent populations. In all tests, longer alleles were associated with higher migratory activity. The consistency of results among different populations and levels of analysis suggests that ADCYAP1 is one of the genes controlling the expression of migratory behaviour. Moreover, the multiple described functions of the gene product indicate that this gene might act at multiple levels modifying the shift between migratory and non-migratory states.     

The ADAM metalloprotease Kuzbanian is crucial for proper heart formation in Drosophila melanogaster

We have screened a collection of EMS mutagenized fly lines in order to identify genes involved in cardiogenesis. In the present work, we have studied a group of alleles exhibiting a hypertrophic heart. Our analysis revealed that the ADAM protein (A Disintegrin And Metalloprotease) Kuzbanian, which is the functional homologue of the vertebrate ADAM10, is crucial for proper heart formation. ADAMs are a family of transmembrane proteins that play a critical role during the proteolytic conversion (shedding) of membrane bound proteins to soluble forms. Enzymes harboring a sheddase function recently became candidates for causing several congenital diseases, like distinct forms of the Alzheimer disease. ADAMs play also a pivotal role during heart formation and vascularisation in vertebrates, therefore mutations in ADAM genes potentially could cause congenital heart defects in humans. In Drosophila, the zygotic loss of an active form of the Kuzbanian protein results in a dramatic excess of cardiomyocytes, accompanied by a loss of pericardial cells. Our data presented herein suggest that Kuzbanian acts during lateral inhibition within the cardiac primordium. Furthermore we discuss a second function of Kuzbanian in heart cell morphogenesis.     

Identification and expression analysis of a heat-shock protein 70 gene in Polycelis sp.

Heat-shock protein 70 (HSP70) is ubiquitously found in a variety of organisms and plays an important role in cytoprotection, environmental monitoring, and disease resistance. In this study, the full-length complementary DNA (cDNA) of hsp70 from planarian Polycelis sp. was first cloned using rapid amplification of cDNA ends (RACE). The expression levels of Pyhsp70 were analyzed in the presence of various stressors by real-time PCR, and its temporal-spatial expression patterns were also examined in both intact and regenerative animals by whole-mount in situ hybridization. The results show that (1) the deduced amino acid sequence of Pyhsp70 includes three typical HSP70 family signature motifs and is highly conserved during evolution; (2) Pyhsp70 expression is induced by prolonged starvation, tissue damage, and ionic liquid but inhibited by high or low temperatures; and (3) Pyhsp70 mRNA is mainly expressed in the head peripheral region and in the regenerating blastema during regeneration. These results suggest that the highly expressed Pyhsp70 gene may contribute to enhance cytoprotection and tolerance against stress-induced molecular damage, and the migration of neoblasts to the wound, which might also be involved in the proliferation and differentiation of neoblasts. Our work provides basic data for the study of stress responses and regenerative mechanism in freshwater planarians.     

Evaluation of the 23S rRNA gene as target for qPCR based quantification of Frankia in soils

The 23S rRNA gene was evaluated as target for the development of Sybr Green-based quantitative PCR (qPCR) for the analysis of nitrogen-fixing members of the genus Frankia or subgroups of these in soil. A qPCR with a primer combination targeting all nitrogen-fixing frankiae (clusters 1, 2 and 3) resulted in numbers similar to those obtained with a previously developed qPCR using nifH gene sequences, both with respect to introduced and indigenous Frankia populations. Primer combinations more specifically targeting three subgroups of the Alnus host infection group (cluster 1) or members of the Elaeagnus host infection group (cluster 3) were specific for introduced strains of the target group, with numbers corresponding to those obtained by quantification of nitrogen-fixing frankiae with both the 23S rRNA and nifH genes as target. Method verification on indigenous Frankia populations in soils, i.e. in depth profiles from four sites at an Alnus glutinosa stand, revealed declining numbers in the depth profiles, with similar abundance of all nitrogen-fixing frankiae independent of 23S rRNA or nifH gene targets, and corresponding numbers of one group of frankiae of the Alnus host infection only, with no detections of frankiae representing the Elaeagnus, Casuarina, or a second subgroup of the Alnus host infection groups.     

Identification and characterization of LIM gene family in Brassica rapa

Background

LIM (Lin-11, Isl-1 and Mec-3 domains) genes have been reported to trigger the formation of actin bundles, a major higher-order cytoskeletal assembly, in higher plants; however, the stress resistance related functions of these genes are still not well known. In this study, we collected 22 LIM genes designated as Brassica rapa LIM (BrLIM) from the Brassica database, analyzed the sequences, compared them with LIM genes of other plants and analyzed their expression after applying biotic and abiotic stresses in Chinese cabbage.

Results

Upon sequence analysis these genes were confirmed as LIM genes and found to have a high degree of homology with LIM genes of other species. These genes showed distinct clusters when compared to other recognized LIM proteins upon phylogenetic analysis. Additionally, organ specific expression of these genes was observed in Chinese cabbage plants, with BrPLIM2a, b, c, BrDAR1, BrPLIM2e, f and g only being expressed in flower buds. Furthermore, the expression of these genes (except for BrDAR1 and BrPLIM2e) was high in the early flowering stages. The remaining genes were expressed in almost all organs tested. All BrDAR genes showed higher expression in flower buds compared to other organs. These organ specific expressions were clearly correlated with the phylogenetic grouping. In addition, BrWLIM2c and BrDAR4 responded to Fusarium oxysporum f. sp. conglutinans infection, while commonly two BrDARs and eight BrLIMs responded to cold, ABA and pH (pH5, pH7 and pH9) stress treatments in Chinese cabbage plants.

Conclusion

Taken together, the results of this study indicate that BrLIM and BrDAR genes may be involved in resistance against biotic and abiotic stresses in Brassica.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-641) contains supplementary material, which is available to authorized users.     

Identification of parasite DNA in common bile duct stones by PCR and DNA sequencing

We attempted to identify parasite DNA in the biliary stones of humans via PCR and DNA sequencing. Genomic DNA was isolated from each of 15 common bile duct (CBD) stones and 5 gallbladder (GB) stones. The patients who had the CBD stones suffered from cholangitis, and the patients with GB stones showed acute cholecystitis, respectively. The 28S and 18S rDNA genes were amplified successfully from 3 and/or 1 common bile duct stone samples, and then cloned and sequenced. The 28S and 18S rDNA sequences were highly conserved among isolates. Identity of the obtained 28S D1 rDNA with that of Clonorchis sinensis was higher than 97.6%, and identity of the 18S rDNA with that of other Ascarididae was 97.9%. Almost no intra-specific variations were detected in the 28S and 18S rDNA with the exception of a few nucleotide variations, i.e., substitution and deletion. These findings suggest that C. sinensis and Ascaris lumbricoides may be related with the biliary stone formation and development.     

Geographic distribution of haplotype diversity at the bovine casein locus

The genetic diversity of the casein locus in cattle was studied on the basis of haplotype analysis. Consideration of recently described genetic variants of the casein genes which to date have not been the subject of diversity studies, allowed the identification of new haplotypes. Genotyping of 30 cattle breeds from four continents revealed a geographically associated distribution of haplotypes, mainly defined by frequencies of alleles at CSN1S1 and CSN3. The genetic diversity within taurine breeds in Europe was found to decrease significantly from the south to the north and from the east to the west. Such geographic patterns of cattle genetic variation at the casein locus may be a result of the domestication process of modern cattle as well as geographically differentiated natural or artificial selection. The comparison of African Bos taurus and Bos indicus breeds allowed the identification of several Bos indicus specific haplotypes (CSN1S1*C-CSN2*A²-CSN3*A_I/CSN3*H) that are not found in pure taurine breeds. The occurrence of such haplotypes in southern European breeds also suggests that an introgression of indicine genes into taurine breeds could have contributed to the distribution of the genetic variation observed.     

Four exons of the serotonin receptor 4 gene are associated with multiple distant branch points

Splicing of vertebrate introns involves recognition of three consensus elements at the 3′ end. The branch point (BP) and polypyrimidine tract (PPT) are usually located within 40 nucleotides (nt) of the 3′ splice site (3′ ss), AG, but can be much more distant. A characteristic of the region between distant BPs (dBPs) and the 3′ ss is the absence of intervening AG dinucleotides, leading to its designation as the “AG exclusion zone” (AGEZ). The human HTR4 gene, which encodes serotonin receptor 4 and has been associated with schizophrenia, bipolar disease, and gastrointestinal disorders, has four exons with extensive AGEZs. We have mapped the BPs for HTR4 exons 3, 4, 5, and g generated by in vitro splicing, and validated them by mutagenesis in exon-trapping vectors. All exons used dBPs up to 273 nt upstream of the exon. Strikingly, exons 4 and 5 used combinations of both distant and conventionally located BPs, suggesting that successful splicing of these exons can occur by distinct pathways. Our results emphasize the importance for single nucleotide polymorphism resequencing projects to take account of potential dBPs, as the extended AGEZs are vulnerable to mutations that could affect splicing itself or regulation of alternative splicing.     

Identification of pseudouridine methyltransferase in Escherichia coli

In ribosomal RNA, modified nucleosides are found in functionally important regions, but their function is obscure. Stem–loop 69 of Escherichia coli 23S rRNA contains three modified nucleosides: pseudouridines at positions 1911 and 1917, and N3 methyl-pseudouridine (m³Ψ) at position 1915. The gene for pseudouridine methyltransferase was previously not known. We identified E. coli protein YbeA as the methyltransferase methylating Ψ1915 in 23S rRNA. The E. coli ybeA gene deletion strain lacks the N3 methylation at position 1915 of 23S rRNA as revealed by primer extension and nucleoside analysis by HPLC. Methylation at position 1915 is restored in the ybeA deletion strain when recombinant YbeA protein is expressed from a plasmid. In addition, we show that purified YbeA protein is able to methylate pseudouridine in vitro using 70S ribosomes but not 50S subunits from the ybeA deletion strain as substrate. Pseudouridine is the preferred substrate as revealed by the inability of YbeA to methylate uridine at position 1915. This shows that YbeA is acting at the final stage during ribosome assembly, probably during translation initiation. Hereby, we propose to rename the YbeA protein to RlmH according to uniform nomenclature of RNA methyltransferases. RlmH belongs to the SPOUT superfamily of methyltransferases. RlmH was found to be well conserved in bacteria, and the gene is present in plant and in several archaeal genomes. RlmH is the first pseudouridine specific methyltransferase identified so far and is likely to be the only one existing in bacteria, as m³Ψ1915 is the only methylated pseudouridine in bacteria described to date.     

Molecular characterization of porcine MMP19 and MMP23B genes and its association with immune traits

MMP19 and MMP23B belong to the Matrix metalloproteases (MMPs) family, which are zinc-binding endopeptidases that are capable of degrading various components of the extracellular matrix. They are thought to play important roles in embryonic development, reproduction and tissue remodeling, as well as in cell proliferation, differentiation, migration, angiogenesis, apoptosis and host defense. However, they are poorly understood in pigs. Here, we obtained the full length coding region sequence and genomic sequence of the porcine MMP19 and MMP23B genes and analyzed their genomic structures. The deduced amino acid sequence shares similar precursor protein domains with human and mouse MMP19 and MMP23B protein, respectively. Using IMpRH panel, MMP19 was mapped to SSC5p12-q11 (closely linked to microsatellite DK) and MMP23B was mapped to SSC8q11-q12 (linked to microsatellite Sw2521). Quantitative real-time PCR showed that MMP19 was abundantly expressed in the liver, while MMP23B was strongly expressed in the ovarian and heart. Furthermore, both genes were all expressed increasingly in prenatal skeletal muscle during development. Three SNPs were detected by sequencing and PCR-RFLP methods, and association analysis indicated that C203T at exon 5 of MMP19 has a significant association with the blood parameters WBC (G/L) and IgG2 (mg/mL) (P<0.05), SNP C131T at exon 3 of MMP23B is significantly associated with the blood parameters HGB (g/L) and MCH (P<0.05), and A150G in exon 4 has no significant association with the economic traits in pigs.     

Drought tolerance through over-expression of the expansin gene TaEXPB23 in transgenic tobacco

Expansins are proteins that are the key regulators of wall extension during plant growth. To investigate the role of TaEXPB23, a wheat expansin gene, we analyzed TaEXPB23 mRNA expression levels in response to water stress in wheat and examined the drought resistance of transgenic tobaccos over-expressing TaEXPB23. We found that the expression of TaEXPB23 corresponded to wheat coleoptile growth and the response to water stress. The results also indicated that the transgenic tobacco lines lost water more slowly than the wild-type (WT) plants under drought stress; their cells could sustain a more integrated structure under water stress than that of WT. Other physiological and biochemical parameters under water stress, such as electrolyte leakage, malondialdehyde (MDA) level, photosynthetic rate, F_v/F_m and ΦPSII, also suggested that the transgenic tobaccos were more drought resistant than WT plants.