League of Fat Bastards Cigars has appointed Eddie Tarazona as the company’s new director of sales for the USA and Americas For more than 20 years, Tarazona ran his own cigar company, Tarazona Cigars. In 2022, he campaigned for a U.S. Congressional seat, representing Florida’s 18th Congressional District, though fell short in the race for the seat “I am honored that Adam and the team have entrusted me with his role not just to share our incredible cigars but to be part of the critical cause to tackle men’s mental health,” said Tarazona in a press release “I’m excited for the opportunity to help grow its impact across the US and the Americas It’s a real pleasure to be joining something bigger than myself.” The company said that Tarazona’s first task will be to build a sales team for the U.S the brand’s vp of business development for the Asia-Pacific Middle East and North Africa regions to further develop the brand in other areas of the world “In building his own brand for the last twenty years he has been there and done it,” said Adam J the founder and ceo of League of Fat Bastards Cigars “I love his work ethic and his aim to always find a win for everybody ‘do it with passion’ and that fits so well with our mission-focussed brand.” League of Fat Bastards Cigars currently has distribution in the USA while distributors in several European countries will be receiving their first shipment in the near future I strive to capture the essence of a cigar and the people behind them in my work – every cigar you light up is the culmination of the work of countless people and often represents generations of struggle and stories it’s about so much more than the cigar – it’s about the story behind it the experience of enjoying the work of artisans and the way that a good cigar can bring people together I’m the public address announcer for the Colorado Rockies and Arizona Diamondbacks during spring training as well as for the Salt River Rafters of the Arizona Fall League and previously the Arizona Rattlers of the Indoor Football League I also work in a number of roles for Major League Baseball I covered the Phoenix and national cigar scene for Examiner.com Utah — A Utah DoorDash driver has now been arrested and faces sexual battery charges after he allegedly touched a Taco Bell employee inappropriately multiple times Angel Tineo Tarazona was arrested on Monday According to court documents obtained by FOX 13 News Saratoga Springs police were called to the Taco Bell at 1300 E The victim called the police to the restaurant to report the incident The victim explained that she gave Tarazona the order he needed when he pointed to the menu behind her she says Tarazona put his hand on her butt and pulled her towards him who denied any inappropriate touching occurred Tarazona explained to police that he had simply been placing his hand on her hip to teach her to dance This story was supported by Tarazona's wife who claimed to have witnessed the interaction and stated that nothing inappropriate occurred officers trespassed Tarazona from the location as they continued their investigation officers received video footage of the incident from a lobby camera inside the Taco Bell Police say the video shows Tarazona pointing to the menu and then placing his hands on the worker and pulling her towards him Angel Tarazona then allegedly went to the bathroom for 2 minutes before calling the victim over again Detectives say he again pointed to the menu and placed his hands on her and pulled her in towards him Tarazona was arrested at his home in Orem Monday and faces a sexual battery charge Utah (KUTV) — A DoorDash driver was arrested after he allegedly touched a Taco Bell employee inappropriately several times as he was picking up an order was arrested on suspicion of sexual battery officers with the Saratoga Springs Police Department responded to Taco Bell on 1300 East State Street Tarazona was still at the scene when officers arrived The caller told oficers that after she gave the order to Tarazona he pointed at the menu and asked a question When she turned around to look at the menu who denied anything happened but said he "put his hand on her hip to teach her to dance," the affidavit states Tarazona's wife claimed she was a witness and that nothing inappropriate had happened officers obtained security footage that showed Tarazona touching the victim inappropriately several times Tarazona was arrested from his home and booked into the Utah County Jail WELLINGTON — Buena Vibra rallied in the second half to defend its title in the U.S. Open Women’s Polo Championship Sunday at National Polo Center's U.S In front of the largest crowd to witness a women's tournament to win its second consecutive women's title La Fe Eastern Hay (Naomi Tachibana Marlough Women's Open awarded $25,000 in prize money from U.S divided between first- ($17,000) and second-place ($8,000) teams Each team also received $2,500 to donate to the polo charity of their choice More: Concord Equity Group rallies to win George Miller Memorial in polo opener La Fe Eastern Hay, led by 10-goaler Hope Arellano, dominated the first half jumping out to leads of 4-1, 6-3 and 7-5 at the half. Arellano scored six goals in the first three chukkers, finishing with eight goals. The momentum shifted to Buena Vibra, which translates to "Good Vibes," in the second half. Buena Vibra scored three consecutive goals in the fourth quarter with Hine tying the game at 8 with 54 seconds left. Hine gave Buena Vibra its first lead of the game early in the fifth chukker when she converted a 30-yard penalty for a 9-8 advantage. Hine finished with eight goals, including four penalty conversions. "I'm so happy, it still hasn't sunk in yet," said Hine, the top scorer with 10 goals. "I'm so happy to be a part of this team. I think we were all a little nervous to play on this incredible field. There was the pressure of being in the finals. In the second half, we went back to playing basics." Buena Vibra outscored La Fe Eastern Hay 3-1 in the final chukker. La Fe Eastern Hay had its share of scoring opportunities but couldn't capitalize. "I don't even know how I am feeling," said 14-year-old Valentina Tarazona, who grew up playing in the Wellington-based Polo Training Foundation. "I can't believe I won. It feels surreal to win again." Her sister Giuliana is 15. "I was a little bit nervous when I got here," said Giuliana, 15. "Definitely in the first chukker I was frozen. It pushed me to get over the nerves when they were leading. We were definitely off in the first half but we came back. It's always been a dream to win with my sister." Cande Fernández Araujo, a 9-goaler, and her horse, Machitos Felpa, swept Most Valuable Player and Best Playing Pony honors. She is rated 10 goals in Argentina. "It's amazing, I've never won the U.S. Open before," Araujo said. "Playing on this field in front of all these people is incredible. It's really a big achievement for me. We just did the little things with force and intelligence to win." The U.S. Open Women's Polo Championship was presented by Brad and Kathy Coors Foundation. The next major women's tournament is the WCT Finals in April at Grand Champions Polo Club. The event honors the legacy of Hall of Famer Sunny Hale, who created the WCT and increased interest in women's polo to what it is today. Women were only allowed to join the USPA in 1972. In the semifinals, Buena Vibra advanced with an 8-6 win over Mint Eco Car Wash and La Fe Eastern Hay advanced with an 11-9 victory over North Star. 06 Apr 2025 10:00:00 GMT?.css-1txiau5-AnswerContainer{color:var(--GlobalColorScheme-Text-secondaryText2);}FC Andorra won 1–0 over SD Tarazona on Sun This is 31 of the Primera Federacion - Group 1 Predicted lineups are available for the match a few days in advance while the actual lineup will be available about an hour ahead of the match The current head to head record for the teams are FC Andorra 1 win(s) Have scored 6 goals in their last 5 matches Who won between FC Andorra and SD Tarazona on Sun 06 Apr 2025 10:00:00 GMT?FC Andorra won 1–0 over SD Tarazona on Sun 06 Apr 2025 10:00:00 GMT.InsightsHave scored 5 goals in their last 5 matches FC Andorra is playing home against SD Tarazona on Sun It wasn’t to be. Despite dominating from start to finish, Barça Atlètic drew 1-1 with SD Tarazona at the Estadi Johan Cruyff this Saturday. Although the visitors went ahead through Llácer, Toni Fernández equalised before the break. The blaugranes couldn’t find a winner in the second half and remain in the relegation zone, now 11 games without a victory with 10 games to go. As expected, Barça looked to control possession, while Tarazona set up defensively and eyed chances on the counter. Breaking their five-man backline proved difficult early on. The first real chance came in the 26th minute when Cédric hit the crossbar after a Dacosta cross. Barça’s misfortune grew when Llácer capitalized on a defensive lapse to make it 1-0 to the visitors. But Toni Fernández responded quickly with a fine left-footed strike from the edge of the box to level the score on 36 minutes, with no further change before half-time. Barça stayed calm and controlled the second half, with Toni and Dacosta providing width and threat. Cédric remained the main danger, having a goal ruled out for a tight offside and narrowly missing another chance. Sergi Milà made several changes—bringing on Trilli, Guille, Jan Virgili, Darvich, and Jofre Torrents—but Barça lacked precision in the final third. Guille had two late chances, and Trilli tested Fuoli with a powerful shot in stoppage time, but the score stayed as it was and it all ended 1-1. Barça Atlètic: Kochen, Farré (Torrents, 79'), Cortés, Sergi Domínguez (Trilli, 66'), Mbacke, Garrido, Soma (Darvich, 79'), Rubén López (Guille Fernández, 66'), Toni Fernández, Raúl Dacosta, Cédric (Virgili, 70').  SD Tarazona: Fuoli, López (Álvarez, 75'), Trilles, Camus, Vadik, Llácer (Chechu Martínez, 40'), Fuentes, Cedeño, Mena, Agüero (Buyla, 75'), Pradera (Areso, 83'). Goals: 0-1, Llácer (min 30); 1-1, Toni Fernández (min 36). Referee: Juncal Moreira. Yellows for Dacosta & Toni Fernández (Barça) and Chechu Martínez, Vadik & Trilles (Tarazona). Your Ads Privacy ChoicesIMDb CRIME 6:01 PM | Updated: Apr 2 FILE - A sign hangs at a Taco Bell on May 23 Declaring a mission to liberate "Taco Tuesday" for all to force Wyoming-based Taco John's to abandon its longstanding claim to the trademark BY KSL TV LEHI — A DoorDash driver is being accused of inappropriately grabbing a Taco Bell employee while asking her a question the incident happened Friday at a Taco Bell in Lehi Saratoga Springs police officers said they received reports from the victim who was touched inappropriately by a DoorDash driver later identified as 59-year-old Angel Tarazona police said that Tarazona and the victim were there and both were separated immediately “The complainant explained after she gave the driver the order (Tarazona) pointed at the menu and asked a question about an item he put his hand on her buttocks and pulled her towards him,” the affidavit stated Tarazona denied inappropriately touching her but admitted to placing his hands on her hips to teach her how to dance Tarazona’s wife claimed to be a witness and said she didn’t see anything inappropriate Police trespassed Tarazona from the Taco Bell and began to further investigate police obtained security footage of when the alleged incident happened “The video shows (Tarazona) points to the menu as if to ask a question he places his hand on her lower back and pulls her body towards him,” the affidavit stated The footage showed Tarazona going into the bathroom for about two minutes and calling for the victim again Tarazona grabbed her again and pulled her toward him Tarazona was located by another police agency at his him Tarazona was booked into the Utah County Jail on suspicion of sexual battery Editor’s Note: The story originally said that Tarazona was arrested on Friday Follow @KSL5TV The dates displayed for an article provide information on when various publication milestones were reached at the journal that has published the article activities on preceding journals at which the article was previously under consideration are not shown (for instance submission All content on this site: Copyright © 2025 Elsevier B.V., its licensors, and contributors. All rights are reserved, including those for text and data mining, AI training, and similar technologies. For all open access content, the relevant licensing terms apply. 15 Feb 2025 16:00:00 GMT?.css-1txiau5-AnswerContainer{color:var(--GlobalColorScheme-Text-secondaryText2);}CD Arenteiro won 2–1 over SD Tarazona on Sat This is 24 of the Primera Federacion - Group 1 The current head to head record for the teams are CD Arenteiro 1 win(s) Haven't scored in their last 2 matches Have scored 5 goals in their last 5 matches Have kept the most clean sheets in the competition (11) Who won between CD Arenteiro and SD Tarazona on Sat 15 Feb 2025 16:00:00 GMT?CD Arenteiro won 2–1 over SD Tarazona on Sat 15 Feb 2025 16:00:00 GMT.InsightsHave scored 5 goals in their last 5 matches CD Arenteiro is playing home against SD Tarazona on Sat Your Ads Privacy ChoicesIMDb, an Amazon company© 1990-2025 by IMDb.com, Inc. Barça Atlètic picked up their first point away from home this season thanks to a goalless draw at Tarazona in a game that could have gone either way with chances for both sides.  The home side were happy to delegate possession to the blaugranes and look to cause problems on the counter attack. A stop start first half saw debut keeper Yaakobishvili make an excellent save for Barça as the home side threatened. The contest was an even one with both teams cancelling each other out.  After the break the game opened up. Unai Hernández almost scored one of the goals of the season but keeper Fuoili saved well. Tarazona had a goal disallowed in a frenetic final period in which Barça also went close to taking the win. Barça Atlètic are next in action in a game against Sestao away from home on Sunday.  Tarazona: Fuoli; Camus, Trilles, Murria, Llacer; Marc Álvarez (Iker Gil 82'), Cedeño (Romero 77'), Buyla, Javi Martín (Pradera 66'); Fuentes (Manu Rico 82') and Cubillas (Areso 66'). Barça Atlètic: Yaakobishvili; Cuenca (Trilli 61'), Olmedo, Cortés, Edu Sánchez; Pau Prim (Pedro Soma 76'), Guillermo (Rubén López 61'), Unai; Toni Fernández (Ureña 61'), Brberá and Iván Cédric (Juan Hernández 87'). 19 Jan 2025 18:30:00 GMT?.css-1txiau5-AnswerContainer{color:var(--GlobalColorScheme-Text-secondaryText2);}SD Amorebieta vs SD Tarazona on Sun This is 20 of the Primera Federacion - Group 1 The current head to head record for the teams are SD Amorebieta 0 win(s) Who won between SD Amorebieta and SD Tarazona on Sun 19 Jan 2025 18:30:00 GMT?SD Amorebieta vs SD Tarazona on Sun 19 Jan 2025 18:30:00 GMT ended in a 2–2 tie.InsightsHave scored 4 goals in their last 5 matches SD Amorebieta is playing home against SD Tarazona on Sun 01 Mar 2025 18:00:00 GMT?.css-1txiau5-AnswerContainer{color:var(--GlobalColorScheme-Text-secondaryText2);}SD Tarazona won 2–1 over Lugo on Sat This is 26 of the Primera Federacion - Group 1 The current head to head record for the teams are Lugo 0 win(s) Have scored 4 goals in their last 5 matches Who won between Lugo and SD Tarazona on Sat 01 Mar 2025 18:00:00 GMT?SD Tarazona won 2–1 over Lugo on Sat 01 Mar 2025 18:00:00 GMT.InsightsHave scored 8 goals in their last 5 matches Lugo is playing home against SD Tarazona on Sat 15 Mar 2025 17:00:00 GMT?.css-1txiau5-AnswerContainer{color:var(--GlobalColorScheme-Text-secondaryText2);}Barca Atletic vs SD Tarazona on Sat This is 28 of the Primera Federacion - Group 1 The current head to head record for the teams are Barca Atletic 2 win(s) Who won between Barca Atletic and SD Tarazona on Sat 15 Mar 2025 17:00:00 GMT?Barca Atletic vs SD Tarazona on Sat 15 Mar 2025 17:00:00 GMT ended in a 1–1 tie.InsightsHave scored 6 goals in their last 5 matches Barca Atletic is playing home against SD Tarazona on Sat 26 Jan 2025 14:30:00 GMT?.css-1txiau5-AnswerContainer{color:var(--GlobalColorScheme-Text-secondaryText2);}SD Tarazona won 1–0 over Zamora on Sun This is 21 of the Primera Federacion - Group 1 The current head to head record for the teams are SD Tarazona 1 win(s) Have scored 7 goals in their last 5 matches Who won between SD Tarazona and Zamora on Sun 26 Jan 2025 14:30:00 GMT?SD Tarazona won 1–0 over Zamora on Sun 26 Jan 2025 14:30:00 GMT.InsightsHave scored 4 goals in their last 5 matches SD Tarazona is playing home against Zamora on Sun Your Ads Privacy ChoicesIMDb, an Amazon company© 1990-2025 by IMDb.com, Inc. 08 Feb 2025 17:00:00 GMT?.css-1txiau5-AnswerContainer{color:var(--GlobalColorScheme-Text-secondaryText2);}Real Sociedad B won 1–0 over SD Tarazona on Sat This is 23 of the Primera Federacion - Group 1 The current head to head record for the teams are SD Tarazona 0 win(s) Have kept the most clean sheets in the competition (10) Have scored 8 goals in their last 5 matches Real Sociedad B have won the previous 3 matches against SD Tarazona Who won between SD Tarazona and Real Sociedad B on Sat 08 Feb 2025 17:00:00 GMT?Real Sociedad B won 1–0 over SD Tarazona on Sat 08 Feb 2025 17:00:00 GMT.InsightsHave scored 5 goals in their last 5 matches SD Tarazona is playing home against Real Sociedad B on Sat ha compartido su trayectoria y vivencias en el programa ‘Eméritos’ de TRECE Con una vida marcada por su labor en Roma tras el Concilio Vaticano II y su episcopado en una de las diócesis más pequeñas de España retiros espirituales y la escritura de su último libro mantiene su compromiso con la oración y el servicio Tras media vida en Roma poniendo en pie los documentos que surgieron del Concilio Vaticano II en el año 2011 pasaría a ser obispo de una de las diócesis más pequeñas de España Pero lejos de una vida retirada se mantiene muy activo  “El Papa quería que entrásemos toda la Iglesia en una sinfonía de oración poniendo el Padrenuestro en el centro de la oración” También encuentra tiempo para visitar su pueblo natal Un pueblo agrícola que destaca por el “espárrago que en esos años fue una fuente de riqueza con tres fábricas Su padre sin embargo era transportista y su madre ama de casa. El mayor de una familia con cuatro hermanos, Eusebio Hernández abandona Cárcar a los doce años para ingresar en el seminario con los Agustinos Recoletos. El obispo emérito de Tarazona vio nacer su vocación a comienzo de los años cincuenta gracias al agustino Jesús Pardo Ojer en el Amazonas y nacido en el mismo pueblo navarro que Eusebio Hernández Jesús Pardo Ojer estaba un día con los niños cerca el río y él estaría rezando y los niños brasileños jugando al fútbol se lanzaron a por el balón y la barca dio la vuelta enseguida y comenzaban a ahogarse los cuatro niños un infarto y murió en la orilla del río Purús La ordenación sacerdotal llegó en un día y año icónico: San Fermín de 1968 Las revueltas de París habían llegado a los seminarios españoles donde incluso los futuros sacerdotes organizaban pequeñas huelgas: “Hicimos una pequeña huelga a un profesor en 1967 Éramos un curso un poco revoltoso y nos dijeron que como pena todos a Filipinas porque entonces se consideraba una penitencia ir a Filipinas” donde los ecos procedentes de París eran más intensos en las universidades: “Viví intensamente estos años en la residencia y compartí la vida con estudiantes me tocó hasta correr delante de los caballos de los grises y meterme debajo de los pupitres porque cuando entraban no había contemplaciones Fueron años de grata experiencia en el que pasamos un poco de miedo los estudios se hacían difíciles porque no era fácil compaginar los estudios con el ambiente que había” En 1975, el obispo emérito de Tarazona fue trasladado a Roma, donde se incorporó a la Congregación para la Vida Consagrada, en el campo de las conferencias que agrupan los institutos religiosos, que en el caso de España es la CONFER donde Eusebio Hernández contribuyó a crear las conferencias episcopales en diferentes países “Cuando llegué habría unas ochenta conferencias episcopales y cuando me fui había 150 Pensar que cuando llegué estaba el Muro de Berlín y no había posibilidad de crear conferencias en los países del Este y me tocó ir a Rusia a crear esta conferencia Y eso de las conferencias me permitió conocer países de Asia la India… Esto nos hacía sensibles a la problemática que iban surgiendo como la Teología de la Liberación la opción por los pobres… estos problemas los vi en primera mano” Eusebio Hernández se ha confesado amante del Concilio Vaticano II: “De Gaulle decía que fue el acontecimiento más importante de la historia de la humanidad salir para ver qué sucede en el mundo para poder responderles” El obispo lo recuerda como años apasionantes pero también de tensiones “al haber cantidad de salidas” de sacerdotes que se secularizaron: “Me tocó cerrar algún Instituto de Teología Intercongregacional en México donde había una Teología de la Liberación fuerte los obispos estaban en desacuerdo y tuve que poner orden en el conflicto” Tras 35 años en Roma, el Papa Benedicto XVI le nombra obispo de Tarazona en 2011 recuerda esta etapa “con un cariño especial” destacando la sencillez y humildad de sus 60.000 feligreses Uno de los recuerdos del obispo fue al final de su ministerio en la diócesis aragonesa, el 2022, cuando estalló la guerra de Ucrania. El Obispado se volcó en la ayuda: “Surgió gente que quería ir a Ucrania a llevarles comestibles ropa… Salieron ocho furgonetas hasta la frontera y se trajeron las furgonetas de vuelta llena de gente” Todos ellos fueron acogidos en el seminario de Tarazona: “Estuvieron en las mejores habitaciones buenas camas… Al principio a alguno le molestó pero les dije que eran los más necesitados Llegamos a tener hasta cien ucranianos en el seminario Era la urgencia más necesaria y había que ayudarles” A Utah DoorDash driver has been arrested and facing sexual battery charges Allegedly he touched a Taco Bell employee inappropriately on multiple accastions In Saratoga Springs the victim had been called the police to the Taco Bell When the police got there both of the parties were in the restaurant The victim states that she was giving Tarazona the order when the victim turned thats when he put his hand on her butt and pulled her towards him he denied it by saying he had been trying to teach her a dance and thats why he placed his hand on her hip Tarazona’s wife claimed she witnessed the interaction and nothing inappropriate occurred officers recived video footage of the incident and it shows that everything the victim said was true Tarazona was arrested at his home in Orem and now faces sexual battery charges Read More Donald Trump has implemented tariffs on several uninhabited islands as a part of an confused effort to establish equal trade The White House published a list of all the countries that have been slapped with “reciprocal tariffs” Social media users were the ones to point out that the list includes uninhabited islands and territory where the only significant settlement was U.S One of the uninhabited islads include Heard and McDonald Islands one of the most remote places on Earth was slapped with 10% “reciprocal tariff.” Along with the British Indian Ocean Terriroty was also slapped with a 10% reciprocal tariff despite the fact the he only inhabitants of the islands are U.S.-U.K Read More a 35 year old man has been arrested after police said he accidentally shot his friend in the groin with an AK-47 has been booked into Salt Lake County Jail as an investigation of aggravated assult resulting in serious injury illegal discharge of a firearm causing injury two counts of domestic violence in the presence off a child and providing alcohol to a minor The 19 year old victim had been at Angell’s house Tuesday night to “hang out” the two “began to drink and smoke weed” thats when “Jason then informed (the victim) that he wanted to show him his new gun.” That is when Angell pulled out the AK-47 The victim staes “he told Jason multiple times to make sure the gun was clear but he did not appear to know how to use the gun in manner that (the victim) was not able to describe pointed the gun at the victim and shot him.” Angell told the investigators that “the gun just went off” and that when she saw the victim with his arms in the air and saying “You just shot me.” Read More  Notifications can be managed in browser preferences. Over 60 people perished in Paiporta when a wave of water rushed down the Poyo canal I would like to be emailed about offers, events and updates from The Independent. Read our Privacy notice The mangled car in which Jorge Tarazona’s three-year-old niece and sister-in-law perished in last week’s catastrophic flooding in Spain now hangs halfway off the ragged edge of a highway His brother managed to survive, clinging to a fence. He and his family had been caught in traffic driving home to Paiporta on Valencia’s southern outskirts, Tarazona said. They had no chance to escape when the tsunami-like wave quickly overflowed the nearby drainage canal and swept away everything in its path. “They did not have time to do anything,” Mr Tarazona told The Associated Press, a week after the 29 October flash floods. “My brother was dragged away and ended up clinging to a fence.” His sister-in-law “could not get out and died with her little girl”. Mr Tarazona had ridden a bike back to the site and taped a note on the car asking for whoever eventually removed the wreck off the side of the highway, to call him. “It all happened so fast,” he said, tears coming to his eyes. “In half an hour the current had carried away the car. There was no time, no time. She managed to send me the location of their car hoping for a rescue. “The next day she was found dead inside,” he said. It’s unclear if the two are included in the official toll of the 217 confirmed dead as fatalities tick up, eight days after the deadliest floods in Spain this century. Paiporta has been labelled by Spanish media as the “ground zero” of the natural disaster that has also left 89 people still missing, while officials say the real figure could be higher. Over 60 people perished in Paiporta when a wave of water rushed down the Poyo canal that cuts through its centre. Frustration over the survivors’ sense of abandonment exploded in Paiporta on Sunday when a crowd greeted Spain’s royals and officials with a barrage of mud and other objects. Prime minister Pedro Sanchez was rushed away and the royal couple had to eventually cancel the visit after speaking to several distraught neighbours amid a chaotic scene. The Civil Guard said they rescued two people who had been trapped in their Paiporta home, almost a miraculous exception among the tragedy. The mayor of Paiporta, a middle-class community of 30,000, pleaded on Tuesday for a “higher authority” to step in and take control of her municipality because the floods had made it impossible to go on. Mayor Maribel Albalat said all the municipal buildings, from the town hall to the local police, had been severely damaged and that many of the local civil servants “are in a state of shock”. “Paiporta is a strong village, but this overwhelms our capacities as a local administration,” she said. The air-throbbing “thup, thup, thup” of the huge, two-propeller Chinook helicopters that have flown overhead with the arrival of the army has added to the post-apocalyptic atmosphere. The destruction, however, went far beyond Paiporta and covered a huge swath of municipalities, above all on the southern flank of Valencia city on the Mediterranean coast. Seventy-eight localities had at least one person die from the floods. Police have expanded their search to the nearby marshes and coastline, where the waters carried some away. The residents, businesses and town councils of the affected localities can apply for financial help from a €10.6bn (£8.8bn) relief package from Spain’s government. The regional Valencia government, which is being slammed for not alerting the populace of the danger in time, has asked the central government in Madrid for €31bn to ensure the recovery. Over a week later, the clean-up goes on to get rid of tons of mud and debris that clog street after street, filling thousands of ground floors, and destroying living rooms and kitchens. Neighbourhoods were left without shops and supermarkets after all their products were ruined. Many houses still don´t have drinking water. An impromptu army of volunteers were the first helpers on the ground, shovelling and sweeping away the sticky brown mire covering everything, and helping to start removing pile after pile of debris that made access to cars impossible in many areas. Authorities eventually mobilised 15,000 soldiers and police reinforcements to help firefighters search for bodies and start extracting thousands of wrecked cars strewn over streets and sunk in canal beds. At every corner, cars are piled on top of one another or smashed into buildings, light poles, trees and bridge overpasses. “There is still so much to do,” said volunteer Juanma Baztan Lopez, who is helping churn through the muck in Catarroja, which borders on Paiporta, in his four-wheel drive. He has helped transport doctors to people in need, deliver essential products, and tow away wrecked cars. “It will take a year to get this back to normal,” he said. Jorge Tarazona attaches a poster to a car in Paiporta, Valencia, Spain 22 Mar 2025 15:00:00 GMT?.css-1txiau5-AnswerContainer{color:var(--GlobalColorScheme-Text-secondaryText2);}SD Tarazona won 2–0 over Sestao on Sat This is 29 of the Primera Federacion - Group 1 The current head to head record for the teams are Sestao 1 win(s) Who won between Sestao and SD Tarazona on Sat 22 Mar 2025 15:00:00 GMT?SD Tarazona won 2–0 over Sestao on Sat 22 Mar 2025 15:00:00 GMT.InsightsHave scored 5 goals in their last 5 matches Sestao is playing home against SD Tarazona on Sat All popular browsers allow zooming in and out by pressing the Ctrl (Cmd in OS X) and + or - keys Or alternatively hold down the Ctrl key and scroll up or down with the mouse hooded individuals arbitrarily detained Carlos Correa a renowned human rights defender and executive director of Venezuelan NGO Espacio Público His detention is solely based on his human rights work similarly to four other defenders currently arbitrarily detained for their activism in Venezuela: Javier Tarazona We demand Nicolás Maduro ensures their immediate and unconditional release and Source: ghanasoccernet.com « Prev Next » Comments (2) Listen to Article Afriyie joined CD Lugo in February on loan Jerry Afriyie during his first game for CD Lugo Black Stars forward Jerry Afriyie made his first appearance for Spanish club CD Lugo in their league game against SD Tarazona The 18-year-old climbed off the bench to star for Lugo despite their 2-1 defeat to Tarazona Afriyie replaced Jon Cabo after the break and immediately made an impact winning a penalty which was converted by Brian Mendoza The visitors raced into a two-goal lead before the hour mark through a brace from Adrián Fuentes who netted his first with a strong header before he added his second from the spot Afriyie joined CD Lugo in February on loan after sealing a deal to Al-Qadsiah from Ghanaian lower-tier side Thoughts FC The Ghana U20 star will spend the rest of the season at Lugo as part of his development before moving to Saudi to connect with his parent club The versatile forward is expected to make the Black Stars team for the upcoming World Cup qualifiers against Chad and Madagascar later in March Watch highlights of Afriyie's debut below: 𝗥𝗘𝗦𝗨𝗠𝗘𝗡 | #LugoTarazona 1-2 | @CDeportivoLugo - @SDTarazona@adryyfuentes (0-1 ⚽ @adryyfuentes (0-2 📺 Partido completo: https://t.co/aMfSPND3Ug#PrimeraFederación | #VersusELearning pic.twitter.com/E7ldxn7WqP WomenAndrea Tarazona and Levante UD Extend Contract Until 2026There are no reactions yet Be the first!Levante UD and goalkeeper Andrea Tarazona (born in Riba-roja 2004) have agreed to extend their contract until June 30 This renewal adds an additional season to the previously signed agreement and officially includes the young goalkeeper as a first-team player This formal change recognizes Tarazona's vital role in achieving the goals of Levante UD's women's first team since her debut The U-20 international for Spain made her first appearance for Levante in the 2022-23 season opener against Alhama de Murcia securing a 3-2 victory for Sánchez Vera's team Tarazona has played 34 official matches across all competitions she recorded the second-best goal-conceding record in the top national league Tarazona joined Levante UD in the summer of 2020 to play for their B team becoming a key figure in the team's project Her development includes regular first-team training marking her as one of the most promising talents to emerge from Levante's academy in recent years SUMMER 2015Won the bronze medal at the Phillips 66 National Championships in the 400 Free Relay with a 3:45.54 ... She finished fourth in the 200 Fly at the Arena Pro Swim Series .. Also made the B Final in the 100 Fly (15th) and the D race in the 50 Fly (28th) .. Had a pair of top-10 finishes at the Los Angeles Invitational placing sixth in the 200 IM championship final (2:18.93) and 10th in the 100 Fly consol race (1:00.90) 2014-2015Swam career-bests in the 100 Fly (52.82) and the 200 IM (1:59.10) at the NCAA Championships .. Set new school records in the 400 Medley Relay (3:34.25) at NCAAs ... Qualified for NCAAs in the 200 Medley Relay while setting a new UCLA record (1:37.17) .. Swam a season-best 2:00.47 in the 200 IM at Pac-12s and won the B Final of the 100 Fly (season-best 52.94) .. Set a season-best 1:55.15 in the 200 Fly A Final to place second .. Earned First-Team CSCAA Scholar All-America honors and Second-Team Pac-12 All-Academic honors .. SUMMER 2014Tarazona swam at the Phillips 66 U.S placing second in the B Final (10th overall) in the 100 Butterfly with a personal-best 59.71 … also had a 10th-place finish and personal-best swim of 2:10.49 in the 200 Fly … part of the fifth-place 400 Freestyle Relay team (3:47.81) and 21st in the 50 Fly (27.40) … also swam at the Los Angeles Invitational finishing seventh in both the 200 Individual Medley (2:19.19) and 100 Fly (1:00.68) 2013-2014Was a CSCAA Honorable Mention All-America after finishing 11th in the 200 Butterfly (1:54.88) at the NCAA Championships … was also a part of the 12th-place 800 Freestyle Relay team (7:06.49) 26th in the 400 Individual Medley (season-best 4:11.74) and 27th in the 100 Fly (personal-best 52.89 4th in UCLA history) … placed second in the 200 Fly at the Pac-12 Championships with a personal-best 1:54.80 which is second in school history … also had three other Top 10 finishes going fifth as part of the 400 Free Relay (3:18.07) and seventh in both the 100 Fly (52.71) and 400 IM (4:14.62) … won six individual events in dual meets including the 200 Fly three times against USC (1:56.50) Arizona State and Arizona … also claimed the 400 IM against ASU and the 100 Fly and 200 Free at San Diego … was a part of the 200 Medley and 400 Medley winning teams against Utah and Washington State … second in the 200 Fly against California Utah and Oregon State/UCSB/FGCU and third versus Stanford … second in the 200 IM against WSU and Arizona and second in the 100 Fly versus WSU and OSU/UCSB/FGCU … finished second in the 200 Fly at the Texas Invitational (1:56.37) also placing 10th in the 400 IM and 100 Fly and 13th in the 200 IM … won the 200 Fly at the UNLV Fall Invitational 12th in the 100 Fly and 17th in the 400 IM … was named Second Team Pac-12 All-Academic and made the Athletic Director’s Honor Roll all three quarters making the finals in the 200 Butterfly (2:12.53 13th) and the 400 Individual Medley (4:52.62 won the Championship final in the 200 IM at the Los Angeles Invitational and was a part of the winning 400 Free Relay squad .. was also second in the 200 Fly and the 400 Medley Relay (4:26.85) and third in the 400 IM and the 100 Fly .. seventh in the 100 Fly (1:02.25) and eighth in the 200 IM (2:19.32) and the 400 IM at the Speedo Grand Challenge in Irvine 2012-2013Earned CSCAA Honorable Mention All-America accolades thanks to a 13th-place finish and school-record breaking swim in the 800 Freestyle Relay (7:06.23) at the NCAA Championships .. also placed 17th in the 200 Butterfly (1:56.61) 22nd in the 400 Individual Medley (4:11.72 third all-time at UCLA) and 38th in the 200 IM (career-best 1:59.31) at the NCAA’s .. 11th in the 400 IM (4:12.50) and 20th in the 200 IM (2:00.76) at the Pac-12 Championships also swimming on the sixth-place 200 Medley Relay (1:39.44) and the 800 Free Relay (7:07.61) teams and setting a personal best in the 100 Fly (53.39) .. out-touching 2012 Olympic bronze medalist Caitlin Leverenz for the second time this season .. San Diego State/Kansas and Oregon State/UC Santa Barbara/Utah and finished second in that event at USC .. finished third in the 200 Fly at the AT&T Winter National Championships (Austin won the 200 Fly at the Arena Invitational in Long Beach also placing third in the 800 Free Relay and fourth in the 200 Medley Relay and the 400 Medley Relay .. won the 400 IM against UC Davis/Washington State .. named to the Athletic Director’s Honor Roll for the Spring and Winter quarters .. also earned Pac-12 All-Academic Honorable Mention accolades SUMMER 2012Competed in the 100/200 Butterfly and 200/400 IM at the U.S 2011-2012Placed 18th in the 200 Fly and 32nd in the 400 IM at the NCAA Championships 11th in the 400 IM (freshman school record 4:12.94) and 12th in the 200 IM at the Pac-12 Championships .. placed second in the 200 Fly and fourth in the 400 IM vs finished second in both the 200 Fly and 100 Fly vs was second in the 400 IM and third in the 200 IM against Stanford .. placed first in the 400 IM and second in both the 200 Fly and 100 Fly against Oregon State/Utah .. finished first in the 200 Breast and was second in the 200 Fly vs placed second in the 400 IM and third in the 200 Fly vs CLAREMONT HIGH SCHOOLCLUB: CLAREMONTMember of the National Junior PanPac team in 2010 where she won the 200m Fly title … All-American in the 100 Fly and 200 IM all four years of high school … Scholastic All-American … First Team all-league … league MVP her senior year … two-time CIF champion in the 100 Fly (‘10/’11) .. also won the 2009 CIF title in the 200 … four-time National Youth Team member … competed on the National Junior World Cup team in 2009 … holds age group records in the 200m and 200-yd fly events (15-16) and IM … competed for the Claremont Club under coach John Ries where she was inducted into the Hall of Fame … daughter of Jennifer and John Tarazona … has two brothers DJ and Dylan … majoring in psychology WELLINGTON — Buena Vibra had to wait a week but it was worth it after winning the rain-delayed U.S. Open Women's Polo Championship Friday at the National Polo Center The youngest team in the tournament rallied in the sixth chukker to defeat 90210 in front of a good crowd to win its Open debut on U.S The final was rained out last Sunday and rescheduled a week later The foursome of 1-goaler Valentina Tarazona 4-goaler Cory Williams and 8-goalers Clara Cassino and Milly Hine was the only team to finish the eight-team tournament presented by Brad and Kathy Coors Foundation Tarazona is the youngest player to win the Women's U.S breaking Hope Arellano's six-year record at age 14 Last year at 12 the Colombian was one of the youngest to compete in the Open as a sub "It means a lot to me to be the youngest player Her being 14 and playing was incredible to me I dreamed about playing the Open since I was a little girl the leading scorer with a team-high eight goals oldest daughter of polo great Adolfo Cambiaso I almost forgot I had to play," Hine said about the rescheduled game "I kept thinking it's a week and then it was like "We met last week for what we thought was going to be the final We met last night again like we do before every game We went out there and showed them what we were capable of and enjoyed it were greeted with pink smoke near the players' tent as they rode off the field after the game Cassino scored back-to-back goals tying the game following up a missed tail shot by Williams and scoring again to give Buena Vibra a 12-11 lead with 1:11 left Hine then scored an insurance goal with a 30-yard penalty conversion It was the first time all four players won a U.S "To win this kind of final is a dream," Cassino said "This is what I have been dreaming my whole life playing this kind of tournament with these teams Winning is the best thing that can happen to a polo player "I think the most important thing are the horses It's not easy to play against that organization We needed a good strategy and to know what we were going to do We needed to be very well-prepared and what we needed to do for this final." 7-goaler Catalina Lavinia and 8-goaler Mia Cambiaso finished the tournament with a 3-2 record and won $8,000 which means "good vibes," knocked off defending champion La Fe (Naomi Tachibana Marlough More: Shariah Harris makes history as first Black woman to play in US Open Women's Polo Championship sponsored by Grand Champions Polo Club President Melissa Ganzi ending Shariah Harris' history-making journey The Work To Ride alum was the first Black woman to compete in the Women's U.S we didn't even know if all of us were going to be playing together "For us to come together so last minute and make it to the semis and to have a real shot at making it to the final I think is definitely something we should be proud of The amount of growth we had from our very first game together to now is crazy NEW: The Effects of the US Foreign Aid Freeze on Freedom House and use of the justice system as a tool for political persecution may amount to crimes against humanity The undersigned organizations wholly reject and express profound concern about the July 2 detention and incarceration of the three human rights defenders (HRDs) from the NGO Fundacion Redes (Fundaredes) as well as activists Rafael Tarazona and Omar García Their detention is the latest event reflecting rapidly escalating political persecution and criminalization of the HRDs in Venezuela who independently monitor and report on the country’s critical human rights situation the Bolivarian Intelligence Service (SEBIN by its Spanish acronym) detained the three activists outside the office of the Prosecutor General of the state of Falcon The activists had been denouncing harassment by SEBIN the previous day the three Fundaredes’ activists were transferred to Caracas and presented before a terrorism tribunal where they were charged with inciting hatred and terrorism—charges frequently used in Venezuela to criminalize HRDs and journalists seeking to defend fundamental freedoms All three activists remained in prison after being denied access to their lawyers their lawyers were denied access to the case filings hampering any possibility of providing an adequate defense we remind Venezuelan authorities that arbitrary detentions and incarceration imposing false charges for political reasons as well as the consistent use of the justice system as a tool for political persecution all amount to possible crimes against humanity as recognized by the United Nations’ Independent Fact-Finding Mission for Venezuela the detention of the Fundaredes’ activists only one day after the publication of the UN Office of the High Commissioner for Human Rights’ (UNOHCHR) report on Venezuela’s feeble compliance with UN recommendations reinforces the Maduro government’s consistent lack of commitment to international cooperation or respect for international standards.  Venezuela should guarantee the basic conditions for HRDs to perform their work freely and safely We demand that Venezuelan authorities immediately release Javier Tarazona and cease to criminalize the critical work of HRDs to document abuses and defend the rights of Venezuelans.  Fundaredes is one of the few organizations that monitor and report on violence and human rights abuses committed by armed groups in Venezuela’s border regions—dangerous and difficult-to-access areas where access to information is scarce Fundaredes has been  documenting the armed conflict between Venezuela’s armed forces and a faction of former FARC dissidents that began in March 2021 in the state of Apure Fundaredes’ activists have consistently faced harassment by Maduro government officials for their work such that the Inter-American Commission on Human Rights (IACHR) issued precautionary measures for Tarazona The conflict in Apure displaced more than 5,000 Venezuelans across the border to the Colombian department of Arauca Fundaredes has faced increasing surveillance and harassment by state security forces for documenting and reporting on grave human rights violations Fundaredes’ documentation on the links between Maduro government officials and the Colombian guerilla groups ELN and FARC dissidents have also fed recent accusations Venezuela’s Prosecutor General Tarek William Saab further justified the activists’ detention arguing that Fundaredes’ reporting is baseless and undermines Venezuela’s security Freedom in the World has been widely used by policymakers Donate today to help us ensure the future of this vital resource.  Eddie Tarazona of Tarazona Cigars has come up short in his bid to be the Republican candidate for Florida’s 18th Congressional District an area that covers a large amount of land to the west and northwest of Lake Okeechobee in central Florida Tarazona billed himself as a conservative leader that is “pro-Constitution pro-freedom and pro-you,” who promised to “take the fight to the establishment and stand for Florida and #AmericaFirst.” Tarazona issued the following statement about the outcome of the primary election to his campaign’s Twitter account: — Eddie Tarazona for Congress (FL-18) (@EddieTarazonaFL) August 24, 2022 was arbitrarily detained on 2 July 2021 after attempting to report harassment from security officers at the Attorney General’s Office in the city of Coro (Western Venezuela) His pre-trial hearing took place on December 16 Javier Tarazona is a prisoner of conscience having been arbitrarily detained for his human rights work Tarazona’s health has seriously deteriorated due to lack of medical treatment We urge the authorities to release him immediately and unconditionally FundaREDES is a Venezuelan human rights organization that promotes and defends human rights in the bordering states of Táchira Its work includes documentation and reporting of human rights violations and human rights abuses by non-state actors in these regions they have focused on documenting the violence in Apure state in what has been claimed to be a conflict with FARC non-demobilized groups Tarazona Cigars next release draws inspiration from the Cigar City The Tarazona La Tampeña will debut in a 6 x 52 toro size that uses a habano criollo wrapper from Jalapa Nicaragua over a Nicaraguan binder and Nicaraguan fillers “pays homage to the women that worked the cigar factories of Tampa during its golden years especially the women in our family that came over from Cuba to Tampa in the late 1800’s due to the cigar trade.” It has an MSRP of $16 per cigar and will be limited to 1,000 boxes of 10 cigars Tarazona says it will begin shipping on June 8 Tarazona says this is the first cigar in its new Discovery Series I am an editor and co-founder of halfwheel.com/Rueda Media I previously co-founded and published TheCigarFeed I have written about the cigar industry for more than a decade covering everything from product launches to regulation to M&A I handle a lot of the behind-the-scenes stuff here at halfwheel wearing sweatshirts year-round and eating gyros Eddie Tarazona of Tarazona Cigars has announced that he is launching a candidacy for a seat in the U.S Tarazona has filed paperwork with the Federal Election Commission and is running as a Republican for Florida’s newly-created 18th Congressional District, an area that covers a large amount of land to the west and northwest of Lake Okeechobee in central Florida. In a Tweet announcing his candidacy Tarazona called himself a conservative leader that is “pro-Constitution pro-freedom and pro-you,” who promised to “take the fight to the establishment and stand for Florida and #AmericaFirst.” He said he plans to launch his website on Monday which will contain more details about his candidacy Tarazona told halfwheel that he will continue to run it during his campaign he will likely be handing off a chunk of the day-to-day duties He said that he feels this is an opportunity to serve his country much like he did during his time in the U.S He added that he thinks having a person from the premium cigar industry in Congress will help the chances of pushing back against FDA regulation and other issues facing the industry as well as provide a perspective as to how a healthy and thriving cigar industry would help the conditions in countries where cigars are produced which in turn could help with international stability and immigration issues Tarazona said that he is the only candidate to have entered the race for the District 18 seat though both federal and state records still show multiple candidates for the seat That appears to be related to a recent redistricting of Florida’s congressional districts and a lag in getting the websites of both the Federal Elections Commission and the Florida Department of State caught up with the changes several of the candidates listed for the District 18 seat have already updated their own websites to show that they are running for other districts following that redrawing of the map It is that redrawing that might also present another challenge for Tarazona as legal challenges have arisen after the state’s legislature recently approved the revised map of its districts a plan that was pushed for and signed into law by Gov Ron DeSantis amidst a significant amount of controversy The governor and Republican legislators have been accused of gerrymandering the districts to reduce the influence of Black voters throughout the state while making the vast majority of districts into Republican strongholds there has been speculation that DeSantis pushed for the changes to help ensure he would win Florida’s electoral votes should he run for president in 2024 as well as moved three districts considered to be highly competitive into ones with a significant partisan lean The net effect has been forecast at a gain of four seats in the House for the GOP including the League of Women Voters of Florida the Black Voters Matter Capacity Building Institute have already filed a lawsuit in state court over the redistricting They say the new map violates the state’s constitution as it dilutes the voting power of minorities while also being drawn to directly benefit one political party A federal lawsuit is also expected to be filed though the generally slow speed at which the court system moves may prevent a decision from being rendered prior to the upcoming elections The groups have also filed a request for a temporary injunction though that is centered around another district in northern Florida The Florida primary election is scheduled for Aug with the winning candidates advancing to the general election on Nov The calendar may have just turned to September but the newest release from Tarazona heading to stores is focused on the end of October The cigar is the Tarazona La Bruja Broomstick The cigar features three Ecuadorian-grown wrappers Underneath those is a Nicaraguan binder and fillers the latter of which come from the country’s Condega and Jalapa regions The wrappers are cut in such a way so as to resemble a broom handle while the exposed fillers are designed to evoke an image of a broom’s bristles The cigar is offered in a 7 1/8 x 57 figurado vitola and comes with an MSRP of $17 cigar or $170 for a box of 10 cigars Tarazona has a follow-up slated for the Christmas season, a 5 x 50 vitola called Misfit Toys that comes in gift-wrapped ten packs featuring art by David Angelo Roman an artist who has done work for the Rick and Morty franchise Tarazona released two brand new limited edition cigars that were named after Eddie Tarazona’s business partner There have been two different sizes announced so far for the Caraballo828 Edición Limitada line so far: Here is what I wrote in my original review a little more than one year ago: then the Tarazona Caraballo828 Silverback Edición Limitada was seemingly custom-built for your enjoyment that aforementioned combination that this blend features so prominently comes at a cost: specifically having negative impacts on the complexity nuance and—albeit to a lesser extent—the balance that is present in the cigars while it shows flashes in the second third the full strength basically waits to hit you virtually all at once in the final third which when combined with the overwhelming earth flavors leads to a very linear profile earth-dominant blends should seek these out immediately but those looking for almost any other flavors will find them in short supply As with the cigars I smoked for my first review this Tarazona Caraballo828 Silverback Edición Limitada is covered in a dark espresso brown wrapper that is sandpaper rough to the touch and features a copious amount of oil While there are just as many very prominent veins running up and down the length as I remember from the samples in my first review this cigar is quite a bit more firm when squeezed than those were The aroma from the wrapper is a combination of bitter dark chocolate while the cold draw brings flavors of coca nibs the Tarazona Caraballo828 Silverback Edición Limitada starts off with a very obvious combination of gritty earth and dark chocolate along with secondary flavors of freshly brewed espresso There is plenty of black pepper to be had on the retrohale albeit not nearly was much as I remember from the first time around nor enough to drown out the distinct pomegranate sweetness that makes itself known early on and never lets up The second half of the Tarazona features a touch less black pepper and spice on both the retrohale and on my tongue Although the main flavors continue to be the same gritty earth and dark chocolate that dominated the first half the flavors somehow seem a bit more rich on the palate and leather flit in and out in various amounts as the cigar burns down while the pomegranate sweetness on the retrohale increases a tad and remains at that constant level until the end of the cigar the Silverback Edición Limitada features an excellent draw throughout after a Dickman cut but the burn needs touching up twice: once in the first third and again at the end of the second third Smoke production was both dense and copious for the entire cigar and the overall strength increased evenly from a very strong mild at the beginning to end a just under the full mark by the end of the cigar When it comes to the imperfect process that makes up the art of aging cigars blends like the Tarazona Caraballo828 Silverback Edición Limitada fascinate me I wondered whether the earth note that dominated the profile melt together with the other flavors enough to provide more balance to the profile While there is still quite a bit of both earth and black pepper in the profile there is also significantly more sweetness on the retrohale and the overall profile is both more balanced and more enjoyable I do think that some more time will probably continue the trend of increasing both the balance and complexity the Silverback Edición Limitada is a much better cigar now than it was when it was released I have worn many hats in my life up to this point: I started out as a photojournalist for the Dallas Morning News and the Fort Worth Star-Telegram then transitioned to photographing weddings—both internationally and in the U.S.—for more than a decade After realizing that there was a need for a cigar website containing better photographs and more in-depth information about each release SmokingStogie quickly became one of the more influential cigar blogs on the internet extremely hard-to-find and expensive cigars and it was one of the predecessors to halfwheel Saint Patrick’s Day is the cue for many cigar companies to release limited edition candela-wrapped cigars and this year Tarazona will join the tradition Its newest release is the Tarazona Corned Beef & Cabbage a 6 x 52 barberpole that uses three different wrappers: Ecuadorian Candela Ecuadorian Connecticut and Mexican San Andrés maduro Underneath is a Nicaraguan binder and Nicaraguan filler There will be just 600 bundles of 15 cigars made and the MSRP is set at $9 per cigar The Tarazona Corned Beef & Cabbage will be shown off at next week’s TPE 2022 trade show and begin shipping to retailers in February will soft launch the Corned Beef & Cabbage before it begins shipping to other stores Metrics details Multi-omic studies combine measurements at different molecular levels to build comprehensive models of cellular systems The success of a multi-omic data analysis strategy depends largely on the adoption of adequate experimental designs and on the quality of the measurements provided by the different omic platforms the field lacks a comparative description of performance parameters across omic technologies and a formulation for experimental design in multi-omic data scenarios we propose a set of harmonized Figures of Merit (FoM) as quality descriptors applicable to different omic data types we formulate the MultiPower method to estimate and assess the optimal sample size in a multi-omics experiment MultiPower supports different experimental settings and includes graphical for experimental design decision-making an algorithm to estimate sample size for machine learning classification problems based on multi-omic data the so-called multi-platform or multi-omics studies are becoming popular the success of a multi-omic project in revealing complex molecular interconnections strongly depends on the quality of the omic measurement and on the synergy between a carefully designed experimental setup and a suitable data integration strategy multi-omic measurements should derive from the same samples and variance distributions similar if the planned approach for data integration relies on correlation networks and analysis expectations are frustrated by underpowered experimental design and the lack of a realistic integration method these tools do not apply to non-sequencing omics (i.e. metabolomics) and are not conceived to compare platform performance nor support multi-omic experimental design choices Figures of Merit (FoM) are performance metrics typically used in analytical chemistry to describe devices and methods and dynamic range; descriptors also applied to omic technologies the definition of each FoM acquires a slightly different specification depending on the omic technology considered and each omic platform possesses different critical FoM comprehensive coverage of the targeted space (i.e. while this is not the case for shotgun proteomics which is strongly biased toward abundant proteins we currently lack both a systematic description of FoM discrepancies across omic assays and a definition of a common performance language to support discussions on the multi-omic experimental design the multi-omics field currently lacks a comparative description of performance metrics across omic technologies and methods to estimate the number of samples required for their multiple applications we propose a formal definition of FoM applicable across several omics and provide a common language to describe the performance of high-throughput methods frequently combined in multi-omic studies We leverage this harmonized quality control vocabulary to develop MultiPower an approach for power calculations in multi-omic experiments applicable to across omics platforms and types of data an R method to obtain the optimal sample size required by ML approaches to achieve a target classification ER together with the MultiPower and MultiML calculations proposed here constitute a framework for quality control and precision in the design of multi-omic experiments a Calibration line for omics measuring the levels of the target features (MS platforms and gene-based sequencing platforms). b Calibration line for omics not measuring concentrations but finding genomic regions, where a biological event occurs (region-based sequencing platforms). The table summarizes critical aspects of each FoM and omic while MS platforms are able to measure hundreds to thousands of metabolites in a single sample Sensitivity in sequencing platforms depends on the number of reads associated with the feature a parameter influenced by sequencing depth Features with an elevated number of reads are more accurately measured and hence smaller relative changes can be detected where the goal is to identify genomic regions where a certain event occurs the above definition is difficult to apply and sensitivity is described in terms of true positive rate or recall the proportion of true sites or regions identified as such given the number of reads in the seq output reproducibility for DNA variant calling is associated with the balance between read coverage at each genome position and the technology sequencing errors or exclude features when repeatedly falling under the LOD LOD is assumed to be zero and data do not contain missing values features with few counts in many samples risk exclusion from downstream analyses The dynamic range of an omic feature indicates the interval of true signal levels that can be measured by the platform, while the linear range represents the interval of true signal levels with a linear relationship between the measured signal value and the true signal value (Fig. 1) These FoM influence the reliability of the quantification value and as detection of the true effect size depends on the width of these ranges linear ranges usually span 3–4 orders of magnitude while dynamic ranges increase to 4–5 orders and can be extended using the isotopic peak of the analytes A combination of analytical methods can increase the dynamic range as different instruments may better capture either high or low concentration metabolites NMR has a high dynamic range and can measure highly abundant metabolites with precision although it is constrained by a high detection limit the dynamic range strongly depends on sequencing depth and values can range from zero counts to up to hundreds of thousands linear range boundaries are difficult to establish and may require the use of calibration RNAs (spike-ins) linear ranges may be feature dependent and affected by sequence GC content which then requires specific normalization while in proteomics selectivity strongly depends on sample complexity which determines whether a given feature is detected or not selectivity relates to competition of fragments to undergo sequencing a highly expressed transcript or a highly accessible chromatin region) may outcompete low-abundant elements this problem is difficult to address by means other than increasing the sequencing depth the application of normalized libraries can partially alleviate selectivity problems; however at the cost of compromising expression level estimations variant detection is compromised at repetitive regions and can be alleviated with strategies that increase read length The coverage of a platform is defined here as the proportion of detected features in the space defined by the type of biomolecule (aka feature space) Targeted MS platforms measure a small subset of compounds with high accuracy; hence the coverage is restricted and lower than for untargeted approaches As sample complexity is frequently much larger than the sampling capacity of current instruments identification is limited to compounds with the highest abundances Repeating sample measurement excluding the features identified in the first run is an efficient strategy to improve coverage although this requires increased instrument runtime and sample amounts Coverage in seq-based methods strongly depends on sequencing depth and can potentially reach the complete feature space associated with each library preparation protocol capture of antisense transcripts imposes strand-specific protocols and microRNAs require specific small RNA protocols coverage also relates to the applied protocol Whole-genome bisulfite sequencing has greater genome-wide coverage of CpGs when compared to Reduced Representation by Bisulfite Sequencing (RRBS) while RRBS and MeDIP provide greater coverage at CpG islands region-based sequencing approaches can cover the whole reference genome with coverage depending on the efficiency of the protocol employed to enrich the targeted regions a Relationship between omic platforms (MS in blue and sequencing in red) MS- and seq-based omics have different key properties that determine their FoM which in turn affect differently to power parameters used as input by MultiPower to compute sample size in multi-omics experiments MultiPower works either with no prior data and with pilot data (red boxes) The blue box contains the parameters to be set by users and purple boxes represent the data types accepted by MultiPower and the rest of power parameters estimated by MultiPower from prior data when available or given by the user MultiPower solves an optimization problem to estimate the optimal sample size and to provide a power study of the experiment (gray boxes) Each color represents a different omic data type b Expected percentage of differentially expressed (DE) features for each omic c Pooled standard deviation (PSD) per gene and per omic for the estimated DE features (pseudo-DE features Boxplots represent the median and interquartile range (IQR) Whiskers depict the minimum and maximum of data without outliers which are the values outside the interval (Q1 − 1.5 × IQR d–g MultiPower results with parameters same sample size in all the omics d Statistical power curves for each omic for sample sizes between 2 and 35 Squared dots indicate power at the optimal sample size (n = 16) e Statistical power curves for each omic when considering different percentiles of PSD Squared dots indicate power for the dispersion used in power calculations for each omic (75th percentile) f Curve relating the initial Cohen’s d to the optimal sample size needed to detect each magnitude of change the red arrows and text highlight the magnitude of change to be detected (Cohen’s d = 1.98 in this case) g Statistical power per omic using the optimal sample size (n = 16) with Cohen’s d = 0.8 and the maximum sample size allowed by the user (n = 4) with Cohen’s d = 1.98 a Given a multi-omic dataset with O different omic types and NO samples per omic and a classification vector Y MultiML obtains the maximum number of common observations Nmax (Nmax ≤ NO) for the combination of omic types specified by the user MultiML then creates subdatasets (ticks) having different number of observations (nt) from the available Nmax and obtains the variables that best explain each tick through LASSO regression Taking provided machine learning algorithm ML and cross-validation method CV the selected variables at each tick are used to predict the classification error rate (ER) on a different random subdatasets of size nt This variable selection-ER prediction process is repeated through m iterations to obtain an average ER and confidence interval for each tick a first-order-smoothed penalized P-spline regression is adjusted to estimate the learning curve that will allow the prediction of classification ERs for sample sizes larger than Nmax b–d MultiML results on the TCGA Glioblastoma data b Predicted sample size for different omics combinations using a target ER of 0.01 Margin of error is given as number of samples c Example of MultiML graphical output with the predictive P-spline for a combination of four omics types An increasing number of observations and omic types were run as input data in MultiML having as ERtarget the ER achieved with complete data for four omics (69 observations) their ER are recovered from the actual data and their deviation from ERtarget is calculated as ΔError The ΔError is plotted against the size of the input data Accuracy in ER predictions increases with the size of input data and the number of omics where ERs rapidly decrease as the number of samples increases to reach a stable classification performance This graph can be used to calculate the number of samples required at different ER levels We concluded that the sample size could be accurately predicted when input data represents ~40–60% of the required sample size the field has yet to address aspects that are essential to understand the complexity of multi-omic analysis such as the definition of performance parameters across omic technologies and the formulation of an experimental design strategy in multi-omic data scenarios we addressed both issues by proposing FoM as a language to compare omic platforms and by providing algorithms to estimate sample sizes in multi-omic experiments aiming at differential features analysis (MultiPower) or at sample classification using ML (MultiML) FoM have traditionally been used in analytical chemistry to describe the performance of instruments and methods These terms can also be intuitively applied to sequencing platforms but we noticed that the meaning and relevance of FoM slightly differ for both types of technologies we explain FoM definitions across omic platforms and discuss which of these metrics are critical to each data type and identification are FoM with critical influence on the number of features comprising the omics dataset which in turn affects the power of the technology to identify features with true signal changes Power diminishes as the number of features increases due to the application of multiple testing corrections to control false positives Reproducibility and dynamic range may also be very different across platforms and these have a direct impact on the within-condition variability and across-condition differences of the study while sequencing depth critically affects many of the described FoM of sequencing platforms the choice of a targeted or untargeted method strongly influences FoM values The number of features detected by the omic platform together with the different measurement variabilities across features and the magnitude of change to be detected represent major components of power calculations for omics data The highly heterogeneous nature of these factors across omics platforms calls for specific methods for power calculations in multi-omic experiments MultiPower solves the optimization problem of obtaining the sample size that minimizes the cost of the multi-omic experiment while ensuring both a required power per omic and a global power MultiPower indicates the sample size required to detect a targeted effect size given a significance threshold MultiPower graphically represents the relationship between sample size and facilitates exploration of alternative experimental design choices Estimates for MultiPower parameters are optimally calculated from pilot data although they can also be manually provided count and binary data to facilitate the integration of omic technologies of different analytical nature Data should have been properly preprocessed and eventual batch effects we showcase MultiPower functionalities using sequencing although the method could be applied to other technologies the optimal sample size can be computed under two different requirements: an equal sample size for all platforms ensuring a common minimal power or different sample size per omic to achieve the same power This is relevant for the choice of downstream statistical analysis Methods that rely on co-variance analysis typically require uniform sample sizes and MultiPower will provide this while revealing the differences in power across data modalities Methods that combine data based on effect estimates allow sample size differences and can benefit from the equally powered effect estimation We illustrate the MultiPower method in three scenarios where different parameterizations are assessed and include both controlled laboratory experiments and cohort data to highlight the general applicability of the method By discussing interpretations of power plots and the factors that contribute to sample size results we provide a means to make informed decisions on experimental design and to control the quality of their integrative analysis The MultiPower approach is not directly applicable to ML methods used for sample classification as in this case the basic parameters of the power calculation—significance threshold and effect size—are not applicable sample size estimation is relevant as multi-omic approaches are frequently used to build classifiers of biological samples feature relationships within the multivariate space are instrumental in ML algorithms The MultiML strategy calculates samples size for multi-omic applications where the classification ER can be used as a measure of performance MultiML is itself a learning algorithm that learns the relationship between sample size and classification error and uses this to estimate the number of observations required to achieve the desired classification performance MultiML has been designed to be flexible for the ML algorithm and to evaluate multiple combinations of omics types in order to identify optimized multi-omic predictors a uniform description of performance parameters across omic technologies and offers computational tools to calculate power and sample size for the diversity of multi-omics applications We anticipate MultiPower and MultiML will be useful resources to boost powered multi-omic studies by the genomics community which we broadly classified into two groups: quantitative or analytical FoM include sensitivity To describe how they apply to omic technologies we distinguish between MS and seq-based platforms we also consider two subgroups: “feature-based” and “region-based” a genome annotation file defines the target features to be quantified the definition of the target feature to be measured (usually genomics regions) is part of the data analysis process Applications of RNA-seq to annotate genomes could be considered a region-based assay We consider here omic assays with a dynamic component regarding the genotype while analytical FoM are defined at the feature level omic platforms by nature measure many features simultaneously and consequently the FoM may not be uniform for all of them we consider FoM globally and discuss how technological or experimental factors affect the FoM of different ranges of features within the same platform The MultiPower R method (Fig. 3b) performs a joint power study that minimizes the cost of a multi-omics experiment while requiring both a minimum power for each omic and an average power for all omics MultiPower calculations are defined for a two groups contrast and implemented in the R package to support the application of the method to single and multiple pairwise comparisons The parameters required to compute power can be estimated from multi-omic available data (pilot data or data from previous studies) or The method considers multiple testing corrections by adjusting the significance level to achieve the indicated FDR MultiPower accepts normally distributed data and optimal sample size (number of replicates or observations per condition) can be computed either requiring the same or allowing different sample sizes for each omic the monetary cost is considered as an additional parameter in the power maximization problem MultiPower can be used to both design a new multi-omic experiment and to assess if an already generated multi-omic dataset provides enough power for statistical analysis MultiPower minimizes the total cost of the multi-omics experiment while ensuring a minimum power per omic (Pi) and a minimum average power for the whole experiment (A). Equation (1) describes the optimization problem to be solved to estimate the optimal number of biological replicates for each omic (xi): ci is the cost of generating a replicate for omic i The effect size is calculated with the fold change (ωi) between both groups as well as the average counts (μi) For binary data with 0/1 or TRUE/FALSE values indicating or if a transcription factor is bound or not the goal is comparing the percentages of 1 or TRUE values between two populations the parameters needed to estimate power are related to the difference between proportions to be detected (the effect size) and the sample size where m is the number of features in a particular omic r1 is the expected number of true detections an assorted range of experimental designs can be found All omic assays may be obtained on the same biological samples or individuals which would result in identical replicates number for all data types this is not always possible due to restrictions in cost or biological material sample size differs among omic types and yet the data are to be analyzed in an integrative fashion MultiPower contemplates these two scenarios Under the first scenario, adding the constraint of equal sample size for all omics (xi = x for all I = 1, …, I) to the optimization problem in Eq. (1) results in a straightforward solution the minimum sample size required to meet the constraint on the minimum power per omic (xi) is calculated and the initial estimation of x is set to x = maxI{xi} the second constraint on the average power is evaluated for x x is increased until the constraint is met the cost per replicate does not influence the optimal solution xopt Under the second scenario, allowing different sample sizes for each omic, the optimization problem in Eq. (1) becomes a nonlinear integer programming problem as the statistical power is a nonlinear function of xi The optimization problem can be transformed into a 0–1 linear integer programming problem by defining the auxiliary variables \(z_n^i\) for each omic i and each possible sample size n from 2 to a fixed maximum value nimax where \(z_n^i = 1\) when the sample size for omic i is n The new linear integer programming problem can be formulated as follows: where fi(n) is the power for sample size n The MultiPower R package implements the described method together with several functionalities to support both the selection of input parameters required by the method and the subsequent interpretation of the results. The R package and user’s manual are freely available at https://github.com/ConesaLab/MultiPower The MultiPower package requires different parameters to compute the optimal sample size As the choice of parameters can be challenging MultiPower can estimate them from pilot or similar existing data The multi-omic pilot data must have two groups and at least two replicates per group The algorithm assumes data are already preprocessed Both normally distributed and count data are accepted We recommend raw count data to be provided for sequencing technologies as MultiPower deals with the sequencing depth bias when count data contains other sources of technical noise we recommend previous transformations to meet normality (e.g. with log or voom transformation) and indicating MultiPower that data are normally distributed We also recommend removing low count features from sequencing data these must be removed or imputed before running MultiPower the power.prop.test() function from stats R package is used MultiPower transforms parameters provided by the user to fit specific arguments required by these functions in such a way that magnitudes are maintained roughly comparable for all data types MultiPower computes the Cohen’s d value for all omic features and applies this value to estimate the set M1 of DE features (those with d > d0) We recommend setting the same Cohen’s d initial value (d0) for all the omics although MultiPower allows different values for each one of them The equivalent to Cohen’s d when comparing proportions in groups A and B is Cohen’s h, which can be defined as h = |φA − φB|, where φi = 2 arcsin √pi, and pi is the proportion of 1 or TRUE values in group i. Classification in Supplementary Table 7 is also valid for Cohen’s h the fold change of pseudo-DE features (ω) is estimated as the percentile P100−k of the fold changes corresponding to pseudo-DE features with a PSD between percentiles Pk−5 and Pk+5 the proportions chosen to estimate power correspond to the P100−k of pseudo-DE features for Cohen’s d and are stored in the delta output parameter of MultiPower As variability is not considered for this data type dispersion power plots are not generated in this case Once the power parameters are obtained for each omic and the user sets the minimum power per omic and the average power for the experiment, the optimization problem in Eq. (4) is completely defined and MultiPower makes use of lpmodeler and Rsymphony R packages to solve it Note that the application of these packages to solve the problem is only needed when the number of replicates for each omic differs MultiPower returns a summary table with the provided and/or estimated power parameters While parameter estimation from previous data is recommended for MultiPower analysis MultiPower requires parameters to be provided by users we provide recommendations for critical MultiPower parameters The average value for the standard deviation per omic feature and condition partially depends on the reproducibility of omics technology Overestimating this parameter guarantees that the sample size fits the power needed but may lead to too large sample sizes a good value for the standard deviation is 1 Value for the expected proportion of DE features per omic should be set according to results seen in similar studies A high percentage is expected for cell differentiation processes or diseases like cancer while small perturbations or other types of diseases may induce fewer changes Typical values for the minimum fold change between conditions and the mean of counts for the DE features when using count data are 2 and 30 but again these values depend on the sequencing depth of the experiment and the magnitude of the expected molecular changes After obtaining the optimal sample size with MultiPower especially when this sample size exceeds the available budget for the experiment: How much reduction of the sample size can we afford without losing too much power If power cannot be decreased but the sample size has to be reduced how will this reduction influence the effect size to be detected Can we remove any omic platform with negligible changes between conditions since this platform imposes a too large optimal sample size To provide answers to these and similar questions, MultiPower returns several diagnostic plots (see Fig. 4d–g for instance) sample size plot shows variations in power as a function of the sample size dispersion plot also displays the power curves but for different dispersion values (PSD in normal data or ϕ parameter in count data) the power for the estimated optimal sample size or the fixed dispersion value is represented by a square dot If the optimal sample size estimated by MultiPower exceeds the available budget researchers may opt for increasing the effect size (given by the Cohen’s d) to be detected and allowing a smaller sample size without modifying the required power the postMultiPower() function can be applied which computes the optimal sample size for different values of Cohen’s d from the initial value set by the user to dmax = mini in I{P90i(d)} where P90i(d) is the 90th percentile of the Cohen’s d values for all the features in omic i This choice ensures sufficient pseudo-DE features to estimate the rest of the parameters needed to compute power Although MultiPower algorithm is essentially defined for a two groups comparison the MultiPower R package supports experimental designs with multiple groups Assuming that a pilot dataset is available the MultiGroupPower() function automatically performs all the possible pairwise comparisons (or those comparisons indicated by the user) and returns both a summary of the power and optimal sample size for each comparison and a numerical and graphical global summary for all the comparisons the global optimal sample size is computed as the maximum of optimal sizes obtained for the individual comparisons the solution given by MultiPower method does not allow a different sample size for each group Users must be aware that different optimal sample size could be obtained in this case each omics data matrix is reduced to its significant variables and analysis is performed maintaining a three-way N × P × O structure and O are the different omics in the study returning either an overall classification ER or a balanced ER calculated on the left-out samples only ER is available as both methods give similar results MultiML also provides different prediction distances for PLS-DA and RF used to assign a category to samples These prediction distances can be defined as a model with H components Given NnewInds new individuals and their corresponding omic data matrix Xnewinds the predicted response variable Ŷnewinds can be computed as follows: where W is a p (variables) × H matrix containing the loading vectors associated with X; D is a p × H matrix containing the regression coefficients of X on its H latent components; and B is an H × n (individuals) matrix containing the regression coefficients of Y on the H latent components associated to X The predicted scores (Tpred) are computed as: The degrees of freedom of this model are the number of observation subsets evaluated (ticks) minus 1 The model is used to predict the sample size required to obtain a given classification error MultiML is a computationally intensive algorithm The function RequiredTimeTest() allows users to estimate the time required to run the full predictive model at the local installation For users who can benefit from parallelization options we have implemented the slurm_creator() R function that creates a sh script to run MultiML calculations in a SLURM cluster The function Previous_CER() allows the utilization of a previous MultiML result with N samples in a new MultiML calculation that expands the size of the prior dataset by M samples thereby significantly accelerating the calculation of the new model The plot also includes the fitted penalized smooth spline model to graphically obtain the sample size for an ERtarget not achievable with the pilot dataset The user can also create a comparative plot of all omic combinations to determine the best contributing data modality to an accurate classification by using Comparative_ERPlot() function Further information on research design is available in the Nature Research Reporting Summary linked to this article The MultiPower and MultiML methods are available at GitHub repository https://github.com/ConesaLab/MultiPower Strategies for integrated analysis of genetic and gene expression variation in cancer: addressing the challenges DNA accessibility and gene expression data to build regulatory maps in an organism A multivariate analysis approach to the integration of proteomic and gene expression data A multiway approach to data integration in systems biology based on Tucker3 and N-PLS A survey of best practices for RNA-seq data analysis ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia esATAC: an easy-to-use systematic pipeline for ATAC-seq data analysis SAAP-RRBS: streamlined analysis and annotation pipeline for reduced representation bisulfite sequencing Using MetaboAnalyst 3.0 for comprehensive metabolomics data analysis Galaxy-M: a Galaxy workflow for processing and analyzing direct infusion and liquid chromatography mass spectrometry-based metabolomics data Experimental design and data-analysis in label-free quantitative LC/MS proteomics: a tutorial with MSqRob Platforms and pipelines for proteomics data analysis and management MethylSig: a whole genome DNA methylation analysis pipeline Andrews S. FASTQC. A Quality Control Tool for High Throughput Sequence Data. http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (2014) Qualimap: evaluating next-generation sequencing alignment data Qualimap 2: advanced multi-sample quality control for high-throughput sequencing data SAMStat: monitoring biases in next generation sequencing data MultiQC: summarize analysis results for multiple tools and samples in a single report Feasibility of sample size calculation for RNA-seq studies Power and sample size calculations for high-throughput sequencing-based experiments Scaling to very very large corpora for natural language disambiguation In Proceedings of the 39th Annual Meeting of the Association for Computational Linguistics 26–33 (Association for Computational Linguistics Predicting sample size required for classification performance Metabolomics: current analytical platforms and methodologies Chang, C.-Y. et al. Protein significance analysis in selected reaction monitoring (SRM) measurements. Mol. Cell. Proteomics 11, M111.014662 https://doi.org/10.1074/mcp.M111.014662 (2012) A two-component model for measurement error in analytical chemistry New figures of merit for comprehensive functional genomics data: the metabolomics case Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry NMR-based Metabolomics P001–P368 (The Royal Society of Chemistry Differential expression in RNA-seq: a matter of depth Evaluation and optimization of metabolome sample preparation methods for Saccharomyces cerevisiae Ultra-high-pressure RPLC hyphenated to an LTQ-Orbitrap Velos reveals a linear relation between peak capacity and number of identified peptides Mass spectrometry-based targeted quantitative proteomics: achieving sensitive and reproducible detection of proteins Strategies for discovery and validation of methylated and hydroxymethylated DNA biomarkers limit of detection and limit of quantitation and dynamic range of Arabidopsis proteome profiling using (15)N metabolic labeling and label-free approaches Kuhn, E. et al. Interlaboratory evaluation of automated, multiplexed peptide immunoaffinity enrichment coupled to multiple reaction monitoring mass spectrometry for quantifying proteins in plasma. Mol. Cell. Proteomics 11, M111.013854 https://doi.org/10.1074/mcp.M111.013854 (2012) Multiple reaction monitoring in mass spectrometry/mass spectrometry for direct analysis of complex mixtures SetupX–a public study design database for metabolomic projects A HUPO test sample study reveals common problems in mass spectrometry-based proteomics Target-decoy search strategy for increased confidence in large-scale protein identifications by mass spectrometry Fragment assignment in the cloud with eXpress-D Mapping and quantifying mammalian transcriptomes by RNA-Seq a comprehensive multi-omics dataset of B-cell differentiation in mouse Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA analysis and interpretation of ‘omics’ data: focus on human endometrium Sample size calculation based on exact test for assessing differential expression analysis in RNA-seq data A direct approach to false discovery rates Sample size for FDR-control in microarray data analysis Estimating the positive false discovery rate under dependence RnaSeqSampleSize: real data based sample size estimation for RNA sequencing Statistical Power Analysis for the Behavioral Sciences (L Quick calculation for sample size while controlling false discovery rate with application to microarray analysis mixOmics: An R package for ‘omics feature selection and multiple data integration An Introduction to Statistical Learning Vol Regularization paths for generalized linear models via coordinate descent Inference using shape-restricted regression splines Download references This work has been funded by FP7 STATegra project agreement 306000 and Spanish MINECO grant BIO2012–40244 work in the Imhof lab has been funded by the (DFG; CIPSM and SFB1064) has been funded by the University of Florida Startup funds Institute for Food and Agricultural Research Munich Center of Integrated Protein Science LMU Munich Biological and Environmental Sciences and Engineering Division Electrical and Mathematical Sciences and Engineering Division King Abdullah University of Science and Technology North-West University (Potchefstroom Campus) and drafted the manuscript; L.B.-N.: implemented and tested the MultiML method and drafted the manuscript; D.G.-C: contributed to FOM analysis and MultiML methods; A.S.: contributed to FOM analysis of proteomics data; A.I.: contributed to FOM analysis of proteomics data; T.H.: discussed FOM analysis of metabolomics data; J.T.: supervised parts of the work; J.A.W.: contributed to FOM definition and analysis; and A.C.: supervised All authors contributed to and approved the final version of the manuscript The authors declare no competing interests Peer review information Nature Communications thanks Timothy Ebbels reviewer(s) for their contribution to the peer review of this work Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations Download citation DOI: https://doi.org/10.1038/s41467-020-16937-8 Anyone you share the following link with will be able to read this content: Sorry, a shareable link is not currently available for this article. Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily. Volume 7 - 2020 | https://doi.org/10.3389/fmars.2020.00569 This article is part of the Research TopicChemically Mediated Interactions Between Marine Macrophytes and MicrobesView all 9 articles The dinoflagellate genus Prorocentrum is globally represented by a wide variety of species found upon benthic and/or epiphytic substrates Many epibenthic Prorocentrum species produce lipophilic polyether toxins some of which act as potent protein phosphatase inhibitors and tumor-promoters associated with Diarrheic Shellfish Poisoning (DSP) Most members of the Prorocentrum lima species complex (PLSC) commonly found in the tropics and sub-tropics are toxigenic Epiphytic and planktonic bacteria co-occur with toxigenic Prorocentrum but reciprocal allelochemical interactions are under-investigated The aim of the present study was to identify the culturable bacteria collected together with isolates of the PLSC from seagrass (Thalassia testudinum) and macroalgae along tropical Atlantic coasts of Mexico and to explore potential species interactions with selected isolates Twenty-one bacterial genera belonging to Proteobacteria and Bacteroidetes were identified by amplification of the 16S rRNA gene marker from nine clonal Prorocentrum cultures with γ-proteobacteria comprising the dominant class A positive correlation was found between the bacterial genera associated with two Prorocentrum clones and the esterified toxin analog DTX1a-D8 but there was no apparent correlation between the other PLSC clones and their associated bacteria with the other five DSP toxins detected No bacteriostatic or allelochemical response was found for cell- and culture medium extracts of five Prorocentrum isolates assayed for bioactivity against Staphylococcus sp was only effective in reducing bacterial load in the initial growth stages but did not yield axenic cultures or lower bacterial cell densities throughout the culture cycle Antibiotic treatment did not impair growth or survival of the dinoflagellate There was no significant correlation between Prorocentrum cell volume or cellular toxin concentration over the entire time-series culture cycle Benthic Prorocentrum and associated bacterial communities comprise highly diverse and characteristic microbiomes upon substrates but this study provides little evidence that allelochemical interactions among Prorocentrum cells and associated bacteria originating from epibenthic substrates play a definable role in growth and toxigenicity The extent to which associated bacteria affect growth and polyketide toxin production by epibenthic Prorocentrum in natural microbiomes or in monoclonal cultures remains a controversial issue the current study aimed to determine the composition and diversity of the bacterial flora associated with unialgal cultures of the PLSC from natural populations and different substrates from seagrass (Thalassia testudinum) and macroalgae from tropical Mexican reef systems The second objective was to identify putative effects of associated bacteria on growth and cellular toxin content and composition throughout a culture cycle under controlled environmental conditions and to evaluate potential allelochemical interactions between Prorocentrum and bacteria in culture was prepared from heat-sterilized seawater stock at salinity 36 Clonal isolates were cultured by incubation at 25 ± 1°C on a 12:12 h light:dark cycle and illumination of 50 μmol photons m–2 s–1 Geographical origin of the monoclonal isolates of the Prorocentrum lima species complex (PLSC) obtained for this study was initially brought into culture and subjected to preliminary analysis including of the associated bacterial community the culture soon exhibited heavy contamination by cyanobacteria and hence was excluded from further analysis of cell morphology Cultures were maintained at 24 ± 1°C on a light:dark cycle of 14:10 h and illumination of 86 μmol photons m–2 s–1 measured with a LI-1000 light sensor (DataLogger LI-COR Specimens of Prorocentrum isolates from stable cultures were analyzed in detail by light and scanning electronic microscopy (SEM) for morphological characteristics. Cells stained with 0.2% Calcofluor White M2R in aqueous solution (Fritz and Triemer, 1985) were studied under epifluorescence microscopy (Axio Scope.A1 Germany) with filter set 18 shift free EX BP 390-420 (excitation) Light micrographs were taken with an Axiocam 506 color digital camera (Zeiss Cultured dinoflagellate cells for SEM were fixed with 2% glutaraldehyde for 90 min Specimens were washed in 1.5 mL distilled water (5°C) and centrifuged at 1200 × g at 5°C for 6 min The wash procedure was repeated four times Samples (1.5 mL each) underwent a graded ethanol dehydration series (10 with centrifugation as for the sample wash at each dehydration step a drop (250 μL) of hexamethyldisilazane (HDMS) air-drying agent was placed upon the sample on aluminum SEM stubs and specimens were gold sputter-coated for 5 min Bacteria for SEM analysis were inoculated into marine broth (Difco United States) and incubated for 6 days under the previously referred conditions Each sample (1.5 mL) was fixed with 100 μL of 25% glutaraldehyde for 45 min Specimens then were washed and dehydrated as described for Prorocentrum cells 1 mL of HDMS:ethanol 1:1 was added to the microtube for 3 min 1 mL of pure HDMS was added for 3 min; 10 μL of sample was pipetted onto round coverslips and gold sputter-coated as above Specimens of dinoflagellates and bacteria were observed with a JEOL JSM6360LV electron microscope equipped with a backscattered electron detector under 8 kV voltage acceleration and at 15 mm working distance Sequence reactions were prepared with big dye (v3.1) (Applied Biosystems United States) following manufacturer’s instructions and processed in a 3730xl DNA Analyzer (Applied Biosystems-Hitachi Sequences were deposited in GenBank with the accession numbers MN853604–MN853648 The software R (R Core Team, 2018) was used to obtain the rarefaction curves to identify the sampling effort of bacterial genera among the isolates A non-metric multidimensional scaling (NMDS) analysis an indirect gradient approach for producing an ordination plot based upon a distance matrix a Bray-Curtis similarity matrix was plotted to detect differences in the distribution of bacterial genera associated with the Prorocentrum isolates excluding Pseudomonas and Microbacterium due to their cosmopolitan character The identified toxins from each isolate were transformed and then fitted onto the NMDS ordination plot Permutation tests (n = 1000) were used to test the significance of vector fits and only significant vectors were depicted (p < 0.1) Dense actively growing Prorocentrum cultures cultivated in plastic Petri plates under the conditions detailed in Section “Sample Collection and Initial Dinoflagellate Culture,” were harvested with a 1 mL micropipette after sub-sampling 2 mL of stirred culture for cell counts in a Sedgewick Rafter cell-counting chamber under a light microscope (Olympus BH2 The remaining cells were transferred into 2 mL cryotubes and centrifuged at 5000 × g at 4°C for 5 min The cell pellets were resuspended in 0.2 μm-filter-sterilized seawater (5°C) and centrifuged again under the same conditions The supernatant was again decanted and the cryotubes were placed in a thermomixer at 100°C for 5 min to inactivate esterase enzymes that could modify the toxin profile Cell pellets were stored frozen (−20°C) for later extraction Cell pellets were resuspended in 500 μL 50% methanol in FastPrep tubes After adding 0.9 g FastPrep lysing matrix D (Thermo Savant cells were homogenized by reciprocal shaking in a FastPrep FP120 (Bio 101 Homogenized samples were centrifuged at 16,000 × g for 15 min at 11°C The supernatant from each sample was transferred into a 0.45 μm pore-size spin-filter (Millipore Ultrafree Germany) and centrifuged at 11,000 × g for 1.5 min at 11°C filtrates were transferred into 2 mL LC-autosampler vials (Agilent Germany) for analysis by liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) The analytical system consisted of an SCIEX- 4000 Q-Trap (Sciex triple quadrupole mass spetrometer equipped with a TurboSpray® interface coupled to an 1100 liquid chromatograph (LC) (Agilent The LC equipment included a solvent reservoir and a temperature-controlled column oven (G1316A) Separation of toxins was achieved following previous protocols (Krock et al., 2008; Nielsen et al., 2013) after injection of 5 μL extract onto a C8 analytical column packed with 3 mm HypercloneTM 3 μm BDS C8 130 Å 50 × 2 mm (Phenomenex The flow rate was 0.2 mL min–1 and gradient elution was performed by eluents A (water Initial conditions were 12 min column equilibration with 5% B followed by a linear gradient to 100% B in 10 min and isocratic elution until 16 min with 100% B The program was then returned to initial conditions until 19 min (total run time: 31 min) Analytical standards of OA (1 ng μL–1) and DTX 2 (DTX2) (500 pg μL–1) from the Institute for Marine Biosciences were used to identify and quantify DSP toxins in extracts of Prorocentrum cells Due to lack of standards for derivatives such as OA diol-ester (OA-D8) and DTX 1-diol ester (DTX1-D8) cell quotas were expressed as OA and DTX1 equivalents considering the following detection limits for OA (47 pg μL–1) Data acquisition and processing were performed with the Analyst Software (Version 1.5 were selected for larger-volume batch cultures based on relative toxin composition These isolates were scaled-up into 2.8 L glass Fernbach flasks and harvested (total volume = 1 L) in active growth phase by centrifugation in 45 mL aliquots at 1600 × g at 10°C for 10 min (Eppendorf 5810R Cell pellets were combined for each isolate into a single tube and freeze-dried The dried pellets were sequentially extracted with hexane and ethanol:water (1:1) after ultrasonication (Sonorex Digitec After the pellet had precipitated at each step the supernatant was transferred into a new scintillation vial and 2 mL of the next solvent was added The supernatant was rotary-evaporated (Laborota 4002 Crude residues were dissolved in their respective solvent (∼500 μL) and stored at −20°C prior to disk application The pooled supernatants for each isolate were passed through a solid phase extraction (SPE) cartridge (Bond Elut PPL The column was first conditioned with methanol the column was washed with Milli Q deionized water to remove unbound compounds and then eluted with 1 mL methanol to release target molecules Cell pellet extracts and SPE-concentrated supernatants were tested by disk diffusion assay (Bauer et al., 1966) for their effect against bacterial growth and methanol] was evaporated on sterile 8 mm paper (two replicates per extract) for ∼30 min at room temperature Cell equivalents were calculated based on the number of harvested cells and the extract volume applied in the assay (PA1: 382.5 Aliquots (60 μL) of overnight-grown (∼108 cells mL–1) bacterial target strains Staphylococcus sp HEL66 were spread over the surface of marine agar plates (Difco United States) with a sterile glass spreader The dried filters containing the extracts were placed symmetrically upon the agar surface (four filters per plate) Paper disks with each solvent served as controls The agar plates were incubated at 20°C for 24 h and the inhibition zones recorded The final 2 L culture (cell concentration: 4.84 × 105 cells L–1) was retained untreated by either cell washing or antibiotics and served as a growth control of the initial conditions from the bulk culture The treated and control cultures were incubated in 5 L glass bottles for acclimation under the same conditions detailed in Section “Sample Collection and Initial Dinoflagellate Culture.” After 4 days the content from each bottle was distributed into 150 × 12 mm plastic Petri dishes after stirring for 20 min for homogenization to yield approximately equal cell numbers among the plates Replicate (n = 5) cultures from each treatment and control were harvested every 3–5 days from Day 0 to 35 The final harvest was extended to the beginning of senescence phase but only two treated replicates were evaluated at Day 47 and no control sample was considered at this point An aliquot of 1.5 mL was pipetted from a well-mixed sample for Prorocentrum cell counts at each time-point for each replicate. Cells were immediately fixed with acidic Lugol’s iodine solution and counted in a Sedgewick-Rafter chamber with an inverted microscope (AxioVert.A1, Carl Zeiss, Oberkochen, Germany) at 100X magnification (according to Reguera et al., 2016) Prorocentrum cell size was measured by the software Capture Express with an HDMI 16MDPX camera (DeltaPix Denmark) and volume calculated by the following equation: Associated bacteria in the well-mixed culture aliquots were fixed with 0.25 mL 20% phosphate-buffered saline (PBS)-formaldehyde and prepared for counting following the fluorescence method of Porter and Feig (1980) The bacterial cells were stained by adding 0.5 mL of preserved sample to 5 mL sterilized water and 25 μL of 4,6-diamidino-2-phenylindole (DAPI with incubation (20 min) at 4°C in the dark 0.2 μm pore size black polycarbonate Nucleopore Track-Etch filter (Sigma–Aldrich Germany) mounted over a 25 mm Whatman GF/C glass microfiber filter (Sigma–Aldrich Filters were placed on a glass slide and the bacteria counted by epifluorescence microscopy (Axioskop 2 plus Statistical analysis was performed with the software Infostat V.2018 (Di Rienzo et al., 2018) with inclusion of mean values and standard deviations An analysis of variance (two-way ANOVA) was followed by a post hoc Tukey’s HSD test for mutually significant differences in interactions among multiple variables; namely and toxin content among the treated (washed cells plus antibiotics and washed cells) and non-treated cultures all these variables were analyzed and compared at each timepoint throughout the 47 day growth experiment All data complied with the assumptions of normality and homoscedasticity Pearson’s correlation coefficient (r) was determined to establish a linear correlation among the same variables Morphometric characteristics (n = 20) of cultured cells of the Prorocentrum lima species complex (PLSC) from the Gulf of Mexico and Caribbean coasts of Mexico observed by epifluorescence microscopy after staining with Calcofluor White M2R Scanning electron microscopy of cultured Prorocentrum showed some dinoflagellate cells with attached bacteria, but most bacteria appeared as free-living coccoid or rod-like forms (sometimes filaments) often clumped in mucilaginous aggregates in the culture medium (Figure 1) Bacteria attached to Prorocentrum cells: (a) PA3; (b) PA19 and isolated from cell aggregates versus culture medium: (c) Thalassospira sp isolated exclusively from cells belonging to PA11; (d) Kocuria sp isolated exclusively from the culture medium of PA11 The phylogenetic tree of bacteria isolated from nine Prorocentrum cell cultures shows the phylogenetic association pattern among four dominant classes: Actinobacteria, Flavobacteriia, α-Proteobacteria, and γ-Proteobacteria (Figure 2) Class γ-Proteobacteria was the most abundant with 11 genera within the orders Alteromonadales Class α-Proteobacteria comprised five genera whereas Actinobacteria presented four genera Genus Euzebyella was the only representative of the class Flavobacteriia Prorocentrum isolate PA17 showed the highest number of bacterial genera (13 of 21 identified genera) whereas PA1 exhibited the lowest abundance Neighbor-joining phylogenetic tree of bacterial isolates from Prorocentrum cell cultures from the Veracruz Reef System and the Mexican Caribbean coast The tree is constructed from 16S rRNA sequences with classification according to the Ribosomal Database Project and with 10,000 bootstrap repetitions (blue dots) The tree scale indicates the evolutionary distance among the isolates The branch color denotes the corresponding bacterial class: Actinobacteria (blue); Flavobacteriia (orange); α-Proteobacteria (red); γ-Proteobacteria (green) the pie charts show the proportion of colonies of each isolate associated with Prorocentrum cell aggregates (D) versus the surrounding culture medium (CM) whereas the colored stacked histogram bars (non-cumulative) indicate the number of colonies of each bacterial species (isolate) found in the nine Prorocentrum cultures Abundance of cultured bacterial genera/species (A) and classes (B) from dinoflagellate cell-aggregates (D) and culture medium (CM) Prorocentrum isolates PA19 and PA26 were clustered mainly by the abundance of Halomonas and Marinobacter while PA2 and PA12 were clustered by Marispirillum and Sphingobium Isolates PA11 and PA17 were clustered primarily by the high abundance of Alteromonas which showed a positive correlation with DTX1a-D8 Other Prorocentrum clones and their associated bacterial isolates were not significantly correlated with the identified DSP toxins Bacterial genera isolated exclusively from specific Prorocentrum isolates and culture compartments: Prorocentrum cell aggregates versus surrounding medium NMDS plot based upon a Bray–Curtis matrix of the distribution of clones of the Prorocentrum lima species complex (PLSC) in relation to the abundance of their associated isolated bacterial genera and DSP toxin composition in culture The scale of the axes is arbitrary as is the orientation of the plot with closeness of the clustered elements indicating the degree of similarity and association The symbol shape denotes the bacterial class; color indicates genera; and size denotes abundance Vector length indicates magnitude and correlation of the toxin with the Prorocentrum clones and their isolated associated bacteria Relative toxin composition (%total concentration) of DSP toxins (OA and DTX analogs) from cultured benthic Prorocentrum isolates from the Veracruz Reef System and the Mexican Caribbean coast The agar disk diffusion assay showed no inhibition zones nor apparent changes in growth of bacteria near the disk area for either Prorocentrum cell pellet extracts or supernatants Solvent controls were similarly negative for growth effects The first assay series did show an apparent reduction in the number of colonies of Staphylococcus sp around the disks for the 50% ethanolic extract of Prorocentrum PA1 cells but subsequent application of extracts in a dilution series (1:10) Since there was no apparent effect on bacterial growth Epifluorescence microscopy after DAPI staining of bacteria in Prorocentrum PA1 cultures showing apparent increase in bacterial cell density over time in all treatments: AB there were no significant differences among treatments (two-way ANOVA p < 0.05 and Tukey’s HSD test Time-series of bacterial cell density from Prorocentrum PA1 cultures under various treatments Error bars indicate mean ± standard deviation (SD) of replicate cell counts Prorocentrum cell volume remained constant through the growth cycle and among treatments with only minor non-significant fluctuations on a daily basis among treatments (two-way ANOVA Time-series of bacterial cell density along the growth curve of Prorocentrum PA1 subjected to alternative treatments but the toxin quota varied significantly during the first 16 days of the experiment Toxin cell quotas in antibiotic-treated (AB) and washed cultures without antibiotics (WT) were substantially lower at Day 4 compared to initial values (Day 0) where the AB culture reached minimum quota (1.58 × 102 fmol cell–1) WT and NT control reached minimum values: 1.36 × 102 and 1.43 × 102 fmol cell–1 During late exponential to stationary growth phase from Day 16 to 47 although AB-cultures reached the highest value (3.36 × 102 fmol cell–1) on Day 25 Variation of total DSP toxin cell quota over time-series growth subjected to alternative treatments Total cell quota (fmol cell– 1) calculated as sum of all detectable DSP toxin analogs: OA Error bars indicate mean ± standard deviation (SD) of replicate measurements These generalities have been defined with respect to dense pelagic blooms in the water column but it is not clear that such growth dynamics apply to consortia of bacteria–dinoflagellates in benthic systems where the “blooms” of dinoflagellates and associated bacteria are more spatially constrained upon and adjacent to the substrates the high differential selection inherent in culture experiments means that the results from co-culture of bacterial assemblages and dinoflagellates isolates cannot be simply extrapolated back to natural benthic systems as representative of in situ bacterial composition and diversity Previous studies on cultured isolates of toxigenic benthic Prorocentrum also noted, as found herein for tropical Atlantic populations, that the microbiota from the bulk culture medium and the phycosphere could be considerably different (Prokic et al., 1998; Perez et al., 2008) Even though five bacterial genera were exclusively isolated from Prorocentrum cell aggregates among dinoflagellate clones from tropical coasts of Mexico the distribution of the identified cultivable bacteria between cell-aggregates and surrounding medium was very similar This is likely attributable to the unique and nutrient-rich microhabitat provided by Prorocentrum cells but genus Roseobacter was not found among bacterial isolates in the current study Isolate PA14 was heavily colonized by cyanobacteria shortly after bacterial isolation which may indicate the association of CFB-member Euzebyella with the senescence of this Prorocentrum in high cell density culture Further experiments are required to determine the algicidal potency of specific bacterial isolates against benthic Prorocentrum clones it is now doubtful that bacteria alone can produce OA analogs they do not possess the complete genetic machinery None of the other Prorocentrum strains and their isolated bacteria showed any correlation with the identified toxins in this study implying that these bacteria do not interfere directly with toxin production or composition in the dinoflagellates This provides further evidence that toxins are produced by Prorocentrum alone and not by their associated extracellular bacteria Other OA-producing Prorocentrum species were also tested under the same conditions but no antimicrobial activity was found; these authors suggested the presence of specific antimicrobial compounds (presumably not DSP toxins) for P Prorocentrum species can either inhibit or stimulate, or have no effect, on bacterial growth, depending on the target bacteria (De Vera et al., 2018) The bioactivity screening in the experiments reported here did not reveal any apparent allelochemical effect on the bacterial growth of Staphylococcus sp HEL66 from either cell-extracts or supernatants from Prorocentrum PA1 Initial results showed a slightly positive growth inhibition but this was not sustained through the dilution series it is possible that higher biomass equivalents in the cell extracts could yield bioactivity against bacteria from Prorocentrum metabolites but not at environmentally realistic exposure concentrations The bacterial isolates for the bioassays originated from the North Sea and were deliberately selected rather than bacterial isolates from the Prorocentrum cultures from Mexico to eliminate any conceivable possibility of resistance or co-evolution from prior exposure Negative evidence of allelochemical interactions of Prorocentrum metabolites antibiotic effect against these two naïve strains does not preclude the possibility that allelochemicals and/or polyketide toxins could be effective growth inhibitors against other potentially competing bacteria in their own native benthic environment and even against other strains of Staphylococcus and Vibrio Further experiments considering extracts prepared from axenic versus non-axenic Prorocentrum cultures would help to resolve the issue of bioactivity against its naturally occurring associated bacteria (i.e. isolated from Prorocentrum cultures or directly from natural assemblages) Disk diffusion assays with purified DSP toxins would demonstrate the specific potency against bacteria in the absence of potentially interfering components in crude extracts This corresponds with the exponential growth phases from PA1 of about 20–25 days under all treatments and this likely interferes with the exposure of these bacteria to the antibiotics The selected antibiotic cocktail and final concentration were chosen with respect to conventional methods to achieve axenic cultures of phytoplankton but were ineffective for dense assemblages of benthic dinoflagellates in culture it is appropriate to be somewhat skeptical of previous references to “axenic” cultures of benthic dinoflagellates where convincing evidence of bacterial absence is not provided The cell toxin quota in strain PA1 decreased at the beginning of the culture cycle, i.e., during acclimation and lag growth phase, but then remained stable through the growth curve, while toxin composition remained stable. Similarly, in a previous study, no significant differences in toxin composition between axenic and non-axenic cultures were found for the closely related toxigenic benthic species P. hoffmannianum (Morton et al., 1994) Further experiments considering N and/or P deprivation and micronutrient and trace element availability would provide a closer estimation of the potential interactions between Prorocentrum and associated bacteria in the natural benthic environment Association of specific bacterial assemblages with the evolution of HABs in the pelagic zone in particular with the high abundance of γ-Proteobacteria may also be represented in the benthic realm for dense aggregations of benthic dinoflagellates it is premature to conclude that the culturable bacterial community yields an appropriate insight into the bacteria–dinoflagellate interactions modulated in natural blooms A metagenomic analysis of the bacterial and microeukaryotic assemblages (including the toxigenic dinoflagellates) is required as the first step for such a comparison High inorganic nutrient loads in culture medium promote optimal growth of cultured benthic dinoflagellates which in turn yield enhanced organic substrates for associated bacterial growth The high component of γ-Proteobacteria in the Prorocentrum cultures may indicate that this microcosm reflects a mature high cell density “bloom” with concomitant selection features The different taxonomic distribution of cultivable bacteria among culture compartments Prorocentrum cell aggregates versus the surrounding medium does suggest that the phycosphere constitutes a unique microenvironment with high selective capacity for bacteria The lack of bioactivity of cell extracts and culture supernatants against selected model bacteria does not rule out the possible production of allelochemical metabolites but does not support arguments for high general potency against bacteria related to phosphatase inhibition by DSP toxins Production of large volume axenic dinoflagellate cultures is essential to unraveling the potential role of bacterial interactions in growth and toxin production This approach would allow to assay the effect of taxonomically identified specific consortia of naturally co-occurring bacteria on axenic cultures of P lima to determine effects on growth and toxin production current studies provide little evidence that extracellular bacteria play a critical role in regulation of DSP toxin production in Prorocentrum or in the modulation of growth under non-nutrient limited conditions in culture The raw data supporting the conclusions of this article will be made available by the authors All authors contributed actively to the preparation of data content and writing of this manuscript and YO conducted the dinoflagellate field sampling LD-R and AC prepared preliminary isolates for culture LD-R and YO performed the morphotaxonomic analysis of Prorocentrum isolates and AC maintained the clonal strains at UNAM MP-Z and SP were responsible for bacterial identification and phylogenetic analysis UT-J and JT conducted the bioactivity assays against bacteria UT-J performed the toxin extractions and data analysis for quantification of DSP toxins under supervision of BK and JT UT-J and AC conceived and designed the time series experiment The coordination of the drafting of the manuscript was led by UT-J Financial contribution for the preparation and publication of this manuscript and associated research activities was provided to AC and JT via the PACES II Research Program (Topic II Coast: WP3) of the Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung The Deutscher Akademischer Austauschdienst (DAAD) made possible the research stay of LD-R at the AWI in order to complete this project on benthic dinoflagellates from Mexico The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest The authors thank Manuel Victoria-Muguira of Dorado Buceo for boat access and logistical support with field sampling and Manuel Rodríguez-Gómez from the Acuario de Veracruz for donating filtered seawater Laura Elena Gómez-Lizárraga from the Scanning Electron Microscopy Laboratory at Instituto de Ciencias del Mar y Limnología provided micrographs and Hugo Pérez-López assisted with SEM sample preparation The help of Annegret Müller for technical assistance in the laboratory and Thomas Max for LC-MS/MS measurements of the DSP toxins at AWI is much appreciated Jennifer Bergemann assisted with the antimicrobial assays within laboratory facilities graciously provided by Tilmann Harder The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2020.00569/full#supplementary-material Potentially toxic epiphytic Prorocentrum (Dinophyceae) species in Greek coastal waters Antimicrobial bioactive compounds from marine algae: a mini review Google Scholar Culturable bacterial flora associated with the dinoflagellate green Noctiluca miliaris during active and declining bloom phases in Northen Arabian Sea Antibiotic suceptibility testing by a standardized single disk method CrossRef Full Text | Google Scholar Toxicity and growth assessments of three thermophilic benthic dinoflagellates (Ostreopsis cf Prorocentrum lima and Coolia monotis) developing in the Southern Mediterranean basin Reproductive compatibility among four global populations of the toxic dinoflagellate Gymnodinium catenatum (Dinophyceae) Master recyclers: features and functions of bacteria associated with phytoplankton blooms Chemical ecology of eukaryotic microalgae in marine ecosystems CrossRef Full Text | Google Scholar Toward automatic reconstruction of a highly resolved Tree of Life Ribosomal database project: data and tools for high throughput rRNA analysis Interactions Between a Dinoflagellate and its Associated Bacterial Microflora: Role of Bacteria in the Toxicity of Prorocentrum lima Ehrenberg (Dodge) Google Scholar The toxic dinoflagellate Prorocentrum lima and its associated bacteria Marine microalgae: promising source for new bioactive compounds Universidad Nacional de Córdoba; Argentina: 2018 Google Scholar Revision of the taxonomy within the genus Prorocentrum CrossRef Full Text | Google Scholar Durán-Riveroll A review on the biodiversity and biogeography of toxigenic benthic marine dinoflagellates of the coasts of Latin America Antimicrobial compounds from eukaryotic microalgae against human pathogens and diseases in aquaculture Google Scholar Epiphytic abundance and toxicity of Prorocentrum lima populations in the Fleet Lagoon CrossRef Full Text | Google Scholar A rapid simple technique utilizaing calcofluor white M2R for the visualization of dinoflagellate thecal plates CrossRef Full Text | Google Scholar Molecular diversity and ecology of microbial plankton PubMed Abstract | CrossRef Full Text | Google Scholar Microbial dynamics during harmful dinoflagellate Ostreopsis cf ovata growth: bacterial succession and viral abundance pattern BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT Google Scholar New insights on cytological and metabolic features of Ostreopsis cf ovata Fukuyo (Dinophyceae): a multidisciplinary approach Taxonomy and phylogeny of the benthic Prorocentrum species (Dinophyceae)-A proposal and review Sulfated diesters of okadaic acid and DTX-1: self-protective precursors of diarrhetic shellfish poisoning (DSP) toxins Polyketides from marine dinoflagellates of the genus Prorocentrum biosynthetic origin and bioactivity of their okadaic acid analogues Media for the culture of oceanic ultraphytoplankton CrossRef Full Text | Google Scholar a polycyclic aromatic hydrocarbon-degrading bacterium isolated from seawater LC-MS-MS aboard ship: tandem mass spectrometry in the search for phycotoxins and novel toxigenic plankton from the North Sea MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets a new marine bacterium isolated from the phycosphere of the toxin-producing dinoflagellate Prorocentrum lima “16S/23S rRNA sequencing,” in Nucleic Acid Techniques in Bacterial Systematics Google Scholar a putative link to diarrhetic shellfish poisoning in Nova Scotia Proceedings of the VIII International Conference on Harmful Algae Wyatt (Vigo: Xunta de Galicia and IOC of UNESCO) Colonization and growth of the toxic dinoflagellate Prorocentrum lima and associated fouling macroalgae on mussels in suspended culture The mechanism of diarrhetic shellfish poisoning toxin production in Prorocentrum spp.: physiological and molecular perspectives A marine algicidal Thalassospira and its active substance against the harmful algal bloom species Karenia mikimotoi molecular phylogeny and okadaic acid production of epibenthic Prorocentrum (Dinophyceae) species from the northern South China Sea Observations on the effects of germanium dioxide on the growth of macro-algae and diatoms Google Scholar Algicidal bacteria in the sea and their impact on algal blooms Google Scholar Microbial community interactions and population dynamics of an algicidal bacterium active against Karenia brevis (Dinophyceae) Google Scholar Environmental effects on the production of okadaic acid from Prorocentrum hoffmannianum Faust I salinity and light intensity on the growth and seasonality of toxic dinoflagellates associated with ciguatera Morphological and biochemical variability of the toxic dinoflagellate Prorocentrum lima isolated from three locations at Heron Island Google Scholar Species boundaries in the toxic dinoflagellate Prorocentrum lima (Dinophyceae based on morphological and phylogenetic characters Antimicrobial activities of polyether compounds of dinoflagellate origins CrossRef Full Text | Google Scholar Morphology and phylogeny of Prorocentrum caipirignum sp a new tropical toxic benthic dinoflagellate toxin composition and pigment content of Prorocentrum lima strains isolated from a coastal lagoon in southern UK pectenotoxin-2 and a novel dinophysistoxin from the DSP-causing marine dinoflagellate Dinophysis acuta - Effects of light Abundance of the benthic dinoflagellate Prorocentrum and the diversity and diarrhetic shellfish toxin production of Prorocentrum lima complex and P Cell cycle and toxin production in the benthic dinoflagellate Prorocentrum lima CrossRef Full Text | Google Scholar Importance of free living and particle asssociated bacteria for the growth of the harmful dinoflagellate Prorocentrum minimum: evidence in culture stages Google Scholar Dynamics of bacterial community structure during blooms of Cochlodinium polykrikoides (Gymnodiniales Changes of bacterial population during the decomposition process of red tide dinoflagellate Cochlodinium polykrikoides in the marine sediment addition of Yellow Loess Google Scholar A survey of epiphytic dinoflagellates from the coastal waters of the island of Hawai’i CrossRef Full Text | Google Scholar Diverse bacterial PKS sequences derived from okadaic acid-producing dinoflagellates CrossRef Full Text | Google Scholar The use of DAPI for identifying aquatic microfloral Google Scholar Bacteria of the genus Roseobacter associated with the toxic dinoflagellate Prorocentrum lima R Core Team (2018) R: A Language and Environment for Statistical Computing Vienna: R Foundation for Statistical Computing Google Scholar Guide for Designing and Implementing a Plan to Monitor Toxin-Producing Microalgae Paris: Intergovernmental Oceanographic Commission Google Scholar The Microbial Community Associated With the Florida Red Tide Dinoflagellate Karenia Brevis: Algicidal and Antagonistic Interactions Google Scholar DNA sequencing with chain-terminating inhibitors “Biological toxins from marine and freshwater microalgae,” in Histamine in Fish and Fishery Products Sañudo-Wilhelmy The role of B vitamins in marine biogeochemistry Dinoflagellatae (Peridineae) in Monographischer Behandlung Google Scholar Isolation of an algicidal bacterium and its effects against the harmful-algal-bloom dinoflagellate Prorocentrum donghaiense (Dinophyceae) Die Naturgeschichte der Arthrodelen Flagellaten. Leipzig: Einleitung und Eklärung der Abbildungen Google Scholar Toxin composition of a Prorocentrum lima strain isolated from the Portuguese coast Effects of different levels of N- and P-deficiency on cell yield protein and carbohydrate dynamics in the benthic dinoflagellate Prorocentrum lima Unbalanced N:P ratios and nutrient stress controlling growth and toxin production of the harmful dinoflagellate Prorocentrum lima (Ehrenberg) Dodge Strategies for culture of “unculturable” bacteria Morphotypes of Prorocentrum lima (Dinophyceae) from Hainan Island South China Sea: morphological and molecular characterization Krock B and Durán-Riveroll LM (2020) Associated Bacteria and Their Effects on Growth and Toxigenicity of the Dinoflagellate Prorocentrum lima Species Complex From Epibenthic Substrates Along Mexican Coasts Copyright © 2020 Tarazona-Janampa, Cembella, Pelayo-Zárate, Pajares, Márquez-Valdelamar, Okolodkov, Tebben, Krock and Durán-Riveroll. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY) distribution or reproduction in other forums is permitted provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited in accordance with accepted academic practice distribution or reproduction is permitted which does not comply with these terms *Correspondence: Allan D. Cembella, YWxsYW4uY2VtYmVsbGFAYXdpLmRl; Lorena M. Durán-Riveroll, bG9yZW5hLmR1cmFuQGF3aS5kZQ==; bGR1cmFuQGNvbmFjeXQubXg= Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher 94% of researchers rate our articles as excellent or goodLearn more about the work of our research integrity team to safeguard the quality of each article we publish Icon DONATE Mauricio Tarazona is the head of Treasury and Trade Solutions (TTS) for Panama with the responsibility to lead the areas of international trade he managed the TTS business for Costa Rica from 2012 to 2015 Tarazona Joined Citi Peru in 2005 as a product manager for trade finance and services covering Peru and Bolivia and since then has taken on new positions with expanded roles In 2010 he became head of the trade and corporate lending business for Peru Tarazona lead the identification and disbursement of major structured trade/supply chain/EAF and lending deals in excess of $1.4 billion financial institutions and public-sector entities Tarazona has over 15 years of banking experience and worked in regional and global financial institutions mainly in the areas of international trade He is currently a member of the board of the ACH Clearing System in Panama (Telered) and throughout the years has participated in business audits around Latin America and business assessments in countries such as Brazil Tarazona holds a BA in economics and an MBA from Universidad de Lima in addition to post-graduate studies in finance and international business © 2025 Americas Society/Council of the Americas