Combining RT-PCR-seq and RNA-seq to catalog all genic elements encoded in the human genome.

Within the ENCODE Consortium, GENCODE aimed to accurately annotate all protein-coding genes, pseudogenes, and noncoding transcribed loci in the human genome through manual curation and computational methods. Annotated transcript structures were assessed, and less well-supported loci were systematically, experimentally validated. Predicted exon-exon junctions were evaluated by RT-PCR amplification followed by highly multiplexed sequencing readout, a method we called RT-PCR-seq. Seventy-nine percent of all assessed junctions are confirmed by this evaluation procedure, demonstrating the high quality of the GENCODE gene set. RT-PCR-seq was also efficient to screen gene models predicted using the Human Body Map (HBM) RNA-seq data. We validated 73% of these predictions, thus confirming 1168 novel genes, mostly noncoding, which will further complement the GENCODE annotation. Our novel experimental validation pipeline is extremely sensitive, far more than unbiased transcriptome profiling through RNA sequencing, which is becoming the norm. For example, exon-exon junctions unique to GENCODE annotated transcripts are five times more likely to be corroborated with our targeted approach than with extensive large human transcriptome profiling. Data sets such as the HBM and ENCODE RNA-seq data fail sampling of low-expressed transcripts. Our RT-PCR-seq targeted approach also has the advantage of identifying novel exons of known genes, as we discovered unannotated exons in ~11% of assessed introns. We thus estimate that at least 18% of known loci have yet-unannotated exons. Our work demonstrates that the cataloging of all of the genic elements encoded in the human genome will necessitate a coordinated effort between unbiased and targeted approaches, like RNA-seq and RT-PCR-seq.

Statistical properties of population differentiation estimators under stepwise mutation in a finite island model.

Microsatellite loci mutate at an extremely high rate and are generally thought to evolve through a stepwise mutation model. Several differentiation statistics taking into account the particular mutation scheme of the microsatellite have been proposed. The most commonly used is R(ST) which is independent of the mutation rate under a generalized stepwise mutation model. F(ST) and R(ST) are commonly reported in the literature, but often differ widely. Here we compare their statistical performances using individual-based simulations of a finite island model. The simulations were run under different levels of gene flow, mutation rates, population number and sizes. In addition to the per locus statistical properties, we compare two ways of combining R(ST) over loci. Our simulations show that even under a strict stepwise mutation model, no statistic is best overall. All estimators suffer to different extents from large bias and variance. While R(ST) better reflects population differentiation in populations characterized by very low gene-exchange, F(ST) gives better estimates in cases of high levels of gene flow. The number of loci sampled (12, 24, or 96) has only a minor effect on the relative performance of the estimators under study. For all estimators there is a striking effect of the number of samples, with the differentiation estimates showing very odd distributions for two samples.

New microsatellite markers developed from an enriched microsatellite common bean library

The objective of this work was to develop new microsatellite markers in common bean. Ninety nine new microsatelitte loci were developed from a microsatellite enriched library for (CT)8 and (GT)8 motifs, from CAL-143 line. The majority of microsatellite sequences (51%) was related to cellular metabolism. The remaining sequences were associated to transcription functions. Only 17.2% of the sequences presented some level of similarity with other plant species genes.

Evidence for morphological and adaptive genetic divergence between lake and stream habitats in European minnows (Phoxinus phoxinus, Cyprinidae).

Natural selection drives local adaptation, potentially even at small temporal and spatial scales. As a result, adaptive genetic and phenotypic divergence can occur among populations living in different habitats. We investigated patterns of differentiation between contrasting lake and stream habitats in the cyprinid fish European minnow (Phoxinus phoxinus) at both the morphological and genomic levels using geometric morphometrics and AFLP markers, respectively. We also used a spatial correlative approach to identify AFLP loci associated with environmental variables representing potential selective forces responsible for adaptation to divergent habitats. Our results identified different morphologies between lakes and streams, with lake fish presenting a deeper body and caudal peduncle compared to stream fish. Body shape variation conformed to a priori predictions concerning biomechanics and swimming performance in lakes vs. streams. Moreover, morphological differentiation was found to be associated with several environmental variables, which could impose selection on body and caudal peduncle shape. We found adaptive genetic divergence between these contrasting habitats in the form of 'outlier' loci (2.9%) whose genetic divergence exceeded neutral expectations. We also detected additional loci (6.6%) not associated with habitat type (lake vs. stream), but contributing to genetic divergence between populations. Specific environmental variables related to trophic dynamics, landscape topography and geography were associated with several neutral and outlier loci. These results provide new insights into the morphological divergence and genetic basis of adaptation to differentiated habitats.

Genetic variation in GIPR influences the glucose and insulin responses to an oral glucose challenge.

Glucose levels 2 h after an oral glucose challenge are a clinical measure of glucose tolerance used in the diagnosis of type 2 diabetes. We report a meta-analysis of nine genome-wide association studies (n = 15,234 nondiabetic individuals) and a follow-up of 29 independent loci (n = 6,958-30,620). We identify variants at the GIPR locus associated with 2-h glucose level (rs10423928, beta (s.e.m.) = 0.09 (0.01) mmol/l per A allele, P = 2.0 x 10(-15)). The GIPR A-allele carriers also showed decreased insulin secretion (n = 22,492; insulinogenic index, P = 1.0 x 10(-17); ratio of insulin to glucose area under the curve, P = 1.3 x 10(-16)) and diminished incretin effect (n = 804; P = 4.3 x 10(-4)). We also identified variants at ADCY5 (rs2877716, P = 4.2 x 10(-16)), VPS13C (rs17271305, P = 4.1 x 10(-8)), GCKR (rs1260326, P = 7.1 x 10(-11)) and TCF7L2 (rs7903146, P = 4.2 x 10(-10)) associated with 2-h glucose. Of the three newly implicated loci (GIPR, ADCY5 and VPS13C), only ADCY5 was found to be associated with type 2 diabetes in collaborating studies (n = 35,869 cases, 89,798 controls, OR = 1.12, 95% CI 1.09-1.15, P = 4.8 x 10(-18)).

Genetic and karyotypic structure of a hybrid zone between the chromosomal races Cordon and Valais in the common shrew, Sorex araneus

Hybrid zones between genetically differentiated populations provide material to study evolutionary processes. Since the discovery of chromosomal races in Sorex araneus, contact zones have attracted attention of scientists. So far, studies on genetic markers in Sorex hybrid zones are missing. The acrocentric chromosomal race Cordon and the highly metacentric race Valais meet and hybridize at Les Houches in the Western Alps. On a transect through the hybrid zone, 273 shrews were caught at 15 localities over 4 years. Karyotype as well as the nuclear protein loci Alb and Pg were analyzed. F-st and F-is values were calculated by F-statistics. An analysis on pooled samples revealed the genetical differences between the hybridizing races as the only cause of population structuring. Genetical markers show dines with very strong frequency shifts at a mountain torrent, but behave differently through the zone. The performance of the torrent in maintaining the hybrid zone, selection against hybrids, possible assortative mating and linkage of the Valais Pg allele to a diagnostic chromosome arm, are discussed.

QTL mapping for protein content in soybean cultivated in two tropical environments

The objectives of this study were to detect quantitative trait loci (QTL) for protein content in soybean grown in two distinct tropical environments and to build a genetic map for protein content. One hundred eighteen soybean recombinant inbred lines (RIL), obtained from a cross between cultivars BARC 8 and Garimpo, were used. The RIL were cultivated in two distinct Brazilian tropical environments: Cascavel county, in Paraná, and Viçosa county, in Minas Gerais (24º57'S, 53º27'W and 20º45'S, 42º52'W, respectively). Sixty-six SSR primer pairs and 65 RAPD primers were polymorphic and segregated at a 1:1 proportion. Thirty poorly saturated linkage groups were obtained, with 90 markers and 41 nonlinked markers. For the lines cultivated in Cascavel, three QTL were mapped in C2, E and N linkage groups, which explained 14.37, 10.31 and 7.34% of the phenotypic variation of protein content, respectively. For the lines cultivated in Viçosa, two QTL were mapped in linkage groups G and #1, which explained 9.51 and 7.34% of the phenotypic variation of protein content. Based on the mean of the two environments, two QTL were identified: one in the linkage group E (9.90%) and other in the group L (7.11%). In order for future studies to consistently detect QTL effects of different environments, genotypes with greater stability should be used.

Efficient targeted transcript discovery via array-based normalization of RACE libraries.

Rapid amplification of cDNA ends (RACE) is a widely used approach for transcript identification. Random clone selection from the RACE mixture, however, is an ineffective sampling strategy if the dynamic range of transcript abundances is large. To improve sampling efficiency of human transcripts, we hybridized the products of the RACE reaction onto tiling arrays and used the detected exons to delineate a series of reverse-transcriptase (RT)-PCRs, through which the original RACE transcript population was segregated into simpler transcript populations. We independently cloned the products and sequenced randomly selected clones. This approach, RACEarray, is superior to direct cloning and sequencing of RACE products because it specifically targets new transcripts and often results in overall normalization of transcript abundance. We show theoretically and experimentally that this strategy leads indeed to efficient sampling of new transcripts, and we investigated multiplexing the strategy by pooling RACE reactions from multiple interrogated loci before hybridization.

Reduced genetic diversity, increased isolation and multiple introductions of invasive giant hogweed in the western Swiss Alps.

Abstract The giant hogweed (Heracleum mantegazzianum) has successfully invaded 19 European countries as well as parts of North America. It has become a problematic species due to its ability to displace native flora and to cause public health hazards. Applying population genetics to species invasion can help reconstruct invasion history and may promote more efficient management practice. We thus analysed levels of genetic variation and population genetic structure of H. mantegazzianum in an invaded area of the western Swiss Alps as well as in its native range (the Caucasus), using eight nuclear microsatellite loci together with plastid DNA markers and sequences. On both nuclear and plastid genomes, native populations exhibited significantly higher levels of genetic diversity compared to invasive populations, confirming an important founder event during the invasion process. Invasive populations were also significantly more differentiated than native populations. Bayesian clustering analysis identified five clusters in the native range that corresponded to geographically and ecologically separated groups. In the invaded range, 10 clusters occurred. Unlike native populations, invasive clusters were characterized by a mosaic pattern in the landscape, possibly caused by anthropogenic dispersal of the species via roads and direct collection for ornamental purposes. Lastly, our analyses revealed four main divergent groups in the western Swiss Alps, likely as a consequence of multiple independent establishments of H. mantegazzianum.

Molecular features of hepatosplenic T-cell lymphoma unravels potential novel therapeutic targets.

The pathogenesis of hepatosplenic T-cell lymphoma (HSTL), a rare entity mostly derived from γδ T cells and usually with a fatal outcome, remains largely unknown. In this study, HSTL samples (7γδ and 2αβ) and the DERL2 HSTL cell line were subjected to combined gene-expression profiling and array-based comparative genomic hybridization. Compared with other T-cell lymphomas, HSTL had a distinct molecular signature irrespective of TCR cell lineage. Compared with peripheral T-cell lymphoma, not otherwise specified and normal γδ T cells, HSTL overexpressed genes encoding NK-cell-associated molecules, oncogenes (FOS and VAV3), the sphingosine-1-phosphatase receptor 5 involved in cell trafficking, and the tyrosine kinase SYK, whereas the tumor-suppressor gene AIM1 (absent in melanoma 1) was among the most down-expressed. We found highly methylated CpG islands of AIM1 in DERL2 cells, and decitabine treatment induced a significant increase in AIM1 transcripts. Syk was present in HSTL cells and DERL2 cells contained phosphorylated Syk and were sensitive to a Syk inhibitor in vitro. Genomic profiles confirmed recurrent isochromosome 7q (n = 6/9) without alterations at the SYK and AIM1 loci. Our results identify a distinct molecular signature for HSTL and highlight oncogenic pathways that offer rationale for exploring new therapeutic options such as Syk inhibitors and demethylating agents.

Molecular characterization of grapevine from Santa Catarina, Brazil, using microsatellite markers

The objective of this work was to characterize the grape germplasm in Santa Catarina, Brazil, using microsatellite DNA markers (simple sequence repeats - SSR). The DNA samples were collected from leaves and shoots of accessions of public and private collections from the counties Urussanga, Nova Trento, Rodeio, São Joaquim, Campos Novos, Videira, and Água Doce. Ten SSR loci (VVS2, VVMD5, VVMD7, VVMD27, VrZAG62, VrZAG79, VVMD25, VVMD28, VVMD31, and VVMD32) were analysed by capillary electrophoresis. Molecular profiling was conducted for 190 grapevines (European, American, and hybrids), and 67 genotypes were obtained. The data were compared with each other and with those from the literature and from online databases, in order to identify varieties and discover cases of synonymy and homonymy. Forty molecular profiles corresponded to known varieties, while 27 genotypes were described for the first time. The existence of typical germplasm composed mainly of American and hybrid varieties is an important finding for local viticulture. Applications of the results rely on quality control and certification at the nursery level. Increasing precision in the characterization of grapevine genotypes may help breeding programs.

In situ and genetic characterization of Gossypium barbadense populations from the states of Pará and Amapá, Brazil

The objective of this work was to characterize the populations of Gossypium barbadense in the states of Amapá and Pará, Brazil. In situ characterization was conducted through interviews with the owners of the plants and environmental observations. Leaf or petal tissue as well as seed samples were collected for genetic characterization by single sequence repeats markers and for storage in germplasm banks, respectively. The plants were maintained in dooryards and used mainly for medical purposes. The genetic analysis showed no heterozygous plants at the loci tested (f = 1), indicating that reproduction occurs mainly through selfing. The total genetic diversity was high (He = 0.39); and a high level of differentiation was observed between cotton plants from the two states (F ST = 0.36). Conventional methods of in situ maintenance of G. barbadense populations are not applicable. The conservation of the genetic variability of populations present in the two states could be achieved through germplasm collection and establishing of ex situ seed banks.

Genetic variants influencing circulating lipid levels and risk of coronary artery disease.

OBJECTIVE: Genetic studies might provide new insights into the biological mechanisms underlying lipid metabolism and risk of CAD. We therefore conducted a genome-wide association study to identify novel genetic determinants of low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and triglycerides. METHODS AND RESULTS: We combined genome-wide association data from 8 studies, comprising up to 17 723 participants with information on circulating lipid concentrations. We did independent replication studies in up to 37 774 participants from 8 populations and also in a population of Indian Asian descent. We also assessed the association between single-nucleotide polymorphisms (SNPs) at lipid loci and risk of CAD in up to 9 633 cases and 38 684 controls. We identified 4 novel genetic loci that showed reproducible associations with lipids (probability values, 1.6×10(-8) to 3.1×10(-10)). These include a potentially functional SNP in the SLC39A8 gene for HDL-C, an SNP near the MYLIP/GMPR and PPP1R3B genes for LDL-C, and at the AFF1 gene for triglycerides. SNPs showing strong statistical association with 1 or more lipid traits at the CELSR2, APOB, APOE-C1-C4-C2 cluster, LPL, ZNF259-APOA5-A4-C3-A1 cluster and TRIB1 loci were also associated with CAD risk (probability values, 1.1×10(-3) to 1.2×10(-9)). CONCLUSIONS: We have identified 4 novel loci associated with circulating lipids. We also show that in addition to those that are largely associated with LDL-C, genetic loci mainly associated with circulating triglycerides and HDL-C are also associated with risk of CAD. These findings potentially provide new insights into the biological mechanisms underlying lipid metabolism and CAD risk.

A novel homozygous R764H mutation in crumbs homolog 1 causes autosomal recessive retinitis pigmentosa.

PURPOSE: Retinitis pigmentosa (RP; MIM 268000) is a hereditary disease characterized by poor night vision and progressive loss of photoreceptors, eventually leading to blindness. This degenerative process primarily affects peripheral vision due to the loss of rods. Autosomal recessive RP (arRP) is clinically and genetically heterogeneous. It has been associated with mutations in different genes, including CRB1 (crumbs homolog 1). The aim of this study was to determine the causative gene in a Tunisian patient with arRP born to non-consanguineous parents. METHODS: Four accessible family members were included. They underwent full ophthalmic examination with best-corrected Snellen visual acuity, fundus photography and fluorescein angiography. Haplotype analysis was used to evaluate homozygosity in the family to 20 arRP loci. All exons and intron-exon junctions of candidate genes not excluded by haplotype analysis were PCR amplified and directly sequenced. RESULTS: The proband was a 43-year-old female patient. Best-corrected visual acuity was 20/63 (right eye) and 20/80 (left eye). Visual loss began during the third decade. Funduscopic examination and fluorescein angiography revealed typical advanced RP changes with bone spicule-like pigment deposits in the posterior pole and the midperiphery along with retinal atrophy, narrowing of the vessels, and waxy optic discs. Haplotype analysis revealed homozygosity with microsatellite markers D1S412 and D1S413 on chromosome 1q31.3. These markers flanked CRB1. Our results excluded linkage of all the other arRP loci/genes tested. Sequencing of the 12 coding exons and splice sites of CRB1 disclosed a homozygous missense mutation in exon 7 at nucleotide c. 2291G>A, resulting in an arginine to histidine substitution (p.R764H). CONCLUSIONS: R764H is a novel mutation associated with CRB1-related arRP. Previously, an R764C mutation was reported. Extending the mutation spectrum of CRB1 with additional families is important for genotype-phenotype correlations and characterization of the scope of mutation.

Diversidad genética de piracanjuba usada en programas de repoblación con marcadores microsatélites

El objetivo de este trabajo fue estimar la diversidad genética de un lote de Brycon orbignyanus usado en programas de repoblación, a través de marcadores microsatélites. Se analizaron muestras de 44 reproductores, de 70 larvas y de 69 alevinos, con la amplificación de cinco loci descritos para Brycon opalinus. El número de alelos, la heterozigosidad observada (Ho) y esperada (He), el índice de Shannon (IS), la diversidad genética de Nei (DGN), el coeficiente de endogamia (Fis), la distancia (DG) e identidad genética (IG), el número efectivo de alelos, el test del equilibrio de Hardy-Weinberg (EHW) y el desequilibrio de ligación fueron calculados. Reproductores y progenie tuvieron un número similar de alelos en los loci evaluados. La Ho media, IS, DGN, DG e IG mostraron que existe menor distancia genética entre parentales y larvas y una disminución de variabilidad genética en los alevinos. Fueron observados desvíos en EHW y desequilibrio de ligación en seis pares de loci. El Fis mostró exceso de heterocigotos en parentales y larvas y déficit de heterocigotos en los alevinos. El lote de reproductores está en proceso de pérdida de alelos y hubo disminución de la variabilidad genética entre la fase de larva y alevino.