Identification of alternate polyadenylation sites and analysis of their tissue distribution using EST data.

Alternate polyadenylation affects a large fraction of higher eucaryote mRNAs, producing mature transcripts with 3' ends of variable length. This variation is poorly represented in the current transcript catalogs derived from whole genome sequences, mostly because such posttranscriptional events are not detectable directly at the DNA level. Alternate polyadenylation of an mRNA is better understood by comparison to EST databases. Comparing ESTs to mRNAs, however, is a difficult task subjected to the pitfalls of internal priming, presence of intron sequences, repeated elements, chimerical ESTs or matches with EST from paralogous genes. We present here a computer program that addresses these problems and displays ESTs matches to a query mRNA sequence to predict alternate polyadenylation and to suggest library-specific forms. The output highlights effective polyadenylation signals, possible sources of artifacts such as A-rich stretches in the mRNA sequences, and allows for a direct visualization of EST libraries using color codes. Statistical biases in the distribution of alternative mRNA forms among EST libraries were systematically sought. About 1450 human and 200 mouse mRNAs displayed such biases, suggesting in each case a tissue- or disease-specific regulation of polyadenylation.

Patterns of variant polyadenylation signal usage in human genes.

The formation of mature mRNAs in vertebrates involves the cleavage and polyadenylation of the pre-mRNA, 10-30 nt downstream of an AAUAAA or AUUAAA signal sequence. The extensive cDNA data now available shows that these hexamers are not strictly conserved. In order to identify variant polyadenylation signals on a large scale, we compared over 8700 human 3' untranslated sequences to 157,775 polyadenylated expressed sequence tags (ESTs), used as markers of actual mRNA 3' ends. About 5600 EST-supported putative mRNA 3' ends were collected and analyzed for significant hexameric sequences. Known polyadenylation signals were found in only 73% of the 3' fragments. Ten single-base variants of the AAUAAA sequence were identified with a highly significant occurrence rate, potentially representing 14.9% of the actual polyadenylation signals. Of the mRNAs, 28.6% displayed two or more polyadenylation sites. In these mRNAs, the poly(A) sites proximal to the coding sequence tend to use variant signals more often, while the 3'-most site tends to use a canonical signal. The average number of ESTs associated with each signal type suggests that variant signals (including the common AUUAAA) are processed less efficiently than the canonical signal and could therefore be selected for regulatory purposes. However, the position of the site in the untranslated region may also play a role in polyadenylation rate.

Size estimate of the alpha beta TCR repertoire of naive mouse splenocytes.

The diversity of the T cell repertoire of mature T splenocytes is generated, in the thymus, by pairing of alpha and beta variable domains of the alpha beta TCR and by the rearrangements of various gene segments encoding these domains. In the periphery, it results from competition between various T cell subpopulations including recent thymic migrants and long-lived T cells. Quantitative data on the actual size of the T cell repertoire are lacking. Using PCR methods and extensive sequencing, we have measured for the first time the size of the TCR-alpha beta repertoire of naive mouse T splenocytes. There are 5-8 x 105 different nucleotide sequences of BV chains in the whole spleen of young adult mice. We have also determined the size of the BV repertoire in a subpopulation of AV2+ T splenocytes, which allows us to provide a minimum estimate of the alpha beta repertoire. We find that the mouse spleen harbors about 2 x 106 clones of about 10 cells each. This figure, although orders of magnitude smaller than the maximum theoretical diversity (estimated up to 1015), is still large enough to maintain a high functional diversity.

Quantitative analysis of the T cell repertoire selected by a single peptide-major histocompatibility complex.

The positive selection of CD4+ T cells requires the expression of major histocompatibility complex (MHC) class II molecules in the thymus, but the role of self-peptides complexed to class II molecules is still a matter of debate. Recently, it was observed that transgenic mice expressing a single peptide-MHC class II complex positively select significant numbers of diverse CD4+ T cells in the thymus. However, the number of selected T cell specificities has not been evaluated so far. Here, we have sequenced 700 junctional complementarity determining regions 3 (CDR3) from T cell receptors (TCRs) carrying Vbeta11-Jbeta1.1 or Vbeta12-Jbeta1.1 rearrangements. We found that a single peptide-MHC class II complex positively selects at least 10(5) different Vbeta rearrangements. Our data yield a first evaluation of the size of the T cell repertoire. In addition, they provide evidence that the single Ealpha52-68-I-Ab complex skews the amino acid frequency in the TCR CDR3 loop of positively selected T cells. A detailed analysis of CDR3 sequences indicates that a fraction of the beta chain repertoire bears the imprint of the selecting self-peptide.