Currently a few tools are capable of detecting genome-wide Copy Number Variations (CNVs) based on sequencing of multiple samples. Although aberrations in mate pair insertion sizes provide additional hints for the CNV detection based on multiple samples, the majority of the current tools rely only on the depth of coverage. Here, we propose a new algorithm (MSeq-CNV) which allows detecting common CNVs across multiple samples. MSeq-CNV applies a mixture density for modeling aberrations in depth of coverage and abnormalities in the mate pair insertion sizes. Each component in this mixture density applies a Binomial distribution for modeling the number of mate pairs with aberration in the insertion size and also a Poisson distribution for emitting the read counts, in each genomic position. MSeq-CNV is applied on simulated data and also on real data of six HapMap individuals with high-coverage sequencing, in 1000 Genomes Project. These individuals include a CEU trio of European ancestry and a YRI trio of Nigerian ethnicity. Ancestry of these individuals is studied by clustering the identified CNVs. MSeq-CNV is also applied for detecting CNVs in two samples with low-coverage sequencing in 1000 Genomes Project and six samples form the Simons Genome Diversity Project.

No abstract

Previous local and national Iranian publications indicate that all Iranian hepatitis B virus (HBV) strains belong to HBV genotype D. The aim of this study was to analyze the evolutionary history of HBV infection in Iran for the first time, based on an intensive phylodynamic study. The evolutionary parameters, time to most recent common ancestor (tMRCA), and the population dynamics of infections were investigated using the Bayesian Monte Carlo Markov chain (BMCMC). The effective sample size (ESS) and sampling convergence were then monitored. After sampling from the posterior distribution of the nucleotide substitution rate and other evolutionary parameters, the point estimations (median) of these parameters were obtained. All Iranian HBV isolates were of genotype D, sub-type ayw2. The origin of HBV is regarded as having evolved first on the eastern border, before moving westward, where Isfahan province then hosted the virus. Afterwards, the virus moved to the south and west of the country. The tMRCA of HBV in Iran was estimated to be around 1894, with a 95% credible interval between the years 1701 and 1957. The effective number of infections increased exponentially from around 1925 to 1960. Conversely, from around 1992 onwards, the effective number of HBV infections has decreased at a very high rate. Phylodynamic inference clearly demonstrates a unique homogenous pattern of HBV genotype D compatible with a steady configuration of the decreased effective number of infections in the population in recent years, possibly due to the implementation of blood donation screening and vaccination programs. Adequate molecular epidemiology databases for HBV are crucial for infection prevention and treatment programs.

The high rate of hepatitis C virus (HCV) infection among transfusion related risk groups such as patients with inherited bleeding disorders highlighting the investigation on prevalent subtypes and their epidemic history among this group. In this study, 166 new HCV NS5B sequences isolated from patients with inherited bleeding disorders together with 29 sequences related to hemophiliacs obtained from a previous study on diversity of HCV in Iran were analyzed. The most prevalent subtype was 1a (65%), followed by 3a (18.7%),1b (14.5%),4(1.2%) and 2k (0.6%). Subtypes 1a and 3a showed exponential expansion during the 20th century. Whereas expansion of 3a started around 20 years earlier than 1a among the study patients, the epidemic growth of 1a revealed a delay of about 10 years compared with that found for this subtype in developed countries. Our results supported the view that the spread of 3a reached the plateau 10 years prior to the screening of blood donors for HCV. Rather, 1a reached the plateau when screening program was implemented. The differences observed in the epidemic behavior of HCV-1a and 3a may be associated with different transmission routes of two subtypes. Indeed, expansion of 1a was more commonly linked to blood transfusion, while 3a was more strongly associated to drug use and specially IDU after 1960. Our findings also showed HCV transmission through blood products has effectively been controlled from late 1990s. In conclusion, the implementation of strategies such as standard surveillance programs and subsiding antiviral treatments seems to be essential to both prevent new HCV infections and to decline the current and future HCV disease among Iranian patients with inherited bleeding disorders.