Silvia Bonfiglio

The Multifaceted Origin of Taurine Cattle Reflected by the Mitochondrial Genome

A Neolithic domestication of taurine cattle in the Fertile Crescent from local aurochsen (Bos primigenius) is generally accepted, but a genetic contribution from European aurochsen has been proposed. Here we performed a survey of a large number of taurine cattle mitochondrial DNA (mtDNA) control regions from numerous European breeds confirming the overall clustering within haplogroups (T1, T2 and T3) of Near Eastern ancestry, but also identifying eight mtDNAs (1.3%) that did not fit in haplogroup T. Sequencing of the entire mitochondrial genome showed that four mtDNAs formed a novel branch (haplogroup R) which, after the deep bifurcation that gave rise to the taurine and zebuine lineages, constitutes the earliest known split in the mtDNA phylogeny of B. primigenius. The remaining four mtDNAs were members of the recently discovered haplogroup Q. Phylogeographic data indicate that R mtDNAs were derived from female European aurochsen, possibly in the Italian Peninsula, and sporadically included in domestic herds. In contrast, the available data suggest that Q mtDNAs and T subclades were involved in the same Neolithic event of domestication in the Near East. Thus, the existence of novel (and rare) taurine haplogroups highlights a multifaceted genetic legacy from distinct B. primigenius populations. Taking into account that the maternally transmitted mtDNA tends to underestimate the extent of gene flow from European aurochsen, the detection of the R mtDNAs in autochthonous breeds, some of which are endangered, identifies an unexpected reservoir of genetic variation that should be carefully preserved.

The Enigmatic Origin of Bovine mtDNA Haplogroup R: Sporadic Interbreeding or an Independent Event of Bos primigenius Domestication in Italy?

Background When domestic taurine cattle diffused from the Fertile Crescent, local wild aurochsen (Bos primigenius) were still numerous. Moreover, aurochsen and introduced cattle often coexisted for millennia, thus providing potential conditions not only for spontaneous interbreeding, but also for pastoralists to create secondary domestication centers involving local aurochs populations. Recent mitochondrial genomes analyses revealed that not all modern taurine mtDNAs belong to the shallow macro-haplogroup T of Near Eastern origin, as demonstrated by the detection of three branches (P, Q and R) radiating prior to the T node in the bovine phylogeny. These uncommon haplogroups represent excellent tools to evaluate if sporadic interbreeding or even additional events of cattle domestication occurred. Methodology The survey of the mitochondrial DNA (mtDNA) control-region variation of 1,747 bovine samples (1,128 new and 619 from previous studies) belonging to 37 European breeds allowed the identification of 16 novel non-T mtDNAs, which after complete genome sequencing were confirmed as members of haplogroups Q and R. These mtDNAs were then integrated in a phylogenetic tree encompassing all available P, Q and R complete mtDNA sequences. Conclusions Phylogenetic analyses of 28 mitochondrial genomes belonging to haplogroups P (N = 2), Q (N = 16) and R (N = 10) together with an extensive survey of all previously published mtDNA datasets revealed major similarities between haplogroups Q and T. Therefore, Q most likely represents an additional minor lineage domesticated in the Near East together with the founders of the T subhaplogroups. Whereas, haplogroup R is found, at least for the moment, only in Italy and nowhere else, either in modern or ancient samples, thus supporting an origin from European aurochsen. Haplogroup R could have been acquired through sporadic interbreeding of wild and domestic animals, but our data do not rule out the possibility of a local and secondary event of B. primigenius domestication in Italy.

Exome sequencing and pathway analysis for identification of genetic variability relevant for bronchopulmonary dysplasia (BPD) in preterm newborns: A pilot study

BACKGROUND Bronchopulmonary dysplasia (BPD) is the most common chronic lung disease in infancy, affecting preterm children with low birth weight. The disease has a multifactorial aetiology with a significant genetic component; until now published association studies have identified several candidate genes but only few of these data has been replicated. In this pilot study, we approached exome sequencing aimed at identifying non-common variants, which are expected to have a stronger phenotypic effect. MATERIALS AND METHODS We performed this study on 26 Italian severely affected BPD preterm unrelated newborns, homogeneously selected from a large prospective cohort. We used an Illumina HiSeq 2000 for sequencing. Data analysis was focussed on genes previously associated to BPD susceptibility and to new candidates in related pathways, highlighted by a prioritization analysis performed using ToppGene Suite. RESULTS By exome sequencing, we identified 3369 novel variants, with a median of 400 variations per sample. The top candidate genes highlighted were NOS2, MMP1, CRP, LBP and the toll-like receptor (TLR) family. All of them have been confirmed with Sanger sequencing. CONCLUSIONS Potential candidate genes have been discovered in this preliminary study; the pathogenic role of identified variants will need to be confirmed with functional and segregation studies and possibly with further methods, able to evaluate the collective influence of rare variants. Moreover, additional candidates will be tested and genetic analysis will be extended to all affected children.

Mitogenomes from Egyptian Cattle Breeds: New Clues on the Origin of Haplogroup Q and the Early Spread of Bos taurus from the Near East

Background 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. Methodology 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. Conclusions 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.

A novel USP9Y polymorphism allowing a rapid and unambiguous classification of Bos taurus Y chromosomes into haplogroups

A new sequence-tagged site (STS) was identified within intron 26 of the bovine USP9Y gene, showing an 81-base pair insertion (g.76439_76440ins81 in sequence with GenBank accession FJ195366) able to distinguish Y2 and Y3 Bos Y haplogroups from Y1. Moreover, four Y3-specific sequence variants allow a distinction from haplogroup Y2. The typing of a Bison bison Y chromosome indicates that the ancestral allele for the USP9Y 81-bp insertion is the short Y1 version. The results from typing the new STS in 1230 cattle Y chromosomes are fully consistent with their classification through standard methods. Thanks to the newly identified STS, it is now possible to assign cattle Y chromosomes to the currently known haplogroups using a single marker.

Characterisation and Validation of Insertions and Deletions in 173 Patient Exomes

Recent advances in genomics technologies have spurred unprecedented efforts in genome and exome re-sequencing aiming to unravel the genetic component of rare and complex disorders. While in rare disorders this allowed the identification of novel causal genes, the missing heritability paradox in complex diseases remains so far elusive. Despite rapid advances of next-generation sequencing, both the technology and the analysis of the data it produces are in its infancy. At present there is abundant knowledge pertaining to the role of rare single nucleotide variants (SNVs) in rare disorders and of common SNVs in common disorders. Although the 1,000 genome project has clearly highlighted the prevalence of rare variants and more complex variants (e.g. insertions, deletions), their role in disease is as yet far from elucidated. We set out to analyse the properties of sequence variants identified in a comprehensive collection of exome re-sequencing studies performed on samples from patients affected by a broad range of complex and rare diseases (N = 173). Given the known potential for Loss of Function (LoF) variants to be false positive, we performed an extensive validation of the common, rare and private LoF variants identified, which indicated that most of the private and rare variants identified were indeed true, while common novel variants had a significantly higher false positive rate. Our results indicated a strong enrichment of very low-frequency insertion/deletion variants, so far under-investigated, which might be difficult to capture with low coverage and imputation approaches and for which most of study designs would be under-powered. These insertions and deletions might play a significant role in disease genetics, contributing specifically to the underlining rare and private variation predicted to be discovered through next generation sequencing.

A CD8α− Subset of CD4+SLAMF7+ Cytotoxic T Cells Is Expanded in Patients With IgG4-Related Disease and Decreases Following Glucocorticoid Treatment

An unconventional population of CD4+ signaling lymphocytic activation molecule family member 7–positive (SLAMF7+) cytotoxic effector memory T (TEM) cells (CD4+ cytotoxic T lymphocytes [CTLs]) has been linked causally to IgG4‐related disease (IgG4‐RD). Glucocorticoids represent the first‐line therapeutic approach in patients with IgG4‐RD, but their mechanism of action in this specific condition remains unknown. We undertook this study to determine the impact of glucocorticoids on CD4+ CTLs in IgG4‐RD.

REST-Governed Gene Expression Profiling in a Neuronal Cell Model Reveals Novel Direct and Indirect Processes of Repression and Up-Regulation.

The role of REST changes in neurons, including the rapid decrease of its level during differentiation and its fluctuations during many mature functions and diseases, is well established. However, identification of many thousand possible REST-target genes, mostly based on indirect criteria, and demonstration of their operative dependence on the repressor have been established for only a relatively small fraction. In the present study, starting from our recently published work, we have expanded the identification of REST-dependent genes, investigated in two clones of the PC12 line, a recognized neuronal cell model, spontaneously expressing different levels of REST: very low as in neurons and much higher as in most non-neural cells. The molecular, structural and functional differences of the two PC12 clones were shown to depend largely on their different REST level and the ensuing variable expression of some dependent genes. Comprehensive RNA-Seq analyses of the 13,700 genes expressed, validated by parallel RT-PCR and western analyses of mRNAs and encoded proteins, identified in the high-REST clone two groups of almost 900 repressed and up-regulated genes. Repression is often due to direct binding of REST to target genes; up-regulation to indirect mechanism(s) mostly mediated by REST repression of repressive transcription factors. Most, but not all, genes governing neurosecretion, excitability, and receptor channel signaling were repressed in the high REST clone. The genes governing expression of non-channel receptors (G protein-coupled and others), although variably affected, were often up-regulated together with the genes of intracellular kinases, small G proteins, cytoskeleton, cell adhesion, and extracellular matrix proteins. Expression of REST-dependent genes governing functions other than those mentioned so far were also identified. The results obtained by the parallel investigation of the two PC12 clones revealed the complexity of the REST molecular and functional role, deciphering new aspects of its participation in neuronal functions. The new findings could be relevant for further investigation and interpretation of physiological processes typical of neurons. Moreover, they could be employed as tools in the study of neuronal diseases recently shown to depend on REST for their development.

Epigenomics of Neural Cells: REST-Induced Down- and Upregulation of Gene Expression in a Two-Clone PC12 Cell Model

Cell epigenomics depends on the marks released by transcription factors operating via the assembly of complexes that induce focal changes of DNA and histone structure. Among these factors is REST, a repressor that, via its strong decrease, governs both neuronal and neural cell differentiation and specificity. REST operation on thousands of possible genes can occur directly or via indirect mechanisms including repression of other factors. In previous studies of gene down- and upregulation, processes had been only partially investigated in neural cells. PC12 are well-known neural cells sharing properties with neurons. In the widely used PC12 populations, low-REST cells coexist with few, spontaneous high-REST PC12 cells. High- and low-REST PC12 clones were employed to investigate the role and the mechanisms of the repressor action. Among 15,500 expressed genes we identified 1,770 target and nontarget, REST-dependent genes. Functionally, these genes were found to operate in many pathways, from synaptic function to extracellular matrix. Mechanistically, downregulated genes were predominantly repressed directly by REST; upregulated genes were mostly governed indirectly. Among other factors, Polycomb complexes cooperated with REST for downregulation, and Smad3 and Myod1 participated in upregulation. In conclusion, we have highlighted that PC12 clones are a useful model to investigate REST, opening opportunities to development of epigenomic investigation.

Whole exome sequencing of independent lung adenocarcinoma, lung squamous cell carcinoma, and malignant peritoneal mesothelioma: A case report

AbstractThe presence of multiple primary tumors (MPT) in a single patient has been identified with an increasing frequency. A critical issue is to establish if the second tumor represents an independent primary cancer or a metastasis. Therefore, the assessment of MPT clonal origin might help understand the disease behavior and improve the management/prognosis of the patient.Herein, we report a 73-year-old male smoker who developed 2 primary lung cancers (adenocarcinoma and squamous cell carcinoma) and a malignant peritoneal mesothelioma (PM).Whole exome sequencing (WES) of the 3 tumors and of germline DNA was performed to determine the clonal origin and identify genetic cancer susceptibility.Both lung cancers were characterized by a high mutational rate with distinct mutational profiles and activation of tumor-specific pathways. Conversely, the PM harbored a relative low number of genetic variants and a novel mutation in the WT1 gene that might be involved in the carcinogenesis of nonasbestos-related mesothelioma. Finally, WES of the germinal DNA displayed several single nucleotide polymorphisms in DNA repair genes likely conferring higher cancer susceptibility.Overall, WES did not disclose any somatic genetic variant shared across the 3 tumors, suggesting their clonal independency; however, the carcinogenic effect of smoke combined with a deficiency in DNA repair genes and the patient advanced age might have been responsible for the MPT development. This case highlights the WES importance to define the clonal origin of MPT and susceptibility to cancer.