Comparison of Statistical Population Reconstruction Using Full and Pooled Adult Age-Class Data

BACKGROUND: Age-at-harvest data are among the most commonly collected, yet neglected, demographic data gathered by wildlife agencies. Statistical population construction techniques can use this information to estimate the abundance of wild populations over wide geographic areas and concurrently estimate recruitment, harvest, and natural survival rates. Although current reconstruction techniques use full age-class data (0.5, 1.5, 2.5, 3.5, … years), it is not always possible to determine an animal's age due to inaccuracy of the methods, expense, and logistics of sample collection. The ability to inventory wild populations would be greatly expanded if pooled adult age-class data (e.g., 0.5, 1.5, 2.5+ years) could be successfully used in statistical population reconstruction. METHODOLOGY/PRINCIPAL FINDINGS: We investigated the performance of statistical population reconstruction models developed to analyze full age-class and pooled adult age-class data. We performed Monte Carlo simulations using a stochastic version of a Leslie matrix model, which generated data over a wide range of abundance levels, harvest rates, and natural survival probabilities, representing medium-to-big game species. Results of full age-class and pooled adult age-class population reconstructions were compared for accuracy and precision. No discernible difference in accuracy was detected, but precision was slightly reduced when using the pooled adult age-class reconstruction. On average, the coefficient of variation (i.e., SE(θ)/θ) increased by 0.059 when the adult age-class data were pooled prior to analyses. The analyses and maximum likelihood model for pooled adult age-class reconstruction are illustrated for a black-tailed deer (Odocoileus hemionus) population in Washington State. CONCLUSIONS/SIGNIFICANCE: Inventorying wild populations is one of the greatest challenges of wildlife agencies. These new statistical population reconstruction models should expand the demographic capabilities of wildlife agencies that have already collected pooled adult age-class data or are seeking a cost-effective method for monitoring the status and trends of our wild resources

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation.

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence

A Novel Human Cytomegalovirus Locus Modulates Cell Type-Specific Outcomes of Infection

Clinical strains of HCMV encode 20 putative ORFs within a region of the genome termed ULb′ that are postulated to encode functions related to persistence or immune evasion. We have previously identified ULb′-encoded pUL138 as necessary, but not sufficient, for HCMV latency in CD34+ hematopoietic progenitor cells (HPCs) infected in vitro. pUL138 is encoded on polycistronic transcripts that also encode 3 additional proteins, pUL133, pUL135, and pUL136, collectively comprising the UL133-UL138 locus. This work represents the first characterization of these proteins and identifies a role for this locus in infection. Similar to pUL138, pUL133, pUL135, and pUL136 are integral membrane proteins that partially co-localized with pUL138 in the Golgi during productive infection in fibroblasts. As expected of ULb′ sequences, the UL133-UL138 locus was dispensable for replication in cultured fibroblasts. In CD34+ HPCs, this locus suppressed viral replication in HPCs, an activity attributable to both pUL133 and pUL138. Strikingly, the UL133-UL138 locus was required for efficient replication in endothelial cells. The association of this locus with three context-dependent phenotypes suggests an exciting role for the UL133-UL138 locus in modulating the outcome of viral infection in different contexts of infection. Differential profiles of protein expression from the UL133-UL138 locus correlated with the cell-type dependent phenotypes associated with this locus. We extended our in vitro findings to analyze viral replication and dissemination in a NOD-scid IL2Rγcnull-humanized mouse model. The UL133-UL138NULL virus exhibited an increased capacity for replication and/or dissemination following stem cell mobilization relative to the wild-type virus, suggesting an important role in viral persistence and spread in the host. As pUL133, pUL135, pUL136, and pUL138 are conserved in virus strains infecting higher order primates, but not lower order mammals, the functions encoded likely represent host-specific viral adaptations

A Comprehensive Genome-Wide Map of Autonomously Replicating Sequences in a Naive Genome

Eukaryotic chromosomes initiate DNA synthesis from multiple replication origins. The machinery that initiates DNA synthesis is highly conserved, but the sites where the replication initiation proteins bind have diverged significantly. Functional comparative genomics is an obvious approach to study the evolution of replication origins. However, to date, the Saccharomyces cerevisiae replication origin map is the only genome map available. Using an iterative approach that combines computational prediction and functional validation, we have generated a high-resolution genome-wide map of DNA replication origins in Kluyveromyces lactis. Unlike other yeasts or metazoans, K. lactis autonomously replicating sequences (KlARSs) contain a 50 bp consensus motif suggestive of a dimeric structure. This motif is necessary and largely sufficient for initiation and was used to dependably identify 145 of the up to 156 non-repetitive intergenic ARSs projected for the K. lactis genome. Though similar in genome sizes, K. lactis has half as many ARSs as its distant relative S. cerevisiae. Comparative genomic analysis shows that ARSs in K. lactis and S. cerevisiae preferentially localize to non-syntenic intergenic regions, linking ARSs with loci of accelerated evolutionary change

Evidence for Sequential and Increasing Activation of Replication Origins along Replication Timing Gradients in the Human Genome

Genome-wide replication timing studies have suggested that mammalian chromosomes consist of megabase-scale domains of coordinated origin firing separated by large originless transition regions. Here, we report a quantitative genome-wide analysis of DNA replication kinetics in several human cell types that contradicts this view. DNA combing in HeLa cells sorted into four temporal compartments of S phase shows that replication origins are spaced at 40 kb intervals and fire as small clusters whose synchrony increases during S phase and that replication fork velocity (mean 0.7 kb/min, maximum 2.0 kb/min) remains constant and narrowly distributed through S phase. However, multi-scale analysis of a genome-wide replication timing profile shows a broad distribution of replication timing gradients with practically no regions larger than 100 kb replicating at less than 2 kb/min. Therefore, HeLa cells lack large regions of unidirectional fork progression. Temporal transition regions are replicated by sequential activation of origins at a rate that increases during S phase and replication timing gradients are set by the delay and the spacing between successive origin firings rather than by the velocity of single forks. Activation of internal origins in a specific temporal transition region is directly demonstrated by DNA combing of the IGH locus in HeLa cells. Analysis of published origin maps in HeLa cells and published replication timing and DNA combing data in several other cell types corroborate these findings, with the interesting exception of embryonic stem cells where regions of unidirectional fork progression seem more abundant. These results can be explained if origins fire independently of each other but under the control of long-range chromatin structure, or if replication forks progressing from early origins stimulate initiation in nearby unreplicated DNA. These findings shed a new light on the replication timing program of mammalian genomes and provide a general model for their replication kinetics