Structural genomic variation as risk factor for idiopathic recurrent miscarriage.

Recurrent miscarriage (RM) is a multifactorial disorder with acknowledged genetic heritability that affects ∼3% of couples aiming at childbirth. As copy number variants (CNVs) have been shown to contribute to reproductive disease susceptibility, we aimed to describe genome-wide profile of CNVs and identify common rearrangements modulating risk to RM. Genome-wide screening of Estonian RM patients and fertile controls identified excessive cumulative burden of CNVs (5.4 and 6.1 Mb per genome) in two RM cases possibly increasing their individual disease risk. Functional profiling of all rearranged genes within RM study group revealed significant enrichment of loci related to innate immunity and immunoregulatory pathways essential for immune tolerance at fetomaternal interface. As a major finding, we report a multicopy duplication (61.6 kb) at 5p13.3 conferring increased maternal risk to RM in Estonia and Denmark (meta-analysis, n = 309/205, odds ratio = 4.82, P = 0.012). Comparison to Estonian population-based cohort (total, n = 1000) confirmed the risk for Estonian female cases (P = 7.9 × 10(-4) ). Datasets of four cohorts from the Database of Genomic Variants (total, n = 5,846 subjects) exhibited similar low duplication prevalence worldwide (0.7%-1.2%) compared to RM cases of this study (6.6%-7.5%). The CNV disrupts PDZD2 and GOLPH3 genes predominantly expressed in placenta and it may represent a novel risk factor for pregnancy complications.

Position-specific scoring matrix and hidden Markov model complement each other for the prediction of conopeptide superfamilies.

Classified into 16 superfamilies, conopeptides are the main component of cone snail venoms that attract growing interest in pharmacology and drug discovery. The conventional approach to assigning a conopeptide to a superfamily is based on a consensus signal peptide of the precursor sequence. While this information is available at the genomic or transcriptomic levels, it is not present in amino acid sequences of mature bioactives generated by proteomic studies. As the number of conopeptide sequences is increasing exponentially with the improvement in sequencing techniques, there is a growing need for automating superfamily elucidation. To face this challenge we have defined distinct models of the signal sequence, propeptide region and mature peptides for each of the superfamilies containing more than 5 members (14 out of 16). These models rely on two robust techniques namely, Position-Specific Scoring Matrices (PSSM, also named generalized profiles) and hidden Markov models (HMM). A total of 50 PSSMs and 47 HMM profiles were generated. We confirm that propeptide and mature regions can be used to efficiently classify conopeptides lacking a signal sequence. Furthermore, the combination of all three-region models demonstrated improvement in the classification rates and results emphasise how PSSM and HMM approaches complement each other for superfamily determination. The 97 models were validated and offer a straightforward method applicable to large sequence datasets.

High-resolution picture of a venom gland transcriptome: case study with the marine snail Conus consors.

Although cone snail venoms have been intensively investigated in the past few decades, little is known about the whole conopeptide and protein content in venom ducts, especially at the transcriptomic level. If most of the previous studies focusing on a limited number of sequences have contributed to a better understanding of conopeptide superfamilies, they did not give access to a complete panorama of a whole venom duct. Additionally, rare transcripts were usually not identified due to sampling effect. This work presents the data and analysis of a large number of sequences obtained from high throughput 454 sequencing technology using venom ducts of Conus consors, an Indo-Pacific living piscivorous cone snail. A total of 213,561 Expressed Sequence Tags (ESTs) with an average read length of 218 base pairs (bp) have been obtained. These reads were assembled into 65,536 contiguous DNA sequences (contigs) then into 5039 clusters. The data revealed 11 conopeptide superfamilies representing a total of 53 new isoforms (full length or nearly full-length sequences). Considerable isoform diversity and major differences in transcription level could be noted between superfamilies. A, O and M superfamilies are the most diverse. The A family isoforms account for more than 70% of the conopeptide cocktail (considering all ESTs before clustering step). In addition to traditional superfamilies and families, minor transcripts including both cysteine free and cysteine-rich peptides could be detected, some of them figuring new clades of conopeptides. Finally, several sets of transcripts corresponding to proteins commonly recruited in venom function could be identified for the first time in cone snail venom duct. This work provides one of the first large-scale EST project for a cone snail venom duct using next-generation sequencing, allowing a detailed overview of the venom duct transcripts. This leads to an expanded definition of the overall cone snail venom duct transcriptomic activity, which goes beyond the cysteine-rich conopeptides. For instance, this study enabled to detect proteins involved in common post-translational maturation and folding, and to reveal compounds classically involved in hemolysis and mechanical penetration of the venom into the prey. Further comparison with proteomic and genomic data will lead to a better understanding of conopeptides diversity and the underlying mechanisms involved in conopeptide evolution.

Identification and analysis of papillomavirus E2 protein binding sites in the human genome.

Papillomavirus E2 protein is required for the replication and maintenance of viral genomes and transcriptional regulation of viral genes. E2 functions through sequence-specific binding to 12-bp DNA motifs-E2 binding sites (E2BS)-in the virus genome. Papillomaviruses are able to establish persistent infection in their host and have developed a long-term relationship with the host cell in order to guarantee the propagation of the virus. In this study, we have analyzed the occurrence and functionality of E2BSs in the human genome. Our computational analysis indicates that most E2BSs in the human genome are found in repetitive DNA regions and have G/C-rich spacer sequences. Using a chromatin immunoprecipitation approach, we show that human papillomavirus type 11 (HPV11) E2 interacts with a subset of cellular E2BSs located in active chromatin regions. Two E2 activities, sequence-specific DNA binding and interaction with cellular Brd4 protein, are important for E2 binding to consensus sites. E2 binding to cellular E2BSs has a moderate or no effect on cellular transcription. We suggest that the preference of HPV E2 proteins for E2BSs with A/T-rich spacers, which are present in the viral genomes and underrepresented in the human genome, ensures E2 binding to specific binding sites in the virus genome and may help to prevent extensive and possibly detrimental changes in cellular transcription in response to the viral protein.

Fluoride-cleavable, fluorescently labelled reversible terminators: synthesis and use in primer extension.

Fluorescent 2'-deoxynucleotides containing a protecting group at the 3'-O-position are reversible terminators that enable array-based DNA sequencing-by-synthesis (SBS) approaches. Herein, we describe the synthesis and full characterisation of four reversible terminators bearing a 3'-blocking moiety and a linker-dye system that is removable under the same fluoride-based treatment. Each nucleotide analogue has a different fluorophore attached to the base through a fluoride-cleavable linker and a 2-cyanoethyl moiety as the 3'-blocking group, which can be removed by using a fluoride treatment as well. Furthermore, we identified a DNA polymerase, namely, RevertAid M-MuLV reverse transcriptase, which can incorporate the four modified reversible terminators. The synthesised nucleotides and the optimised DNA polymerase were used on CodeLink slides spotted with hairpin oligonucleotides to demonstrate their potential in a cyclic reversible terminating approach.

Haplotypes in the human Foxo1a and Foxo3a genes; impact on disease and mortality at old age.

Recently, the Daf-16 gene has been shown to regulate the lifespan of nematodes and flies. In mammals, the Daf-16 homologues are forkhead (FOXO) transcription factors, of which specific functions have been identified for Foxo1a and Foxo3a. Despite that, their influence on human age-related trajectories and lifespan is unknown. Here, we analysed the effect of genetic variance in Foxo1a and Foxo3a on metabolic profile, age-related diseases, fertility, fecundity and mortality. This study was carried out in the prospective population-based Leiden 85-plus Study, which includes 1245 participants, aged 85 years or more. The mean follow-up time was 4.4 years. Haplotype analyses of Foxo1a revealed that carriers of haplotype 3 'TCA' have higher HbA1c levels (P=0.025) and a 1.14-fold higher all-cause mortality risk (P=0.021). This increase in mortality was attributable to death from diabetes, for which a 2.43-fold increase was observed (P=0.025). The analyses with Foxo3a haplotypes revealed no differences in metabolic profile, fertility or fecundity. However, increased risks of stroke were observed for Foxo3a block-A haplotype 2 'GAGC' (P=0.007) and haplotype 4 'AAAT' (P=0.014) carriers. In addition, the haplotype 2 'GAGC' carriers had a 1.13-fold increased risk for all-cause mortality (P=0.036) and 1.19-fold increased risk for cardiovascular mortality (P=0.052). In conclusion, this study shows that genetic variation in evolutionarily conserved Foxo1a and Foxo3a genes influences lifespan in our study population.

Highly expressed proteins have an increased frequency of alanine in the second amino acid position.

BACKGROUND: Although the sequence requirements for translation initiation regions have been frequently analysed, usually the highly expressed genes are not treated as a separate dataset. RESULTS: To investigate this, we analysed the mRNA regions downstream of initiation codons in nine bacteria, three archaea and three unicellular eukaryotes, comparing the dataset of highly expressed genes to the dataset of all genes. In addition to the detailed analysis of the nucleotide and codon frequencies we compared the N-termini of highly expressed proteins to the N-termini of all proteins coded in the genome. CONCLUSION: The most conserved pattern was observed at the amino acid level: strong alanine over-representation was observed at the second amino acid position of highly expressed proteins. This pattern is well conserved in all three domains of life.