Search

Abstract While many countries employed digital contact tracing to contain the spread of SARS-CoV-2, the contribution of cospace-time interaction (i.e., individuals who shared the same space and time) to transmission and to super-spreading in the real world has seldom been systematically studied due to the lack of systematic sampling and testing of contacts. To address this issue, we utilized data from 2230 cases and 220,878 contacts with detailed epidemiological information during the Omicron outbreak in Beijing in 2022. We observed that contact number per day of tracing for individuals in dwelling, workplace, cospace-time interactions, and community settings could be described by gamma distribution with distinct parameters. Our findings revealed that 38% of traced transmissions occurred through cospace-time interactions whilst control measures were in place. However, using a mathematical model to incorporate contacts in different locations, we found that without control measures, cospace-time interactions contributed to only 11% (95%CI: 10%–12%) of transmissions and the super-spreading risk for this setting was 4% (95%CI: 3%–5%), both the lowest among all settings studied. These results suggest that public health measures should be optimized to achieve a balance between the benefits of digital contact tracing for cospace-time interactions and the challenges posed by contact tracing within the same setting.

Abstract African swine fever virus (ASFV) causes a devastating disease affecting domestic and wild pigs. ASF was first introduced in Sardinia in 1978 and until 2019 only genotype I isolates were identified. A remarkable genetic stability of Sardinian ASFV isolates was described, nevertheless in 2019 two wild boar isolates with a sustained genomic deletion (4342 base pairs) were identified (7303WB/19, 7212WB/19). In this study, we therefore performed in vitro experiments with monocyte-derived macrophages (moMФ) to unravel the phenotypic characteristics of these deleted viruses. Both 7303WB/19 and 7212WB/19 presented a lower growth kinetic in moMФ compared to virulent Sardinian 26544/OG10, using either a high (1) or a low (0.01) multiplicity of infection (MOI). In addition, flow cytometric analysis showed that both 7303WB/19 and 7212WB/19 presented lower intracellular levels of both early and late ASFV proteins. We subsequently investigated whether deleted virus variants were previously circulating in wild boars in Sardinia. In the four years preceding the last genotype I isolation (February 2015–January 2019), other eight wild boar isolates were collected, all belonging to p72 genotype I, B602L subgroup X, but none of them presented a sustained genomic deletion. Overall, we observed the deleted virus isolates in Sardinia only in 2019, at the end of a strong eradication campaign, and our data suggest that it might possess an attenuated phenotype in vivo. A better understanding of ASFV evolution in endemic territories might contribute to development of effective control measures against ASF.

In agriculture, new and sustainable strategies are increasingly demanded to integrate the traditional management of viral diseases based on the use of virus-free propagation materials and resistant or tolerant cultivars and on the control of insect vectors. Among the possible Integrated Pest Management (IPM) approaches, organic biostimulants have shown promising results in enhancing plant tolerance to virus infections by improving plant fitness and productivity and modulating metabolic functions. In this study, the combination of two organic biostimulants, Alert D-Max and Resil EVO Q, composed of seaweed and alfalfa extracts, enzymatic hydrolysates, and micronized zeolite, was applied on the leaves and roots of zucchini squashes, both healthy and infected by zucchini yellow mosaic virus (ZYMV). Four applications were scheduled based on ZYMV inoculation timing, and plant vegetative and reproductive parameters were recorded along with the virus titre and symptom severity. The modulation of the expression of specific genes potentially involved in pattern-triggered immunity (PTI), systemic acquired resistance (SAR), and oxidative stress defence pathways was also investigated. Besides increasing the general fitness of the healthy plants, the biostimulants significantly improved the production of flowers and fruits of the infected plants, with a potential positive impact on their productivity. The repeated biostimulant applications also led to a one-tenth reduction in ZYMV titre over time and induced a progressive slowdown of symptom severity. Genes associated with SAR and PTI were up-regulated after biostimulant applications, suggesting the biostimulant-based priming of plant defence mechanisms. Due to the observed beneficial effects, the tested biostimulant mix can be an effective component of the IPM of cucurbit crops, acting as a sustainable practice for enhancing plant fitness and tolerance to potyviruses.

Canine distemper virus (CDV) is the etiological agent of a highly prevalent viral infectious disease of carnivores, which could seriously lead to a threat to the conservation of the affected species worldwide [...]

Conservation, characterization and exploitation of agrobiodiversity are key factors to guarantee food security and face future challenges such as climate changes. These issues are the subject of a series of international agreements, such as the Convention of Biological Diversity, with its Nagoya Protocol, and the International Treaty on Plant Genetic Resources for Food and Agriculture (ITPGRFA) adopted in 2001 and entered into force in 2004. Italy ratified the Treaty in 2004 and instituted a long-lasting program, RGV-FAO, to implement it. CREA is one of the three organizations involved in the RGV-FAO Program, together with the National Research Council (CNR) and Reti Semi Rurali. CREA maintains a total of 40,186 accessions including cereals, vegetables, fruits, forages, industrial crops, forest and woody crops, medicinal and aromatic plants, and their wild relatives. Accessions are conserved using different ex situ conservation systems (seeds, in vivo plants, vegetative organs and in vitro plantlets), and characterized using genetic, morpho-phenological and/or biochemical methods. Herein, we will present the CREA long-lasting program RGV-FAO with some examples of the use of plant genetic resources in breeding programs, including molecular approaches. Some critical issues related to access and benefit sharing in PGRFA, such as the Nagoya Protocol and the Digital Sequence Information, will be discussed, highlighting their potential impact on food security and on the advancement of knowledge.

Porcine parvoviruses (PPVs) are among the most important agents of reproductive failure in swine worldwide. PPVs comprise eight genetically different species ascribed to four genera: Protoparvovirus (PPV1, PPV8), Tetraparvovirus (PPV2-3), Copiparvovirus (PPV4-6), and Chaphamaparvovirus (PPV7). In 2016, PPV7 was firstly detected in the USA and afterwards in Europe, Asia, and South America. Recently, it was also identified in Italy in pig farms with reproductive failure. This study aimed to evaluate the circulation of PPV7 in domestic and wild pigs in Sardinia, Italy. In addition, its coinfection with Porcine Circovirus 2 (PCV2) and 3 (PCV3) was analysed, and PPV7 Italian strains were molecularly characterised. PPV7 was detected in domestic pigs and, for the first time, wild pigs in Italy. The PPV7 viral genome was detected in 20.59% of domestic and wild pig samples. PPV7 detection was significantly lower in domestic pigs, with higher PCV2/PCV3 co-infection rates observed in PPV7-positive than in PPV7-negative domestic pigs. Molecular characterisation of the NS1 gene showed a very high frequency of recombination that could presumably promote virus spreading.

[This corrects the article DOI: 10.1371/journal.ppat.1011461.].

The NHS COVID-19 digital contact tracing app was designed to notify people of potential exposure to SARS-CoV-2. Here, the authors summarise the uptake and engagement with the app in its first year, and estimate its epidemiological impact in terms of numbers of cases, hospitalisations, and deaths averted.

In this study, we evaluated the impact of viral variant, in addition to other variables, on within-host viral burden, by analysing cycle threshold (Ct) values derived from nose and throat swabs, collected as part of the UK COVID-19 Infection Survey. Because viral burden distributions determined from community survey data can be biased due to the impact of variant epidemiology on the time-since-infection of samples, we developed a method to explicitly adjust observed Ct value distributions to account for the expected bias. By analysing the adjusted Ct values using partial least squares regression, we found that among unvaccinated individuals with no known prior exposure, viral burden was 44% lower among Alpha variant infections, compared to those with the predecessor strain, B.1.177. Vaccination reduced viral burden by 67%, and among vaccinated individuals, viral burden was 286% higher among Delta variant, compared to Alpha variant, infections. In addition, viral burden increased by 17% for every 10-year age increment of the infected individual. In summary, within-host viral burden increases with age, is reduced by vaccination, and is influenced by the interplay of vaccination status and viral variant.

The standard neutral model of molecular evolution has traditionally been used as the null model for population genomics. We gathered a collection of 45 genome-wide site frequency spectra from a diverse set of species, most of which display an excess of low and high frequency variants compared to the expectation of the standard neutral model, resulting in U-shaped spectra. We show that multiple merger coalescent models often provide a better fit to these observations than the standard Kingman coalescent. Hence, in many circumstances these under-utilized models may serve as the more appropriate reference for genomic analyses. We further discuss the underlying evolutionary processes that may result in the widespread U-shape of frequency spectra.

Testing capacity continues to limit detection of COVID-19 infections and impacts reliability of mortality estimates. Here, the authors develop a statistical model to estimate COVID-19 attributable deaths using all-cause mortality data from Iran and estimate that around half of these deaths have been reported.

Abstract Cell lines allow studying various biological processes that may not be easily tractable in whole organisms. Here, we have established the first male-specific cell line from the African malaria mosquito, Anopheles gambiae. The cells, named AgMM and derived from the sex-sorted neonate larvae, were able to undergo spontaneous contractions for a number of passages following establishment, indicating their myoblast origin. Comparison of their transcriptome to the transcriptome of an A. gambiae-derived Sua5.1 hemocyte cells revealed distinguishing molecular signatures of each cell line, including numerous muscle-related genes that were highly and uniquely expressed in the AgMM cells. Moreover, the AgMM cells express the primary sex determiner gene Yob and support male sex determination and dosage compensation pathways. Therefore, the AgMM cell line represents a valuable tool for molecular and biochemical in vitro studies of these male-specific processes. In a broader context, a rich transcriptomic data set generated in this study contributes to a better understanding of transcribed regions of the A. gambiae genome and sheds light on the biology of both cell types, facilitating their anticipated use for various cell-based assays.

Hop (Humulus lupulus L.) is a minor ingredient in the beer production but has a strong influence on the beer quality due to the high chemical complexity of the cones used in brewing. One of the major factors that can severely affect the chemical composition of the hop cones and their marketability is the presence of viral infections in the plant. Amongst the five major hop viruses, three belong to the Carlavirus genus: hop mosaic virus (HpMV), hop latent virus (HpLV), and American hop latent virus (AHLV). The occurrence of carlaviruses on hop germplasm in Italy was firstly recorded in 2017 but, in that context, a generic detection was only performed and no information on the infecting Carlavirus species was provided. To fill this gap, 51 hop samples previously found infected by carlaviruses were analysed by RT-PCR employing primer pairs specific for the coat protein (CP) of HpMV, HpLV and AHLV, respectively. HpLV resulted largely prevalent as it was detected in 96.1% of tested samples whereas for HpMV and AHLV an infection rate of 5.9% and 3.9% was recorded, respectively. CP nucleotide sequences from 13 selected virus isolates were obtained and analysed; moreover, the complete genome sequence of 7 isolates was obtained by using high throughput sequencing (HTS). Phylogenetic analysis showed close relationships among isolates from different geographical origin, including European and non-European countries, according to the worldwide movement of hop germplasm due to global trade. This is the first report of HpMV, HpLV and AHLV on hop germplasm in Italy.

Summary: Insights into the evolution of non-model organisms are limited by the lack of reference genomes of high accuracy, completeness, and contiguity. Here, we present a chromosome-level, karyotype-validated reference genome and pangenome for the barn swallow (Hirundo rustica). We complement these resources with a reference-free multialignment of the reference genome with other bird genomes and with the most comprehensive catalog of genetic markers for the barn swallow. We identify potentially conserved and accelerated genes using the multialignment and estimate genome-wide linkage disequilibrium using the catalog. We use the pangenome to infer core and accessory genes and to detect variants using it as a reference. Overall, these resources will foster population genomics studies in the barn swallow, enable detection of candidate genes in comparative genomics studies, and help reduce bias toward a single reference genome.

Among the cucurbit-infecting viruses, watermelon mosaic virus (WMV) and zucchini yellow mosaic virus (ZYMV) (Potyvirus: Potyviridae) are responsible for severe symptoms on cucumber, melon, watermelon, and zucchini cultivations worldwide. In this study, reverse transcription real-time PCR (real-time RT-PCR) and droplet-digital PCR (RT-ddPCR) assays targeting the coat protein (CP) genes of WMV and ZYMV were developed and validated according to the international standards of plant pest diagnosis (EPPO PM 7/98 (5)). First, the diagnostic performance of WMV-CP and ZYMV-CP real-time RT-PCRs was evaluated, and the assays displayed an analytical sensitivity of 10−5 and 10−3, respectively. The tests also showed an optimal repeatability, reproducibility and analytical specificity, and were reliable for the virus detection in naturally infected samples and across a wide range of cucurbit hosts. Based on these results, the real-time RT-PCR reactions were adapted to set up RT-ddPCR assays. These were the first RT-ddPCR assays aiming at the detection and quantification of WMV and ZYMV and showed a high sensitivity, being able to detect until 9 and 8 copies/µL of WMV or ZYMV, respectively. The RT-ddPCRs allowed the direct estimation of the virus concentrations and opened to a broad range of applications in disease management, such as the evaluation of partial resistance in breeding processes, identification of antagonistic/synergistic events, and studies on the implementation of natural compounds in the integrated management strategies.

Abstract Background During the Neolithic expansion, cattle accompanied humans and spread from their domestication centres to colonize the ancient world. In addition, European cattle occasionally intermingled with both indicine cattle and local aurochs resulting in an exclusive pattern of genetic diversity. Among the most ancient European cattle are breeds that belong to the so-called Podolian trunk, the history of which is still not well established. Here, we used genome-wide single nucleotide polymorphism (SNP) data on 806 individuals belonging to 36 breeds to reconstruct the origin and diversification of Podolian cattle and to provide a reliable scenario of the European colonization, through an approximate Bayesian computation random forest (ABC-RF) approach. Results Our results indicate that European Podolian cattle display higher values of genetic diversity indices than both African taurine and Asian indicine breeds. Clustering analyses show that Podolian breeds share close genomic relationships, which suggests a likely common genetic ancestry. Among the simulated and tested scenarios of the colonization of Europe from taurine cattle, the greatest support was obtained for the model assuming at least two waves of diffusion. Time estimates are in line with an early migration from the domestication centre of non-Podolian taurine breeds followed by a secondary migration of Podolian breeds. The best fitting model also suggests that the Italian Podolian breeds are the result of admixture between different genomic pools. Conclusions This comprehensive dataset that includes most of the autochthonous cattle breeds belonging to the so-called Podolian trunk allowed us not only to shed light onto the origin and diversification of this group of cattle, but also to gain new insights into the diffusion of European cattle. The most well-supported scenario of colonization points to two main waves of migrations: with one that occurred alongside with the Neolithic human expansion and gave rise to the non-Podolian taurine breeds, and a more recent one that favoured the diffusion of European Podolian. In this process, we highlight the importance of both the Mediterranean and Danube routes in promoting European cattle colonization. Moreover, we identified admixture as a driver of diversification in Italy, which could represent a melting pot for Podolian cattle.

The tiger mosquito (Aedes albopictus) is one of the most invasive species in the world and a competent vector for numerous arboviruses, thus the study and monitoring of its fast worldwide spread is crucial for global public health. The small extra-nuclear and maternally-inherited mitochondrial DNA represents a key tool for reconstructing phylogenetic and phylogeographic relationships within a species, especially when analyzed at the mitogenome level. Here the mitogenome variation of 76 tiger mosquitoes, 37 of which new and collected from both wild adventive populations and laboratory strains, was investigated. This analysis significantly improved the global mtDNA phylogeny of Ae. albopictus, uncovering new branches and sub-branches within haplogroup A1, the one involved in its recent worldwide spread. Our phylogeographic approach shows that the current distribution of tiger mosquito mitogenome variation has been strongly affected by clonal and sub-clonal founder events, sometimes involving wide geographic areas, even across continents, thus shedding light on the Asian sources of worldwide adventive populations. In particular, different starting points for the two major clades within A1 are suggested, with A1a spreading mainly along temperate areas from Japanese and Chinese sources, and A1b arising and mainly diffusing in tropical areas from a South Asian source.

In the last decades, the interest in biological activity of natural compounds has been growing. In plant protection, essential oils have been reported to exhibit antiviral, antimycotic, and antiparasitic activities, and are regarded as promising for the formulation of safe antimicrobial agents. Attention has also been focused on hydrosols, the by-products of hydro-distillation of essential oils. Their production is easy, fast, and cheap, and they seem to arise less concern for human health than essential oils. Plant viruses represent a major concern for agricultural crops since no treatment compound is available for virus control. This work was aimed at evaluating the antiphytoviral effectiveness of treatments with three essential oils and corresponding hydrosols extracted from Origanum vulgare, Thymus vulgaris, and Rosmarinus officinalis on Cucurbita pepo plants infected by zucchini yellow mosaic virus or tomato leaf curl New Delhi virus. Treatments were applied either concurrently or after virus inoculation to ascertain an inhibition or curative activity, respectively. Symptoms were observed and samplings were performed weekly. Virus titer and expression levels of phenylalanine ammonia lyase gene (PAL) were measured on treated and untreated infected plants by real-time PCR. PAL gene plays an important role in plant defense response as it is involved in tolerance/resistance to phytopathogens. Results indicated that treatments were effective against tomato leaf curl New Delhi virus whether applied simultaneously with the inoculation or after. A major inhibition was observed with O. vulgare essential oil and hydrosol, resulting in 10–4-fold decrease of virus titer 3 weeks after treatment. Curative activity gave maximum results with all three essential oils and T. vulgaris and R. officinalis hydrosols, recording from 10–2-fold decrease to virus not detected 4 weeks after treatment. An induction of PAL gene expression was recorded at 12 d.p.i. and then was restored to the levels of untreated control. This allows to hypothesize an early plant defense response to virus infection, possibly boosted by treatments. Plant extracts’ composition was characterized by gas chromatography-mass spectrometry. Phenols were largely main components of O. vulgare and T. vulgaris extracts (carvacrol and thymol, respectively), while extracts from R. officinalis were based on monoterpene hydrocarbons (essential oil) and oxygenated monoterpenes (hydrosol).

In recent years, natural compounds have gained attention in many fields due to their wide-range biological activity. In particular, essential oils and their associated hydrosols are being screened to control plant pests, exerting antiviral, antimycotic and antiparasitic actions. They are more quickly and cheaply produced and are generally considered safer for the environment and non-target organisms than conventional pesticides. In this study, we report the evaluation of the biological activity of two essential oils and their corresponding hydrosols obtained from Mentha suaveolens and Foeniculum vulgare in the control of zucchini yellow mosaic virus and its vector, Aphis gossypii, in Cucurbita pepo plants. The control of the virus was ascertained with treatments applied either concurrently with or after virus infection; choice tests were performed to verify repellency activity against the aphid vector. The results indicated that treatments could decrease virus titer as measured using real-time RT-PCR, while the experiments on the vector showed that the compounds effectively repelled aphids. The extracts were also chemically characterized using gas chromatography–mass spectrometry. Mentha suaveolens and Foeniculum vulgare hydrosol extracts mainly comprised fenchone and decanenitrile, respectively, while essential oils analysis returned a more complex composition, as expected.

Understanding how geography and human mobility shape the patterns and spread of infectious diseases such as COVID-19 is key to control future epidemics. An interesting example is provided by the second wave of the COVID-19 epidemic in Europe, which was facilitated by the intense movement of tourists around the Mediterranean coast in summer 2020. The Italian island of Sardinia is a major tourist destination and is widely believed to be the origin of the second Italian wave. In this study, we characterize the genetic variation among SARS-CoV-2 strains circulating in northern Sardinia during the first and second Italian waves using both Illumina and Oxford Nanopore Technologies Next Generation Sequencing methods. Most viruses were placed into a single clade, implying that despite substantial virus inflow, most outbreaks did not spread widely. The second epidemic wave on the island was actually driven by local transmission of a single B.1.177 subclade. Phylogeographic analyses further suggest that those viral strains circulating on the island were not a relevant source for the second epidemic wave in Italy. This result, however, does not rule out the possibility of intense mixing and transmission of the virus among tourists as a major contributor to the second Italian wave.

Introduction: Many countries with an early outbreak of SARS-CoV-2 struggled to gauge the size and start date of the epidemic mainly due to limited testing capacities and a large proportion of undetected asymptomatic and mild infections. Iran was among the first countries with a major outbreak outside China. Methods: We constructed a globally representative sample of 802 genomes, including 46 samples from patients inside or with a travel history to Iran. We then performed a phylogenetic analysis to identify clades related to samples from Iran and estimated the start of the epidemic and early doubling times in cases. We leveraged air travel data from 36 exported cases of COVID-19 to estimate the point-prevalence and the basic reproductive number across the country. We also analysed the province-level all-cause mortality data during winter and spring 2020 to estimate under-reporting of COVID-19-related deaths. Finally, we use this information in an SEIR model to reconstruct the early outbreak dynamics and assess the effectiveness of intervention measures in Iran. Results: By identifying the most basal clade that contained genomes from Iran, our phylogenetic analysis showed that the age of the root is placed on 2019-12-21 (95 % HPD: 2019-09-07 – 2020-02-14). This date coincides with our estimated epidemic start date on 2019-12-25 (95 %CI: 2019-12-11 – 2020-02-24) based air travel data from exported cases with an early doubling time of 4.0 (95 %CI: 1.4–6.7) days in cases. Our analysis of all-cause mortality showed 21.9 (95 % CI: 16.7–27.2) thousand excess deaths by the end of summer. Our model forecasted the second epidemic peak and suggested that by 2020-08-31 a total of 15.0 (95 %CI: 4.9–25.0) million individuals recovered from the disease across the country. Conclusion: These findings have profound implications for assessing the stage of the epidemic in Iran despite significant levels of under-reporting. Moreover, the results shed light on the dynamics of SARS-CoV-2 transmissions in Iran and central Asia. They also suggest that in the absence of border screening, there is a high risk of introduction from travellers from areas with active outbreaks. Finally, they show both that well-informed epidemic models are able to forecast episodes of resurgence following a relaxation of interventions, and that NPIs are key to controlling ongoing epidemics.

African swine fever virus (ASFV) is the etiological agent of a lethal disease of domestic pigs and wild boars. ASF threatens the pig industry worldwide due to the lack of a licensed vaccine or treatment. The disease has been endemic for more than 40 years in Sardinia (Italy), but an intense campaign pushed it close to eradication; virus circulation was last detected in wild boars in 2019. In this study, we present a genomic analysis of two ASFV strains isolated in Sardinia from two wild boars during the 2019 hunting season. Both isolates presented a deletion of 4342 base pairs near the 5′ end of the genome, encompassing the genes MGF 360-6L, X69R, and MGF 300-1L. The phylogenetic evidence suggests that the deletion recently originated within the Sardinia ecosystem and that it is most likely the result of a non-allelic homologous recombination driven by a microhomology present in most Sardinian ASFV genomes. These results represent a striking example of a genomic feature promoting the rapid evolution of structural variations and plasticity in the ASFV genome. They also raise interesting questions about the functions of the deleted genes and the potential link between the evolutionary timing of the deletion appearance and the eradication campaign.

SARS-CoV-2 has spread across the world, causing high mortality and unprecedented restrictions on social and economic activity. Policymakers are assessing how best to navigate through the ongoing epidemic, with computational models being used to predict the spread of infection and assess the impact of public health measures. Here, we present OpenABM-Covid19: an agent-based simulation of the epidemic including detailed age-stratification and realistic social networks. By default the model is parameterised to UK demographics and calibrated to the UK epidemic, however, it can easily be re-parameterised for other countries. OpenABM-Covid19 can evaluate non-pharmaceutical interventions, including both manual and digital contact tracing, and vaccination programmes. It can simulate a population of 1 million people in seconds per day, allowing parameter sweeps and formal statistical model-based inference. The code is open-source and has been developed by teams both inside and outside academia, with an emphasis on formal testing, documentation, modularity and transparency. A key feature of OpenABM-Covid19 are its Python and R interfaces, which has allowed scientists and policymakers to simulate dynamic packages of interventions and help compare options to suppress the COVID-19 epidemic.

Inversions are a little-studied type of genomic variation that could contribute to phenotypic traits. Here the authors characterize 45 common polymorphic inversions in human populations and investigate their evolutionary and functional impact.

Abstract Background Recent genome studies of modern and ancient samples have proposed that Native Americans derive from a subset of the Eurasian gene pool carried to America by an ancestral Beringian population, from which two well-differentiated components originated and subsequently mixed in different proportion during their spread in the Americas. To assess the timing, places of origin and extent of admixture between these components, we performed an analysis of the Y-chromosome haplogroup Q, which is the only Pan-American haplogroup and accounts for virtually all Native American Y chromosomes in Mesoamerica and South America. Results Our analyses of 1.5 Mb of 152 Y chromosomes, 34 re-sequenced in this work, support a “coastal and inland routes scenario” for the first entrance of modern humans in North America. We show a major phase of male population growth in the Americas after 15 thousand years ago (kya), followed by a period of constant population size from 8 to 3 kya, after which a secondary sign of growth was registered. The estimated dates of the first expansion in Mesoamerica and the Isthmo-Colombian Area, mainly revealed by haplogroup Q-Z780, suggest an entrance in South America prior to 15 kya. During the global constant population size phase, local South American hints of growth were registered by different Q-M848 sub-clades. These expansion events, which started during the Holocene with the improvement of climatic conditions, can be ascribed to multiple cultural changes rather than a steady population growth and a single cohesive culture diffusion as it occurred in Europe. Conclusions We established and dated a detailed haplogroup Q phylogeny that provides new insights into the geographic distribution of its Eurasian and American branches in modern and ancient samples.

Summary: Background: In May 2020, the UK National Health Service (NHS) Test and Trace programme was launched in England in response to the COVID-19 pandemic. The programme was first rolled out on the Isle of Wight and included version 1 of the NHS contact tracing app. The aim of the study was to make a preliminary assessment of the epidemiological impact of the Test and Trace programme using publicly available data. Methods: We used COVID-19 daily case data from Public Health England to infer incidence of new infections and estimate the reproduction number (R) for each of the 150 Upper-Tier Local Authorities (UTLAs) in England and nationally, before and after the launch of the Test and Trace programme on the Isle of Wight. We used Bayesian and maximum-likelihood methods to estimate R and compared the Isle of Wight with other UTLAs using a synthetic control method. Findings: We observed significant decreases in incidence and R on the Isle of Wight immediately after the launch of the Test and Trace programme. The Isle of Wight had a marked reduction in R, from 1·3 before the Test and Trace programme to 0·5 after by one of our measures, and went from having the third highest R before the Test and Trace programme, to the twelfth lowest afterwards compared with other UTLAs. Interpretation: Our results show that the epidemic on the Isle of Wight was controlled quickly and effectively after the launch of Test and Trace. Our findings highlight the need for further research to determine the causes of the reduction in the spread of the disease, as these could be translated into local and national non-pharmaceutical intervention strategies in the period before a treatment or vaccination for COVID-19 becomes available. Funding: Li Ka Shing Foundation and UK Economic and Social Research Council.

The use of organic substances in integrated pest management can contribute to human- and environment-safe crop production. In the present work, a combination of organic biostimulants (Fullcrhum Alert and BioVeg 500) and an inorganic corroborant (Clinogold, zeolite) was tested for the effects on the plant response to the quarantine pest tomato leaf curl New Delhi virus (ToLCNDV). Biostimulants were applied to healthy and infected greenhouse-grown zucchini plants, and the vegetative parameters and viral titer were evaluated. Although no antiviral effects were observed in terms of both virus replication and symptom expression, these biostimulants were shown to influence plant fitness. A significant increase in biomass and in leaf, flower, and fruit production was induced in both healthy and infected plants. Biostimulants also enhanced the production of metabolites commonly involved in plant response to virus infection, such as carbohydrates, phenylpropanoids and free amino acids. These results encourage new field trials to evaluate the actual productivity of infected plants after treatments and the possible application of organic biostimulants in agriculture.

[This corrects the article DOI: 10.1371/journal.ppat.1008235.].

Although recombination is known to occur in foot-and-mouth disease virus (FMDV), it is considered only a minor determinant of virus sequence diversity. Analysis at phylogenetic scales shows inter-serotypic recombination events are rare, whereby recombination occurs almost exclusively in non-structural proteins. In this study we have estimated recombination rates within a natural host in an experimental setting. African buffaloes were inoculated with a SAT-1 FMDV strain containing two major viral sub-populations differing in their capsid sequence. This population structure enabled the detection of extensive within-host recombination in the genomic region coding for structural proteins and allowed recombination rates between the two sub-populations to be estimated. Quite surprisingly, the effective recombination rate in VP1 during the acute infection phase turns out to be about 0.1 per base per year, i.e. comparable to the mutation/substitution rate. Using a high-resolution map of effective within-host recombination in the capsid-coding region, we identified a linkage disequilibrium pattern in VP1 that is consistent with a mosaic structure with two main genetic blocks. Positive epistatic interactions between co-evolved variants appear to be present both within and between blocks. These interactions are due to intra-host selection both at the RNA and protein level. Overall our findings show that during FMDV co-infections by closely related strains, capsid-coding genes recombine within the host at a much higher rate than expected, despite the presence of strong constraints dictated by the capsid structure. Although these intra-host results are not immediately translatable to a phylogenetic setting, recombination and epistasis must play a major and so far underappreciated role in the molecular evolution of the virus at all scales.

Tobamovirus species represent a threat to solanaceous crops worldwide, due to their extreme stability and because they are seed borne. In particular, recent outbreaks of tomato brown rugose fruit virus in tomato and pepper crops led to the establishment of prompt control measures, and the need for reliable diagnosis was urged. Another member of the genus, tomato mottle mosaic virus, has recently gained attention due to reports in different continents and its common features with tomato brown rugose fruit virus. In this study, a new real-time RT-PCR detection system was developed for tomato brown rugose fruit virus and tomato mottle mosaic virus on tomato leaves and seeds using TaqMan chemistry. This test was designed to detect tomato mottle mosaic virus by amplifying the movement protein gene in a duplex assay with the tomato brown rugose fruit virus target on the CP-3’NTR region, which was previously validated as a single assay. The performance of this test was evaluated, displaying analytical sensitivity 10−5–10−6-fold dilution for seeds and leaves, respectively, and good analytical specificity, repeatability, and reproducibility. Using the newly developed and validated test, tomato brown rugose fruit virus detection was 100% concordant with previously performed analyses on 106 official samples collected in 2021 from different continents.

In 2020, a test performance study (TPS) for the specific detection of tomato brown rugose fruit virus (ToBRFV) was organized in the frame of the H2020 Valitest project. Since no validated tests were available, all the protocols reported in the literature were at first screened, performing preliminary studies in accordance with the EPPO standard PM 7/98 (4). Five molecular tests, two conventional RT-PCR and three real-time RT-PCR were found to be suitable and were included in the TPS. Thirty-four laboratories from 18 countries worldwide took part in TPS, receiving a panel of 22 blind samples. The panel consisted of sap belonging to symptomatic or asymptomatic leaves of Solanum lycopersicum and Capsicum annuum. The results returned by each laboratory were analyzed and diagnostic parameters were assessed for each test: reproducibility, repeatability, analytical sensitivity, diagnostic sensitivity and diagnostic specificity. All the evaluated tests resulted in being reliable in detecting ToBRFV and were included in an EPPO Standard PM 7/146—Diagnostics.

Background: Due to its central and strategic position in Europe and in the Mediterranean Basin, the Italian Peninsula played a pivotal role in the first peopling of the European continent and has been a crossroad of peoples and cultures since then. Aim: This study aims to gain more information on the genetic structure of modern Italian populations and to shed light on the migration/expansion events that led to their formation. Subjects and methods: High resolution Y-chromosome variation analysis in 817 unrelated males from 10 informative areas of Italy was performed. Haplogroup frequencies and microsatellite haplotypes were used, together with available data from the literature, to evaluate Mediterranean and European inputs and date their arrivals. Results: Fifty-three distinct Y-chromosome lineages were identified. Their distribution is in general agreement with geography, southern populations being more differentiated than northern ones. Conclusions: A complex genetic structure reflecting the multifaceted peopling pattern of the Peninsula emerged: southern populations show high similarity with those from the Middle East and Southern Balkans, while those from Northern Italy are close to populations of North-Western Europe and the Northern Balkans. Interestingly, the population of Volterra, an ancient town of Etruscan origin in Tuscany, displays a unique Y-chromosomal genetic structure.

African swine fever virus (ASFV) is the etiological agent of the devastating disease African swine fever (ASF), for which there is currently no licensed vaccine or treatment available. ASF is defined as one of the most serious animal diseases identified to date, due to its global spread in regions of Africa, Europe and Asia, causing massive economic losses. On the Italian island of Sardinia, the disease has been endemic since 1978, although the last control measures put in place achieved a significant reduction in ASF, and the virus has been absent from circulation since April 2019. Like many large DNA viruses, ASFV mutates at a relatively slow rate. However, the limited availability of whole-genome sequences from spatial-localized outbreaks makes it difficult to explore the small-scale genetic structure of these ASFV outbreaks. It is also unclear if the genetic variability within outbreaks can be captured in a handful of sequences, or if larger sequencing efforts can improve phylogenetic reconstruction and evolutionary or epidemiological inference. The aim of this study was to investigate the phylogenetic patterns of ASFV outbreaks between 1978 and 2018 in Sardinia, in order to characterize the epidemiological dynamics of the viral strains circulating in this Mediterranean island. To reach this goal, 58 new whole genomes of ASFV isolates were obtained, which represents the largest ASFV whole-genome sequencing effort to date. We provided a complete description of the genomic diversity of ASFV in terms of nucleotide mutations and small and large indels among the isolates collected during the outbreaks. The new sequences capture more than twice the genomic and phylogenetic diversity of all the previously published Sardinian sequences. The extra genomic diversity increases the resolution of the phylogenetic reconstruction, enabling us to dissect, for the first time, the genetic substructure of the outbreak. We found multiple ASFV subclusters within the phylogeny of the Sardinian epidemic, some of which coexisted in space and time.

African swine fever virus (ASFV) is one of the pathogens of highest concern worldwide. Despite different virus lineages co-circulating in several areas, dual infections in the same animal have been rarely observed, suggesting that ASF superinfections are infrequent events. Here we present the first genome-wide detection and analysis of two intragenotype dual ASFV infections. The dual infections have been detected in a hunted wild boar and in a pig carcass, both infected by ASFV genotype I in Sardinia in 1984 and 2018, respectively. We characterize the genetic differences between the two sequences, their intra-host frequency, and their phylogenetic relationship among fully sequenced ASFV strains from Sardinia. Both dual infections involve pairs of closely related but different viruses that were circulating in Sardinia in the same period. The results imply that dual ASFV infections or similar ASFV strains are more common than expected, especially in ASF endemic areas, albeit difficult to detect.

High-throughput sequencing (HTS) technologies and bioinformatic analyses are of growing interest to be used as a routine diagnostic tool in the field of plant viruses. The reliability of HTS workflows from sample preparation to data analysis and results interpretation for plant virus detection and identification must be evaluated (verified and validated) to approve this tool for diagnostics. Many different extraction methods, library preparation protocols, and sequence and bioinformatic pipelines are available for virus sequence detection. To assess the performance of plant virology diagnostic laboratories in using the HTS of ribosomal RNA depleted total RNA (ribodepleted totRNA) as a diagnostic tool, we carried out an interlaboratory comparison study in which eight participants were required to use the same samples, (RNA) extraction kit, ribosomal RNA depletion kit, and commercial sequencing provider, but also their own bioinformatics pipeline, for analysis. The accuracy of virus detection ranged from 65% to 100%. The false-positive detection rate was very low and was related to the misinterpretation of results as well as to possible cross-contaminations in the lab or sequencing provider. The bioinformatic pipeline used by each laboratory influenced the correct detection of the viruses of this study. The main difficulty was the detection of a novel virus as its sequence was not available in a publicly accessible database at the time. The raw data were reanalysed using Virtool to assess its ability for virus detection. All virus sequences were detected using Virtool in the different pools. This study revealed that the ribodepletion target enrichment for sample preparation is a reliable approach for the detection of plant viruses with different genomes. A significant level of virology expertise is needed to correctly interpret the results. It is also important to improve and complete the reference data.

The expression of accessory non-structural proteins V and W in Newcastle disease virus (NDV) infections depends on RNA editing. These proteins are derived from frameshifts of the sequence coding for the P protein via co-transcriptional insertion of one or two guanines in the mRNA. However, a larger number of guanines can be inserted with lower frequencies. We analysed data from deep RNA sequencing of samples from in vitro and in vivo NDV infections to uncover the patterns of mRNA editing in NDV. The distribution of insertions is well described by a simple Markov model of polymerase stuttering, providing strong quantitative confirmation of the molecular process hypothesised by Kolakofsky and collaborators three decades ago. Our results suggest that the probability that the NDV polymerase would stutter is about 0.45 initially, and 0.3 for further subsequent insertions. The latter probability is approximately independent of the number of previous insertions, the host cell, and viral strain. However, in LaSota infections, we also observe deviations from the predicted V/W ratio of about 3:1 according to this model, which could be attributed to deviations from this stuttering model or to further mechanisms downregulating the abundance of W protein.

Bilateria are the predominant clade of animals on Earth. Despite having evolved a wide variety of body plans and developmental modes, they are characterized by common morphological traits. By default, researchers have tried to link clade-specific genes to these traits, thus distinguishing bilaterians from non-bilaterians, by their gene content. Here we argue that it is rather biological processes that unite Bilateria and set them apart from their non-bilaterian sisters, with a less complex body morphology. To test this hypothesis, we compared proteomes of bilaterian and non-bilaterian species in an elaborate computational pipeline, aiming to search for a set of bilaterian-specific genes. Despite the limited confidence in their bilaterian specificity, we nevertheless detected Bilateria-specific functional and developmental patterns in the sub-set of genes conserved in distantly related Bilateria. Using a novel multi-species GO-enrichment method, we determined the functional repertoire of genes that are widely conserved among Bilateria. Analyzing expression profiles in three very distantly related model species—D. melanogaster, D. rerio and C. elegans—we find characteristic peaks at comparable stages of development and a delayed onset of expression in embryos. In particular, the expression of the conserved genes appears to peak at the phylotypic stage of different bilaterian phyla. In summary, our study illustrate how development connects distantly related Bilateria after millions of years of divergence, pointing to processes potentially separating them from non-bilaterians. We argue that evolutionary biologists should return from a purely gene-centric view of evolution and place more focus on analyzing and defining conserved developmental processes and periods.

The Site Frequency Spectrum (SFS) and the heterozygosity of allelic variants are among the most important summary statistics for population genetic analysis of diploid organisms. We discuss the generalization of these statistics to populations of autopolyploid organisms in terms of the joint Site Frequency/Dosage Spectrum and its expected value for autopolyploid populations that follow the standard neutral model. Based on these results, we present estimators of nucleotide variability from High-Throughput Sequencing (HTS) data of autopolyploids and discuss potential issues related to sequencing errors and variant calling. We use these estimators to generalize Tajima's D and other SFS-based neutrality tests to HTS data from autopolyploid organisms. Finally, we discuss how these approaches fail when the number of individuals is small. In fact, in autopolyploids there are many possible deviations from the Hardy–Weinberg equilibrium, each reflected in a different shape of the individual dosage distribution. The SFS from small samples is often dominated by the shape of these deviations of the dosage distribution from its Hardy–Weinberg expectations.

Over the past 15 years, 300 out of 6000 breeds of all farm animal species identified by the Food and Agriculture Organization of the United Nations (FAO) have gone extinct. Among cattle, many Podolian breeds are seriously endangered in various European areas. Podolian cattle include a group of very ancient European breeds, phenotypically close to the aurochs ancestors (Bos primigenius). The aim of the present study was to assess the genetic diversity of Podolian breeds and to reconstruct their origin.The mitochondrial DNA (mtDNA) control-regions of 18 Podolian breeds have been phylogenetically assessed. Nine non-Podolian breeds have been also included for comparison.The overall analysis clearly highlights some peculiarities in the mtDNA gene pool of some Podolian breeds. In particular, a principal component analysis point to a genetic proximity between five breeds (Chianina, Marchigiana, Maremmana, Podolica Italiana and Romagnola) reared in Central Italy and the Turkish Grey. We here propose the suggestive hypothesis of a dual ancestral contribution to the present gene pool of Podolian breeds, one deriving from Eastern European cattle; the other arising from the arrival of Middle Eastern cattle into Central Italy through a different route, perhaps by sea, ferried by Etruscan boats. The historical migration of Podolian cattle from North Eastern Europe towards Italy has not cancelled the mtDNA footprints of this previous ancient migration.

Several studies have investigated the differentiation of mitochondrial DNA in Eurasian, African and American cattle as well as archaeological bovine material. A global survey of these studies shows that haplogroup distributions are more stable in time than in space. All major migrations of cattle have shifted the haplogroup distributions considerably with a reduction of the number of haplogroups and/or an expansion of haplotypes that are rare or absent in the ancestral populations. The most extreme case is the almost exclusive colonization of Africa by the T1 haplogroup, which is rare in Southwest Asian cattle. In contrast, ancient samples invariably show continuity with present-day cattle from the same location. These findings indicate strong maternal founder effects followed by limited maternal gene flow when new territories are colonized. However, effects of adaptation to new environments may also play a role.

Recombination is one of the determinants of genetic diversity in the foot-and-mouth disease virus (FMDV). FMDV sequences have a mosaic structure caused by extensive intra- and inter-serotype recombination, with the exception of the capsid-encoding region. While these genome-wide patterns of broad-scale recombination are well studied, not much is known about the patterns of recombination that may exist within infected hosts. In addition, detection of recombination among viruses evolving at the within-host level is challenging due to the similarity of the sequences and the limitations in differentiating recombination from point mutations. Here, we present the first analysis of recombination events between closely related FMDV sequences occurring within buffalo hosts. The detection of these events was made possible by the occurrence of co-infection of two viral swarms with about 1% nucleotide divergence. We found more than 15 recombination events, unequally distributed across eight samples from different animals. The distribution of these events along the FMDV genome was neither uniform nor related to the phylogenetic distribution of recombination breakpoints, suggesting a mismatch between within-host evolutionary pressures and long-term selection for infectivity and transmissibility.

Doppel, a prion-like protein, is a GPI-membrane anchored protein generally not expressed in the Central Nervous System (CNS) of different mammalian species, including human. Nevertheless, in astrocytomas, a particular kind of glial tumors, the doppel encoding gene (PRND) is over-expressed and the corresponding protein product (Dpl) is ectopically localized in the cytoplasm of the tumor cells. In this study we have analysed the sub-cellular localization of Dpl using double-immunofluorescence staining and confocal microscopy examinations in two astrocytoma-derived human cell lines (IPDDC-A2 and D384-MG). Our results confirmed that Dpl is localized in the cytoplasm of the astrocytoma cells and indicated that it is mostly associated with Lamp-1 and Limp-2 positive lysosomal vesicles and, marginally, to the Golgi apparatus and other cellular organelles. Noticeably, none of the examined tumor cells showed a membrane-Dpl localization. The membrane-associated Dpl expression was restored after the transfection of the astrocytoma cells with mutated Dpl-expression vectors in its glycosylation sites. Additionally, Dpl showed altered expression and traffic using the acidotropic agent ammonium chloride, leading to the accumulation of Dpl in nascent exocytic vesicles. Altogether, these results indicated that in the astrocytic tumor cells Dpl has an altered biosynthetic trafficking, likely derived from abnormal post-translational processes: these modifications do not permit the localization of Dpl in correspondence of the plasma membrane and lead to its intracellular accumulation in the lysosomes. In these proteolytic compartments, the astrocytic tumor cells might provide to the degradation of the excess of a potentially cytotoxic Dpl product.

Genetic studies support the scenario that Bos taurus domestication occurred in the Near East during the Neolithic transition about 10 thousand years (ky) ago, with the likely exception of a minor secondary event in Italy. However, despite the proven effectiveness of whole mitochondrial genome data in providing valuable information concerning the origin of taurine cattle, until now no population surveys have been carried out at the level of mitogenomes in local breeds from the Near East or surrounding areas. Egypt is in close geographic and cultural proximity to the Near East, in particular the Nile Delta region, and was one of the first neighboring areas to adopt the Neolithic package. Thus, a survey of mitogenome variation of autochthonous taurine breeds from the Nile Delta region might provide new insights on the early spread of cattle rearing outside the Near East.Using Illumina high-throughput sequencing we characterized the mitogenomes from two cattle breeds, Menofi (N = 17) and Domiaty (N = 14), from the Nile Delta region. Phylogenetic and Bayesian analyses were subsequently performed.Phylogenetic analyses of the 31 mitogenomes confirmed the prevalence of haplogroup T1, similar to most African cattle breeds, but showed also high frequencies for haplogroups T2, T3 and Q1, and an extremely high haplotype diversity, while Bayesian skyline plots pointed to a main episode of population growth ~12.5 ky ago. Comparisons of Nile Delta mitogenomes with those from other geographic areas revealed that (i) most Egyptian mtDNAs are probably direct local derivatives from the founder domestic herds which first arrived from the Near East and the extent of gene flow from and towards the Nile Delta region was limited after the initial founding event(s); (ii) haplogroup Q1 was among these founders, thus proving that it underwent domestication in the Near East together with the founders of the T clades.

Geologically, Panama belongs to the Central American land-bridge between North and South America crossed by Homo sapiens >14 ka ago. Archaeologically, it belongs to a wider Isthmo-Colombian Area. Today, seven indigenous ethnic groups account for 12.3% of Panama's population. Five speak Chibchan languages and are characterized by low genetic diversity and a high level of differentiation. In addition, no evidence of differential structuring between maternally and paternally inherited genes has been reported in isthmian Chibchan cultural groups. Recent data have shown that 83% of the Panamanian general population harbour mitochondrial DNAs (mtDNAs) of Native American ancestry. Considering differential male/female mortality at European contact and multiple degrees of geographical and genetic isolation over the subsequent five centuries, the Y-chromosome Native American component is expected to vary across different geographic regions and communities in Panama. To address this issue, we investigated Y-chromosome variation in 408 modern males from the nine provinces of Panama and one indigenous territory (the comarca of Kuna Yala). In contrast to mtDNA data, the Y-chromosome Native American component (haplogroup Q) exceeds 50% only in three populations facing the Caribbean Sea: the comarca of Kuna Yala and Bocas del Toro province where Chibchan languages are spoken by the majority, and the province of Colón where many Kuna and people of mixed indigenous-African-and-European descent live. Elsewhere the Old World component is dominant and mostly represented by western Eurasian haplogroups, which signal the strong male genetic impact of invaders. Sub-Saharan African input accounts for 5.9% of male haplotypes. This reflects the consequences of the colonial Atlantic slave trade and more recent influxes of West Indians of African heritage. Overall, our findings reveal a local evolution of the male Native American ancestral gene pool, and a strong but geographically differentiated unidirectional sex bias in the formation of local modern Panamanian populations.

Phytoplasmas are cell wall-less microorganisms associated with plant diseases worldwide. Many important food, vegetable and fruits crops as well as ornamental plants can be severely affected by these pathogens, with significant economic impacts. Phytoplasma diseases of ornamentals have been described worldwide in a wide range of plant genera, and 11 different 16Sr groups have been identified. In Italy, many ornamental plant species belonging to several botanical families have been found to be infected by phytoplasmas, classified into the ribosomal groups 16SrI, 16SrII, 16SrV and 16SrXII. During a survey carried out in commercial gardens in Rome, some marguerite daisy (Argyranthemum frutescens) plants showing symptoms of phytoplasma-like disease, were collected and submitted to molecular analyses. Cloning and sequencing of the portion of the 16S rRNA gene followed by BLAST analysis, real and virtual restriction fragment length polymorphism anlaysis with AluI and RsaI, allowed assignment of the detected phytoplasma to the 16SrIX-C group (Picris echioides yellows, PEY).

Citrus tristeza virus (CTV) outbreaks have been reported in the main citrus-growing regions of Italy in the past 10 years. In some areas where eradication efforts failed to suppress spread, high CTV incidence is now observed. Recently, potentially severe CTV strains were detected in Calabria (southern Italy), one of the major citrus-growing area. As a result, investigations of the virulence and molecular features of CTV populations spreading in this region were undertaken. Virus was detected by enzyme-linked immunosorbent assay (ELISA) using a broad spectrum polyclonal antiserum, and was differentiated into potential virulent categories with the severe-strain discriminating monoclonal antibody MCA13. Isolate genotyping was conducted using reverse-transcription polymerase chain reaction (RT-PCR) with multiple molecular markers (MMM), single-strand conformation polymorphism (SSCP) analysis of the amplicons from the genes coding for the coat protein (CP) p25 and the non-structural p20 protein as well as sequence analysis. Based on the serological reactivity, the isolates were differentiated in two distinct serogroups: MCA13-reactive and MCA13 non-reactive. Similarly, based on the molecular profile, the isolates were grouped in two genetically distinct phylogenetic clusters, and associated either with a T30-like or with a T3-like genotype. These data were related to the results of biological indexing on standard indicator plants, which distinguished isolates causing mild or severe seedling yellow reactions. The study has demonstrated the presence of MCA13-reactive isolates associated with a T3-like genotype and causing severe seedling yellows in sour orange, grapefruit and lemon seedlings, and stem pitting in Mexican lime.

The coalescent with recombination is a fundamental model to describe the genealogical history of DNA sequence samples from recombining organisms. Considering recombination as a process which acts along genomes and which creates sequence segments with shared ancestry, we study the influence of single recombination events upon tree characteristics of the coalescent. We focus on properties such as tree height and tree balance and quantify analytically the changes in these quantities incurred by recombination in terms of probability distributions. We find that changes in tree topology are often relatively mild under conditions of neutral evolution, while changes in tree height are on average quite large. Our results add to a quantitative understanding of the spatial coalescent and provide the neutral reference to which the impact by other evolutionary scenarios, for instance tree distortion by selective sweeps, can be compared.

BACKGROUND: Most genetic studies on modern cattle have established a common origin for all taurine breeds in the Near East, during the Neolithic transition about 10 thousand years (ka) ago. Yet, the possibility of independent and/or secondary domestication events is still debated and is fostered by the finding of rare mitochondrial DNA (mtDNA) haplogroups like P, Q and R. Haplogroup T1, because of its geographic distribution, has been the subject of several investigations pointing to a possible independent domestication event in Africa and suggesting a genetic contribution of African cattle to the formation of Iberian and Creole cattle. Whole mitochondrial genome sequence analysis, with its proven effectiveness in improving the resolution of phylogeographic studies, is the most appropriate tool to investigate the origin and structure of haplogroup T1. METHODOLOGY: A survey of >2200 bovine mtDNA control regions representing 28 breeds (15 European, 10 African, 3 American) identified 281 subjects belonging to haplogroup T1. Fifty-four were selected for whole mtDNA genome sequencing, and combined with ten T1 complete sequences from previous studies into the most detailed T1 phylogenetic tree available to date. CONCLUSIONS: Phylogenetic analysis of the 64 T1 mitochondrial complete genomes revealed six distinct sub-haplogroups (T1a-T1f). Our data support the overall scenario of a Near Eastern origin of the T1 sub-haplogroups from as much as eight founding T1 haplotypes. However, the possibility that one sub-haplogroup (T1d) arose in North Africa, in domesticated stocks, shortly after their arrival from the Near East, can not be ruled out. Finally, the previously identified "African-derived American" (AA) haplotype turned out to be a sub-clade of T1c (T1c1a1). This haplotype was found here for the first time in Africa (Egypt), indicating that it probably originated in North Africa, reached the Iberian Peninsula and sailed to America, with the first European settlers.