Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Fractional hepatic de novo lipogenesis in healthy subjects during near-continuous oral nutrition and bed rest: a comparison with published data in artificially fed, critically ill patients.

BACKGROUND AND AIMS: In critically ill patients, fractional hepatic de novo lipogenesis increases in proportion to carbohydrate administration during isoenergetic nutrition. In this study, we sought to determine whether this increase may be the consequence of continuous enteral nutrition and bed rest. We, therefore, measured fractional hepatic de novo lipogenesis in a group of 12 healthy subjects during near-continuous oral feeding (hourly isoenergetic meals with a liquid formula containing 55% carbohydrate). In eight subjects, near-continuous enteral nutrition and bed rest were applied over a 10 h period. In the other four subjects, it was extended to 34 h. Fractional hepatic de novo lipogenesis was measured by infusing(13) C-labeled acetate and monitoring VLDL-(13)C palmitate enrichment with mass isotopomer distribution analysis. Fractional hepatic de novo lipogenesis was 3.2% (range 1.5-7.5%) in the eight subjects after 10 h of near continuous nutrition and 1.6% (range 1.3-2.0%) in the four subjects after 34 h of near-continuous nutrition and bed rest. This indicates that continuous nutrition and physical inactivity do not increase hepatic de novo lipogenesis. Fractional hepatic de novo lipogenesis previously reported in critically ill patients under similar nutritional conditions (9.3%) (range 5.3-15.8%) was markedly higher than in healthy subjects (P<0.001). These data from healthy subjects indicate that fractional hepatic de novo lipogenesis is increased in critically ill patients.

Carbohydrate metabolism and de novo lipogenesis in human obesity.

Respiratory exchange was measured during 14 consecutive hours in six lean and six obese individuals after ingestion of 500 g of dextrin maltose to investigate and compare their capacity for net de novo lipogenesis. After ingestion of the carbohydrate load, metabolic rates rose similarly in both groups but fell earlier and more rapidly in the obese. RQs also rose rapidly and remained in the range of 0.95 to 1.00 for approximately 8 h in both groups. During this time, RQ exceeded 1.00 for only short periods of time with the result that 4 +/- 1 g and 5 +/- 3 g (NS) of fat were synthesized via de novo lipogenesis in excess of concomitant fat oxidation in the lean and obese subjects, respectively. Results demonstrate that net de novo lipid synthesis from an unusually large carbohydrate load is not greater in obese than in lean individuals.

De novo mutation in the KCNQ1 gene causal to Jervell and Lange-Nielsen syndrome.

Jervell and Lange-Nielsen syndrome (JLNS) is an autosomal recessive disorder, clinically characterized by severe cardiac arrhythmias [due to prolonged QTc interval in electrocardiogram (ECG)] and bilateral sensory neural deafness. Molecular defects causal to JLNS are either homozygous or compound heterozygous mutations, predominantly in the KCNQ1 gene and occasionally in the KCNE1 gene. As the molecular defect is bi-allelic, JLNS patients inherit one pathogenic mutation causal to the disorder from each parent. In this report, we show for the first time that such a disorder could also occur due to a spontaneous de novo mutation in the affected individual, not inherited from the parent, which makes this case unique unlike the previously reported JLNS cases.

Development and testing of novel chloroplast microsatellite markers for Lolium perenne and other grasses (Poaceae) from de novo sequencing and in silico sequences

Twelve primers to amplify microsatellite markers from the chloroplast genome of Lolium perenne were designed and optimized using de novo sequencing and in silico sequences. With one exception, each locus was polymorphic with a range from two to nine alleles in L. perenne. The newly developed primer pairs cross-amplified in different species of Lolium and in 50 other grass species representing nine grass subfamilies.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Early occurrence of lung adenocarcinoma and breast cancer after radiotherapy of a chest wall sarcoma in a patient with a de novo germline mutation in TP53.

We report a 26-year-old female patient who was diagnosed within 4 years with chest sarcoma, lung adenocarcinoma, and breast cancer. While her family history was unremarkable, DNA sequencing of TP53 revealed a germline de novo non-sense mutation in exon 6 p.Arg213X. One year later, she further developed a contralateral ductal carcinoma in situ, and 18 months later a jaw osteosarcoma. This case illustrates the therapeutic pitfalls in the care of a young cancer patient with TP53 de novo germline mutations and the complications related to her first-line therapy. Suggestion is made to use the less stringent Chompret criteria for germline TP53 mutation screening. Our observation underlines the possibly negative effect of radiotherapy in generating second tumors in patients with a TP53 mutation. We also present a review of six previously reported cases, comparing their cancer phenotypes with those generally produced by TP53 mutations.

Recurrent de novo mitochondrial DNA mutations in respiratory chain deficiency.

Starting from a cohort of 50 NADH-oxidoreductase (complex I) deficient patients, we carried out the systematic sequence analysis of all mitochondrially encoded complex I subunits (ND1 to ND6 and ND4L) in affected tissues. This approach yielded the unexpectedly high rate of 20% mutation identification in our series. Recurrent heteroplasmic mutations included two hitherto unreported (T10158C and T14487C) and three previously reported mutations (T10191C, T12706C and A13514G) in children with Leigh or Leigh-like encephalopathy. The recurrent mutations consistently involved T-->C transitions (p<10(-4)). This study supports the view that an efficient molecular screening should be based on an accurate identification of respiratory chain enzyme deficiency.

Molecular characterization of 39 de novo sSMC: contribution to prognosis and genetic counselling, a prospective study.

Small supernumerary marker chromosomes (sSMCs) are structurally abnormal chromosomes that cannot be characterized by karyotype. In many prenatal cases of de novo sSMC, the outcome of pregnancy is difficult to predict because the euchromatin content is unclear. This study aimed to determine the presence or absence of euchromatin material of 39 de novo prenatally ascertained sSMC by array-comparative genomic hybridization (array-CGH) or single nucleotide polymorphism (SNP) array. Cases were prospectively ascertained from the study of 65,000 prenatal samples [0.060%; 95% confidence interval (CI), 0.042-0.082]. Array-CGH showed that 22 markers were derived from non-acrocentric markers (56.4%) and 7 from acrocentic markers (18%). The 10 additional cases remained unidentified (25.6%), but 7 of 10 could be further identified using fluorescence in situ hybridization; 69% of de novo sSMC contained euchromatin material, 95.4% of which for non-acrocentric markers. Some sSMC containing euchromatin had a normal phenotype (31% for non-acrocentric and 75% for acrocentric markers). Statistical differences between normal and abnormal phenotypes were shown for the size of the euchromatin material (more or less than 1 Mb, p = 0.0006) and number of genes (more or less than 10, p = 0.0009). This study is the largest to date and shows the utility of array-CGH or SNP array in the detection and characterization of de novo sSMC in a prenatal context.

Garnet-chloritoid-kyanite metapelites from the Raspas Complex (SW Ecuador): a key eclogite-facies assemblage

The Raspas Complex (Ecuador) contains one of the few eclogitic bodies in the northern Andes. It consists of metaperidotites, eclogites, and metapelites. The latter display three assemblages: (i) garnet + chloritoid + kyanite, (ii) garnet + chloritoid and (iii) garnet + chlorite, in all cases with quartz and muscovite in addition. The growth of these assemblages was coeval with the main ductile deformation, and was followed by minor reequilibration (chlorite growth in garnet + chloritoid samples and chloritoid + quartz aggregates replacing garnet and kyanite in garnet + chloritoid + kyanite samples). Detailed microprobe analyses show increasing magnesian compositions for garnet (from core to rim) and chloritoid (inclusions within garnet compared to matrix grains) in kyanite-bearing samples. The above data are interpreted in the framework of the KFMASH system. Reaction progress along the divariant reaction Cld = Grt + Ky explains the change in chemistry of coexisting phases. The divariant Grt-Cld-Ky assemblage has a narrow stability field, and the P-T conditions are estimated at about 20 kbar, 550-600degreesC. Decompression, recorded by chloritoid-quartz pseudomorphs of garnet, probably occurred as temperature decreased.

Genome reconstruction and comparative genomics of pneumocystis jirovecii

Pneumocystis jirovecii is a fungus belonging to a basal lineage of the Ascomycotina, the Taphrinomycotina subphylum. It is a parasite specific to humans that dwells primarily in the lung and can cause severe pneumonia in individuals with debilitated immune system. Despite its clinical importance, many aspects of its biology remain poorly understood, at least in part because of the lack of a continuous in vitro cultivation system. The present thesis consists in the genome reconstruction and comparative genomics of P. jirovecii. It is made of three parts: (i) the de novo sequencing of P. jirovecii genome starting from a single broncho- alveolar lavage fluid of a single patient (ii) the de novo sequencing of the genome of the plant pathogen Taphrina deformans, a fungus closely related to P. jirovecii, and (iii) the genome scale comparison of P. jirovecii to other Taphrinomycotina members. Enrichment in P. jirovecii cells by immuno-precipitation, whole DNA random amplification, two complementary high throughput DNA sequencing methods, and in silico sorting and assembly of sequences were used for the de novo reconstruction of P. jirovecii genome from the microbiota of a single clinical specimen. An iterative ad hoc pipeline as well as numerical simulations was used to recover P. jirovecii sequences while purging out contaminants and assembly or amplification chimeras. This strategy produced a 8.1 Mb assembly, which encodes 3,898 genes. Homology searches, mapping on biochemical pathways atlases, and manual validations revealed that this genome lacks (i) most of the enzymes dedicated to the amino acids biosyntheses, and (ii) most virulence factors observed in other fungi, e.g. the glyoxylate shunt pathway and specific peptidases involved in the degradation of the host cell membrane. The same analyses applied to the available genomic sequences from Pneumocystis carinii the species infecting rats and Pneumocystis murina the species infecting mice revealed the same deficiencies. The genome sequencing of T. deformans yielded a 13 Mb assembly, which encodes 5,735 genes. T. deformans possesses enzymes involved plant cell wall degradation, secondary metabolism, the glyoxylate cycle, detoxification, sterol biosynthesis, as well as the biosyntheses of plant hormones such as abscisic acid or indole-3-acetic acid. T. deformans also harbors gene subsets that have counterparts in plant saprophytes or pathogens, which is consistent with its alternate saprophytic and pathogenic lifestyles. Mating genes were also identified. The homothallism of this fungus suggests a mating-type switching mechanism. Comparative analyses indicated that 81% of P. jirovecii genes are shared with eight other Taphrinomycotina members, including T. deformans, P. carinii and P. murina. These genes are mostly involved in housekeeping activities. The genes specific to the Pneumocystis genus represent 8%, and are involved in RNA metabolism and signaling. The signaling is known to be crucial for interaction of Pneumocystis spp with their environment. Eleven percent are unique to P. jirovecii and encode mostly proteins of unknown function. These genes in conjunction with other ones (e.g. the major surface glycoproteins) might govern the interaction of P. jirovecii with its human host cells, and potentially be responsible of the host specificity. P. jirovecii exhibits a reduced genome in size with a low GC content, and most probably scavenges vital compounds such as amino acids and cholesterol from human lungs. Consistently, its genome encodes a large set of transporters (ca. 22% of its genes), which may play a pivotal role in the acquisition of these compounds. All these features are generally observed in obligate parasite of various kingdoms (bacteria, protozoa, fungi). Moreover, epidemiological studies failed to evidence a free-living form of the fungus and Pneumocystis spp were shown to co-evolved with their hosts. Given also the lack of virulence factors, our observations strongly suggest that P. jirovecii is an obligate parasite specialized in the colonization of human lungs, and which causes disease only in individuals with compromised immune system. The same conclusion is most likely true for all other Pneumocystis spp in their respective mammalian host. - Pneumocystis jirovecii est un champignon appartenant à ine branche basale des Ascomycotina, le sous-embranchement des Taphrinomycotina. C'est un parasite spécifique aux humains qui réside principalement dans les poumons, et qui peut causer des pneumonies sévères chez des individus ayant un système immunitaire déficient. En dépit de son importance clinique, de nombreux aspects de sa biologie demeurent,largement méconnus, au moins en partie à cause de l'absence d'un système de culture in vitro continu. Cette thèse traite de la reconstruction du génome et de la génomique comparative de P. jirovecii. Elle comporte trois parties: (i) le séquençage de novo du génome de P. jirovecii à partir d'un lavage broncho-alvéolaire provenant d'un seul patient, (ii) le séquençage de novo du génome d'un champignon pathogène de plante Taphrina deformans qui est phylogénétiquement proche de P. jirovecii, et (iii) la comparaison du génome de P. jirovecii à celui d'autres membres du sous-embranchement des Taphrinomycotina. Un enrichissement en cellules de P. jirovecii par immuno-précipitation, une amplification aléatoire des molécules d'ADN, deux méthodes complémentaires de séquençage à haut débit, un tri in silico et un assemblage des séquences ont été utilisés pour reconstruire de novo le génome de P. jirovecii à partir du microbiote d'un seul échantillon clinique. Un pipeline spécifique ainsi que des simulations numériques ont été utilisés pour récupérer les séquences de P. jirovecii tout en éliminant les séquences contaminants et les chimères d'amplification ou d'assemblage. Cette stratégie a produit un assemblage de 8.1 Mb, qui contient 3898 gènes. Les recherches d'homologies, de cartographie des voies métaboliques et des validations manuelles ont révélé que ce génome est dépourvu (i) de la plupart des enzymes dédiées à la biosynthèse des acides aminés, et (ii) de la plupart des facteurs de virulence observés chez d'autres champignons, par exemple, le cycle du glyoxylate ainsi que des peptidases spécifiques impliquées dans la dégradation de la membrane de la cellule hôte. Les analyses appliquées aux données génomiques disponibles de Pneumocystis carinii, l'espèce infectant les rats, et de Pneumocystis murina, l'espèce infectant les souris, ont révélé les mêmes déficiences. Le séquençage du génome de T. deformans a généré un assemblage de 13.3 Mb qui contient 5735 gènes. T. deformans possède les gènes codant pour les enzymes impliquées dans la dégradation des parois cellulaires des plantes, le métabolisme secondaire, le cycle du glyoxylate, la détoxification, la biosynthèse des stérols ainsi que la biosynthèse d'hormones de plantes telles que l'acide abscissique ou l'acide indole 3-acétique. T. deformans possède également des sous-ensembles de gènes présents exclusivement chez des saprophytes ou des pathogènes de plantes, ce qui est consistent avec son mode de vie alternatif saprophyte et pathogène. Des gènes impliqués dans la conjugaison ont été identifiés. L'homothallisme de ce champignon suggère mécanisme de permutation du type conjuguant. Les analyses comparatives ont démontré que 81% des gènes de P. jirovecii sont présent chez les autres membres du sous-embranchement des Taphrinomycotina. Ces gènes sont essentiellement impliqués dans le métabolisme basai. Les gènes spécifiques au genre Pneumocystis représentent 8%, et sont impliqués dans le métabolisme de l'ARN et la signalisation. La signalisation est connue pour être cruciale pour l'interaction des espèces de Pneumocystis avec leur environnement. Les gènes propres à P. jirovecii représentent 11% et codent en majorité pour des protéines dont la fonction est inconnue. Ces gènes en conjonction avec d'autres (par exemple, les glycoprotéines de surface), pourraient être déterminants dans l'interaction de P. jirovecii avec les cellules de l'hôte humain, et être potentiellement responsable de la spécificité d'hôte. P. jirovecii possède un génome de taille réduite à faible pourcentage en GC et récupère très probablement des composés vitaux comme les acides aminés et le cholestérol à partir des poumons humains. De manière consistante, son génome code pour de nombreux transporteurs (22% de ses gènes), qui pourraient jouer un rôle essentiel dans l'acquisition de ces composés. Ces caractéristiques sont généralement observées chez les parasites obligatoires de plusieurs règnes (bactéries, protozoaires, champignons). De plus, les études épidémiologiques n'ont pas réussi à prouver l'existence d'ime forme vivant librement du champignon. Etant donné également l'absence de facteurs de virulence, nos observations suggèrent que P. jirovecii est un parasite obligatoire spécialisé dans la colonisation des poumons humains, ne causant une maladie que chez des individus ayant un système immunitaire compromis. La même conclusion est très probablement applicable à toutes les autres espèces de Pneumocystis dans leur hôte mammifère respectif.

De novo characterization of the Timema cristinae transcriptome facilitates marker discovery and inference of genetic divergence.

Adaptation to different ecological environments can promote speciation. Although numerous examples of such 'ecological speciation' now exist, the genomic basis of the process, and the role of gene flow in it, remains less understood. This is, at least in part, because systems that are well characterized in terms of their ecology often lack genomic resources. In this study, we characterize the transcriptome of Timema cristinae stick insects, a system that has been researched intensively in terms of ecological speciation, but for which genomic resources have not been previously developed. Specifically, we obtained >1 million 454 sequencing reads that assembled into 84,937 contigs representing approximately 18,282 unique genes and tens of thousands of potential molecular markers. Second, as an illustration of their utility, we used these genomic resources to assess multilocus genetic divergence within both an ecotype pair and a species pair of Timema stick insects. The results suggest variable levels of genetic divergence and gene flow among taxon pairs and genes and illustrate a first step towards future genomic work in Timema.

Visualization and quality assessment of de novo genome assemblies.

Recent technological progress has greatly facilitated de novo genome sequencing. However, de novo assemblies consist in many pieces of contiguous sequence (contigs) arranged in thousands of scaffolds instead of small numbers of chromosomes. Confirming and improving the quality of such assemblies is critical for subsequent analysis. We present a method to evaluate genome scaffolding by aligning independently obtained transcriptome sequences to the genome and visually summarizing the alignments using the Cytoscape software. Applying this method to the genome of the red fire ant Solenopsis invicta allowed us to identify inconsistencies in 7%, confirm contig order in 20% and extend 16% of scaffolds.Scripts that generate tables for visualization in Cytoscape from FASTA sequence and scaffolding information files are publicly available at https://github.com/ksanao/TGNet.