1
|
Zink FA, Tembrock LR, Timm AE, Gilligan TM. Ultra-deep sequencing of 45S rDNA to discern intragenomic diversity in three Chrysodeixis species for molecular identification. Sci Rep 2023; 13:13017. [PMID: 37563256 PMCID: PMC10415407 DOI: 10.1038/s41598-023-39673-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2023] [Accepted: 07/28/2023] [Indexed: 08/12/2023] Open
Abstract
Species identification is necessary to prevent introductions of exotic plant pests through global trade. Many of these pests are understudied and lack publicly available DNA sequence data on which rapid molecular identification methods can be based. One such lineage is the genus Chrysodeixis, which includes three species of potential concern for United States trade initiatives: C. includens, C. chalcites, and C. eriosoma. Here we describe a method to generate robust 45S rDNA profiles using long read sequencing in order to clarify evolutionary relationships and develop a real-time PCR identification technique. Such an identification tool will be useful in rapidly differentiating between Chrysodeixis species of quarantine concern where traditional morphological identification methods are insufficient. Molecular methods such as this greatly reduce the time spent identifying each specimen, allow for detection of eDNA, vastly increase throughput, and increase the probability of detection. The methods presented here will be generally adaptable to any understudied lepidopteran taxa that necessitates a molecular diagnostic assay and, with adjustment or testing of the primers, could be applied to any group for which development of rDNA profiles in a benchtop setting would prove useful.
Collapse
Affiliation(s)
- Frida A Zink
- Department of Agricultural Biology, Colorado State University, Fort Collins, CO, USA
| | - Luke R Tembrock
- Department of Agricultural Biology, Colorado State University, Fort Collins, CO, USA.
| | - Alicia E Timm
- Department of Agricultural Biology, Colorado State University, Fort Collins, CO, USA
| | - Todd M Gilligan
- Pest Identification Technology Laboratory, USDA-APHIS-PPQ-Science and Technology, Fort Collins, CO, USA
| |
Collapse
|
2
|
Lack of ITS sequence homogenization in Erysimum species (Brassicaceae) with different ploidy levels. Sci Rep 2022; 12:16907. [PMID: 36207443 PMCID: PMC9546898 DOI: 10.1038/s41598-022-20194-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2022] [Accepted: 09/09/2022] [Indexed: 11/29/2022] Open
Abstract
The internal transcribed spacers (ITS) exhibit concerted evolution by the fast homogenization of these sequences at the intragenomic level. However, the rate and extension of this process are unclear and might be conditioned by the number and divergence of the different ITS copies. In some cases, such as hybrid species and polyploids, ITS sequence homogenization appears incomplete, resulting in multiple haplotypes within the same organism. Here, we studied the dynamics of concerted evolution in 85 individuals of seven plant species of the genus Erysimum (Brassicaceae) with multiple ploidy levels. We estimated the rate of concerted evolution and the degree of sequence homogenization separately for ITS1 and ITS2 and whether these varied with ploidy. Our results showed incomplete sequence homogenization, especially for polyploid samples, indicating a lack of concerted evolution in these taxa. Homogenization was usually higher in ITS2 than in ITS1, suggesting that concerted evolution operates more efficiently on the former. Furthermore, the hybrid origin of several species appears to contribute to the maintenance of high haplotype diversity, regardless of the level of ploidy. These findings indicate that sequence homogenization of ITS is a dynamic and complex process that might result in varying intra- and inter-genomic diversity levels.
Collapse
|
3
|
Comparative Analyses of 35 Complete Chloroplast Genomes from the Genus Dalbergia (Fabaceae) and the Identification of DNA Barcodes for Tracking Illegal Logging and Counterfeit Rosewood. FORESTS 2022. [DOI: 10.3390/f13040626] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The genus Dalbergia contains more than 200 species, several of which are trees that produce traditional medicines and extremely high-value timber commonly referred to as rosewood. Due to the rarity of these species in the wild, the high value of the timber, and a growing international illicit trade, CITES (Convention on International Trade in Endangered Species of Wild Fauna and Flora) has listed the entire genus in appendix II and the species Dalbergia nigra in appendix I because species in this genus are considered at risk of extinction. Given this, and the fact that species or even genus level determination is nearly impossible from cut timber morphology, alternative molecular methods are needed to identify and track intercepted rosewood. To better identify rosewood using molecular methods, we sequenced and assembled eight chloroplast genomes including D. nigra as well as conducted comparative analyses with all other available chloroplast genomes in Dalbergia and closely related lineages. From these analyses, numerous repeats including simple sequence repeats (SSR) and conserved nucleotide polymorphisms unique to subclades within the genus were detected. From phylogenetic analysis based on the CDS from 77 chloroplast genes, the groups Siam rosewood and scented rosewood resolved as monophyletic, supporting the morphological traits used to delimit these species. In addition, several instances of paraphyly and polyphyly resulting from mismatches between taxonomic determinations and phylogenetic tree topology were identified. Ultimately, the highly variable regions in the chloroplast genomes will provide useful plastid markers for further studies regarding the identification, phylogeny, and population genetics of Dalbergia species, including those frequently intercepted in illegal trade.
Collapse
|
4
|
Milani D, Ruiz-Ruano FJ, Camacho JPM, Cabral-de-Mello DC. Out of patterns, the euchromatic B chromosome of the grasshopper Abracris flavolineata is not enriched in high-copy repeats. Heredity (Edinb) 2021; 127:475-483. [PMID: 34482369 PMCID: PMC8551250 DOI: 10.1038/s41437-021-00470-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2021] [Revised: 08/24/2021] [Accepted: 08/25/2021] [Indexed: 02/07/2023] Open
Abstract
In addition to the normal set of standard (A) chromosomes, some eukaryote species harbor supernumerary (B) chromosomes. In most cases, B chromosomes show differential condensation with respect to A chromosomes and display dark C-bands of heterochromatin, and some of them are highly enriched in repetitive DNA. Here we perform a comprehensive NGS (next-generation sequencing) analysis of the repeatome in the grasshopper Abracris flavolineata aimed at uncovering the molecular composition and origin of its B chromosome. Our results have revealed that this B chromosome shows a DNA repeat content highly similar to the DNA repeat content observed for euchromatic (non-C-banded) regions of A chromosomes. Moreover, this B chromosome shows little enrichment for high-copy repeats, with only a few elements showing overabundance in B-carrying individuals compared to the 0B individuals. Consequently, the few satellite DNAs (satDNAs) mapping on the B chromosome were mostly restricted to its centromeric and telomeric regions, and they displayed much smaller bands than those observed on the A chromosomes. Our data support the intraspecific origin of the B chromosome from the longest autosome by misdivision, isochromosome formation, and additional restructuring, with accumulation of specific repeats in one or both B chromosome arms, yielding a submetacentric B. Finally, the absence of B-specific satDNAs, which are frequent in other species, along with its euchromatic nature, suggest that this B chromosome arose recently and might still be starting a heterochromatinization process. On this basis, it could be a good model to investigate the initial steps of B chromosome evolution.
Collapse
Affiliation(s)
- Diogo Milani
- Departamento de Biologia Geral e Aplicada, UNESP - Univ Estadual Paulista, Instituto de Biociências/IB, Rio Claro, São Paulo, Brazil
| | - Francisco J Ruiz-Ruano
- Evolutionary Biology Centre, Department of Organismal Biology - Systematic Biology, Uppsala University, Uppsala, Sweden
- Norwich Research Park, School of Biological Sciences, University of East Anglia, Norwich, UK
| | - Juan Pedro M Camacho
- Departamento de Genética, Facultad de Ciencias, UGR - Univ de Granada, Granada, Spain
| | - Diogo C Cabral-de-Mello
- Departamento de Biologia Geral e Aplicada, UNESP - Univ Estadual Paulista, Instituto de Biociências/IB, Rio Claro, São Paulo, Brazil.
| |
Collapse
|
5
|
Assay Optimization Can Equalize the Sensitivity of Real-Time PCR with ddPCR for Detection of Helicoverpa armigera (Lepidoptera: Noctuidae) in Bulk Samples. INSECTS 2021; 12:insects12100885. [PMID: 34680654 PMCID: PMC8538000 DOI: 10.3390/insects12100885] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/12/2021] [Revised: 09/14/2021] [Accepted: 09/24/2021] [Indexed: 11/17/2022]
Abstract
Simple Summary Invasive species are a constant threat to agriculture throughout the world against which early detection is one of the primary defenses. The Old World bollworm is one of the most important invasive agricultural pests in the world. While historically absent from the Americas, this species was first found in South America in 2013 and poses an ongoing threat of spreading into North America. Surveys are conducted each year, which result in hundreds or thousands of traps that must be screened for this species. Unfortunately, the most common non-target is the native corn earworm, which is nearly identical morphologically to the Old World bollworm and cannot be easily separated. Molecular methods have been developed to screen these trap samples, but the required equipment is expensive and not commonly available. This study details improvements to current molecular methods that will allow for screening of bulk trap samples using standard laboratory instruments and protocols. The ability to perform these methods in nearly any molecular biology lab will greatly enhance our ability to detect and exclude this important pest. Abstract Helicoverpa armigera (Hübner) is one of the most important agricultural pests in the world. This historically Old World species was first reported in Brazil in 2013 and has since spread throughout much of South America and into the Caribbean. Throughout North America, H. armigera surveys are ongoing to detect any incursions. Each trap is capable of capturing hundreds of native Helicoverpa zea (Boddie). The two species cannot be separated without genitalic dissection or molecular methods. A ddPCR assay is currently used to screen large trap samples, but this equipment is relatively uncommon and expensive. Here, we optimized a newly designed assay for accurate and repeatable detection of H. armigera in bulk samples across both ddPCR and less costly, and more common, real-time PCR methods. Improvements over previously designed assays were sought through multiple means. Our results suggest bulk real-time PCR assays can be improved through changes in DNA extraction and purification, so that real-time PCR can be substituted for ddPCR in screening projects. While ddPCR remains a more sensitive method for detection of H. armigera in bulk samples, the improvements in assay design, DNA extraction, and purification presented here also enhance assay performance over previous protocols.
Collapse
|
6
|
Martí E, Milani D, Bardella VB, Albuquerque L, Song H, Palacios-Gimenez OM, Cabral-de-Mello DC. Cytogenomic analysis unveils mixed molecular evolution and recurrent chromosomal rearrangements shaping the multigene families on Schistocerca grasshopper genomes. Evolution 2021; 75:2027-2041. [PMID: 34155627 DOI: 10.1111/evo.14287] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2021] [Revised: 05/11/2021] [Accepted: 05/26/2021] [Indexed: 11/26/2022]
Abstract
Multigene families are essential components of eukaryotic genomes and play key roles either structurally and functionally. Their modes of evolution remain elusive even in the era of genomics, because multiple multigene family sequences coexist in genomes, particularly in large repetitive genomes. Here, we investigate how the multigene families 18S rDNA, U2 snDNA, and H3 histone evolved in 10 species of Schistocerca grasshoppers with very large and repeat-enriched genomes. Using sequenced genomes and fluorescence in situ hybridization mapping, we find substantial differences between species, including the number of chromosomal clusters, changes in sequence abundance and nucleotide composition, pseudogenization, and association with transposable elements (TEs). The intragenomic analysis of Schistocerca gregaria using long-read sequencing and genome assembly unveils conservation for H3 histone and recurrent pseudogenization for 18S rDNA and U2 snDNA, likely promoted by association with TEs and sequence truncation. Remarkably, TEs were frequently associated with truncated copies, were also among the most abundant in the genome, and revealed signatures of recent activity. Our findings suggest a combined effect of concerted and birth-and-death models driving the evolution of multigene families in Schistocerca over the last 8 million years, and the occurrence of intra- and interchromosomal rearrangements shaping their chromosomal distribution. Despite the conserved karyotype in Schistocerca, our analysis highlights the extensive reorganization of repetitive DNAs in Schistocerca, contributing to the advance of comparative genomics for this important grasshopper genus.
Collapse
Affiliation(s)
- Emiliano Martí
- Departamento de Biologia Geral e Aplicada, UNESP - Univ Estadual Paulista, Instituto de Biociências/IB, Rio Claro, 13506-900, Brazil
| | - Diogo Milani
- Departamento de Biologia Geral e Aplicada, UNESP - Univ Estadual Paulista, Instituto de Biociências/IB, Rio Claro, 13506-900, Brazil
| | - Vanessa B Bardella
- Departamento de Biologia Geral e Aplicada, UNESP - Univ Estadual Paulista, Instituto de Biociências/IB, Rio Claro, 13506-900, Brazil
| | - Lucas Albuquerque
- Departamento de Biologia Geral e Aplicada, UNESP - Univ Estadual Paulista, Instituto de Biociências/IB, Rio Claro, 13506-900, Brazil
| | - Hojun Song
- Department of Entomology, Texas A&M University, College Station, Texas, 77843
| | - Octavio M Palacios-Gimenez
- Department of Organismal Biology - Systematic Biology, Evolutionary Biology Centre, Uppsala University, Uppsala, SE-75236, Sweden.,Population Ecology Group, Institute of Ecology and Evolution, Friedrich Schiller University Jena, Jena, DE-07743, Germany
| | - Diogo C Cabral-de-Mello
- Departamento de Biologia Geral e Aplicada, UNESP - Univ Estadual Paulista, Instituto de Biociências/IB, Rio Claro, 13506-900, Brazil
| |
Collapse
|
7
|
Satellite DNA Is an Inseparable Fellow Traveler of B Chromosomes. PROGRESS IN MOLECULAR AND SUBCELLULAR BIOLOGY 2021; 60:85-102. [PMID: 34386873 DOI: 10.1007/978-3-030-74889-0_4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Next-Generation Sequencing (NGS) has revealed that B chromosomes in several species are enriched in repetitive DNA, mostly satellite DNA (satDNA). This raises the question of whether satDNA is important to B chromosomes for functional reasons or else its abundance on Bs is simply a consequence of properties of B chromosomes such as their dispensability and late replication. Here we review current knowledge in this respect and contextualize it within the frame of practical difficulties to perform this kind of research, the most important being the absence of good full genome sequencing for B-carrying species, which is an essential requisite to ascertain the intragenomic origin of B chromosomes. Our review analysis on 16 species revealed that 38% of them showed B-specific satDNAs whereas only one of them (6%) carried an inter-specifically originated B chromosome. This shows that B-specific satDNA families can eventually evolve in intraspecifically arisen B chromosomes. Finally, the possibility of satDNA accumulation on B chromosomes for functional reasons is exemplified by B chromosomes in rye, as they contain B-specific satDNAs which are transcribed and occupy chromosome locations where they might facilitate the kind of drive shown by this B chromosome during pollen grain mitosis.
Collapse
|
8
|
Grant JR, Pilotte N, Williams SA. A Case for Using Genomics and a Bioinformatics Pipeline to Develop Sensitive and Species-Specific PCR-Based Diagnostics for Soil-Transmitted Helminths. Front Genet 2019; 10:883. [PMID: 31608116 PMCID: PMC6768101 DOI: 10.3389/fgene.2019.00883] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2019] [Accepted: 08/21/2019] [Indexed: 12/22/2022] Open
Abstract
The balance of expense and ease of use vs. specificity and sensitivity in diagnostic assays for helminth disease is an important consideration, with expense and ease often winning out in endemic areas where funds and sophisticated equipment may be scarce. In this review, we argue that molecular diagnostics, specifically new assays that have been developed with the aid of next-generation sequence data and robust bioinformatic tools, more than make up for their expense with the benefit of a clear and precise assessment of the situation on the ground. Elimination efforts associated with the London Declaration and the World Health Organization (WHO) 2020 Roadmap have resulted in areas of low disease incidence and reduced infection burdens. An accurate assessment of infection levels is critical for determining where and when the programs can be successfully ended. Thus, more sensitive assays are needed in locations where elimination efforts are approaching a successful conclusion. Although microscopy or more general PCR targets have a role to play, they can mislead and cause study results to be confounded. Hyper-specific qPCR assays enable a more definitive assessment of the situation in the field, as well as of shifting dynamics and emerging diseases.
Collapse
Affiliation(s)
- Jessica R. Grant
- Department of Biological Sciences, Smith College, Northampton, MA, United States
| | - Nils Pilotte
- Department of Biological Sciences, Smith College, Northampton, MA, United States
- Molecular and Cellular Biology Program, University of Massachusetts, Amherst, MA, United States
| | - Steven A. Williams
- Department of Biological Sciences, Smith College, Northampton, MA, United States
- Molecular and Cellular Biology Program, University of Massachusetts, Amherst, MA, United States
| |
Collapse
|
9
|
Milani D, Bardella VB, Ferretti ABSM, Palacios-Gimenez OM, Melo ADS, Moura RC, Loreto V, Song H, Cabral-de-Mello DC. Satellite DNAs Unveil Clues about the Ancestry and Composition of B Chromosomes in Three Grasshopper Species. Genes (Basel) 2018; 9:genes9110523. [PMID: 30373193 PMCID: PMC6265867 DOI: 10.3390/genes9110523] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2018] [Revised: 10/12/2018] [Accepted: 10/21/2018] [Indexed: 11/16/2022] Open
Abstract
Supernumerary (B) chromosomes are dispensable genomic elements occurring frequently among grasshoppers. Most B chromosomes are enriched with repetitive DNAs, including satellite DNAs (satDNAs) that could be implicated in their evolution. Although studied in some species, the specific ancestry of B chromosomes is difficult to ascertain and it was determined in only a few examples. Here we used bioinformatics and cytogenetics to characterize the composition and putative ancestry of B chromosomes in three grasshopper species, Rhammatocerus brasiliensis, Schistocerca rubiginosa, and Xyleus discoideus angulatus. Using the RepeatExplorer pipeline we searched for the most abundant satDNAs in Illumina sequenced reads, and then we generated probes used in fluorescent in situ hybridization (FISH) to determine chromosomal position. We used this information to infer ancestry and the events that likely occurred at the origin of B chromosomes. We found twelve, nine, and eighteen satDNA families in the genomes of R. brasiliensis, S. rubiginosa, and X. d. angulatus, respectively. Some satDNAs revealed clustered organization on A and B chromosomes varying in number of sites and position along chromosomes. We did not find specific satDNA occurring in the B chromosome. The satDNAs shared among A and B chromosomes support the idea of putative intraspecific ancestry from small autosomes in the three species, i.e., pair S11 in R. brasiliensis, pair S9 in S. rubiginosa, and pair S10 in X. d. angulatus. The possibility of involvement of other chromosomal pairs in B chromosome origin is also hypothesized. Finally, we discussed particular aspects in composition, origin, and evolution of the B chromosome for each species.
Collapse
Affiliation(s)
- Diogo Milani
- Instituto de Biociências/IB, Departamento de Biologia, UNESP-Universidade Estadual Paulista, Rio Claro, São Paulo 01049-010, Brazil.
| | - Vanessa B Bardella
- Instituto de Biociências/IB, Departamento de Biologia, UNESP-Universidade Estadual Paulista, Rio Claro, São Paulo 01049-010, Brazil.
| | - Ana B S M Ferretti
- Instituto de Biociências/IB, Departamento de Biologia, UNESP-Universidade Estadual Paulista, Rio Claro, São Paulo 01049-010, Brazil.
| | - Octavio M Palacios-Gimenez
- Instituto de Biociências/IB, Departamento de Biologia, UNESP-Universidade Estadual Paulista, Rio Claro, São Paulo 01049-010, Brazil.
- Department of Evolutionary Biology, Evolutionary Biology Center, Uppsala University, 75236 Uppsala, Sweden.
| | - Adriana de S Melo
- Instituto de Ciências Biológicas, Laboratório de Biodiversidade e Genética de Insetos, UPE-Universidade de Pernambuco, Recife 50100-130, Pernambuco, Brazil.
| | - Rita C Moura
- Instituto de Ciências Biológicas, Laboratório de Biodiversidade e Genética de Insetos, UPE-Universidade de Pernambuco, Recife 50100-130, Pernambuco, Brazil.
| | - Vilma Loreto
- Centro de Biociências/CB, Departamento de Genética, UFPE-Universidade Federal de Pernambuco, Recife 50670-901, Pernambuco, Brazil.
| | - Hojun Song
- Department of Entomology, Texas A&M University, 2475 TAMU, College Station, TX 77843-2475, USA.
| | - Diogo C Cabral-de-Mello
- Instituto de Biociências/IB, Departamento de Biologia, UNESP-Universidade Estadual Paulista, Rio Claro, São Paulo 01049-010, Brazil.
| |
Collapse
|
10
|
Jetybayev IY, Bugrov AG, Dzuybenko VV, Rubtsov NB. B Chromosomes in Grasshoppers: Different Origins and Pathways to the Modern B s. Genes (Basel) 2018; 9:genes9100509. [PMID: 30340429 PMCID: PMC6209997 DOI: 10.3390/genes9100509] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2018] [Revised: 09/30/2018] [Accepted: 10/17/2018] [Indexed: 01/05/2023] Open
Abstract
B chromosomes (Bs) were described in most taxa of eukaryotes and in around 11.9% of studied Orthopteran species. In some grasshopper species, their evolution has led to many B chromosome morphotypes. We studied the Bs in nine species (Nocaracris tardus, Nocaracris cyanipes, Aeropus sibiricus, Chorthippus jacobsoni, Chorthippus apricarius, Bryodema gebleri, Asiotmethis heptapotamicus songoricus, Podisma sapporensis, and Eyprepocnemis plorans), analyzing their possible origin and further development. The studied Bs consisted of C-positive or C-positive and C-negative regions. Analyzing new data and considering current hypotheses, we suggest that Bs in grasshoppers could arise through different mechanisms and from different chromosomes of the main set. We gave our special attention to the Bs with C-negative regions and suggest a new hypothesis of B chromosome formation from large or medium autosomes. This hypothesis includes dissemination of repetitive sequences and development of intercalary heterochromatic blocks in euchromatic chromosome arm followed by deletion of euchromatic regions located between them. The hypothesis is based on the findings of the Eyprepocnemis plorans specimens with autosome containing numerous intercalary repeat clusters, analysis of C-positive Bs in Eyprepocnemis plorans and Podisma sapporensis containing intercalary and terminal C-negative regions, and development of heterochromatic neo-Y chromosome in some Pamphagidae grasshoppers.
Collapse
Affiliation(s)
- Ilyas Yerkinovich Jetybayev
- The Federal Research Center Institute of Cytology and Genetics, Russian Academy of Sciences, Siberian Branch, Lavrentjev Ave., 10, 630090 Novosibirsk, Russia.
- Institute of Systematics and Ecology of Animals, Russian Academy of Sciences, Siberian Branch, Frunze str. 11, 630091 Novosibirsk, Russia.
| | - Alexander Gennadievich Bugrov
- Institute of Systematics and Ecology of Animals, Russian Academy of Sciences, Siberian Branch, Frunze str. 11, 630091 Novosibirsk, Russia.
- Novosibirsk State University, Pirogov str., 2, 630090 Novosibirsk, Russia.
| | | | - Nikolay Borisovich Rubtsov
- The Federal Research Center Institute of Cytology and Genetics, Russian Academy of Sciences, Siberian Branch, Lavrentjev Ave., 10, 630090 Novosibirsk, Russia.
- Novosibirsk State University, Pirogov str., 2, 630090 Novosibirsk, Russia.
| |
Collapse
|
11
|
Ruiz-Ruano FJ, Cabrero J, López-León MD, Sánchez A, Camacho JPM. Quantitative sequence characterization for repetitive DNA content in the supernumerary chromosome of the migratory locust. Chromosoma 2017; 127:45-57. [PMID: 28868580 DOI: 10.1007/s00412-017-0644-7] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2017] [Revised: 08/23/2017] [Accepted: 08/24/2017] [Indexed: 12/23/2022]
Abstract
Repetitive DNA is a major component in most eukaryotic genomes but is ignored in most genome sequencing projects. Here, we report the quantitative composition in repetitive DNA for a supernumerary (B) chromosome, in the migratory locust (Locusta migratoria), by Illumina sequencing of genomic DNA from B-carrying and B-lacking individuals and DNA obtained from a microdissected B chromosome, as well as the physical mapping of some elements. B chromosome DNA of 94.9% was repetitive, in high contrast with the 64.1% of standard (A) chromosomes. B chromosomes are enriched in satellite DNA (satDNA) (65.2% of B-DNA), with a single satellite (LmiSat02-176) comprising 55% of the B. Six satDNAs were visualized by FISH on the B chromosome, and the only A chromosome carrying all these satellites was autosome 9, pointing to this chromosome, along with autosome 8 (which shares histone genes with the B) as putative ancestors of the B chromosome. We found several transposable elements (TEs) showing nucleotidic variation specific to B-carrying individuals, which was also present in B-carrying transcriptomes. Remarkably, an interstitial region of the B chromosome included a 17 kb chimera composed of 29 different TEs, suggesting reiterative TE insertion in this B chromosome region.
Collapse
Affiliation(s)
- Francisco J Ruiz-Ruano
- Departamento de Genética, Facultad de Ciencias, Universidad de Granada, Avda. Fuentenueva s/n, 18071, Granada, Spain.
| | - Josefa Cabrero
- Departamento de Genética, Facultad de Ciencias, Universidad de Granada, Avda. Fuentenueva s/n, 18071, Granada, Spain
| | - María Dolores López-León
- Departamento de Genética, Facultad de Ciencias, Universidad de Granada, Avda. Fuentenueva s/n, 18071, Granada, Spain
| | - Antonio Sánchez
- Departamento de Biología Experimental, Universidad de Jaén, Jaén, Spain
| | - Juan Pedro M Camacho
- Departamento de Genética, Facultad de Ciencias, Universidad de Granada, Avda. Fuentenueva s/n, 18071, Granada, Spain
| |
Collapse
|
12
|
Milani D, Ramos É, Loreto V, Martí DA, Cardoso AL, de Moraes KCM, Martins C, Cabral-de-Mello DC. The satellite DNA AflaSAT-1 in the A and B chromosomes of the grasshopper Abracris flavolineata. BMC Genet 2017; 18:81. [PMID: 28851268 PMCID: PMC5575873 DOI: 10.1186/s12863-017-0548-9] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2017] [Accepted: 08/22/2017] [Indexed: 01/03/2023] Open
Abstract
BACKGROUND Satellite DNAs (satDNAs) are organized in repetitions directly contiguous to one another, forming long arrays and composing a large portion of eukaryote genomes. These sequences evolve according to the concerted evolution model, and homogenization of repeats is observed at the intragenomic level. Satellite DNAs are the primary component of heterochromatin, located primarily in centromeres and telomeres. Moreover, satDNA enrichment in specific chromosomes has been observed, such as in B chromosomes, that can provide clues about composition, origin and evolution of this chromosome. In this study, we isolated and characterized a satDNA in A and B chromosomes of Abracris flavolineata by integrating cytogenetic, molecular and genomics approaches at intra- and inter-population levels, with the aim to understand the evolution of satDNA and composition of B chromosomes. RESULTS AflaSAT-1 satDNA was shared with other species and in A. flavolineata, was associated with another satDNA, AflaSAT-2. Chromosomal mapping revealed centromeric blocks variable in size in almost all chromosomes (except pair 11) of A complement for both satDNAs, whereas for B chromosome, only a small centromeric signal occurred. In distinct populations, variable number of AflaSAT-1 chromosomal sites correlated with variability in copy number. Instead of such variability, low sequence diversity was observed in A complement, but monomers from B chromosome were more variable, presenting also exclusive mutations. AflaSAT-1 was transcribed in five tissues of adults in distinct life cycle phases. CONCLUSIONS The sharing of AflaSAT-1 with other species is consistent with the library hypothesis and indicates common origin in a common ancestor; however, AflaSAT-1 was highly amplified in the genome of A. flavolineata. At the population level, homogenization of repeats in distinct populations was documented, but dynamic expansion or elimination of repeats was also observed. Concerning the B chromosome, our data provided new information on the composition in A. flavolineata. Together with previous results, the sequences of heterochromatic nature were not likely highly amplified in the entire B chromosome. Finally, the constitutive transcriptional activity suggests a possible unknown functional role, which should be further investigated.
Collapse
Affiliation(s)
- Diogo Milani
- Departamento de Biologia, UNESP - Univ Estadual Paulista, Instituto de Biociências/IB, Rio Claro, São Paulo CEP 13506-900 Brazil
| | - Érica Ramos
- Departamento de Morfologia, UNESP - Univ Estadual Paulista, Instituto de Biociências/IB, Botucatu, São Paulo Brazil
| | - Vilma Loreto
- Departamento de Genética, UFPE - Univ Federal de Pernambuco, Centro de Biociências/CB, Recife, Pernambuco Brazil
| | | | - Adauto Lima Cardoso
- Departamento de Morfologia, UNESP - Univ Estadual Paulista, Instituto de Biociências/IB, Botucatu, São Paulo Brazil
| | | | - Cesar Martins
- Departamento de Morfologia, UNESP - Univ Estadual Paulista, Instituto de Biociências/IB, Botucatu, São Paulo Brazil
| | - Diogo Cavalcanti Cabral-de-Mello
- Departamento de Biologia, UNESP - Univ Estadual Paulista, Instituto de Biociências/IB, Rio Claro, São Paulo CEP 13506-900 Brazil
| |
Collapse
|
13
|
Ruiz-Estévez M, Ruiz-Ruano FJ, Cabrero J, Bakkali M, Perfectti F, López-León MD, Camacho JPM. Non-random expression of ribosomal DNA units in a grasshopper showing high intragenomic variation for the ITS2 region. INSECT MOLECULAR BIOLOGY 2015; 24:319-330. [PMID: 25565136 DOI: 10.1111/imb.12158] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
We analyse intragenomic variation of the ITS2 internal transcribed spacer of ribosomal DNA (rDNA) in the grasshopper Eyprepocnemis plorans, by means of tagged PCR 454 amplicon sequencing performed on both genomic DNA (gDNA) and RNA-derived complementary DNA (cDNA), using part of the ITS2 flanking coding regions (5.8S and 28S rDNA) as an internal control for sequencing errors. Six different ITS2 haplotypes (i.e. variants for at least one nucleotide in the complete ITS2 sequence) were found in a single population, one of them (Hap4) being specific to a supernumerary (B) chromosome. The analysis of both gDNA and cDNA from the same individuals provided an estimate of the expression efficiency of the different haplotypes. We found random expression (i.e. about similar recovery in gDNA and cDNA) for three haplotypes (Hap1, Hap2 and Hap5), but significant underexpression for three others (Hap3, Hap4 and Hap6). Hap4 was the most extremely underexpressed and, remarkably, it showed the lowest sequence conservation for the flanking 5.8-28S coding regions in the gDNA reads but the highest conservation (100%) in the cDNA ones, suggesting the preferential expression of mutation-free rDNA units carrying this ITS2 haplotype. These results indicate that the ITS2 region of rDNA is far from complete homogenization in this species, and that the different rDNA units are not expressed at random, with some of them being severely downregulated.
Collapse
Affiliation(s)
- M Ruiz-Estévez
- Departamento de Genética, Universidad de Granada, Granada, Spain
| | | | | | | | | | | | | |
Collapse
|
14
|
Montiel EE, Ruiz-Ruano FJ, Cabrero J, Marchal JA, Sánchez A, Perfectti F, López-León MD, Camacho JPM. Intragenomic distribution of RTE retroelements suggests intrachromosomal movement. Chromosome Res 2015; 23:211-23. [PMID: 25605325 DOI: 10.1007/s10577-014-9461-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2014] [Revised: 11/25/2014] [Accepted: 12/18/2014] [Indexed: 11/25/2022]
Abstract
Much is known about the abundance of transposable elements (TEs) in eukaryotic genomes, but much is still unknown on their behaviour within cells. We employ here a combination of cytological, molecular and genomic approaches providing information on the intragenomic distribution and behaviour of non-long terminal repeat (LTR) retrotransposon-like elements (RTE). We microdissected every chromosome in a single first meiotic metaphase cell of the grasshopper Eyprepocnemis plorans and polymerase chain reaction (PCR) amplified a fragment of the RTE reverse transcriptase gene with specific primers. PCR products were cloned and 139 clones were sequenced. Analysis of molecular variance (AMOVA) showed significant intragenomic structure for these elements, with 4.6 % of molecular variance being found between chromosomes. A maximum likelihood tree built with the RTE sequences revealed the frequent presence of two or more elements showing very high similarity and being located on the same chromosome, thus suggesting intrachromosome movement. The 454 pyrosequencing of genomic DNA gave strong support to the microdissection results and provided evidence for the existence of 5' truncated elements. Our results thus indicate a tendency of RTE elements to reinsert into the same chromosome from where they were transcribed, which could be achieved if retrotranscription and insertion takes place immediately after transcription.
Collapse
Affiliation(s)
- Eugenia E Montiel
- Departamento de Genética, Facultad de Ciencias, Universidad de Granada, Granada, 18071, Spain,
| | | | | | | | | | | | | | | |
Collapse
|
15
|
Vittorazzi SE, Lourenço LB, Recco-Pimentel SM. Long-time evolution and highly dynamic satellite DNA in leptodactylid and hylodid frogs. BMC Genet 2014; 15:111. [PMID: 25316286 PMCID: PMC4201667 DOI: 10.1186/s12863-014-0111-x] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2014] [Accepted: 10/06/2014] [Indexed: 12/18/2022] Open
Abstract
BACKGROUND Satellite DNA sequences are the most abundant components of heterochromatin and are repeated in tandem hundreds to thousands of times in the genome. However, the number of repeats of a specific satellite family can vary even between the genomes of related species or populations. The PcP190 satellite DNA family was identified in the genome of the leptodactylid frog Physalaemus cuvieri, which showed to be derived most likely from the 5S rDNA in an ancestral species. In this study, we investigate the presence of the PcP190 satellite DNA in several P. cuvieri populations and in four closely related species at the chromosomal and molecular level. Furthermore, we investigate the occurrence of this satellite DNA in the genomes of P. marmoratus as well as in representative species of the leptodactylid genus Leptodactylus (L. latrans) and the hylodid family (Crossodactylus gaudichaudii), all with the aim of investigating if the PcP190 satellite DNA presents or not a restricted distribution. RESULTS The PcP190 satellite DNA was detected in all the analyzed species. Some of them exhibited particular sequence differences, allowing the identification of species-specific groups of sequences, but in other species, the sequences were more conserved. However, in a general analysis, conserved and variable domains have been recognized within the PcP190 monomer. The chromosomal analysis performed on P. cuvieri populations and closely related species revealed high variability of the satellite DNA amount and its chromosomal location, which has always been coincident with regions of centromeric/pericentromeric heterochromatin. CONCLUSION The PcP190 satellite DNA was found in representatives of two families, Leptodactylidae and Hylodidae, indicating that these sequences are widely distributed and conserved in these frogs. There is a pattern of non-random variation within the repeating units, indicating interplay between stochastic events and selective pressure along the PcP190 sequences. Karyotypic differences involving the PcP190 satellite DNA prove to be highly dynamic on the chromosomes of the Physalaemus and its differential accumulation has contributed to the differentiation process of the Z and W sex chromosomes in P. ephippifer.
Collapse
Affiliation(s)
- Stenio Eder Vittorazzi
- Departamento de Biologia Estrutural e Funcional, Instituto de Biologia, Universidade Estadual de Campinas, 13083-863 Campinas, São Paulo, Brazil.
| | - Luciana Bolsoni Lourenço
- Departamento de Biologia Estrutural e Funcional, Instituto de Biologia, Universidade Estadual de Campinas, 13083-863 Campinas, São Paulo, Brazil.
| | - Shirlei Maria Recco-Pimentel
- Departamento de Biologia Estrutural e Funcional, Instituto de Biologia, Universidade Estadual de Campinas, 13083-863 Campinas, São Paulo, Brazil.
| |
Collapse
|
16
|
Pansonato-Alves JC, Serrano ÉA, Utsunomia R, Camacho JPM, da Costa Silva GJ, Vicari MR, Artoni RF, Oliveira C, Foresti F. Single origin of sex chromosomes and multiple origins of B chromosomes in fish genus Characidium. PLoS One 2014; 9:e107169. [PMID: 25226580 PMCID: PMC4165761 DOI: 10.1371/journal.pone.0107169] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2014] [Accepted: 08/08/2014] [Indexed: 12/17/2022] Open
Abstract
Chromosome painting with DNA probes obtained from supernumerary (B) and sex chromosomes in three species of fish genus Characidium (C. gomesi, C. pterostictum and C. oiticicai) showed a close resemblance in repetitive DNA content between B and sex chromosomes in C. gomesi and C. pterostictum. This suggests an intraspecific origin for B chromosomes in these two species, probably deriving from sex chromosomes. In C. oiticicai, however, a DNA probe obtained from its B chromosome hybridized with the B but not with the A chromosomes, suggesting that the B chromosome in this species could have arisen interspecifically, although this hypothesis needs further investigation. A molecular phylogenetic analysis performed on nine Characidium species, with two mtDNA genes, showed that the presence of heteromorphic sex chromosomes in these species is a derived condition, and that their origin could have been unique, a conclusion also supported by interspecific chromosome painting with a CgW probe derived from the W chromosome in C. gomesi. Summing up, our results indicate that whereas heteromorphic sex chromosomes in the genus Characidium appear to have had a common and unique origin, B chromosomes may have had independent origins in different species. Our results also show that molecular phylogenetic analysis is an excellent complement for cytogenetic studies by unveiling the direction of evolutionary chromosome changes.
Collapse
Affiliation(s)
- José Carlos Pansonato-Alves
- Universidade Estadual Paulista (UNESP), Instituto de Biociências/IB, Departamento de Morfologia, Botucatu, São Paulo, Brazil
- * E-mail:
| | - Érica Alves Serrano
- Universidade Estadual Paulista (UNESP), Instituto de Biociências/IB, Departamento de Morfologia, Botucatu, São Paulo, Brazil
| | - Ricardo Utsunomia
- Universidade Estadual Paulista (UNESP), Instituto de Biociências/IB, Departamento de Morfologia, Botucatu, São Paulo, Brazil
| | | | - Guilherme José da Costa Silva
- Universidade Estadual Paulista (UNESP), Instituto de Biociências/IB, Departamento de Morfologia, Botucatu, São Paulo, Brazil
| | - Marcelo Ricardo Vicari
- Universidade Estadual de Ponta Grossa, Departamento de Biologia Estrutural, Molecular e Genética, Ponta Grossa, Paraná, Brazil
| | - Roberto Ferreira Artoni
- Universidade Estadual de Ponta Grossa, Departamento de Biologia Estrutural, Molecular e Genética, Ponta Grossa, Paraná, Brazil
| | - Cláudio Oliveira
- Universidade Estadual Paulista (UNESP), Instituto de Biociências/IB, Departamento de Morfologia, Botucatu, São Paulo, Brazil
| | - Fausto Foresti
- Universidade Estadual Paulista (UNESP), Instituto de Biociências/IB, Departamento de Morfologia, Botucatu, São Paulo, Brazil
| |
Collapse
|
17
|
Ruiz-Estévez M, Cabrero J, Camacho JPM, López-León MD. B chromosomes in the grasshopper Eyprepocnemis plorans are present in all body parts analyzed and show extensive variation for rDNA copy number. Cytogenet Genome Res 2014; 143:268-74. [PMID: 25138924 DOI: 10.1159/000365797] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/26/2014] [Indexed: 11/19/2022] Open
Abstract
B chromosomes in the grasshopper Eyprepocnemis plorans are considered to be mitotically stable, because all meiotic (primary spermatocytes and oocytes) or mitotic (embryos, ovarioles, and gastric caecum) cells analyzed within the same individual show the same B chromosome number. Nothing is known, however, about body parts with somatic tissues with no mitotic activity in adult individuals, constituting the immense majority of their body. Therefore, we investigated whether B chromosomes are present in 8 non-mitotically active somatic body parts from both sexes in addition to ovarioles and testes by PCR analysis of 2 B-specific molecular markers. We also elucidated the number of B chromosomes that an individual carried through quantifying the B-located rDNA copy number by qPCR. Our results indicated the amplification of both B-specific markers in all analyzed body parts. However, we found high variation between males for the estimated number of rDNA units in the B chromosomes. These results demonstrate the presence of B chromosomes in all body parts from the same individual and suggest a high variation in the rDNA content of the B chromosomes carried by different individuals from the same population, presumably due to unequal crossovers during meiosis.
Collapse
Affiliation(s)
- Mercedes Ruiz-Estévez
- Departamento de Genética, Facultad de Ciencias, Universidad de Granada, Granada, Spain
| | | | | | | |
Collapse
|
18
|
Ruiz-Estévez M, Badisco L, Broeck JV, Perfectti F, López-León MD, Cabrero J, Camacho JPM. B chromosomes showing active ribosomal RNA genes contribute insignificant amounts of rRNA in the grasshopper Eyprepocnemis plorans. Mol Genet Genomics 2014; 289:1209-16. [PMID: 24997085 DOI: 10.1007/s00438-014-0880-y] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2014] [Accepted: 06/14/2014] [Indexed: 11/25/2022]
Abstract
The genetic inertness of supernumerary (B) chromosomes has recently been called into question after finding several cases of gene activity on them. The grasshopper Eyprepocnemis plorans harbors B chromosomes containing large amounts of ribosomal DNA (rDNA) units, some of which are eventually active, but the amount of rRNA transcripts contributed by B chromosomes, compared to those of the standard (A) chromosomes, is unknown. Here, we address this question by means of quantitative PCR (qPCR) for two different ITS2 amplicons, one coming from rDNA units located in both A and B chromosomes (ITS2(A+B)) and the other being specific to B chromosomes (ITS2(B)). We analyzed six body parts in nine males showing rDNA expression in their B chromosomes in the testis. Amplification of the ITS2(B) amplicon was successful in RNA extracted from all six body parts analyzed, but showed relative quantification (RQ) values four orders of magnitude lower than those obtained for the ITS(A+B) amplicon. RQ values differed significantly between body parts for the two amplicons, with testis, accessory gland and wing muscle showing threefold higher values than head, gastric cecum and hind leg. We conclude that the level of B-specific rDNA expression is extremely low even in individuals where B chromosome rDNA is not completely silenced. Bearing in mind that B chromosomes carry the largest rDNA cluster in the E. plorans genome, we also infer that the relative contribution of B chromosome rRNA genes to ribosome biogenesis is insignificant, at least in the body parts analyzed.
Collapse
Affiliation(s)
- Mercedes Ruiz-Estévez
- Departamento de Genética, Facultad de Ciencias, Universidad de Granada, 18071, Granada, Spain
| | | | | | | | | | | | | |
Collapse
|
19
|
Silva DMZDA, Pansonato-Alves JC, Utsunomia R, Araya-Jaime C, Ruiz-Ruano FJ, Daniel SN, Hashimoto DT, Oliveira C, Camacho JPM, Porto-Foresti F, Foresti F. Delimiting the origin of a B chromosome by FISH mapping, chromosome painting and DNA sequence analysis in Astyanax paranae (Teleostei, Characiformes). PLoS One 2014; 9:e94896. [PMID: 24736529 PMCID: PMC3988084 DOI: 10.1371/journal.pone.0094896] [Citation(s) in RCA: 77] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2014] [Accepted: 03/19/2014] [Indexed: 12/19/2022] Open
Abstract
Supernumerary (B) chromosomes have been shown to contain a wide variety of repetitive sequences. For this reason, fluorescent in situ hybridisation (FISH) is a useful tool for ascertaining the origin of these genomic elements, especially when combined with painting from microdissected B chromosomes. In order to investigate the origin of B chromosomes in the fish species Astyanax paranae, these two approaches were used along with PCR amplification of specific DNA sequences obtained from the B chromosomes and its comparison with those residing in the A chromosomes. Remarkably, chromosome painting with the one-arm metacentric B chromosome probe showed hybridization signals on entire B chromosome, while FISH mapping revealed the presence of H1 histone and 18S rDNA genes symmetrically placed in both arms of the B chromosome. These results support the hypothesis that the B chromosome of A. paranae is an isochromosome. Additionally, the chromosome pairs Nos. 2 or 23 are considered the possible B chromosome ancestors since both contain syntenic H1 and 18S rRNA sequences. The analysis of DNA sequence fragments of the histone and rRNA genes obtained from the microdissected B chromosomes showed high similarity with those obtained from 0B individuals, which supports the intraspecific origin of B chromosomes in A. paranae. Finally, the population hereby analysed showed a female-biased B chromosome presence suggesting that B chromosomes in this species could influence sex determinism.
Collapse
Affiliation(s)
- Duílio M. Z. de A. Silva
- Departamento de Morfologia, Instituto de Biociências, Universidade Estadual Paulista, Distrito de Rubião Junior, Botucatu, São Paulo, Brazil
| | - José Carlos Pansonato-Alves
- Departamento de Morfologia, Instituto de Biociências, Universidade Estadual Paulista, Distrito de Rubião Junior, Botucatu, São Paulo, Brazil
| | - Ricardo Utsunomia
- Departamento de Morfologia, Instituto de Biociências, Universidade Estadual Paulista, Distrito de Rubião Junior, Botucatu, São Paulo, Brazil
| | - Cristian Araya-Jaime
- Departamento de Morfologia, Instituto de Biociências, Universidade Estadual Paulista, Distrito de Rubião Junior, Botucatu, São Paulo, Brazil
| | | | - Sandro Natal Daniel
- Departamento de Ciências Biológicas, Faculdade de Ciências, Universidade Estadual Paulista, Campus de Bauru, Bauru, São Paulo, Brazil
| | - Diogo Teruo Hashimoto
- CAUNESP, Universidade Estadual Paulista, Campus Jaboticabal, Jaboticabal, São Paulo, Brazil
| | - Cláudio Oliveira
- Departamento de Morfologia, Instituto de Biociências, Universidade Estadual Paulista, Distrito de Rubião Junior, Botucatu, São Paulo, Brazil
| | | | - Fábio Porto-Foresti
- Departamento de Ciências Biológicas, Faculdade de Ciências, Universidade Estadual Paulista, Campus de Bauru, Bauru, São Paulo, Brazil
| | - Fausto Foresti
- Departamento de Morfologia, Instituto de Biociências, Universidade Estadual Paulista, Distrito de Rubião Junior, Botucatu, São Paulo, Brazil
| |
Collapse
|